Научная статья на тему 'Осесимметричная задача Лемба для среды Коссера'

Осесимметричная задача Лемба для среды Коссера Текст научной статьи по специальности «Математика»

CC BY
148
20
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
СРЕДА КОССЕРА / ПОВЕРХНОСТНЫЕ ФУНКЦИИ ВЛИЯНИЯ / МЕТОД МАЛОГО ПАРАМЕТРА / ИНТЕГРАЛЬНЫЕ ПРЕОБРАЗОВАНИЯ ЛАПЛАСА И ХАНКЕЛЯ / СВЯЗЬ ПЛОСКОЙ И ОСЕСИММЕТРИЧНОЙ ЗАДАЧ / COSSERAT MEDIUM / SURFACE INflUENCE FUNCTIONS / LAPLACE AND HANKEL INTEGRAL TRANSFORMATIONS / THESMALL PARAMETER METHOD / CONNECTION BETWEEN PLANE AND AXISYMMETRIC PROBLEMS

Аннотация научной статьи по математике, автор научной работы — Чан Ле Тхай, Тарлаковский Д. В.

В статье рассматривается упругое однородное изотропное полупространство, заполненное средойКоссера. На границ е полупространства заданы нормальные давления. В начальный момент време-ни и на бесконечности возмущения отсутствуют. С учетом осевой симметрии разрешающая системауравнений включает в себя три гиперболических уравнения относительно скалярного потенциала иненулевых компонент векторного потенциала и вектора поворота. Решение задачи ищется в видеобобщенных сверток заданного давления с соответствующими поверхностными функциями влияния.Для построения последних применяются преобразования Ханкеля по радиусу и Лапласа по вре-мени. Используется разложение в степенные ряды по малому параметру, характеризующему связьволн сдвига и вращения. Найдены изображения первых двух коэффициентов этих рядов. Соответ-ствующие оригиналы определяются связью плоской и осесимметричной задач. Приведены п римерырасчетов регулярных составляющих функций влияния зернистого композита из алюминиево й дробив эпоксидной матрице

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Чан Ле Тхай, Тарлаковский Д. В.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Axisymmetric Problem Lemba for the Cosserat Medium

The article deals with elastic homogeneous isotropic half-space filled with the Cosserat medium. At the initialinstant of time and at infinity, there are no perturbations. At the boundary of the half-space, normal pressuresare given. All the components of the stress-strain state are supposed to be limited. A cylindrical coordinatesystem is used with a n axis directed inward into the half-space. With allowance for axial symmetry, theresolving system of equations includes three hyperbolic equations with respect to the scalar potential andthe non-zero com pone nts of the vector potential and the rotation vector. The components of displacementvectors, rotation angle, stress tensors and stress moments are related to the potentials by known relationships.The solution of the problem is sought in the form of generalized convolutions of a given pressure with thecorresponding surface influence function s. To construct the latter, Hankel transformations along the radiusand Laplace transformations are applied i n time. We use the expansion in power series for a small parametercharacterizing the connection between the shear and rotation waves. The images of the fir st two coefficientsof t hese series are found. The corresponding originals are determined by the connection between the planeand axisymmetric problems. Examples of calculations of the regular components of the influence of a granularcomposite from an aluminum shot in a n epoxy matrix are given

Текст научной работы на тему «Осесимметричная задача Лемба для среды Коссера»

УДК 539.3

ОСЕСИММЕТРИЧНАЯ ЗАДАЧА ЛЕМБА ДЛЯ СРЕДЫ КОССЕРА

Чан Ле Тхай, Д. В. Тарлаковский

Чан Ле Тхай, аспирант кафедры сопротивления материалов, динамики и прочности машин, Московский авиационный институт (национальный исследовательский университет), Россия, 125993, Москва, Волоколамское шоссе, 4, [email protected]

Тарлаковский Дмитрий Валентинович, доктор физико-математических наук, профессор, заведующий лабораторией динамических испытаний НИИ механики, Московский государственный университет имени М. В. Ломоносова, Россия, 119192, Москва, Мичуринский просп., 1; заведующий кафедрой сопротивления материалов, динамики и прочности машин, Московский авиационный институт (национальный исследовательский университет), Россия, 125993, Москва, Волоколамское шоссе, 4, [email protected]

