Распространение термоупругих волн в среде с учетом релаксации теплового потока

Бабенков, Михаил Борисович

Распространение термоупругих волн в среде с учетом релаксации теплового потока тема автореферата и диссертации по механике, 01.02.04 ВАК РФ

Бабенков, Михаил Борисович АВТОР

кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ

Санкт-Петербург МЕСТО ЗАЩИТЫ

2013 ГОД ЗАЩИТЫ

01.02.04 КОД ВАК РФ

Диссертация по механике на тему «Распространение термоупругих волн в среде с учетом релаксации теплового потока»

Автореферат диссертации на тему "Распространение термоупругих волн в среде с учетом релаксации теплового потока"

На правах рукописи

Бабепков Михаил Борисович

Распространение термоупругих волн в

среде с учетом релаксации теплового

потока

01.02.04 — Мехаиика деформируемого твердого тела

Автореферат £ ^ ^

диссертации на соискание ученой степени кандидата физико-математических паук

Санкт-Петербург 2013

005540522

Работа выполнена в Федеральном Государственном Бюджетном Учреждение Науки Институт Проблем Машиноведения Российской Академии Наук

Научный руководитель:

Официальные оппоненты:

Ведущая организация:

доктор физико-математических наук Иванова Елена Александровна

Муратиков Кирилл Львович доктор физико-математических наук, ФГБУН ФТИ им.Иоффе РАН, вед. науч. сотр. Отделения твердотельной электропики

Еремеев Виктор Анатольевич доктор физико-математических паук, доцент Южный научный центр РАН, зав. лаб. Механики активных материалов

ФГБУН Институт Проблем Машиностроения Российской Академии Наук

Защита состоится 19 декабря 2013 г. в 14— часов на заседании диссертационного совета Д 002.075.01 при Институте проблем машиноведения РАН по адресу: 199178, Санкт-Петербург, В.О., Большой пр., 61.

С диссертацией можно ознакомиться в ОНТИ ИПМаш РАН.

Автореферат разослан « » ноября 2013 г.

Ученый секретарь диссертационного совета д.т.н., профессор

В.В. Дубаренко

Общая характеристика работы

Актуальность темы исследования. Исследование волнового переноса тепла актуально для многих развивающихся технологий. Внутренние источиики ианоразмерного масштабного уровня (например, размер современного транзистора составляет всего несколько нанометров) вызывают более интенсивный рост температуры, чем предсказывает классическая теория, что повышает требования к будущим системам охлаждения. Температурные эффекты вносят значительные изменения в механические свойства нанопла-стин, так как чем тоньше пластина, тем она чувствительнее к изменениям температуры. Исследования процессов волнового переноса тепла могут быть полезны для изучения термомеханических характеристик объектов микро-и нано- масштабного уровня: тонких пластин и стержней, используемых в микро- и пало- электромеханических устройствах (MEMS и NEMS).

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

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

Обобщенный закон Фурье был предложен в работах авторов: С. Cattaneo, Р. Vcrnottc, A.B. Лыкова. В отличии от классического закона, он учитывает инерционность процессов тсплопереноса, которая характеризуется постоянной релаксации теплового потока.

Экспериментальным определением величин релаксации теплового потока занимались: R.H. Matsunaga, I. Santos, W. Kaminski, К. Mitra et. al., А. Grabmann, F. Peters, H. Herwig, K. Bockert, W. Roctzela et. al.

На основании обобщенного закона Фурье в работе Е.Б. Попова впервые получено гиперболическое уравнение связанной термоупругости, в линейном приближении учитывающее конечность скорости распространения тепла. Более общая модель термоупругости была предложена в работе авторов Н. Lord, A. Shulman. Тсрмоупругость Лорда-Шульмана (LS) является предметом изучения в данной работе. Обзор по теоретическим основам неклассической тер-

моупругости и теплопроводности представлен в статьях и книгах авторов: А.Г. Шашков, В.А. Бубнов, С.Ю. Яновский, В.А. Ковалев, Ю.Н. Радаев, Д. Жоу, X. Касас-Баскес, Дж. Лебон, J. Ignaczak, D.S. Chandrasekharaiah.

Исследованием плоских гармонических воли, распространяющихся в термоупругой среде с релаксацией теплового потока занимались: А.Н. Nayfeh, S. Nemat-Nasser, Ю.К. Энгельбрехт, P. Puri, И.М. Штер, Ф.В. Семерак, Ц. Иванов, Ю.К. Энгельбрехт, Р.Х. Швед, A.A. Лопатьев. Результаты некоторых из перечисленных исследований приводятся в книгах авторов: А.Г. Шашков, В.А. Бубнов, С.Ю. Яновский, А.Д. Коваленко, Я.С. Подстри-гач, Ю.М. Коляно.

В исследование гиперболической задачи теплопроводности значительный вклад внесли: А.Г. Шашков, В.А. Бубнов, С.Ю. Яновский, В.А. Ку-динов, Д. Жоу, X. Касас-Басксс, Дж. Лсбон, K.L. Baumeister, T.D. Hamill, M.N. Özisik, M. Lewandowska, D. Zhang, B.S. Yilbas.

Следующие авторы исследовали задачи распространения термоупругих волн LS типа в среде, подверженной воздействию теплового импульса: F.R. Norwood, W.E. Warren, Г.А. Кильчинская, А.Н. Nayfeh, S. Neinat-Nasser, М. Baila, P.M. Jordan, P. Puri, I.A. Abdallah, N. Sarkar, E.F. Hcnain, J.C. Strikwerda, A.M. Scott, H.M. Youssef, A.S. Al-Felali.

Экспериментальным исследованием термоупругих напряжений при неравновесных процессах теплообмена занимались: Н.В. Вовпенко, Б.А. Зимин, Ю.В. Судьенков, K.V. Polctkin, G.G. Gurzadyan, J. Shang, V. Kulish, O.B. Wright, V.E. Guscv и др.

Целью данной работы является изучение поведения термоупругой среды, описываемой уравнениями Лорда-Шульмапа, под воздействием периодических и импульсных возмущений.

Задачи исследования:

1. Исследовать распространение термоупругих волн в широком диапазоне значений постоянной релаксации теплового потока. Данная задача представляется актуальной, поскольку теоретическая оценка постоянной релаксации в металлах, согласно фонопной теории, составляет несколько пикосекунд (10~12с). В то время как экспериментальные данные дают разброс результатов, отличающихся от теоретической оценки на несколько порядков: от 10-8с до 10~ис.

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

3. Сравнить две модели импульсного лазерного воздействия на среду. Первая модель: среда нагревается источниками тепла, заданными на границе, вторая модель: среда нагревается источниками тепла распределенными в объеме.

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

Научную новизну работы представляют следующие положения, выносимые на защиту:

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

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

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

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

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

6. Установлено, что при скоростях воздействия, при которых необходим учет гиперболичности, принципиально важным является учет связанности.

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

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

Методы исследования. Основные результаты работы получены аналитически с использованием широко известных методов математической физики: метода Эйлера, метода функций Грина, метода преобразования Лапласа и метода Гринберга.

Достоверность и апробация результатов. Достоверность изложенных в работе результатов обусловлена строгостью формулировок задач математической физики, использованием фундаментальных принципов механики, а также сравнением с результатами, полученными ранее. Основные результаты работы доложены на: XLI, XL и XXXVIII International Summer School Conference «Advanced Problems in Mechanics» (Россия, С.-Петербург, 2013, 2012 и 2010 гг.), International Conference «Days on Diffraction» (Россия, С.-Петербург, 2013 и 2012 гг.), VI Поляховские чтения (Россия, С.-Петербург, 2012), 2nd International Conference on Material Modelling (Франция, Париж, 2011), на Санкт-Петербургском Семинаре по Вычислительной и Теоретической Акустике Научного Совета РАН по Акустике (руководитель проф. Д.П. Коузов, С.-Петербург, 2009, 2010, 2011 и 2012 гг.), па Городском семинаре по механике (руководитель чл.-корр. РАН, проф. Д.А. Индейцев, С.-Петербург, 2013 г.)

Публикации. Основные результаты по теме диссертации изложены в 10 печатных изданиях, 3 из которых изданы в журналах, рекомендованных ВАК, 7 — в сборниках тезисов докладов.

Структура и объем диссертации. Диссертация состоит из введения, трех глав, заключения и приложения. Полный объем диссертации 102 страницы текста с 26 рисунками и 3 таблицами. Список литературы содержит 93 наименований.

Содержание работы.

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

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

Для моделирования лазерного воздействия можно задать распределение источников тепла на границе или в объеме среды. Для непрозрачных сред задание распределения источников тепла па границе является обоснованным приближением, поскольку глубина проникновения лазерного излучения в материал (толщина скин-слоя) будет небольшой, для металлов порядка 10"° — 10~7м при частотах лазерного излучения от 10ГГЦ до 1ТГц (инфракрасный спектр). Но данная глубина проникновения может оказаться существенной при рассмотрении быстро протекающих процессов. В таком случае необходимо задавать распределение источников тепла в объеме и рассматривать неоднородную задачу теплопроводности. Подобное усложнение задачи должно привести к появлению новых эффектов. Система дифференциальных уравнений1, позволяющая найти распределение температуры при заданных начальных и граничных условиях, состоит из уравнения баланса энергии и обобщенного закона Фурье (закон Каттансо-Вернотте):

Здесь h — вектор теплового потока, V — оператор пабла, q — мощность источников тепла, распределенных в объеме.

Рассмотрим распространение тепловых волн в полупространстве вдоль координаты s. Пусть q(s,t) = 0 и тепловой поток задан на границе полупространства: h|s=o = ho5(t), где /io — множитель, добавленный для сохранения размерности, S(t) — функция Дирака, которая аппроксимирует временной профиль короткого лазерного импульса. Тогда распределение температуры при однородных начальных условиях будет иметь вид2, представленный на рис. 2, а.

1Левапов Е.И., Сотский E.H. Некоторые свойства процессов теплопереноса в неподвижной среде с учетом постоянной |к;лаксации теплового потока // Инженерно-Физический Журнал. 1986. Т.50. N.6. С.1017-1024

2Новиков H.A. Гиперболическое уравнение теплопроводности. Решение прямых и обратных задач для полуограшпеппого стержня // Инженерно-Физический Журнал. 1977. Т.35. N.4. С.734-740

pCvT = - V • h + pq h + rh = -A VT

(1) (2)

I? s

0 I

—

Рис. 1: Слой, находящийся под воздействием лазерного импульса

Зададим распределение источников тепла в объеме неограниченного слоя толщины I (рис. 1). Пусть интенсивность лазерного импульса поглощается средой по закону Бугера, тогда мощность источников тепла можно представить в виде:

д(а,*) = 7о*(*-0)ехр(-/?в) (3)

Здесь /0 - интенсивность света па входе в поглощающий слой, /3 — коэффициент поглощения вещества. На границах слоя заданы условия теплоизоляции: /г|8=о = 0; к\,¡=1 = 0. При однородных начальных условиях, распределение температуры будет иметь вид3, представленный па рис. 2,6. Если устремить значение (5 к бесконечности, то получится распределение температуры, показанное на рис. 2,а.

Рис. 2: Температурные профили при разных значениях времени порядка постоянной релаксации т. График а — тепловой поток задай па границе. График б — распределение источников тепла задано в объеме. Штриховой линией показана огибающая.

На волновом фронте решения задачи с тепловым воздействием на границе значение температуры обращается в бесконечность (рис. 2,а). В задаче с распределенными источниками тепла, мощность которых зависит от времени как ¿-функция, вместо температурного скачка на фронте наблюдается точка излома (рис. 2,6). В задаче с тепловым воздействием па границе график температуры обрывается на волновом фронте, тогда как в задаче с объемным воздействием значение температуры убывает до нуля постепенно. Проведенный сравнительный анализ позволяет сделать вывод, что задача с объемными источниками тепла значительно лучше описывает процесс тсплопсрепоса в слое, находящимся под воздействием лазерного импульса.

^Lewandowska M. Hyperbolic heat conduction in the semi-infinite body with a timc-dcpendcnt laser heat source // Heat and Mass Transfer. 2001. N.37. C.333-342

Результаты сравнения решений граничной и неоднородной задач гиперболической теплопроводности приводятся в Таблице 1.

Таблица 1: Сравнение задач нагрева среды на границе и с тонком слое вблизи границы

Вид воздействия Т на фронте Т за точкой излома

Ч,=о = 5(г - 0) оо 0

4) = 106(г - 0)ехр(—/За') конечна т^О

/ф=0 = Я(4 - 0) конечна 0

<?(М) =/0Я(4 - 0)охр(-/?з) конечна

Рассмотрим процесс волнового переноса тепла при однородных граничных условиях первого рода на температуру и однородных начальных условиях:

П=о = °; = т1=о = ^ = 0-, (4)

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

Решение задачи (1), (2), (3), (4), полученное методом функций Грина имеет вид:

°° п(е"-(-1)")8тН

5П (п2тг2 + /Х272) Бп соэ

Г1 , 1 Г1 1

~Ог) ^ 2 + /Л БШ lJ'nt 2

