Thermal analysis of a rotating micropolar medium using a two-temperature micropolar thermoelastic model with higher-order time derivatives Текст научной статьи по специальности «Физика»

УДК 539.3

Термический анализ микрополярной среды при вращении на основе двухтемпературной микрополярной модели термоупругости с высшими производными по времени

A.E. Abouelregal1'2, R. Alanazi1, A.H. Sofiyev3'4'5, H.M. Sedighi6

1 Университет Джуфа, Аль-Курайят, 75911, Саудовская Аравия

2 Университет Мансура, Мансура, 35516, Египет 3 Университет им. Сулеймана Демиреля, Испарта, 32260, Турция 4 Стамбульский коммерческий университет, Стамбул, 34445, Турция 5 Азербайджанский государственный экономический университет, Баку, 1001, Азербайджан 6 Университет Шахида Чамрана в Ахвазе, Ахваз, 61357-43337, Иран

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

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

DOI 10.55652/1683-805X_2023_26_1_15

Thermal analysis of a rotating micropolar medium using a two-temperature micropolar thermoelastic model with higher-order time derivatives

A.E. Abouelregal1,2, R. Alanazi3, A.H. Sofiyev4,5,6, and H.M. Sedighi7,8

1 Department of Mathematics, College of Science and Arts, Jouf University, Al-Qurayat, 75911, Saudi Arabia 2 Department of Mathematics, Faculty of Science, Mansoura University, Mansoura, 35516, Egypt

3 Department of Computer Science, College of Science and Arts, Jouf University, Al-Qurayat, 75911, Saudi Arabia 4 Department of Civil Engineering, Engineering Faculty, Suleyman Demirel University, Isparta, 32260, Turkey 5 Information Technology Research and Application Center, Istanbul Commerce University, Beyoglu, Istanbul, 34445, Turkey

6 Scientific Research Center for Composition Materials, Azerbaijan State Economic University, Baku, 1001, Azerbaijan 7 Faculty of Engineering, Mechanical Engineering Department, Shahid Chamran University of Ahvaz, Ahvaz, 61357-43337, Iran 8 Drilling Center of Excellence and Research Center, Shahid Chamran University of Ahvaz, Ahvaz, 61357-43337, Iran

In this work, the propagation of planar waves in a homogeneous micropolar thermoelastic medium is studied while the entire body rotates with a uniform angular speed. The coordinate system of the rotating medium is assumed to be stationary, and therefore the kinematic equations have two additional terms, namely, the gravitational and the Coriolis accelerations. The problem is addressed based on the two-temperature ther-moelastic model with higher-order time derivatives and dual-phase lag, which can explain the effect of microscopic features in nonsimple materials. With certain boundary conditions and the normal mode approach, the variations in temperature, displacement, microrotation, and thermal stresses induced by heating are derived. In the absence of rotation and two-temperature factor, comparison is made with the results of classical thermoelastic models.

Keywords: micropolar, two-temperature thermoelasticity, rotating microstructures, phase lag, higherorder time derivatives

© Abouelregal A.E., Alanazi R., Sofiyev A.H., Sedighi H.M., 2023

1. Введение

Распространение термомеханических волн в упругих средах имеет теоретическое и практическое значение во многих научных областях. В последние десятилетия возрос интерес к изучению поведения волн в твердых телах при различных условиях, например при вращении [1]. Например, в [2] теоретически и экспериментально исследовано распространение волн от источника тепла в изотропной упругой среде, вращающейся с постоянной скоростью. Для получения полной информации о скорости распространения волн в упругой среде необходимо изучение плоских волн. Сложной проблемой является описание распространения волн, вызванных локальным возмущением, и уменьшения амплитуды волн по мере удаления от источника возмущения. В настоящем исследовании решается задача о распространении поверхностной волны во вращающейся упругой среде. Следует отметить, что поверхностные волны возникают в медленно вращающихся средах, где скорость поверхностных волн выше. Например, волны Рэлея в среде ведут себя подобно сейсмическим волнам, распространяющимся по поверхности Земли [3]. В последние годы распространение упругих, а также тепловых волн во вращающейся среде активно изучается на основе разных моделей с учетом различных условий [410].

В 1837 г. Ж.-М. Дюамель разработал связанную теорию термоупругости с учетом влияния тепловых полей на деформацию. Ж.-Б. Био [11] ввел связанную теорию термоупругости с основными уравнениями, полученными из закона Фурье, т.е. на основе термодинамики необратимых процессов. В соответствии с данной концепцией, система определяющих уравнений (уравнение теплопроводности) имеет параболический тип. Другими словами, тепловые волны распространяются в среде с бесконечными скоростями, что противоречит природе физических явлений.

Рядом авторов были предложены модели термоупругости, которые позволяют получить гиперболические уравнения теплопроводности с конечной скоростью тепловых волн. Это обобщенная теория Лорда-Шульмана [12], модель Грина-Линдсея, в которой определяющие соотношения между напряжениями и энтропией обобщаются за счет учета двух альтернативных коэффициентов релаксации [13], модели Грина-Нагди I, II и III типа, которые позволяют применять теорию термоупругости для однородных материалов [14-

16]. Дальнейшим расширением теории термоупругости является модель с двухфазным запаздыванием Цзоу [17-19]. За счет включения двухфазного запаздывания в вектор теплового потока и в градиент температуры было получено основное уравнение для описания замедленного тепло- и массопереноса в различных материалах, а также особенностей микроструктуры, таких как элект-рон-фононное взаимодействие и рассеяние фоно-нов [17-19].

В [20-24] разработаны математические модели, описывающие процесс теплообмена в упругих телах с использованием высших производных по времени и фазового запаздывания. Эти модели являются расширением механических основ теорий теплообмена Грина-Нагди [15, 16], модели с трехфазным запаздыванием Чоудхури [25], а также моделей Цзоу [17-19]. Разработанные модели высоких порядков с фазовым запаздыванием обеспечивают теоретическую основу для изучения процессов теплопереноса с учетом различных особенностей микроструктуры, позволяя прогнозировать тепловое поведение конструкций на основе многомасштабной модели.

В работах [26-28] представлена двухтемпера-турная теория термоупругости. В данной модели неравенство Клаузиуса-Дюгема преобразовано путем введения двух разных температур: связанной с проводимостью и термодинамической. Температура проводимости обусловлена тепловым процессом, а термодинамическая температура связана с диссипацией механической энергии при взаимодействии частиц и упругой матрицы. Считается, что разность данных температур пропорциональна количеству подведенного тепла, а в отсутствие подвода тепла они равны. Разница между температурами отсутствует в простых веществах и, наоборот, присутствует в сложных материалах. В [29] обсуждается существование, структурная устойчивость, конвергенция и пространственное поведение решений двухтемпературной модели термоупругости. Модификациями рассматриваемой теории являются обобщенная теория термоупругости с одним коэффициентом релаксации [30], двухтемпературная модель термоупругости с производными по времени высших порядков и тремя различными фазовыми запаздываниями [20]. В [31] на основе метода пространства состояний получена модель одномерных уравнений для обобщенной двухтемпературной теории магнитотермоупругости в идеальной электропроводной среде с двумя временами релакса-