В статье рассматривается упругое однородное изотропное полупространство, заполненное средой Коссера. На границе полупространства заданы нормальные давления. В начальный момент времени и на бесконечности возмущения отсутствуют. С учетом осевой симметрии разрешающая система уравнений включает в себя три гиперболических уравнения относительно скалярного потенциала и ненулевых компонент векторного потенциала и вектора поворота. Решение задачи ищется в виде обобщенных сверток заданного давления с соответствующими поверхностными функциями влияния. Для построения последних применяются преобразования Ханкеля по радиусу и Лапласа по времени. Используется разложение в степенные ряды по малому параметру, характеризующему связь волн сдвига и вращения. Найдены изображения первых двух коэффициентов этих рядов. Соответствующие оригиналы определяются связью плоской и осесимметричной задач. Приведены примеры расчетов регулярных составляющих функций влияния зернистого композита из алюминиевой дроби в эпоксидной матрице.

Ключевые слова: среда Коссера, поверхностные функции влияния, метод малого параметра, интегральные преобразования Лапласа и Ханкеля, связь плоской и осесимметричной задач.

DOI: Г^/Ш^Л 0.18500/1816-9791-2018-18-4-496-506 ВВЕДЕНИЕ

Общая теория моментной упругости впервые была разработана братьями Коссера (Е. СоББега!, Р. СоББега!) в 1909 г. [1]. В ней в отличие от классической теории упругости движение среды описывается не только вектором перемещения, но также вектором поворота. Краткие обзоры публикаций по тематике упругих сред с микроструктурой представлены в работах [2-7].

В настоящий момент моментная теория упругости привлекает внимание многих исследователей. Необходимо, прежде всего, выделить работы [4-6,8-13], в которых исследованы статические и динамические задачи для таких сред. В [11] рассмотрена динамическая связанная осесимметричная задача микрополярной теории термоупругости для изотропных слоя, полупространства или пространства. А в работах [4-6] получены решения аналогичных нестационарных осесимметричных задач для однородной изотропной среды псевдокоссера.

В [12] дана постановка и проведены аналитические исследования задач о действии нестационарной поверхностной нагрузки на моментную упругую полуплоскость и

построены соответствующие нестационарные поверхностные функции влияния для плоской задачи Лемба. В работе [13] исследовано нестационарное осесимметрич-ное движение упругого моментного полупространства под действием нестационарных нормальных поверхностных перемещений. В силу осесимметричного характера задачи наряду с преобразованием Лапласа по времени применяется двойное преобразование Фурье (или преобразование Ханкеля) по пространственным переменным.

1. ПОСТАНОВКА ЗАДАЧИ

Рассматривается упругое однородное изотропное полупространство, заполненное средой Коссера [14]. Уравнения ее осесимметричного движения при отсутствии массовых сил и моментов аналогично [13] записываем относительно скалярного потенциала р перемещений и ненулевых компонент ф и и векторов перемещений и угла поворота в цилиндрической системе координат Ог$г (г ^ 0, —п < $ ^ п, г ^ 0) с направленной в глубь полупространства осью Ог (точками обозначены производные по времени):

ф

Ш = Ар ф = (7- + а) Аф + 2аи — (7-2 + а) ,

(1)

и

72"2 Аи — 2ав Аф — 4ави — г

-2 / -2

и

^ 0.

При этом ненулевые компоненты и и т вектора перемещений связаны с потенциалами так:

дш дф дш 1 д (гф

и = 7Г — ТТ, т = ^ +

дг дг дг г дг

(2)

а нетривиальные физические компоненты тензоров напряжений и моментов ^ ({£, с} = {г, г}) определяются следующим образом:

ди (дт и\

°гг = дг + V дг + г /

ди

и

= о--

дг г

и (ди дт г V дг дг

дт (ди и

^ = + к — + -дг V дг г

дт ди

^ = (7-2 + а) ^ + (г-2 — а) — + 2аи,

дт ( -2 ч ди

дг + — а ъ

2 \ ди ( 2 ч дт - + а)-*: + (7-2 — а) — — 2аи,

дг

и ди

=---ь ,

г дг

ди дг'

ди дг

(3)

В соотношениях (1)-(3) использованы безразмерные величины (при одинаковом начертании они обозначены штрихами, которые в них и далее опущены):

и

т

и = I' т = I' г = I'

г с

= I' т = т

а « = 1+2^'

7 + £

11 =

722 =

с

2

а

а=

1 + 2^'

с=

1 + 2^ Р

2 _ М 2 _ 7 + £

С2 = Р' С3 = ~

/ Ш ф

Ш= 12' ф= 12'

П =

к =

7 — £

7 + £ 1 + 2^

в =

р12

3 '

= 1 — 27-2.

Здесь Ь — время; Ь — некоторый характерный линейный размер; 1, м — упругие постоянные Ламе; а, в, 7, £ — физические параметры моментной среды; р — ее

2

г

2

2

1

1

плотность; J — мера инерции среды при вращении (плотность момента инерции); С1, с2 и с3 — скорости волн растяжения-сжатия, сдвига и кручения соответственно. Соотношения между этими скоростями могут быть различными. Далее будем считать имеют место справедливые для композита в виде алюминиевой дроби в эпоксидной матрице [2] неравенства с2 < с1 < с3, из которых вытекает следующее соотношение:

72 < 1 < 71.

Отметим, что при а = 0 второе и третье уравнения в (1) становятся независимыми.

Считаем, что в начальный момент времени возмущения отсутствуют:

ф |т=0 = ф |т=0 = ф |т=0 = ф |т=0 = с |т=0 = О |т=0 = 0. (4)

Все искомые функции предполагаются ограниченными, а на граничной плоскости подобно классической задаче Лемба [14-18] выполняются такие условия:

^г |г=0 = 0, |г=0 = 5 (х, у) 5 (т), ^ \z=o = 0, (5)

где х, у — декартовы координаты в плоскости г = 0; 5 (£) — дельта-функция Дирака [15].

2. ИЗОБРАЖЕНИЯ ПЕРЕМЕЩЕНИЙ И УГЛА ПОВОРОТА

К начально-краевой задаче (1)-(5) применяем преобразования Лапласа по времени т и Ханкеля по радиусу г нулевого порядка для ф, т и первого порядка для ф, и, и (значки «Ь» и «Н» указывают на соответствующие изображения; 5 и д — параметры этих преобразований). Тогда уравнения (1), равенства (2), (3) и граничные условия (5) переходят в следующие соотношения:

^ - Р

дх-

ИЬ ,-

- к02 (д,5) фИЬ = 0 к-0 (9,5) = ^д2 + 7о52 (70 = 1), Ие^ > 0,

Т2 + а - [ (7Г2 + а) д2 + *2] ФИЬ + 2аиИЬ = 0, (6)

7-2 ¥ - - (7-°92 + 52 + 4ав) сИЬ + 2авд2фИЬ = 0;

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

иИЬ = -дфИЬ - ^, тИЬ = дфИЬ + дфИЬ, аИЬ = д|ИЬ + кдиИЬ,

дг дг ^^ дг

аИгЬ = (7Г° + а) диИЬ - - а) дтИЬ - 2асИЬ, ц*Ь = ^ИЬ

(7)

д^ ---- ,

1

¿=0

оИЬ\^ = 0. |г=0 = 2п мИЬ |^=0 =0 (8)

Общие решения уравнений в (6) с учетом их ограниченности имеют вид [13]:

2

фИЬ (д, г, 5) = С0 (д, 5) Е0 (д, г, 5), фИЬ (д, г, 5) = ^ С, (д, 5) Е, (д, г, 5),

-(„,*) = ^ Л _ .е^ . (9)

(д, г, 5) = — ^ Т (д, 5) С, (д, 5) Е1 (д, г, 5),

1=1

Ej (д, г, 5) = егк>, Т, (д, 5) = (7-2 + а) [д2 - к? (д. *)] +

498 Научный отдел

где Сз- (я, в) — постоянные интегрирования; а к 1,2 (я, в) — корни биквадратного уравнения, приведенного в [13].

Подстановка (9) в (7) и использование условий (8) позволяет определить постоянные интегрирования:

С (я в) 7х2к2 ^ в) К1 ^ в) С (я в) Т12Яко (я,в) Кз (я,в)

С0 (Я,в) —--о п /-ч-5 С 1 (Я, в) — —■

" в 2,, , (10)

2пя2 (я, в) ' у (я, в)

^ , \ 72яко (я, в) К (я, в)

С2 (я5 в) —--

пЯ2 (я, в) где

Я2 (я, в) — -4я2ко (я, в) КА (я, в) - к4 (я, в) К (я, в), К (я, в) — К (я, в) + Кз (я, в), К (я, в) — к1 (я, в) Т1 (я, в), Кз (я, в) — -к2 (я, в) Т2 (я, в), К4 (я, в) — -к1 (я, в) к2 (я, в) Я1 (я, в), Я1 (я, в) — Т1 (я, в) - Т2 (я, в).

Поскольку вид корней к1)2 (я, в) не позволяет находить оригиналы аналитически, то аналогично [13] используем разложения в степенные ряды по малому параметру а, а именно, ограничиваясь линейным приближением, перемещения и угол поворота записываем так:

и (г, г, т) — и0 (г, г, т) + и1 (г, г, т) а, т (г, г, т) — т0 (г, г, т) + т1 (г, г, т) а, (ц) ш (г, г, т) — ш0 (г, г, т) + ш1 (г, г, т) а.

Отметим, что первые члены равенств (11) есть решения классической задачи Лемба для упругой среды. При этом ш0 = 0. Функции же и0 и т0 найдены, например, в [18]. Поэтому далее они не рассматриваются.

Соответствующие (11) равенства для к1)2 (я, в) получены в [13]:

кг (я, в) — к01 (я, в) + акц (я, в) (1 — 1,2),

4 2 22

/--Т4 в2 2^7'

к01 (я, в) — у я2 + Т2в2, кц (я, в) — -----, кц (я, в) —

2к01 (я, в)' ' к02 (я, в)' Т (я, в) — Тц (я, в) а2, Т2 (я, в) — Т20 (я, в) + аТц (я, в), (12)

Т12 (я, в) — -^^, Т20 (я, в) — 7Г2 (7? - 722) в2, 72 - 7! )

Тц (я, в) — -7Г2722 (40 + 7^) .

Использование (7), (9), (10) и (12) позволяет найти изображения искомых коэффициентов в (11). Например, для нормального перемещения они имеют вид

2

т^ (я,г,в) —£ т^ (я, г, в), (13)

з=0

где

тНЬ (я г в) — 76 в2я2к2 (я ,в) к22 (я ,в) егк0(*,-)* тпь ( в) — 0 т10 (я, ^ в) — пк01 (я, в) Я2 (я, в) е , т12 (я, ^ в) — 0,

^ (я,г,в) — - 76в2я2к0 (я,в) Г4я2к0 (я,в) - ^ егк0!

2пк01 (я, в) Я (я, в)

Я(я, в)

Я (я, в) — к34 (я, в) - 4я2к0 (я, в) к01 (я, в), к32 (я, в) — 2я2 + 72в2.

3. ОПРЕДЕЛЕНИЕ ОРИГИНАЛОВ

Далее рассматриваем нормальное перемещение на поверхности полупространства. Соответствующие изображения согласно (13) имеют вид

^ (9,«) = (9,0, в) + (д, 0, в) = (9, в). (14)

Здесь — изображение Лапласа по времени т и Фурье по декартовой координате х (им соответствует верхний индекс Р) соответствующей составляющей решения плоской задачи Лемба для моментно-упругой полуплоскости. Оно получено в [12] и записывается так:

^ (9,в) = /У . (15)

ко 1 (9, в) Я2 (9, в)

Его оригинал также найден в указанной работе (Н (т) — единичная функция Хевисайда):

(16)

(х,т) = £ (х,т) Н (т - 7з |х|),

0 (х,т) = 1^^^, ^1 (х,т) = ^^^, П ут2 - х2 П л/т2 - 72х2

Здесь

Я« (х2, т2) = (т2 ^2), Ро (х2, т2) = (2т2 - 712х2)2 (т2 - х2) ,

Р1 (х2, т2) = (7?х2 - 2т2)4 + 16т4 (т2 - х2) (т2 - 7?х2) , а Р3 (х,т) — многочлен, имеющий следующий вид [15,18]:

Рз (х, т) = 76х3 - 874х2т + 87? (2 + к) хт2 - 8(1 + к) т3 = 716 (х - 4т) Р2 (х, т),

Рз (4, 1) =0, Р2 (х,т)= х2 - 2а2хт + Ь2т2, а2 = 42 - Ся, Ь2 = 161