(5)

где 1? - ТСУ//р; с = 8у/рС„/А; /х = Д; 7 = (}^\/{рСу)] 5„ =

у/4п27г2г — ц2\ Н(£) - функция Хевисайда. Сходимость данного ряда доказана по признаку Дирихле в приложении к диссертационной работе. Распределение температуры (5) показано на рис. 3. При аппроксимации временного профиля гладкой функцией, например, распределением Гаусса, скачки

па рис. 3 сглаживаются и уменьшаются пиковые значения температуры.

При граничных условиях первого рода, в отличии от условий теплоизоляции, можно наблюдать области, где температура опускается ниже начальной (рис. 3,а). Время их существования оценивается величиной порядка постоянной релаксации теплового потока. Результаты сравнения задач гиперболической теплопроводности, рассмотренных в работе, при разных граничных условиях и длительности воздействия представлены в Таблицах 2, 3. Величина Т, отложенная на графиках по вертикальной оси, не является абсолютным значением температуры и равна отклонению от некоторого начального значения Ту. Предполагается, что значения температуры То далеки от абсолютного нуля. В работе установлено, что области охлаждения

4Морс Ф.М., Фстбах Г. Методы теоретической физики Т.2. М.: ИИЛ, 1960. 886 с.

Рис. 3: Профили распределения температуры Т в слое в зависимости от координаты 5, построенные при нескольких значениях Ц а — случай г > т„, б — случай т < г*. Штриховой линией показаны огибающие.

Таблица 2: Сравнение задач теплопроводности. Короткое объемное воздействие: (¡(я, 4) = 1В6Ц - 0)ехр(—/Зз)____

Теплопроводность ГУ Т > 0 Экстремум в точке излома

Гиперболическая Л|»=о,1 = 0 всегда т < т.

Гиперболическая т |а=о,, = 0 т < т» есть

Классическая Л|5=о ,1 = 0 всегда нет

Классическая Г|8=0,1 = 0 всегда нет

Таблица 3: Сравнение задач теплопроводности. Длительное объемное воздействие: г/(.ч, ¿) = /0Я(« - 0)ехр(—/Зз)_

Теплопроводность ГУ Т > 0 Экстремум в точке излома

Гиперболическая Л|а=о,1 = 0 всегда есть

Гиперболическая Г|5=0,; = 0 всегда нет

Классическая 5=0,1 = 0 всегда нет

Классическая Т\3=0,1 = 0 всегда нет

могут обнаруживаться или не обнаруживаться в зависимости от соотношения между постоянной релаксации теплового потока среды т. и величиной т* = 1/(472), характеризующей глубину проникновения лазера в материал. Если г > т», тогда Т принимает отрицательные значения (рис. 3,а), если т < т„, то Т принимает только положительные значения вблизи облучаемой границы слоя, как показано на рис. 3,6.

Величина т„ зависит от частоты лазерного излучения и в пределах от ЮГГц до 1ТГц изменяется в диапазоне от Ю-9 до 10~ис для типичных металлов (алюминий, золото, медь, серебро). Рассмотрим специальные случаи решения задач гиперболической теплопроводности, которые могут быть применены для сравнения с результатами экспериментального измерения температуры. На рис. 4,а приведено решение (5) при ¡3 = 0. В этом случае мощность лазерного излучения не убывает с увеличением расстояния от границы.

1-- Т

2т

4т

6т

Рис. 4: График а — температурные профили при разных значениях времени порядка постоянной релаксации т в равномерно нагретом слое, остывающем при постоянной температуре на границе. График б — температура на границе теплоизолированного слоя, облучаемого коротким лазерным импульсом. Штриховой линией показано классическое решение при т = 0.

Примером может служить тонкий теплопроводящпй стержень, равномерно по всей длине нагретый лазерным импульсом. Так как на границах поддерживается постоянная температура, стсржепь постепенно остывает. При малых временах (графики на рис. 4,а вверху) заметны точки излома. На рис. 4,б представлен график зависимости Т(£) на границе слоя .ч = /, где проще всего проводить измерения. Наличие минимумов и максимумов па графике гиперболической теплопроводности (показаны сплошной линией) обусловлено перемещением тепловой волны внутри слоя.

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

Рассмотрим одномерную задачу о распространении объемных термо-уиругих воли в направлении координаты

Где Т — отклонение от начального значения температуры, р — плотность материала, е„ — удельная теплоемкость при постоянном объеме, Л — коэффициент теплопроводности, т — время релаксации теплового потока, а — объемный коэффициент теплового расширения, К — изотермический модуль

(6)

(К + 46/3) е" - аКТ" = рё

(7)

объемного сжатия, То — температура, при которой производились измерения констант, е — след тензора деформаций, G — модуль сдвига, точкой обозначается производная по времени, штрихом обозначена производная по координате. Для построения дисперсионных соотношений будем искать решение системы уравнений (6), (7) в виде экспонент, затухающих по координате:

ф, t) = Bekse~iwt\ T(s, t) = Dekse~iut; k = -7 + iö\ (8)

где 7 — коэффициент затухания, 5 — волновое число, со — частота. Представим дисперсионные соотношения в параметрическом виде:

V J у/х2 + 4 у2(х) — х

где использованы следующие обозначения:

ш(х) = Л3(1 - Л4 + x)/{Al{ 1 + А + 2х)2у(х)); (10)

. . 1 {1-А4 + х)(ААа-{Х + 1){Х + А)) У = 2 у-Л + Л4 + а:-; (П)

Л = Л2/Л1; А\ = р/(К + 4С/3); Л2 = гЛ3;

А =К I а2К*Т° \ ~А=_1_ (12)

3 Л + К + 4С/3) ' 4 1 + рс„(К + 4С/3)/(а2К2То)'

На основании (9) и (10) получим выражение фазовой скорости Cf.

С^х) = ш(х)/5{х) = у/2/у!Ах (у~х2 + Ау2{х) - ж) (13)

Групповая скорость Сд находится по определению:

(14)

В зависимости от того, выполняется ли условие Л < 1 — 2Л4 или Л > 1 — 2Л4, тепловые и акустические ветви дисперсионных соотношений демонстрируют различное поведение. Данные условия могут быть сформулированы относительно постоянной релаксации теплового потока т. В первом варианте получим: т < то, во втором варианте: т > тц, где значение релаксации го выражается через термомеханические параметры среды:

Ар (рву (К + 4С/3) — а2К2Т0) Т° (ре» (К + 46/3) + а2К2Т0)2 1 '

Значение то прямо пропорционально коэффициенту теплопроводности А и для реальных материалов лежит в пределах оценок величины релаксации

теплового потока, встречающихся в литературе. Для классической термоупругости существуют аналогичные условия: Л4 < 1/2 и А4 > 1/2, их можно сформулировать относительно коэффициента теплового расширения а: а < а0 и а > а0, где а0 находится по формуле: а0 = у/рсУ (К + 4С/3) /Т0/К. Для реальных веществ значения коэффициента теплового расширения приблизительно в 1.5-2 раза меньше ао для твердых веществ и жидкостей, и на треть меньше для газов.

Приведем интервалы значений х, при которых исследуемые функции (9)-(14) будут вещественными и положительными. При т < ть интервал (.т2, хл) соответствует акустическим ветвям, а интервал а^) соответствует тепловым ветвям; при г > ть интервал (ж2, х,4) соответствует тепловым ветвям, а интервал х^) соответствует акустическим ветвям. Значения границ интервалов зависят от термомеханичсских параметров:

XX = (-1 - А + у/1 - 2А + А2 + 4ЛЛ4) /2; х3 =-1 + А4 (16)

х2= (-1-А - л/1-2 Л + Л2 + 4 ЛЛ4) /2; ж4 = -Л - Л4 (17)

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

График у (со) имеет следующие горизонтальные асимптоты (рис. 5,а,в):

ъ =__• 72 =_+ (18)

2ау/2(Ах + А2 - а)' 2а^/2{Ах + Л2 + а)'

где константы а и Ь имеют вид: а = у/{А\ — Л2)2 + 4Л! Л2Л4 и Ь = Л2 + Лх(2Л4 — 1). В классической термоупругости возможна только одна асимптота: 72 = ЛзЛ4/(21/Л7), которая соответствует тепловой ветви при Л4 > 1/2 и акустической при Л4 < 1/2.

График ¿(а;) (рис. 5,б,г) имеет две наклонные асимптоты, выходящие из начала координат, тангенсы углов наклона которых определяются формулами:

= л/{А1+А2-а)/2; Ьёср2 = л/{А1+А2 + а)/2; (19)

где при т < то величина (рх равна углу наклона асимптоты тепловых ветвей, аср2 — углу наклона асимптоты акустических ветвей. При г > ть — наоборот. В классической термоупругости возможна только одна асимптота, тангенс угла наклона которой равен у/Ах это значение соответствует тепловой ветви при Л4 > 1/2 и акустической при Л4 < 1/2.

Асимптоты фазовых (13) и групповых (14) скоростей совпадают (рис.

С):

= у/г/у/Ах +А2-а; С2 = у/2/у/Ах + А2 +а] (20)

Рис. 5: Графики а, в — зависимость коэффициента затухания от частоты; б, г — зависимость волнового числа от частоты; графики а, б построены при г < ть, графики в,г построены при г > ть; кривые 1 - тепловые ветви, 2 - акустические ветви.

Рис. 6: Фазовые (тонкие линиями) и групповые (показаны толстыми линиями) скорости в зависимости от частоты; графики а построены при т < То, графики б - при г > то; тепловые ветви выходят из начала координат.

где при т < т0 величина С\ равна асимптотическому значению тепловых ветвей фазовых и групповых скоростей, С2 — асимптотическому значению акустических ветвей. При т > То — наоборот.

Рассмотрим частоты, соответствующие экстремальным значениям групповых скоростей. Обозначим их (для классической термоупругости) и И, (для гиперболической термоупругости). При т < То экстремумами групповой скорости являются минимум акустической ветви и максимум тепловой ветви, при т > то - минимум тепловой ветви и максимум акустической (рис. 6). При некотором соотношении термомеханических параметров материала, экстремумы могут не наблюдаться; На частоте Л, и при т = То гра-

Рис. 7: а - графики зависимости волнового числа от частоты, б — графики зависимости коэффициента затухания от частоты, построенные в окрестности ÍÍ* при г = го; 1, 3 - тепловые ветви, 2, 4 ~ акустические ветви.

фики 1 и 2 иа рис. 7,а пересекаются; графики 3 и 4 па рис. 7,6 подходят друг к другу, не имея общих точек. Для параметров меди, взятой в качестве примера типичного металла, при постоянной релаксации, равной значению То ~ 5.3 х 10~12с, вычисленному по формуле (15), частота Í!, составляет примерно 700ГГц. На данной частоте длина акустических волн на два порядка больше (2тг/5 те 3.73 х 10-8м) постоянной решетки меди (3.61 х 1О~10м), что не противоречит возможности применения континуальной модели в диапазоне частот близких к ílt. Величина то имеет тот же порядок, что и оценка релаксации теплового потока из фононной теории5: трн и X/(CV(K + 4G/3)). Для тсрмомеханичсских параметров меди оценка постоянной релаксации, вычисленная по данной формуле равна г,,/, ~ 5.6 х 10_12с.