ции. В [32] представлена обобщенная теория термоупругости с двумя температурами и двухфазным запаздыванием. В последние десятилетия активно исследуются различные аспекты двухтем-пературных теорий термоупругости [33-36].

Микрополярная теория упругости, также известная как теория упругости Коссера или микрополярная теория механики сплошной среды, учитывает, что при вращении локальной точки, наряду с перемещением согласно классической теории упругости, на единице площади действуют силовые и моментные напряжения. В общепринятой теории упругости под напряжением подразумевается силовое напряжение ввиду отсутствия других его видов. Концепция моментного напряжения восходит к ранней работе Фойгта по теории упругости. Недавно разработанные теории моментных напряжений используют весь спектр возможностей механики сплошной среды. Поскольку значение коэффициента интенсивности напряжений снижается вблизи отверстий и трещин, обобщенные концепции сплошной среды, такие как теория упругости Коссера, применимы для описания характеристик материала [37].

В последние годы микрополярная теория упругости Эрингена [38] привлекает большое внимание ученых, поскольку в отличие от классической теории позволяет подробно описывать деформационное поведение твердых тел. Микрополярная теория особенно эффективна при изучении материалов, состоящих из вытянутых молекул, испытывающих микроповорот под действием объемных и поверхностных моментов, которые вызваны действием вектора силы на поверхности тела. Микрополярный континуум представляет собой совокупность связанных частиц, которые принимают форму твердых тел малого размера и способны вращаться, наряду с обычным поступательным движением [37].

Линейное описание микрополярной термоупругости было получено при учете влияния температуры в концепции микрополярного континуума. Микрополярная теория термоупругости основана на включении термоупругих свойств в понятие микрополярности [39, 40]. Основные уравнения математической модели микрополярной термоупругости представлены в работе [41].

С использованием модели Грина-Линдсея были построены уравнения обобщенной микрополярной теории термоупругости [42]. Уравнение энергетического баланса и теорема единственнос-

ти для анизотропных материалов получены на основе микрополярной модели термоупругости, в которой определяющие переменные зависят от характеристик теплового потока [43]. Основы трехфазной микрополярной теории термоупругости сформулированы в работе [44]. Авторы установили вариационный принцип для линейного анизотропного и неоднородного термоупругого тела, доказав теоремы взаимности и единственности. В [45] представлена математическая модель магнитного микрополярного термоупругого полупространства с температурно-зависимыми свойствами материала под действием механической полосовой нагрузки. Ряд исследований посвящен детальной разработке микрополярной теории термоупругости [46-53].

В настоящей работе обсуждается модифицированная микрополярная модель термоупругости с высшими производными по времени и двухфазным запаздыванием. Наряду с учетом влияния микроструктуры в процессе теплообмена, предлагаемая модель учитывает также влияние макроструктуры в предположении, что электрон-фо-нонные отклики приводят к задержке роста температуры кристаллической решетки на макро-масштабном уровне. При дальнейшем построении модели осуществляется разложение в ряд Тейлора выражения закона Фурье и устанавливается корреляция между двумя температурами при сохранении условий фазового запаздывания вплоть до допустимых членов более высокого порядка.

Предлагаемая модель с фазовым запаздыванием является расширением двухтемпературной модели термоупругости с одним временем релаксации [30] и двухфазным запаздыванием [32]. В предыдущих работах установлено, что двухтем-пературные модели термоупругости более применимы к решению практических задач физики. По мере накопления практических и теоретических знаний обобщенные теории термоупругости с высшими производными по времени смогут найти применение в современных приложениях.

На основе предложенной модели изучено распространение тепловой волны в изотропном микрополярном полупространстве, вращающемся с постоянной угловой скоростью, внешняя поверхность которого является свободной. Получены выражения для термодинамической температуры, температуры проводимости, микроповорота, перемещений и термических напряжений, а также исследованы распределения физических переменных.

2. Математическая модель и основные уравнения

Введем уравнения, описывающие бесконечную микрополярную среду в рамках двухтемпературной модели теплопроводности с фазовым запаздыванием и высшими производными по времени. В качестве определяющих уравнений используем основные размерно-зависимые уравнения Эрингена [40] в отсутствие объемных сил и моментов, источников тепла и уравновешенной внешней силы [54-56]:

(ц + а)вгу ji -y9ôj, (1)

h(u + P)<j +(и-р)шг,, (2)

сту = askk5у -

m=s®k,ksy

=—(u

2

U j).

(3)

Здесь ац — компоненты тензора напряжений; и, — компоненты вектора перемещений; югк — компоненты вектора микроповорота; тц — компоненты тензора моментных напряжений; X, ц, а, и, в и в — константы материала; 5ц — дельта-функция Кронекера; вц — компоненты тензора малых деформаций; юг- — компоненты поворота; 0 = Т - Т0 — изменение температуры выше начальной температуры Т0; у — постоянная материала, определяемая выражением у = (3Х + 2ц + а)ах.

Уравнения движения и баланса количества движения имеют вид

°ц,, =Ри,, (4)

SijkСjk ■

+ mji, j = P/®,i, i, j, k = —, 2,3,

(5)

где вук — символ перестановки; р — плотность; ц — микроинерция.

Используя определяющие уравнения (1), (2),

исключим напряжения о., движения (4), (5) и получим

mi,

из уравнений

(^ + a)V2 -P—

dt2

(А, + ц, - а)V(V • u) + 2aVx ю = yV9, (6)

(U + P)V2 -pj

dt2

ю

+ (s + u-P)V(V- ю) + 2aVx u = 4аю = 0. Уравнение энергии запишем в виде

PCE ^ + УГо ^Т (divu) = -V • q + pÖ,

dt dt

(7)

(8)

где СЕ — удельная теплоемкость; q — вектор теплового потока. Классический закон Фурье выразим как

q(x, х) ^куе^, х), (9)

где x — радиус-вектор, а K — коэффициент теплопроводности .

Сложное вещество, согласно работам [27, 28], представляет собой материал, в котором напряжение, энергия, энтропия, тепловой поток и термодинамическая температура в конкретный момент времени зависят от величин градиента деформации, температуры проводимости и ее градиента до данного момента времени. В [57] обычный закон Фурье (9) был преобразван в соответствии с двухтемпературной концепцией к виду

q(x, t) = - kv9(x, t ), (10)

где K — теплопроводность; ф — температура проводимости, измеренная при постоянной начальной температуре ф0, которая удовлетворяет соотношению

ф-0 = aV29, (11)

где a > 0 — двухтемпературный коэффициент.

В работах [32, 58] предложена модель теплопроводности с двухфазным запаздыванием, объединенная с двухтемпературной теорией за счет учета микроструктурных эффектов в процессе теплопередачи:

q(x, t + т9) = -KVq>(x, t + Тф), (12)

где Tq и тф — фазовые запаздывания теплового потока и градиента температуры проводимости соответственно.

Выполним разложение в ряд Тейлора для обеих частей уравнения (12) отдельно по фазовым запаздываниям Tq и тф до членов достаточно высоких порядков m и n соответственно [20, 21]:

( ,,, -Г \ i _r ^r Л

m T d' Г=1 r ! dtr

q = -K

1

n T d' rT1 r ! dtr

Vф. (13)

Объединив (13) с уравнением энергии (8), получим модифицированное уравнение для теплопроводности высокого порядка с двумя температурами:

( ~ -Г \ ( ,,, _Г ~г

K

1

V 2ф =

1

h h. _dl

h r ! dtr

n T dr r=1 r ! dtr

PCE ^ + yT0 dt(div u)-PÖ |.

(14)

Ограничения на термомеханические параметры изотропного объекта удовлетворяют следующим неравенствам:

3А, + 2ц> 0, ц> 0,

ц + а > 0, у + в>0, а, у, 8>0.

Авторы [59-61] выявили определенные ограничения на использование более высоких порядков т и п. Например, при т > 5 получается неустойчивая система, которая не может описывать реальные физические свойства. В любом случае, когда порядок аппроксимации меньше или равен четырем, необходимо проверять соответствие второму закону термодинамики.

Будем считать, что однородная микрополярная среда равномерно вращается с угловой скоростью П = Пп, где вектор п указывает направление оси вращения. Под действием вращения возникают центробежное ускорение П х (П х и) и ускорение Кориолиса 2П х и . Тогда уравнения динамики (6) и баланса количества движения (7) в условиях микрополярности и вращения принимают вид

2 д2 (ц + а)У2 -р—

дг2

и + (Я + ц-а)У(У-и)

-2аУхю = уУ9 + рПх(Пхи) + 2рПхи, (16)

2 д2 (и + Р)У2 -рр

дг

ю + (в + и-р)У(У-ю) 2аУ х и = 4аю + ррП х ю.

(17)

2.1. Частные случаи

На основе полученной системы уравнений можно построить различные одно- и двухтемпе-ратурные модели термоупругости с микрополярным эффектом и без него.

2.1.1. Случай 1

При отсутствии микрополярного эффекта (е = и = в = 0) могут быть получены следующие варианты одно- и двухтемпературных моделей термоупругости:

- классическая модель термоупругости (СТЕ) [11] при гд = Тф = 0, a = 0, 0 = ф;

- обобщенная модель термоупругости Лорда-Шульмана (Ь8) [12] при Тд = т0 > 0, a = 0, 0 = ф, Тд, тф ^ 0 и т = 1;

- модель термоупругости с двухфазным запаздыванием (БРЬ) [18, 19] при a = 0, тя, тф > 0, п = 1 и т = 2;

- двухтемпературная модель термоупругости с двухфазным запаздыванием (2БРЬ) [32] при a > 0, Т0 = Тф, т = п = 1;

- обобщенная двухтемпературная модель с одним временем релаксации (2Ь8) [30] при a > 0, Тф = 0, Тд = Т0 > 0, т = п = 1;

- обобщенная двухтемпературная модель термоупругости с двухфазным запаздыванием с производными высокого порядка (2НБРЬ) [20] при

a > 0, Тд, Тф > 0, п, т > 1.

2.1.2. Случай 2

С учетом микрополярности можно получить такие же одно- или двухтемпературные модели термоупругости, как и в предыдущем случае. Далее такие модели будем обозначать добавлением буквы М в название модели (2МЬ8, 2МБРЬ, 2МНБРЬ, МТТЬ8 и МНТТЕ).

3. Постановка задачи

Рассмотрим полупространство из однородного микрополярного материала. Изначально среда не деформируется и не сжимается, находясь при постоянной температуре Т0. Граница среды при у = 0 свободна от усилий и подвергается воздействию источника тепла, мощность которого уменьшается со временем и который влияет на узкую область вокруг оси х. Задача решается в декартовой системе координат (х, у, z), начало которой находится на верхней поверхности у = 0, а ось у направлена внутрь полупространства. Распространение плоской волны в плоскости задается одинаковым смещением частиц вдоль линии, параллельной оси z. В результате все рассматриваемые поля являются только функциями переменных х, у и г и не зависят от координаты z. Исходя из этих предположений, запишем выражения для компонент вектора перемещений и микроповорота: их = (х, у, г), Ыу = у(х, у, г), и2 = 0,

ю = (0, ю, 0).

Тогда, объемное расширение e в плоскости ху описывается уравнением

„г ды ду

e = Чк =т- + Х'. (19)

дх ду

Предполагается, что среда вращается с равномерной угловой скоростью П = (0, 0, П). Соответственно, уравнения движения (16) и баланса количества движения (17) для двумерной задачи сводятся к

дe 2 (Я + ц-а)--+(ц + а)У и

дх

д9 „ дю д 2ы

+ = -рП2ы, (20) дх ду дг

дe 2

(Я + ц-а)— + (ц + а) У2у

ду

(18)

де 2

- у--2<

су

(и + Р)У2ю + 2а

до дх д 2у

ду ды ^

дх ду у

.д 2 о д*2

(21)

. (22)

Уравнение теплопроводности высшего порядка (14) с двумя температурами имеет вид

К

( п Хг дг > гТ1 Г! д*г

V 2ф =

( т хг дг > гТ1 Г! д*г

ЧРСЕ | + УТ0 ^

(23)

Определяющие соотношения (1), (2) в плоскости ху записываются как

О х О,

уу

Я+2ц+а X у X Х+2ц + а у X X у

О^,

ху =

ух

ц + а ц-а ц-а ц + а

= («-Р)

до

-2а 2а

до/дх до/ду

ды/дх ду/ ду

е

ду/ дх ды/ду о

т„

= (и-Р) — = (и-р)ту2,

ф-е = аV2ф, V2 = — + —

дх2 ду2

(24)

(25)

(26)

(27)

(28)

Приведенные выше основные уравнения можно использовать для постановки любых граничных условий. В настоящей работе рассматривается микрополярная термоупругая полуобласть со свободной поверхностью под действием зависящего от времени источника тепла, который воздействует на узкую область вокруг оси абсцисс. Зададим соответствующие механические граничные условия

О уу(х y, *) = О ху(x, y, *) = 0 при у = 0,

ху

туг (х, у, *) = 0 при у = 0

(29)

(30)

и тепловое граничное условие

ф(х, у, *) = ф0И(Ь-1 х |)ехр(-Ь*) при у = 0, (31) где ф0 — константа; Н(.) — ступенчатая функция Хевисайда. Данное соотношение указывает на то,

что ширина области воздействия приложенной термомеханической нагрузки не равна нулю только в направлении оси х.

4. Методика решения

Для сокращения числа переменных в системе уравнений и обобщения задачи определим следующие безразмерные величины:

{х', у', ы', V'} = ^{х, у, ы, V}, {е', ф'} = -Ь{е, ф},

ц

Рс1

Ои =

о = -

V 1е>

О/

а + ц

ао

П

ти =

= СГ Г ^, Те, П

ацту

(32)

сх(а + ц)(и + Р)

2 X + 2ц + а к

- , ц=-

а + ц р рСЕ

Опустив штрихи и используя безразмерные переменные, перепишем уравнения движения и баланса количества движения (20)-(22) и модифицированное уравнение (23):

до ду

2 д 2ы

з2 д2у д2ы п2 де

Р2—7+Р2—-Р2—+2—

1 2 дхду ду2

дх2

дх

Р2 д2ы Р2П2ы = Р1^*2-Р2Пы,

Р2 +р2 -р2 2 ^

" '2 дхду дх2 су дх

(33)

ду2

V 2о + ^

(

=Р2 ду -Р2п2у,

(ду ды ^ дх Су

п т дг

Г=1 г! д*'

\

V 2ф =

д 2о

^ 2о = g3 Ц^ =

( Г \

т т д

1+ И

Г=1 Г! д*'

(34)

(35)

де дв

д*

кк

д*

(36)

где

Р2 =

^Р! = Р? - I, «1 = 2

2а2ц2

а + ц

= 4ац

«2 =

С> + Р)

, «3 =

с2(а + ц)(и + Р)

Р/2 «= Щ2?0

(37)

и + Р

с

При этом безразмерные определяющие уравнения (24)-(28) принимают вид

а хх р2 ^2 Р2 "

а уу = Р12 Р2

_Сгг _ _52 Р2 _

а ху " 1 -2"

а ух _ Л 1 2

Си/дх сУ/ ду 0

ду/ дх ди/ду ю

гу

= 54 дю/дх

4 дю/ду

ф-0 = а1 V ф,

(38)

(39)

(40)

(41)

где

«2 =-

X

ц-а и-р С1 53 = ^-, 54 =-^, а1 = -2а. (42)

' и3

а + ц а + ц и + р п

Ввиду большого количества дифференциальных уравнений, помимо переменных поля, введем потенциальные функции Ф(х, у, х) и Т(х, у, х), которые определяются соотношениями

u = VФ + VxT, Т = (0,0,Т). (43)

На основе полученных уравнений выразим компоненты смещений и(х, у, х) и у(х, у, х):

дФ дТ и =--1--

дх ду

дФ дТ

ду дх

Подставляя (44) и (45) в (33)-(36), получим:

(п2 СФ дТ^| п2 С0 „ди Р2-+--р2 — + 2—

дх ду ] дх ду

у = -

(44)

(45)

2

= Р?

С

Р2

=р?

дт-п2

дх2

дФ дТ

ду дх

С д 2

^г-п2

дх2

л

дФ дТ^ -+-

дх ду

-р. 59-2аю

ду дх

(46)

V

дФ дТ

Л

( С2 > ^2 -gз - g2 7

( п ХГ дГ ^

Г=1 г! дхг

V 2ф =

Су дх

ю = glV 2Т,

( г

т х С

1 +

Г=1 Г! дхг

(47)

(48)

Л

^д0 д2 V2Ф^ -+ g—г V 2Ф

V

дх дх2

(49)

Упростим уравнения (46), (47):

С д 2 > V2 -А_-п2

дх2

V2 — рЦ Дг -р?п2

дх2

Ф = 0, (50)

Т = -2ю. (51)

5. Решение для собственных мод

Для решения уравнений (46)-(51) примем, что

(0, ф, и, V, сгу, Ф, тг], Т, ю)

= (0, ф, и, V, а у, Ф, тг], Т, ю)(у) ехр(ф + Пх), (52)

где О — комплексная постоянная угловая частота; £ — волновое число вдоль оси х. Переменные 0, ф, и, V, а у, Ф, Т и ю являются неопределенными амплитудными функциями только по переменной у и не зависят от времени х и координаты х.

Подставляя (52) в (46)-(51), получим следующие уравнения:

(В2 -б1)Ф = 0, (53)

(В2 -С2)ф = 84^(В2 -С2)Ф), (54)

(55)

ф-0 = аДВ2 -С2)ф,

(в2 -е2)т = -2ю,

(56)

где

(В2-Е3)ю = gl(В2-С 2)Т, (57)

в=Е1 =С2 +^2-п2,

Е2=с 2+р2^2 -р?п2,

Е3 =С2 + g2 + gзQ2, Е4 =-1

(58)

& х1

V

п Огхф

4 = 1 , Ц= 1 +У-ф.

а , г! ф 1 Г!

г=1 ' • г=1 ' •

Исключая 0, ф, Ф из (53)-(55) и ю, Т из (56) и (57), получим:

(в4 -пВ2 +У1){0, ф, Ф}(у) = 0, (59)

(В4 -п2В2 +у2){ю, Т}(у) = 0, (60)

где

П =

С 2 + Е1 + Е4[а1 (Е 2 + gС 2) + ЦС 2 + 1)( g + 1)]

а1Е4( g +1) +1 '

П2 =Е2 +Е3 -gl,

С2 +Е1 +Е4(а1С2 + 1)(Е1 + gС 2) (61)

У1 =

а^( g +1) +1

У2 =Е2е3 - 2glС2.

Решения уравнений (59), (60) удовлетворяют критериям непротиворечивости, при которых функции 9, ф, u, v, а-, Ф, ¥ и ю стремятся к нулю при стремлении y к бесконечности, а выражения для них записываются в виде

Ф( У) = Z Ae-X уУ,

1=1

2

¥ (У) = Z Be

-Ц уУ

j=1

2

9( У) = Z Aje"'уУ,

j=1

2

Ф( У) = Z Aj e— уУ,

j=1

2

ю( у) = Z B'f^,

(62)

(63)

(64)

(65)

(66)

j=i

где А/, А/ А"/, В/ и В/ (г = 1, 2, 3) — параметры, зависящие от £ и Кроме того, параметры X/, ц (г = 1, 2) являются положительными корнями следующих уравнений:

X4 -ц^2 +У! = 0, (67)

ц4 -ц2ц2 + у2 = 0. (68)

Подставляя (62)-(66) в (53)-(57), имеем: Aj=(" 2 -si) Aj, a;=-

2 X j — s1

2 _ \ л ЛИ j 1 A

1 — ai("2 —С2) j'

(69)

= ——S2)Bj , j =1,2.

Подставляя (69) в (64)-(66), получим:

9=Z(X2 —si) Aje-уУ,

9=Z

j=1

i jW1*,

j=ii — ai(X2 —С2) j

1 2

2 M

^ = (Ц —£2)B7e^ уУ .

(70)

(71)

(72)

Для определения компонент смещения ы и у подставим (52) в (44), (45), а затем на основе уравнений (62) и (63) получим:

ы = X]А]в~Х/У - £ X}В}в~ЬУ, (73) /=1 /=1

у = -£ XJAJe-X]УX }В }е~цу. (74) /=1 /=1 Используя метод решения для собственных мод, запишем выражения для компонент тепловых и моментных напряжений О/ и т- :