72 2 71ся

(ся — скорость волн Релея).

При этом для дискриминанта (к) многочлена Р2 имеют место следующие соотношения:

ся < 7Г 1, ^2 (к) > 0 (0 ^ к < к, « 0, 357), В2 (к,) = 0, В2 (к) < 0 (к, < к < 1). Далее для определенности ограничимся случаем к, < к.

Равенства (14) позволяют найти оригинал функции СНЦ с помощью утверждений о связи преобразований Фурье и Ханкеля [15,18]. При этом учитываем, что для этой функции использовалось преобразование Ханкеля нулевого порядка. Кроме того, принимаем во внимание, что функция ^ (х, т) — четная функция по х. В результате получаем (интегралы понимаются в смысле регуляризованных значений):

сю 1

1 /• х^1 (х, т) ¿ж 72 т , ,

»1 (г т) = -П I (х2 - Г2)3/2 = - V ^ ^ (Г, т) Н (т - ^Г),

(17)

1

Y-2

1у (г, т) = ^ [ ^) ^ , (*) = УЛ,

П г2(1 - Ц*) (т2* - г2) т2* - г2

Аю = 8, Лп = 1.

С помощью разложения рациональной части подынтегральной функции на элементарные дроби в последнем интеграле выделяем регулярную и сингулярную части:

'у (г-т) = ^ [К* (г,т) + К* (г,т)], (18)

Yi

где

Y-2

1 Г ф,- (z,r2/т2)

K,,. (г,т)- — j /т ) dz,

ПТ /т^(1 - YjJz) (z - Г2/т2) ф, (z, £) = ГF (z) - F (£) F (z) - F (cR) 1 1 - £

— £ z - c:

R

(£ - cR )2

1 - cR F (z) - Fj (cR) - F'j (cR) (z - cR)

£ - cR (z - cR)2

_ (z, 1) ,

F (z) - wm);

Kjs (г, т) - hji (т) (г2 - т2cR) ;3/2 + hj2 (г, т) (г2 - т2cR) +5/2, hji (т) - - (l1 - cR2 )з/2 ^ Fj (cR) +2 (1 - Y2cR) F'j (cR)] ,

hj2 (г,т)- -т(3 - cR) т2-2г2 Fj (cR) .

V1 - Y'cR

Достаточно просто доказывается, что функция Ф, (z, £) доопределяется непрерывным образом:

ф, (z, cR) - - Fj (z) - Fj (4) -Fj(cR)(z - cR), Ф, (cR, cR) - - 2Fj' (cR).

(z - cR) 2

Выясним, как ведет себя функция G0wi (г, т) на границах носителя ее составляющих. Прежде всего, отметим, что носителями слагаемых в (17) и функций K,s (г, т) являются множества т ^ y,г и г ^ cRт соответственно. При этом функции K,r (г, т) в (18) непрерывные и для них справедливы равенства

lim (г, т) - lim (г, т) - 0. (19)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

т ^Yj r+0 т

Функция Kjs (г, т) при т - y,г принимает конечное значение:

Ks (г, Yj г) - - П22 {(1 - cR) (1 - y,2 cR) F'j (cR) +

+F (4)[ (2 - 4) y? -1}, (20)

а на фронте волны Релея коэффициенты hji и hj2 при ее сингулярных составляющих определяются так:

hj 1 (г/сл) = - - CR2};аз/2 [y?Fj (cRR) +2 (1 - 72cR) Fj (cR)] ,

2c41 - Yj2 4

3(1 _ c2 ) r3 (21)

hj 2 (r,r/CR ) = - 3(1 . Cg) r Fj (4).

24 у1 - T?4

Окончательно с учетом неравенств 1 < y1 < c-1 из (17)-(21) получаем:

Gowi (r, t) = 0 (т < r), lim Gowi (r, т) = -^Kos (r, r), lim Gowi (r, т) = 0,

T ^r+0 nYl T

lim Gowi (r, т) =--10 [Kor (r, Yir) + Kos (r, Yir)] ,

t^yi r-o nYl

lim Gowi (r, т) = lim Gowi (r, т)--^Kb (r, Yir),

t^yi r+o t^yi r-o nYi

lim Gowi (r, т) =--8 [AioKor (r, r/cß) + AnKir (r, r/cß)] .