Рассмотрим задачу о распространении продольных термоупругих волн в направлении координаты s. Будем искать установившееся решение системы дифференциальных уравнений (6), (7) для полупространства в виде экспонент, затухающих по координате:

u(s,t) =Cekse~iut-, T(s,t) = Dekse~iiJt-, к = -7 + ió; (21)

Волны распространяются в направлении, перпендикулярном границе полупространства. Граничные условия имеют вид:

Г|5=0 = 0; u|s=0 = w0sin(wí); (22)

Установившееся решение задачи (6), (7), (22) имеет вид:

u(s, t) = UQe~lhS (ai cos (wt — <5/,s) + аг sin — $hs)) +

uoe~7aS (a3 cos (cut — óas) + a4 sin (uit — óas)) (23)

T(s, t) = uüe~lhS (bi eos (ut — 6hs) + í>2 sin (cut — S^s)) —

uoe"laS (bi eos (ujt — 5as) + 62 sin (tot — 5as)) (24)

5Chcstcr, M. Sccond sound in solids // Physical Reviews. 1963. N.131. C.2013-2015

Вещественные константы зависят только от параметров среды и частоты колебаний:

а4 = "1т(л7^); 61 = ее (л7^) ; б2 = 1т(л7^);

д __Как _

1,2 " (4/3в + К) к^ + рш2'

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

Рис. 8: На графике а толстыми линиями показаны амплитуды перемещений в гиперболической термоупругости, тонкими линиями показаны амплитуды перемещений в классической термоупругости, на графике б амплитуды температуры в гиперболической термоупругости показаны топкими линиями, а амплитуды в классической термоупругости показаны толстыми линиями.

С ростом частоты усиливается затухание акустических волн, при этом максимальное значение амплитуд тепловых волн увеличивается, а расстояние, па которое опи проникают в среду, уменьшается.

При частотах меньших Г2», основной вклад в значения амплитуд (23) и (24) дают квазиупругие составляющие (вторые слагаемые в выражениях для и и Т). Следовательно, поведение решения описывается в большей степени акустической ветвыо дисперсионных соотношений. При частотах, больших Я*, вклад квазитепловых составляющих становится существенным и па графиках решения становится заметна модуляция амплитуд термоупругих волн (рис. 8,а,б).

На графиках (рис. 8), полученных с использованием гиперболической теории термоупругости (толстые линии), видно, что акустические и тепловые волны затухают за небольшое число периодов. В классической теории тсрмоупругости (тонкие линии), затухание практически не заметно.

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

Рассмотрим задачу о нагревании бесконечного слоя коротким лазерным импульсом в связанной и несвязанной постановке. Длительность реального импульса, который можно считать коротким в термоупругости LS типа, определяется как величина, обратная частоте, на которой дисперсионные кривые имеют значения близкие к асимптотическим. Например для алюминия при т = 10_11с акустические и тепловые ветви коэффициента затухания и фазовых скоростей выходят па асимптотические значения при частоте порядка 10/т = 1012Гц что соответствует длительности импульса г/10 = 10"12с. Неоднородные уравнения связанной гиперболической термоупругости LS типа относительно температуры и напряжений с учетом распределенных в объеме источников тепла выглядят так:

о" -о- D{t = 0 (25)

Т" - (à + D3à) - (т + D3T) = -D2e->s (â(D4t - 0) + £>3 j(D4i - 0)) (26) где введены обозначения:

s = stX^3p(4G + 3K)/d; t = 3i*A p/d; a = atd/(3/3KT0);

D1 = 3f32K2T0/d] D2 = 3I0Xp2(AG+3K)/d2- D3 = rd/(3Ap); D4 = 3Ap/d;

d = AGpCv + 3К (pCv + /32KT0) ;

Здесь s» — координата, tt — время, <т„ — напряжения на площадке, перпендикулярной оси st.

Для получения решения задачи термоупругости в несвязанной постановке, распределение температуры, найденное из задачи теплопроводности, подставляется в уравнение движения (25), которое затем решается совместно с начальными и граничными условиями как неоднородное волновое уравнение.

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

В ходе исследования установлено, что учет постоянной релаксации теплового потока вносит наиболее заметный вклад в задачу с однородными граничными условиями первого рода на температуру и напряжения. Временные профили температуры, напряжений и перемещений, соответствующие данной задаче, показаны толстыми линиями на рис. 9,10,а,б,в. Тонкими линиями показаны классические решения при г = 0 (графики классического решения быстро затухают, поэтому для них выбран более крупный масштаб по оси ординат).

Экстремумы, наблюдаемые на графиках, показанных толстыми линиями, соответствуют квазитепловым и квазиакустичсским составляющим термоупругой волны и перемещаются с разными скоростями слева направо в сторону от облучаемой лазерным импульсом границы. Скорость квазитепловых составляющих в несвязанной термоупругости (локальные экстремумы на рис. 9,5,0 слева и на рис. 10,б,в справа) равна скорости тепла в среде: Си = у*л/(рС„т). Скорость квазиакустичсских составляющих в несвязанной термоупругости (локальные экстремумы на рис. 9,б,в справа и на рис. 10,5,6 слева) находится по формуле для скорости упругих волн в среде: с0 = л/р/(К + 4(3/3). Коэффициент затухания квазитспловых составляющих равен = л/рС„/(Хт)/2. Квазиакустические составляющие термоупругой волны в несвязанной постановке задачи не затухают.

Рассмотрим задачу термоупругости для бесконечного слоя, нагреваемого коротким лазерным импульсом по закону Бугера в связанной постановке. Граничные и начальные условия на напряжения и температуру однородны:

Решение задачи (25), (26), (27), полученное по методу Гринберга, принимает вид:

Т'\а=о = Т'\3=1 = 0; ег|4=0 = а\5=г = 0;

Г'|8=0 = Т'\я=1 ~ 0; ст'|5=о = сг'|а=г = 0; Т|»=0 = Т\8=1 - 0; <т|4.=о = <т\„=1 = 0;

Т|,=о = Ти = 0; <т'|5=0 = а%=1 = 0;

(27)

Выражения для трансформант Тп{1) и сгп{£) имеют вид: Т =

П -

Р2СпЩ^ уч (-1)* (Д? + Л„) (Д3д< + 1) ехр(Ш0

ААгСпЩ*) А (-1)*Д2 (РзЪ + 1) ехр(Щ)

(1 - ох) ад и ш Пи № - Ъ)