а« = Z «1 j 1 + Z Pi }B}^'y, (75) 1=1 1=1

а yy = Z « 2 у Aje+ Z P2 JBJe-^jy, (76)

j =1

1=1

аzz = Z «зу A/"'уУ + Z Рз jBje—1, (77) 1=1 1=1

а Xy = Z «4 + Z P4 jBje~^y, (78)

у=1

у=1

а yX = Z «5 jAje~X уУ + Z P5 jBje—1, (79)

1=1 j=1

iC 2 2 7=1'

mXz =^Z (Ц —S2) B,e—1, (80)

= 2 Z Цj (ц2 —s2)Bje~^y, (81) 2 j =1

где

«1 j = (62 — P?)X 2 +P?s1 —С 2P?, p1 j = ^Ц j (62 — 6Д

«2 j =P12S1 —С P2 j = i-СЦ j (Pl2 —62), «3 j = —P12(X57 — S1) — C 262, Рз j = —i^A,

«4 j =—2iCX j, P4 j =C2 —S2 + Ц (1 + 63),

(82)

«5 j =—iCXj (1 + 63), P5 j =63C +S2. Используя уравнение (52), сформулируем граничные условия (29)-(31):

а y* (x, y, t) = 0, а Ху (x, y, t) = 0,

m*z (x, y, t) = 0 at y = 0,

(83)

ф(x, y, t) = Ф0 exp(—[iCx + (Q + b)t])

= f (t, C, Q) at y = 0. (84)

Использование граничных условий (83), (84) приводит к

Z «2 jAj + Z Р2 jBj = 0, (85)

1=1 1=1

Z «4 jAj + Z P4 A- = 0, (86)

1=1 1=1

Z Ц у- (ц2 — S2) Bj = 0, (87)

j=1

1 j S1C 2) Ai = f (t, C, Q), (88) i=11 — a1(X- — C )