T^r/cR+o nYl

Следовательно, график этой функции при фиксированных г имеет в точках т = г и т = 71г разрывы первого рода, а в точке т = г/ся левую вертикальную асимптоту.

Отметим, что решение осесимметричной задачи Лемба для упругой среды имеет степенную особенность порядка -3/2 на фронте волны Релея г = сят. Для мо-ментной же среды, как вытекает из формул (17) и (18), появляется дополнительная особенность порядка -5/2. Кроме того, обращает на себя внимание тот факт, что рассмотренное первое приближение для нормального перемещения не содержит фронта волны кручения.

4. ПРИМЕРЫ РАСЧЕТОВ

Полагаем, что материалом, заполняющим полупространство, является зернистый композит из алюминиевой дроби в эпоксидной матрице со следующими физическими характеристиками [2]:

Л = 7.59 ГПа; ^ = 1.89 ГПа; а = 7.45 МПа; 7 + е = 2.64 кН; 3 = 0.429 • 10-3кг/м.

В качестве характерного линейного размера принимаем Ь = 1м. При этом безразмерные параметры таковы:

71 =2,45; 72 = 0, 92; а = 0, 66 • 10-3; в = 5,1 • 106; сЯ = 0,148.

На рисунке представлены зависимости функций влияния от времени т для трех радиусов: г = 0.2, г = 0.4 и г = 0.6.

G^xio"

10 О -10

.............rr............. .............J.L............ r:............ i

II I T

\o. ,/^X1. I.---- 2

\ 7 Г - - I / 1 1 \ i .........)...:

/.........П .........1. л 1 .........!..;

r r 1 1 ' i

1 1 ...........i. .1. i ..........1 ■

Зависимости функции Gowi от времени т: сплошная кривая соответствует r = 0.2, пунктирная — r = 0.4, а штрихпунктирная — r = 0.6 The function Gowi depends on time т: the solid curve corresponds r = 0.2, dashed — r = 0.4, dot-dash — r = 0.6

Благодарности. Работа выполнена при финансовой поддержке РФФИ (проект № 18-08-00471).

Библиографический список

1. Cosserat E., Cosserat F. Theorie des corps deformables. P. i Librairie Scientifique A. Hermann et Fils, 1909. 226 p.

2. Ерофеев В. И. Волновые процессы в твердых телах с микроструктурой. М. i Изд-во МГУ, 1999. 328 c.

3. Кулеш М. А., Шардаков И. Н. Построение и анализ некоторых точных аналитических решений двумерных упругих задач в рамках континуума Коссера // Вестн. ПГТУ. Математическое моделирование систем и процессов. 2001. № 9. С. 187-201.

4. Лай Тхань Туан, Тарлаковский Д. В. Распространение нестационарных кинематических возмущений от сферической полости в псевдоконтинууме Коссера // Механика композиционных материалов и конструкций. 2011. T. 17, № 2. С. 184-195.

5. Лай Тхань Туан, Тарлаковский Д. В. Распространение нестационарных осе-симметричных возмущений от поверхности шара, заполненного псевдоупругой средой Коссера // Электронный журн. «Труды МАИ». 2012. № 53. URLi httpi//trudymai.ru/published.php?ID=29267/ (дата обращения! 19.04.2012).

6. Лай Тхань Туан, Тарлаковский Д. В. Дифракция нестационарных волн на сферической полости в псевдоконтинууме Коссера // Радиоэлектроника. Наносистемы. Информационные технологии. 2013. T. 5, № 1. С. 119-125.

7. Пальмов В. А. Основные уравнения теории несимметричной упругости // ПММ. 1964. Т. 28, вып. 6. С. 1117-1120.

8. Белоносов С. М. Моментная теория упругости: (Статика). Владивосток i Дальнаука, 1993. 148 с.

9. Бытев В. О., Слезко И. В. Решение задач асимметричной упругости // Математическое и информационное моделирование i сб. науч. тр. Тюмень i Вектор Бук, 2008. Вып. 10. С. 27-32.

10. Атоян А. А., Саркисян С. О. Динамическая теория микрополярных упругих тонких пластин // Экол. вестн. науч. центров ЧЭС. 2004. № 1. С. 18-29.