где для краткости записи введены обозначения:

Сп = ъ1п (1 - (—1)пе_7г) /(72/2 + тг2тг2); Ап = тт2п2/12 Дч, /?,1 — корни уравнения четвертой степени: (Г>1 - 1) Я3Д4 + (А - 1) л3 - (£>з + 1) АПД2 — ХпЯ — А2 = О

(30)

Рис. 9: Профили распределения температуры Т, напряжения а и перемещения и в слое в зависимости от координаты в в момент времени порядка г при т > Тд. Графики а, б, в соответствуют задачам термоупругости в несвязанной постановке, графики г, д, е — задачам в связанной постановке. Тонкими линиями показаны графики классической термоупругости, толстыми линиями показаны графики гиперболической термоупругости.

На рис. 9,10,г,д,е представлены решения (28) связанной задачи термоупругости. В результате учета связанности в постановке задачи, на графике температуры (рис. 9,10,г) возникает второй минимум, соответствующий ква-зиакустичсской составляющей термоупругой волны. В зависимости от параметров задачи в области отрицательных значений Т могут находиться два минимума, один минимум или ни одного.

Скорости квазитепловых и квазиакустических составляющих термоупругой волны в связанной задаче термоупругости находятся по формулам (20), полученным в первой главе:

Сг = у/2/у/А1+А2-а; С2 = ^/х/лТТлТ+а;

Всегда выполняется соотношение: С\ > С2. Значения скорости упругих волн са и скорости тепла с^ лежат в интервале (С2, С1). При коэффициенте теплового расширения равном нулю а = 0, выполняются соотношения: С\ = са и С2 = с/,, если г > То (см. формулу (15)), и если г < то, то наоборот: С\ = Сд и С2 = са. Расхождение скорости упругих волн са со скоростью квазиунругой составляющей термоупругой волны, а также скорости тепловых волн с/, со скоростью квазитепловой составляющей термоупругой волны не превышает 2% (в расчетах использованы термомеханические параметры алюминия, взятого в качестве примера типичного металла; значение постоянной релаксации принято равным т = 10_11с).

Если т > Го, то скорость распространения квазитепловой составляющей термоупругой волны, равная С2, будет меньше скорости распространения квазиакустичсской составляющей С\ (рис. 9,г,д,е). Если т < то, то скорость квазитепловой составляющей термоупругой волны С\ будет больше скорости квазиакустической составляющей термоупругой волны С"2 (рис. 10,г,д,е).

В случае классической термоупругости, если т = 0, то значение С\ обращается в бесконечность, а значение С2 становится равным скорости упругих волн в среде са-

Коэффициенты затухания квазитепловых и квазиакустичсских составляющих термоупругой волны в связанной задаче термоупругости находятся по формулам (18):

= Аз (а - Ь) = Л3(а + Ь)

71 2ау/2(А1+А2 - а)' 72 2ах/2(А1 + А2 + а)

При г < т0 коэффициент 71 является характеристикой затухания квазитермических составляющих термоупругой волны, а 72 — квазиакустических; при т > То наоборот. Если положить коэффициент теплового расширения а равным нулю, то будут выполняться следующие соотношения: 71 = 7д и 72 = 0 при т < То или 71 = 0 и 72 = -ун при т > То-

Графики температуры, построенные с использованием связанной гиперболической термоупругости (рис. 9,10,г) отличаются от графиков, построенных с использованием гиперболической теплопроводности (рис. 9,10, а) скоростями распространения волн, характеристиками их затухания, минимальным и максимальным значением Т, количеством температурных пиков. Все перечисленные отличия вносит учет механических слагаемых в уравнение теплопроводности.

г д е

Рис. 10: Профили распределения температуры Т, напряжения а и перемещения и в слое в зависимости от координаты я в момент времени порядка т при т < го. Графики а, б, в соответствуют задачам термоупругости в несвязанной постановке, графики г, д, е — задачам в связанной постановке. Тонкими линиями показаны графики классической тсрмоупругости, толстыми линиями показаны графики гиперболической термоупругости.

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

В заключении приведен список основных результатов и сформулированы основные выводы диссертационной работы. В приложении приведено доказательство сходимости ряда (5).

1. Получены выражения зависимостей волнового числа, коэффициента затухания, фазовых и групповых скоростей от частоты в параметрической форме. Данные соотношения оказываются удобными для их аналитического исследования: нахождения асимптот и особых точек графиков функций.

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

3. Исследовано поведение дисперсионных соотношений в широком диапазоне значений термомеханических параметров среды. Обнаружены два

Основные результаты работы

варианта поведения дисперсионных кривых в зависимости от времени релаксации теплового потока т. Установлено, что существует такое Тр, что при т < Tq и при т > То характер дисперсионных соотношений качественно отличается. Найдено выражение для То через термомеханичкес-кие параметры среды.

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

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

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

Публикации автора по теме диссертации в изданиях, включенных в Перечень ВАК

1. Бабенков, М.Б. Анализ дисперсионных соотношений связанной задачи термоупругости с учетом релаксации теплового потока / М.Б. Бабенков // Прикл. Мех. Тех. Физ.—2011.—Т.52.— №6—С. 112-121.

2. Бабенков, М.Б. Анализ распространения гармонических возмущений в термоупругой среде с релаксацией теплового потока / М.Б. Бабенков // Прикл. Мех. Тех. Физ.-2013.-Т.54,- №2.-С. 12(5-137.

3. Babenkov, M.B. Analysis of the wave propagation processes in heat transfer problem of the hyperbolic type/ Babenkov, M.B., Ivanova, E.A. // Continuum Mech. Thermodyn.-2013.-doi:10.1007/s00161-013-0315-8

Подписано в печать 12.11.2013 Формат 60x90/16 Бумага офсетная. Усл. печ. л. 1,5 Тираж 100 экз. Заказ 560