2 2

(89)

Z« yjAi + ZP yiA3 = 0 1=1 1=1

Рис. 1. Изменение динамической температуры 0 (а) и температуры проводимости ф (б) для различных значений параметра скорости вращения П (цветной в онлайн-версии)

2 2 Уа ухА + Ур ухА3 = а

3=1 3=1

1 2 3 -У еА = 0, 21=1 3 3

У Е3А3 = /,

3=1

Еу =Ху(X.2 -Е2), 3 = 3,4,

Е3 =

х5 -Е1

. 2 ,3 = 1,2. 3 1 -а^2 -С2)

(90)

(91)

(92)

(93)

Решив приведенную выше систему уравнений, зададим константы А3-, Б3- (/' = 1, 2) и получим выражения для различных физических полей.

6. Численные результаты и обсуждение

В разделе обсуждаются результаты численного моделирования физических полей и используемая математическая модель. Для удобства сравнения исследование проведено для кристаллического магния. В расчетах использовались стандартная температура Т0 = 298 К и следующие значения физических параметров [54]:

X = 9.4 • 1010 Н/м2, | = 4.5 • 1010 Н/м2, а = 0.5 • 1010 Н/м2, и + в = 0.779 • 10-5 Н/м2, К = 1.7 • 102 Дж/(мс - К), Се = 1.04 • 103 Дж/(кг • К), р = 1.74 • 103 кг/м3,3 = 0.2 • 10-19 м2, у = 0.779 • 10-9 Н.

Рис. 2. Изменение нормальных перемещений и (а) и касательных перемещений V (б) для различных значений параметра скорости вращения П (цветной в онлайн-версии)

Рис. 3. Изменение нормальных напряжений о^ (а), ауу (б) и касательных напряжений оух (в) для различных значений параметра скорости вращения П (цветной в онлайн-версии)

Примем, что О является действительной величиной (О = 2) при малых временах. Остальные константы в численных расчетах имеют следующие значения [56]: £ = 2, Ь = 2, а1 = 0.2, ф0 = 1, Ь = 1.

Расчеты выполняются для безразмерного времени х=0.2 по поверхности х = 0.1 на интервале 0< у < 3, при этом действительная составляющая амплитуды переменных поля отложена по вертикальной оси. Распределения температуры проводимости и динамической температуры (ф, 0), термических напряжений (охх, оуу, оух), компонент вектора перемещений (и, V), микроповорота ю и компонент моментных напряжений туг вдоль оси у представлены в графическом виде при различных эффективных параметрах, таких как коэффициент а1, параметр поворота, а также производные высшего порядка по времени т и п. Также проведено сравнение различных моделей с учетом и без учета микрополярных эффектов.

Из рисунков видно, что двухтемпературный коэффициент а1 и высшие производные по времени т, п оказывают существенное влияние на все исследуемые поля. В зависимости от значений производных высшего порядка по времени и коэффициента расхождения а1, тепловые и механические волны могут достигать устойчивого состояния. Учет влияния этих параметров позволяет разграничить вклад термодинамической температуры и температуры проводимости.

6.1. Влияние вращения на распределение параметров

Рассмотрим изменение безразмерной температуры проводимости ф и динамической температуры 0, термических напряжений охх, оуу, оух, компонент вектора перемещения и, V, микроповорота ю и компонент моментных напряжений туг вдоль оси у при различных значениях параметра вращения П. Для вращающейся среды задаются значения П = 1, 1.2, 1.4 и 1.6. В случае отсутствия вращения П = 0. На рис. 1-4 показаны результаты расчетов влияния поворота на различные переменные теплофизического поля, полученные с использованием двухтемпературной микрополярной модели термоупругости с двухфазным запаздыванием и высшими производными по времени (2МЫБРЬ).

На рис. 1 представлены зависимости температуры проводимости ф и динамической температуры 0 для заданных значений параметра вращения П. Из рисунков видно, что параметр вращения влияет на распределение температур. Так, значения динамической температуры 0 при вращении меньше, чем при отсутствии вращения. Также наблюдается слабое влияние параметра П на температуру проводимости ф. Можно сделать вывод,

Рис. 4. Изменение микроповорота ю (а) и моментных напряжений тух (б) для различных значений параметра скорости вращения П (цветной в онлайн-версии)

что обе температуры ф и 9 уменьшаются с увеличением расстояния по оси у при любых П и температура проводимости ф удовлетворяет граничным условиям. Другим важным выводом является то, что тепловая волна имеет ограниченную скорость распространения в среде, и в этом предложенная модель 2МЫБРЬ согласуется с обобщенной теорией термоупругости.

Изменение компонент перемещения ы и у вдоль оси у показано на рис. 2 для различных значений параметра вращения П. Видно, что смещения уменьшаются с увеличением расстояния у. Кроме того, значения тепловых смещений больше при более высоких значениях параметра П. Из рисунков следует, что деформация среды зависит не только от типа граничных условий задачи, но и от скорости вращения среды.

На рис. 3, а представлены кривые нормальных напряжений <зхх для заданных значений параметра вращения П. Из рисунка видно, что в зависимости от угловой скорости П напряжения ахх на поверхности имеют разные начальные значения, которые быстро увеличиваются, пока не достигнут максимальных значений в точке у = 1.2. При дальнейшем увеличении расстояния у кривые напряжения постепенно снижаются и покидают область термической турбулентности. Это указывает на то, что учет вращения в основных уравнениях критически влияет на распределение термических напряжений вхх. С ростом скорости вращения П вблизи поверхности у = 0 (в диапазоне 0 < у < 1.2) напряжения падают, а затем возрастают при у больше 1.2.

На рис. 3, б, в показано влияние параметра вращения П на изменение термических напряжений

вуу и вух соответственно. Напряжения имеют нулевые значения в точке у = 0, что удовлетворяет граничным условиям задачи. Кривые напряжений демонстрируют схожее поведение: выходя из начала координат, они быстро убывают до минимальных значений, а затем снова постепенно возрастают до нулевых значений по мере удаления от поверхности. Выявлено, что с ростом скорости вращения упругого тела амплитуда нормальных напряжений оуу уменьшается, а амплитуда касательных напряжений вух увеличивается.

На рис. 4, а показано изменение значения микроповорота ю вдоль оси у в зависимости от угловой скорости П. Во всех случаях микроповорот ю имеет отрицательные значения на поверхности, а затем постепенно увеличивается до нулевых значений по мере удаления от поверхности. Параметр вращения оказыват заметное влияние на микроповорот ю, который уменьшается с увеличением угловой скорости П.

На рис. 4, б построены кривые моментных напряжений туг при разных значениях скорости вращения П. Поведение кривых на диаграмме одинаково. В соответствии с граничными условиями задачи на поверхности у = 0 моментные напряжения туг равны нулю. Затем они постепенно увеличиваются до максимальных значений при у = 1.2 и постепенно снижаются до нулевых значений по мере удаления от источника тепла и переходной зоны. Кроме того, угловая скорость вращения оказывает существенное влияние на распределение моментных напряжений туг, что выражается в снижении пиковых значений с увеличением параметра П.

— 2МЬБ модель

— 2МБРЬ модель

— 2НМБРЬ модель

2 Н Э Р Ь модель

— 1Н Э Р Ь модель

\а

Рис. 5. Распределение динамической температуры 0 (а) и

термоупругости (цветной в онлайн-версии)

6.2. Сравнение различных микрополярных моделей термоупругости

Рассмотрим изменение значений переменных поля с расстоянием у, расчитанных с помощью двухтемпературной микрополярной модели термоупругости с одним временем релаксации (2МЬБ), двухтемпературной микрополярной модели термоупругости с двухфазным запаздыванием (2МБРЬ), а также одно- и двухтемпературных термоупругих моделей с двухфазным запаздыванием и высшими производными по времени (2МЫБРЬ, 2ЫБРЬ, 1ЫБРЬ). Отметим, что построение двухтемпературных моделей термоупругости возможно только без учета влияния микрополярности: а = в = и = в =3 = 0. Для однотемпе-ратурных моделей а1 = 0.

Рис. 6. Распределение нормальных перемещений и (а) и термоупругости (цветной в онлайн-версии)

температуры проводимости ф (б) в различных моделях

На рис. 5, а показано влияние микрополярности, высших производных по времени т, п и коэффициента а1 на изменение температуры вдоль оси у. Видно, что в направлении распространения волны вдоль оси у динамическая температура постепенно снижается. Кроме того, в рассматриваемых моделях температура быстро достигает максимальных значений («пиковая температура») в результате уменьшения тепловой нагрузки. Затем в пределах определенной области температура постепенно уменьшается до нулевых значений. Приведенные данные подтверждают, что в моделях без учета микрополярности значения температуры намного выше, чем в микрополярных моделях. Сравнение результатов, полученных с помощью моделей 2МЫБРЬ и 2МБРЬ, показывает,

касательных перемещений V (б) в различных моделях

что наличие высших производных т и п приводит к ограничению распространения тепла в среде.

Другим важным наблюдением является влияние двухтемпературного параметра а! на распределение температуры. Температурная кривая в модели 2ЫБРЬ оказалась намного ниже температурной кривой в модели 1ЫБРЬ.

Как показано на рис. 5, б, все случаи с постоянным значением температуры проводимости ф на свободной поверхности среды находятся в полном соответствии с тепловыми граничными условиями. В обобщенной модели начальная температура проводимости имеет самое высокое значение (за счет диффузии тепловой границы) и постепенно снижается в направлении тепловой волны практически до нуля. По сравнению с двух-температурными моделями температура проводимости ф в модели 1ЫБР без коэффициента расхождения а1 имеет меньшие значения. На рисунке наблюдается расхождение кривых, построенных по моделям с учетом и без учета микрополярности. Кроме того, производные высшего порядка т и п существенно снижают температуру теплопроводности.

На рис. 6 показано распределение нормальных и поперечных перемещений и и V, рассчитанных в рамках различных теорий термоупругости. Согласно полученным результатам, наличие или отсутствие двухтемпературного параметра а1 оказывает сильное влияние на распределение перемещений. Производные высшего порядка по времени т и п также оказывают заметное влияние на распределение полей перемещений и и V. Кривые смещений имеют одинаковый характер в трех микрополярных моделях (2МЬБ, 2МБРЬ, 2МЫБРЬ) и существенно отличаются от кривых в двух других моделях без учета микрополярного эффекта (2ЫБРЬ, 1ЫБРЬ).

Изменение нормальных и касательных напряжений ахх, оу,, оух проиллюстрировано на рис. 7 в зависимости от двухтемпературного параметра а1 и производных высшего порядка т и п для моделей с учетом и без учета микрополярности. Из рисунков видно, что двухтемпературный параметр а1 оказывает существенное влияние на все переменные поля. Влияние параметра а1 и производных высшего порядка т, п заключается в снижении уровня нормальных и касательных напряжений. В зависимости от значений разности температур, а также высших производных, волны приближаются к стационарному состоянию.

/ — 2МЬБ модель

/ — 2МБРЬ модель

/ — 2НМБРЬ модель

/ 2 Н Э Р Ь модель1

/ — ИГОРЬ модель [а

1.0 1.4 1.8 2.2 2.6 у

0 1 2 3 4 ><

Рис. 7. Распределение нормальных напряжений ахх (а), Суу (б) и касательных напряжений сух (в) в различных моделях термоупругости (цветной в онлайн-вер-сии)

Кривые на рис. 7, а, построеные по моделям с учетом и без учета микрополярности, указывают на заметное влияние выбора модели на распределение нормального напряжения ахх. Из рис. 7, б, в

Рис. 8. Распределение микроповорота ю (а) и моментных

(цветной в онлайн-версии)

видно, что абсолютные значения напряжения ауу в микрополярных теориях (2MLS, 2MDPL, 2MHDPL) намного меньше значений в моделях 2HDPL и 1HDPL, не учитывающих микрополярные эффекты. По кривым, представленым на рис. 7, б, в для разных теорий, можно сделать вывод, что напряжения ауу и аух удовлетворяют механическим граничным условиям на свободной поверхности среды.

На рис. 8 показано изменение микроповорота ю и касательных моментных напряжений myz для трех микрополярных теорий (2MLS, 2MDPL, 2MHDPL) и двух теорий без учета микрополярного эффекта (2HDPL, 1HDPL). Наличие безразмерного коэффициента расхождения a1, а также производных высшего порядка по времени m и n оказывает существенное влияние на поведение кривых.

Из рис. 8 видно, что микрополярные модели (2MLS, 2MDPL, 2MHDPL) дают более высокие абсолютные значения микроповорота ю по сравнению с моделями 2HDPL и 1HDPL и наблюдается сильное расхождение решений для микрополярных моделей с решениями для моделей 2HDPL и 1HDPL. Кроме того, все кривые мо-ментных напряжений исходят из начала координат, что удовлетворяет граничным условиям и подтверждает точность расчетов.

7. Заключение

В работе представлена двухтемпературная термоупругая модель, разработанная в рамках обобщенной микрополярной теории термоупругости. Первая температура обусловлена механическим процессом (динамическая температура), а вторая

является следствием теплового процесса (температура проводимости). Уравнение теплопроводности включает частные производные высшего порядка. Отсутствие учета коэффициента расхождения позволяет получить различные теории, являющиеся частными случаями предлагаемой модели. Методом анализа собственных мод изучено поведение термических напряжений, температур, перемещений и микроповорота для однородного изотопного вращающегося микрополярного термоупругого материала.

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

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

Частные производные высшего порядка оказывают значительное влияние на распределение переменных поля. Используя данные производные, можно ограничить распространение тепловых и механических волн. Таким образом, модифицированная модель с фазовым запаздыванием и производными высшего порядка является более эффек-

тивной по сравнению с предыдущими обобщенными теориями и лучше подходит для изучения распространения тепловых волн.

Учет микрополярности оказывает существенное влияние на теплофизические поля, приводя к уменьшению их значений, за исключением поперечного напряжения и температуры проводимости.

Работа выполнена при финансовой поддержке Университета Джуфа в рамках гранта № DSR-2021-03-0383. Авторы выражают благодарность Колледжу науки и искусств в Аль-Курайяте за техническую поддержку.

Литература

1. Subba Rao V., Nigam S.D. Wave propagation in rotating elastic media // Mathematica. - 1964. - P. 29-38.

2. Durukan Y., Shevelko M., Peregudov A., Popkova E., Shev-chenko S. The effect of a rotating medium on bulk acoustic wave polarization: From theoretical considerations to perspective angular motion sensor design // Sensors. - 2020. -V. 20(9). - P. 2487.

3. Khan A., Islam S., Khan M. Rayleigh waves in a rotating orthotropic medium // World Appl. Sci. J. - 2012. -V. 16(1). - P. 73-75.

4. Jafari S. Elastic limit angular velocity and acceleration investigation in non-uniform rotating disk under time-dependent mechanical loading // J. Appl. Comput. Mech. - 2022. -V. 8(3). - P. 791-808. - https://doi.org/10.22055/jacm.2020. 32914.2099

5. Abouelregal A.E., Atta D., Sedighi H.M. Vibrational behavior of thermoelastic rotating nanobeams with variable thermal properties based on memory-dependent derivative of heat conduction model // Arch. Appl. Mech. - 2022. -https://doi.org/10.1007/s00419-022-02110-8

6. Sladek J., Sladek V., Repka M. The heat conduction in nano-sized structures // Phys. Mesomech. - 2021. - V. 24. -No. 5. - P. 611-617. - https://doi.org/10.1134/S1029959921 05012X

7. Abouelregal A.E., Sedighi H.M., Malikan M., Eremeyev V.A. Nonlocalized thermal behavior of rotating micromachined beams under dynamic and thermodynamic loads // ZAMM-J. Appl. Math. Mech. - 2022. - V. 102(4). - P. e202100310.

8. Abouelregal A.E., Abo-Dahab S.M. A two-dimensional problem of a mode-I crack in a rotating fibre-reinforced isotropic thermoelastic medium under dual-phase-lag model // Sadha-na. - 2018. - V. 43. - P. 13.

9. Abouelregal A.E., Mohammad-Sedighi H., Faghidian S.A., Shirazi A.H. Temperature-dependent physical characteristics of the rotating nonlocal nanobeams subject to a varying heat source and a dynamic load // Facta Univ. Ser. Mech. Eng. -2021. - V. 19(4). - P. 633-656.

10. Abouelregal A.E., Marin M., Askar S. Thermo-optical mechanical waves in a rotating solid semiconductor sphere using the improved Green-Naghdi III model // Mathematics. -2021. - V. 9(22). - P. 2902.

11. Biot M. Thermoelasticity and irreversible thermodynamics // J. Appl. Phys. - 1956. - V. 27. - P. 240-253.

12. LordH.W., Shulman Y.H. A generalized dynamical theory of thermoelasticity // J. Mech. Phys. Solids. - 1967. - V. 15. -P. 299-309.

13. Green A.E., Naghdi P.M. A re-examination of the basic results of thermomechanics // Proc. Math. Phys. Sci. - 1991. -V. 432. - P. 171-194.

14. Green A.E., Lindsay K.A. Thermoelasticity // J. Elasticity. -1972. - V. 2. - P. 1-7.

15. Green A.E., Naghdi P.M. On undamped heat waves in an elastic solid // J. Therm. Stress. - 1992. - V. 15. - P. 252264.

16. Green A.E., Naghdi P.M. Thermoelasticity without energy dissipation // J. Elasticity. - 1994. - V. 31. - P. 189-208.

17. Tzou D.Y. A unified filed approach for heat conduction from macro to macroscales // ASME J. Heat Transfer. - 1995. -V. 117. - P. 8-16.

18. Tzou D.Y. The generalized lagging response in small-scale and high-rate heating // Int. J. Heat Mass Transfer. - 1995. -V. 38. - P. 3231-3234.

19. Tzou D.Y. Experimental support for the lagging behavior in heat propagation // J. Thermophys. Heat Transfer. - 1995. -V. 9. - P. 686-693.

20. Abouelregal A.E. Two-temperature thermoelastic model without energy dissipation including higher order time-derivatives and two phase-lags // Mater. Res. Express. - 2019. -V. 6. - P. 116535.

21. Abouelregal A.E. On Green and Naghdi thermoelasticity model without energy dissipation with higher order time differential and phase-lags // J. Appl. Comput. Mech. - 2020. -V. 6(3). - P. 445-456.

22. Abouelregal A.E. A novel generalized thermoelasticity with higher-order time-derivatives and three-phase lags // Multi-discipl. Model. Mater. Struct. - 2019. - V. 16. - P. 689-711.

23. Abouelregal A.E. A novel model of nonlocal thermoelas-ticity with time derivatives of higher order // Math. Meth. Appl. Sci. - 2020. - V. 43. - P. 6746-6760.

24. Abouelregal A.E. Three-phase-lag thermoelastic heat conduction model with higher-order time-fractional derivatives // Indian J. Phys. - 2020. - V. 94. - P. 1949-1963.

25. Choudhuri S.R. On a thermoelastic three-phase-lag model // J. Therm. Stress. - 2020. - V. 30. - P. 231-238.

26. Chen P.J., Gurtin M.E. On a theory of heat conduction involving two temperatures // Zeitschrift Angewandte Math. Phys. - 1968. - V. 19. - P. 614-627.

27. Chen P.J., Williams W.O. A note on non-simple heat conduction // Zeitschrift Angewandte Math. Phys. - 1968. -V. 19. - P. 969-970.

28. Chen P.J., Gurtin M.E., Williams W.O. On the thermodynamics of non-simple elastic materials with two temperatures // Zeitschrift Angewandte Math. Phys. - 1969. - V. 20 -P. 107-112.

29. Quintanilla R. On existence, structural stability, convergence and spatial behavior in thermoelasticity with two temperatures // Acta Mech. - 2004. - V. 168. - P. 61-73.

30. Youssef H. Theory of two-temperature-generalized thermoelasticity // IMA J. Appl. Math. - 2006. - V. 71. - P. 383-390.

31. Ezzat M.A., El-Karamany A.S. Two temperature theory in generalized magneto thermoelasticity with two relaxation times // Meccanica. - 2011. - V. 46. - P. 785-794.

32. Mukhopadhyay S., Prasad R., Kumar R. On the theory of two-temperature thermoelasticity with two phase-lags // J. Therm. Stress. - 2011. - V. 34. - P. 352-365.

33. Kumar R., Prasad R., Kumar R. Thermoelastic interactions on hyperbolic two-temperature generalized thermoelasticity in an infinite medium with a cylindrical cavity // Eur. J. Mech. A. Solids. - 2020. - V. 82. - P. 104007.

34. Fernández J.R., Quintanilla R. Uniqueness and exponential instability in a new two-temperature thermoelastic theory // AIMS Mathematics. - 2021. - V. 6(6). - P. 5440-5451.

35. Sarkar N., Mondal S. Two-dimensional problem of two-temperature generalized thermoelasticity using memory-dependent heat transfer: An integral transform approach // Indian J. Phys. - 2020. - V. 94(12). - P. 1965-1974.

36. Hobiny A., Abbas I., Marin M. The influences of the hyperbolic two-temperatures theory on waves propagation in a semiconductor material containing spherical cavity // Mathematics. - 2022. - V. 10(1). - P. 121.

37. Hassanpour S., Heppler G.R. Micropolar elasticity theory: A survey of linear isotropic equations, representative notations, and experimental investigations // Math. Mech. Solids. -2017. - V. 22(2). - P. 224-242.

38. Eringen A.C. Linear theory of micropolar elasticity // J. Appl. Math. Mech. - 1966. - V. 15. - P. 909-923.

39. Nowacki W. Theory of Asymmetric Elasticity. - Oxford: Pergamon, 2004.

40. Eringen A.C. Foundations of Micropolar Thermoelasticity, International Centre for Mechanical Science, Udine Course and Lectures 23. - Berlin: Springer-Verlag, 1970.

41. Tauchert T.R., Claus W.D., Jr., Ariman T. The linear theory of micropolar thermoelasticity // Int. J. Eng. Sci. - 1968. -V. 6. - P. 37-47.

42. Dost S., Tabarrok B. Generalized micropolar thermoelasticity // Int. J. Eng. Sci. - 1978. - V. 16. - P. 173.

43. Chandrasekhariah. Heat flux dependent micropolar elasticity // Int. J. Eng. Sci. - 1986. - V. 24. - P. 1389-1395.

44. El-Karamany A.S., EzzatM.A. On the three-phase-lag linear micropolar thermoelasticity theory // Eur. J. Mech. A. Solids. - 2013. - V. 40. - P. 198-208.

45. Alharbi A.M., Said S.M., Abd-Elaziz E.M., Othman M.I.A. Mathematical model for a magneto-thermoelastic micropolar medium with temperature-dependent material moduli under the effect of mechanical strip load // Acta Mech. - 2021. -V. 232(6). - P. 2331-2346.

46. Chirila A., Agarwal R.P., Marin M. Proving uniqueness for the solution of the problem of homogeneous and anisotropic micropolar thermoelasticity // Bound. Value Probl. - 2017. -V. 3. - P. 2017.

47. Sharma H., Kumari S., Kumar A. Study of micropolar thermo-elasticity // Adv. Appl. Math. Sci. - 2020. - V. 19. -P. 929-941.

48. Hilal M.I.M., Abd-Elaziz E.M., Hanoura S.A. Reflection of plane waves in magneto-micropolar thermoelastic medium

with voids and one relaxation time due to gravity and two-temperature theory // Indian J. Phys. - 2021. - V. 95. -P. 915-924.

49. Kumar R., Prasad R., Kumar R. Thermoelastic interactions on hyperbolic two-temperature generalized thermoelasticity in an infinite medium with a cylindrical cavity // Eur. J. Mech. A. Solids. - 2020. - V. 8. - P. 104007.

50. Lianngenga R., Singh S.S. Reflection of coupled dilatational and shear waves in the generalized micropolar thermoelastic materials // J. Vibr. Control. - 2020. - V. 26. - P. 19481955.

51. Marin M., Chirila A., Ochsner A., Vlase S. About finite energy solutions in thermoelasticity of micropolar bodies with voids // Bound. Value Probl. - 2019. - V. 89. - P. 2019.

52. Abouelregal A.E., Zenkour A.M. Two-temperature thermo-elastic surface waves in micropolar thermoelastic media via dual-phase-lag model // Adv. Aircraft Spacecraft Sci. -2017. - V. 4. - P. 711-727.

53. Guesmia A., MuNoz Rivera J.E., Sepulveda CortEs M.A., Vera VillagrAn O. Well-posedness and stability of a generalized micropolar thermoelastic body with infinite memory // Quarterly J. Mathematics. - 2021. - V. 72. - P. 1495-1515.

54. Kumar R., Abbas I.A. Deformation due to thermal source in micropolar thermoelastic media with thermal and conductive temperatures // J. Comput. Theor. Nanosci. - 2013. -V. 10. - P. 2241-2247.

55. Shaw S., Mukhopadhyay B. Moving heat source response in micropolar half-space with two-temperature theory // Continuum Mech. Thermodyn. - 2013. - V. 25. - P. 523-535.

56. Ezzat M.A., Awad E.S. Constitutive relations, uniqueness of solution and thermal shock application in the linear theory of micropolar generalized thermoelasticity involving two temperatures // J. Therm. Stress. - 2010. - V. 33. - P. 226-250.

57. Quintanilla R. Exponential stability and uniqueness in ther-moelasticity with two temperatures // Dyn. Contin. Discrete Impuls. Syst. A. Math. Anal. - 2004. - V. 11. - P. 57-68.

58. Quintanilla R. A well posed problem for the dual-phase-lag heat conduction // J. Therm. Stress. - 2008. - V. 31. -P. 260-269.

59. Chirita S. On the time differential dual-phase-lag thermoelastic model // Meccanica. - 2017. - V. 52. - P. 349-361.

60. Chirita S., Ciarletta M., Tibullo V. On the thermomechanic consistency of the time differential dual-phase-lag models of heat conduction // Int. J. Heat Mass Transfer. - 2017. -V. 114. - P. 277-285.

61. Chirita S., Ciarletta M., Tibullo V. On the wave propagation in the time differential dual-phase-lag thermoelastic model // Proc. R. Soc. A. - 2015. - V. 471. - P. 20150400.

Поступила в редакцию 03.03.2022 г., после доработки 29.05.2022 г., принята к публикации 30.05.2022 г.

Сведения об авторах

Ahmed E. Abouelregal, Prof., Jouf University, Saudi Arabia; Mansoura University, Egypt, [email protected] Rayan Alanazi, Assist. Prof., Jouf University, Saudi Arabia, [email protected]

Abdullah H. Sofiyev, Prof., Suleyman Demirel University, Istanbul Commerce University, Turkey; Azerbaijan State Economic University, Azerbaijan, [email protected], [email protected]

Hamid M. Sedighi, Ph.D., Assoc. Prof., Shahid Chamran University of Ahvaz, Iran, [email protected], [email protected]