11. Hirdeshwar S. Saxena, Ranjit S. Dhaliwal. Eigenvalue approach to axially-symmetric coupled micropolar thermoelasticity // Bull. Pol. Acad. Sci. Techn. Sci. 1990. Vol. 38, № 1. P. 7-18.

12. Суворов Е. М, Тарлаковский Д. В., Федотенков Г. В. Плоская задача об ударе твердого тела по полупространству, моделируемому средой Коссера // ПММ. 2012. Т. 76, Вып. 5. С. 850-859.

13. Чан Ле Тхай, Тарлаковский Д. В. Нестационарное осесимметричное движение упругого моментного полупространства под действием нестационарных нормальных поверхностных перемещений // Учен. зап. Казан. ун-та. Сер. Физ.-матем. науки. 2017. Т. 159, кн. 2. С. 231-245.

14. Новацкий В. Теория упругости. М. : Мир, 1975. 872 с.

15. Горшков А. Г., Медведский А. Л., Рабинский Л. Н, Тарлаковский Д. В. Волны в сплошных средах. М. : Физматлит, 2004. 472 с.

16. Слепян Л. И. Нестационарные упругие волны. Л. : Судостроение, 1972. 351 с.

17. Поручиков В. Б. Методы динамической теории упругости. М. : Наука, 1986. 328 с.

18. Горшков А. Г., Тарлаковский Д. В. Динамические контактные задачи с подвижными границами. М. : Физматлит, 1995. 352 с.

Образец для цитирования:

Чан Ле Тхай, Тарлаковский Д. В. Осесимметричная задача Лемба для среды Коссера // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2018. Т. 18, вып. 4. С. 496-506. DOI: https://doi.org/10.18500/1816-9791-2018-18-4-496-506

Axisymmetric Problem Lemba for the Cosserat Medium

Tran Le Thai, D. V. Tarlakovskii

Tran Le Thai, https://orcid.org/0000-0002-8733-6150, Moscow Aviation Institute (National Research University), 4 Volokolamskoye Shosse, Moscow 125993, Russia, [email protected] Dmitrii V. Tarlakovskii, https://orcid.org/0000-0002-5694-9253, Institute of Mechanics Lomonosov Moscow University, 1 Michurinsky Ave., Moscow 119192, Russia; Moscow Aviation Institute (National Research University), 4 Volokolamskoye Shosse, Moscow 125993, Russia, [email protected]

The article deals with elastic homogeneous isotropic half-space filled with the Cosserat medium. At the initial instant of time and at infinity, there are no perturbations. At the boundary of the half-space, normal pressures are given. All the components of the stress-strain state are supposed to be limited. A cylindrical coordinate system is used with an axis directed inward into the half-space. With allowance for axial symmetry, the resolving system of equations includes three hyperbolic equations with respect to the scalar potential and the non-zero components of the vector potential and the rotation vector. The components of displacement vectors, rotation angle, stress tensors and stress moments are related to the potentials by known relationships. The solution of the problem is sought in the form of generalized convolutions of a given pressure with the corresponding surface influence functions. To construct the latter, Hankel transformations along the radius and Laplace transformations are applied in time. We use the expansion in power series for a small parameter characterizing the connection between the shear and rotation waves. The images of the first two coefficients of these series are found. The corresponding originals are determined by the connection between the plane and axisymmetric problems. Examples of calculations of the regular components of the influence of a granular composite from an aluminum shot in an epoxy matrix are given.

Key words: Cosserat medium, surface influence functions, Laplace and Hankel integral transformations, the small parameter method, connection between plane and axisymmetric problems.

Acknowledgements: This work was supported by the Russian Foundation for Basic Research (project no. 18-08-00471).

References

1. Cosserat E., Cosserat F. Theorie des corps deformables. Paris, Librairie Scientifique A. Hermann et Fils, 1909. 226 p. (in French).

2. Erofeev V. I. Volnovye processy v tverdyh telah s mikrostrukturoj [Wave Processes in Solids with Microstructure]. Moscow, Moscow Univ. Press, 1999. 328 p. (in Russian).

3. Kulesh M. A., Shardakov I. N. Construction and analysis of some exact analytic solutions of two-dimensional elastic problems within the Cosserat continuum. Vestn. PGTU. Matem-aticheskoe modelirovanie sistem i processov [Bulletin of Perm State Technical University. Mathematical modeling of systems and processes], 2001, no. 9, pp. 187-201 (in Russian).