Отпечатано в типографии «Адмирал» 199178, Санкт-Петербург, В.О., 7-я линия, д. 84 А

Текст научной работы диссертации и автореферата по механике, кандидата физико-математических наук, Бабенков, Михаил Борисович, Санкт-Петербург

Федеральное Государственное Бюджетное Учреждение Науки Институт Проблем Машиноведения Российской Академии Наук

Бабенков Михаил Борисович

Распространение термоупругих волн в среде с учетом релаксации теплового потока

Специальность 01.02.04 — «Механика деформируемого твердого тела»

Диссертация на соискание учёной степени кандидата физико-математических наук

Научный руководитель: доктор физико-математических наук

Иванова Е.А.

Санкт-Петербург - 2013

Содержание

Введение .................................. 4

1 Анализ волновых процессов в задачах теплопроводности гиперболического типа .........................12

1.1 Уравнение теплопроводности гиперболического типа ....... 12

1.2 Волновые процессы в полупространстве, возникающие при тепловом воздействии на границе...................... 15

1.2.1 Тепловой поток на границе зависит от времени как 6-функция Дирака........................ 15

1.2.2 Тепловой поток на границе зависит от времени как функция Хевисайда......................... 17

1.3 Волновые процессы в слое, возникающие под воздействием внутренних источников тепла....................... 19

1.3.1 Мощность внутренних источников зависит от времени как ¿-функция Дирака....................... 20

1.3.2 Мощность внутренних источников зависит от времени как функция Хевисайда...................... 32

1.4 Заключение............................... 38

2 Анализ дисперсионных соотношений в связанной задачи термоупругости гиперболического типа................42

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

2.2 Анализ дисперсионных соотношений..................................43

2.3 Фазовая и групповая скорость в термоупругой среде ..............54

2.4 Распространение плоских гармонических волн в термоупругом

полупространстве........................................................60

2.4.1 Анализ графиков термических и акустических волн .... 63

2.5 Заключение..............................................................66

3 Анализ волновых процессов в задачах термоупругости гиперболического типа ...........................68

3.1 Уравнение движения несвязанной термоупругости..................68

3.2 Перемещения в слое, возникающие под воздействием внутренних источников тепла в несвязанной термоупругости....................70

3.2.1 Границы слоя закреплены......................................70

3.2.2 Границы слоя свободны от нагрузок..........................73

3.3 Температура и перемещения в слое, возникающие под воздействием короткого лазерного импульса в связанной задаче термоупругости ......................................................................76

3.4 Заключение..............................................................81

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

Литература.................................85

Список рисунков .............................99

Список таблиц...............................100

А Приложение ..............................101

Введение

Актуальность .темы исследования. Исследование волнового переноса тепла актуально для многих развивающихся технологий. Внутренние источники наноразмерного масштабного уровня (например, размер современного транзистора составляет всего несколько нанометров) вызывают более интенсивный рост температуры, чем предсказывает классическая теория, что повышает требования к будущим системам охлаждения. Температурные эффекты вносят значительные изменения в механические свойства нанопластин, так как чем тоньше пластина, тем она чувствительнее к изменениям температуры. Исследования процессов волнового переноса тепла могут быть полезны для изучения термомеханических характеристик объектов микро- и нано- масштабного уровня: тонких пластин и стержней, используемых в микро- и нано- электромеханических устройствах (MEMS и NEMS).

ют гиперболическое уравнение теплопроводности на основе обобщенного закона Фурье.

Обзор литературы. Обобщенный закон Фурье был предложен в работах авторов: С. Cattaneo [1], P. Vernotte [2], A.B. Лыкова [3]. В отличии от классического закона, он учитывает инерционность процессов теплопереноса, которая характеризуется постоянной релаксации теплового потока.

Экспериментальным определением величин релаксации теплового потока занимались: R.H. Matsunaga, I. Santos [4], W. Kaminski [5], К. Mitra et. al. [6], A. Grabmann, F. Peters [7], H. Herwig, K. Beckert [8], W. Roetzela et. al. [9]

Известны другие, более сложные модели теплопроводности, имеющие в линейном приближении конечную скорость распространения тепловых возмущений [10-12]: модели с тепловой памятью [13], учитывающие историю нагревания тела; модель [23], учитывающая времена запаздывания теплового потока и градиента температуры и т.д. Обзор моделей теплопроводности представлен в [10].

Связанная термоупругость, в которой вместо классического закона теплопроводности Фурье используется обобщенный закон Фурье (закон Каттанео-Вернотте) была впервые предложена в работе Е.Б. Попова [15], позже была предложена более общая модель связанной гиперболической термоупругости: Н. Lord, A. Shulman [16]. Термоупругость Лорда-Шульмана (LS) является предметом изучения в данной работе. На сегодняшний день существует ряд других теорий термоупругости, учитывающих в линейном приближении конечность скорости распространения тепла: термоупругость Грина-Линдси (GL), Грина-Нагди второго типа (GNII) и третьего типа (GNIII), среды с тепловой памятью. Автором D.S. Chandrasekharaiah выполнен широкий обзор по данной тематике [11,12]. В статье J. Ignaczak [17] подробно обсуждается применимость гиперболической модели к описанию процессов теплопереноса. Обзор экспериментальных работ приведен в книге А.Г. Шашков, В.А. Бубнов, С.Ю. Яновский [20].

Приведем список авторов, занимавшихся исследованием плоских стационарных гармонических волн, распространяющихся полупространстве (LS мо-

дель термоупругости): A.H. Nayfeh, S. Nemat-Nasser [19], Ю.К. Энгельбрехт [85], P. Puri [25], И.М. Штер [26], Ф.В. Семерак [27], Ц. Иванов, Ю.К. Энгельбрехт [28], Р.Х. Швец, A.A. Лопатьев [29]. В перечисленных работах получены и проанализированы асимптотические выражения для фазовых скоростей и коэффициентов затухания квазиупругой и квазитепловой волн при больших и малых значениях частот. Установлено, что фазовые скорости и коэффициенты затухания с ростом частоты выходят на асимптоту, в отличии от классической термоупругости, где они неограниченно возрастают. V.K. Agarwal [30] обобщает результаты P. Puri [25] на термоупругость GL типа. Результаты некоторых из перечисленных исследований приводятся в книгах: А.Г. Шашков, В.А. Бубнов, С.Ю. Яновский [20], А.Д. Коваленко [31], Я.С. Подстригач, Ю.М. Коляно [32]. Исследованием дисперсионных соотношений в рамках классической связанной термоупругости (СТЕ) занимались Р. Chadwick, I.N. Sneddon [33].

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