4. Lai Thanh Tuan, Tarlakovsky D. V. Propagation of non-stationary kinematic perturbations from a spherical cavity in the Cosserat pseudocontinuum. Mekh. Kompoz. Mater. Konstr. [Mechanics of composite materials and structures], 2011, vol. 17, no. 2, pp. 184-195 (in Russian).

5. Lai Thanh Tuan, Tarlakovsky D. V. Propagation of nonstationary axisymmetric perturbations from the surface of a sphere filled with a Cousser pseudoelastic medium. Online journal "Trudy MAI", 2012, no. 53. (in Russian). Available at: http://trudymai.ru/published.php?ID=29267/ (accessed 19 April 2012).

6. Lai Thanh Tuan, Tarlakovsky D. V. Diffraction of waves by a spherical cavity in the Cosserat pseudo-continuum. Radioelectron., Nanosist., Inf. Tekhnol., 2013, vol. 5, no. 1, pp. 119-125 (in Russian).

7. Palmov V. A. Basic equations of the theory of asymmetric elasticity. Applied Mathematics and Mechanics, 1964. vol. 28, iss. 6. pp. 1117-1120 (in Russian).

8. Belonosov S. M. Momentnaya teoriya uprugosti: (Statika) [Moment theory of elasticity: (Statics)]. Vladivostok, Dal'nauka Publ., 1993. 148 p. (in Russian).

9. Bytev V. O., Slezko I. V. Solution of asymmetric elasticity problems. Mathematical and Informational Modeling : Collection of scientific papers. Tyumen : Vector Beech, 2008, iss. 10, pp. 27-32 (in Russian).

10. Atoyan A. A., Sarkisyan S. O. Dynamic theory of micropolar elastic thin plates. Ekol. vestn. nauch. tsentrov CHES, 2004, no. 1, pp. 18-29 (in Russian).

11. Hirdeshwar S. Saxena, Ranjit S. Dhaliwal. Eigenvalue approach to axially-symmetric coupled micropolar thermoelasticity. Bull. Pol. Acad. Sci. Techn. Sci., 1990, vol. 38, no. 1, pp. 7-18.

12. Suvorov E. M., Tarlakovskii D. V., Fedotenkov G. V. The plane problem of the impact of a rigid body on a half-space modelled by a Cosserat medium. J. Appl. Math. Mech., 2012, vol. 76, no. 5, pp. 511-518. DOI: https://doi.org/10.1016/jjappmathmech.2012.11.015

13. Tran Le Thai, Tarlakovskii D. V. Nonstationary axisymmetric motion of an elastic momentum semi-space under non-stationary normal surface movements. Uchenye Zapiski Kazanskogo Universiteta. Seriya Fiziko-Matematicheskie Nauki, 2017, vol. 159, no. 2, pp. 231-245 (in Russian).

14. Novatsky V. Teoriya uprugosti [Theory of Elasticity]. Moscow, Mir, 1975. 872 p. (in Russian).

15. Gorshkov A. G., Medvedskii A. L., Rabinskii L. N., Tarlakovskii D. V. Volny v sploshnyh sredah [Waves in Continuous Media]. Moscow, Fizmatlit, 2004. 472 p. (in Russian).

16. Slepian L. I. Nestacionarnye uprugie volny [Non-stationary elastic waves]. Leningrad, Sudostroenie Publ., 1972. 351 p. (in Russian).

17. Poruchikov V. B. Metody dinamicheskoj teorii uprugosti [Methods of the Dynamic Theory of Elasticity]. Moscow, Nauka, 1986. 328 p. (in Russian).

18. Gorshkov A. G, Tarlakovskii D. V. Dinamicheskie kontaktnye zadachi s podvizhnymi granicami [Dynamic contact problems with moving boundaries]. Moscow, Fizmatlit, 1995. 352 p. (in Russian).

Cite this article as:

Tran Le Thai, Tarlakovskii D. V. Axisymmetric Problem Lemba for the Cosserat Medium. Izv. Saratov Univ. (N.S.), Ser. Math. Mech. Inform., 2018, vol. 18, iss. 4, pp. 496-506 (in Russian). DOI: https://doi.org/10.18500/1816-9791-2018-18-4-496-506

i Надоели баннеры? Вы всегда можете отключить рекламу.