Приведем список авторов работ в которых на основании гиперболического уравнения теплопроводности (LS типа) исследуется распределение температуры в полупрозрачной среде, находящейся под воздействием лазерного импульса. Лазерный импульс моделируется источниками тепла распределенными в объеме, рассматривается одномерный случай: K.L. Baumeister, T.D. Hamill [34] В. Vick, M.N. Özisik [35], D.W. Tang, N. Araki [36,37]. В следующих работах для описания пространственного распределения источников тепла используется закон Бугера: M. Lewandowska [79], D. Zhang et al. [38], T.T. Lam, E. Fong [39], H. Al-Qahtani, B.S. Yilbas [40], B.S. Yilbas, A.Y. Al-Dweik [41]. Во всех вышеперечисленных работах использованы граничные условия термоизоляции.

В книге А.Г. Шашков, В.А. Бубнов, С.Ю. Яновский [20] приводится системно-структурный анализ L-изображений температуры в задачах теплопроводности как с тепловым воздействием на границе, так и с распределенными в объеме источниками тепла при различных граничных условиях: первого рода, термоизоляции, обмена теплом с окружающей средой по закону Ньютона,

теплового контакта. В книге В.А. Ковалева, Ю.Н. Радаева [43] представлено исследование гиперболической термоупругости Грина-Нагди (GNI/CTE, GNII, GNIII) с позиций теории поля. В книге Д. Жоу, X. Касас-Баскес, Дж. Лебон [10] приведено общее решение неоднородного уравнения Каттанео-Вернотте, сделан обзор работ, посвященных задачам гиперболической теплопроводности и теоретическим основам неклассического теплопереноса. В книге В.А. Кудинов, И.В. Кудинов [42] подробно исследовано решение задачи гиперболической теплопроводности для бесконечного слоя при симметричных граничных условиях первого рода, построены графики решения. В книге Ф.М. Морс, Г. Фешбах [44] рассматривается нестационарное нагревание пластины внутренними распределенными источниками тепла с однородными граничными условиями первого рода в рамках классической модели теплопроводности. Как отмечают авторы, решение данной задачи находит применение для описания нестационарного распределения температуры в металлической полосе, нагреваемой переменным током и погруженной в среду с большой теплопроводностью.

Следующие авторы использовали гиперболическую термоупругость LS типа для исследования температурных напряжений, возникающих в полупространстве вследствие внезапного нагревания границы (одномерный случай): F.R. Norwood, W.E. Warren [45], Г.А. Кильчинская [46], А.Н. Nayfeh, S. Nemat-Nasser [47], M. Baila [48], P.M. Jordan, P. Puri [49], I.A. Abdallah [50], N. Sarkar [51], E.F. Henain et al. [52].

Ниже приведен список авторов, рассматривавших задачи термоупругости LS типа в полупространстве с распределенными в объеме источниками (одномерный случай). Для описания пространственного распределения источников тепла в каждой работе используется закон Бугера. J.C. Strikwerda, A.M. Scott [53] используют в своем исследовании граничные условия термоизоляции, граница полупространства свободна от механических нагрузок. Решение для напряжений и температуры получено в приближении для малых значений постоянной релаксации теплового потока, построена асимптотика для больших и малых времен. X. Wang, X. Xu [54,55] в отличии от предыдущих исследователей, получают решение для широкого диапазона значений времени и постоян-

ной релаксации теплового потока. Н.М. Youssef, A.S. Al-Felali [56] рассматривая ту же задачу, ставят граничные условия первого рода на температуру, граница полупространства в рассматриваемой ими задаче закреплена.

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

Экспериментальным исследованием термоупругих напряжений при неравновесных процессах теплообмена занимались: Н.В. Вовненко, Б.А. Зимин, Ю.В. Судьенков [57-60], K.V. Poletkin, G.G. Gurzadyan, et al. [61], O.B. Wright et al. [62,63] и др.

Целью данной работы является изучение поведения термоупругой среды, описываемой уравнениями Лорда-Шульмана, под воздействием периодических и импульсных возмущений. Задачи исследования:

1. Исследовать распространение термоупругих волн в широком диапазоне значений постоянной релаксации теплового потока. Данная задача представляется актуальной, поскольку теоретическая оценка постоянной релаксации в металлах, согласно фононной теории, составляет несколько пикосекунд (10-12с). В то время как экспериментальные данные дают разброс результатов, отличающихся от теоретической оценки на несколько порядков: от 10_8с до 10~ис.

Научную новизну работы представляют следующие положения, выносимые на защиту:

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

5. Обнаружено, что если постоянная релаксации теплового потока меньше установленного значения, то квазитермические составляющие распро-

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

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

таты работы доложены на: XLI, XL и XXXVIII International Summer School Conference «Advanced Problems in Mechanics» (Россия, С.-Петербург, 2013, 2012 и 2010 гг.), International Conference «Days on Diffraction» (Россия, С.-Петербург, 2013 и 2012 гг.), VI Поляховские чтения (Россия, С.-Петербург, 2012), 2nd International Conference on Material Modelling (Франция, Париж, 2011), на Санкт-Петербургском Семинаре по Вычислительной и Теоретической Акустике Научного Совета РАН по Акустике (руководитель проф. Д.П. Коузов, С.Петербург, 2009, 2010, 2011 и 2012 гг.), на Городском семинаре по механике (руководитель чл.-корр. РАН, проф. Д.А. Индейцев, С.-Петербург, 2013 г.)

Публикации. Основные результаты по теме диссертации изложены в 10 печатных изданиях, 3 из которых опубликованы в журналах, рекомендованных ВАК [64-66], 7 — в сборниках тезисов докладов.

Структура и объем дис