ПРИКЛАДНАЯ МАТЕМАТИКА И МЕХАНИКА
Том 78. Вып. 2, 2014
УДК 539.319:536.21
© 2014 г. В. И. Горбачев
ИНТЕГРАЛЬНЫЕ ФОРМУЛЫ В СВЯЗАННОЙ ЗАДАЧЕ ТЕРМОУПРУГОСТИ НЕОДНОРОДНОГО ТЕЛА.
ПРИМЕНЕНИЕ В МЕХАНИКЕ КОМПОЗИТОВ
Рассматривается связанная нестационарная задача термоупругости для неоднородного тела, описываемая системой из четырех дифференциальных уравнений второго порядка в частных производных с переменными по координатам коэффициентами. Наряду с этой исходной задачей рассматривается такая же задача для однородного тела той же формы (сопутствующая задача). Получены интегральные формулы, позволяющие выразить перемещения и температуру в исходной задаче через перемещения и температуру в сопутствующей задаче. Интегральные формулы используются для представления решения исходной задачи в виде рядов по всевозможным производным от решения сопутствующей задачи. Записана система рекуррентных задач для коэффициентов этих рядов. Найдены выражения для коэффициентов сопутствующей задачи (эффективные коэффициенты) и приведены постановки специальных краевых задач, из решения которых находятся конкретные выражения для эффективных коэффициентов термоупругости. Доказана теорема о том, что эффективные коэффициенты удовлетворяют физико-механическим ограничениям, накладываемым на термоупругие константы реальных тел. Рассмотрен случай неоднородного по толщине слоя, и для него получены явные аналитические выражения всех эффективных коэффициентов термоупругости. Подробно рассмотрен случай периодической зависимости коэффициентов термоупругости от координат.
Основы термоупругости однородных изотропных тел и решения многих классических задач изложены в книгах [1—4]. Рассмотрены задачи термоупругости для тел с непрерывной неоднородностью и для кусочно-однородных тел [5]. Н.С. Бахваловым [6—8] был предложен асимптотический метод осреднения для уравнений с периодическими коэффициентами, позволяющий в некоторых случаях получать точные решения, пригодные и в случае непериодической зависимости коэффициентов от координат [9, 10]. Проводилось осреднение связанной нестационарной задачи термоупругости для регулярных композитов [11—14].
Интегральные представления впервые применялись [15—17] для осреднения статической и динамической задач неоднородной упругости.
1. Постановка связанной задачи термоупругости. Система связанных динамических уравнений линейной термоупругости для неоднородного анизотропного тела [1, 2, 4] состоит из уравнений движения
° I], ] + х = Р и г
(1.1)
соотношений Дюгамеля—Неймана
(1.2)
соотношений Коши для деформаций
1
2
&тп ^тпрдЩр, д п (ит, п + ип, т) (1.3)
Здесь ит — компоненты вектора перемещения точки тела, Етп — компоненты тензора малых деформаций, а у — компоненты тензора напряжений, Сутп — компоненты изотермического тензора модулей упругости, ву — компоненты изотермического тензора тепловых напряжений, 9 = Т(х, ?) — Т0 — разность температур текущего и начального состояний (отклонение от начального состояния), р — плотность материала, Х(х1, х2, х3, ?) — компоненты вектора объемных нагрузок.
Тепловая часть уравнений термоупругости состоит из уравнения баланса энтропии
ТБ = - #. + ъ (1.4)
закона теплопроводности Фурье для медленных термодинамических процессов
= -\ijTj = -\nfij (1.5) и из выражения для энтропии при |9/Т0| 1
Б = вцЦ + С | (1.6)
Т0
где q¡ — компоненты вектора плотности теплового потока, ^(хъ х2, х3, ^ — плотность объемных источников или стоков тепла, Т — текущая температура в точке тела, — плотность энтропии, коэффициенты Ху образуют симметричный тензор второго ранга, называемый тензором теплопроводности, через Се обозначена удельная теплоемкость при постоянной деформации.
В уравнения (1.1)—(1.6) входят пять типов термомеханических характеристик материала С щ, ву, Ху, р, Се. В общем случае анизотропии их 35 штук и все они предполагаются интегрируемыми функциями координат х1, х2, х3.
Уравнения (1.1)—(1.6) можно записать в виде связанной системы линейных уравнений термоупругости относительно компонент вектора перемещений и отклонения температуры от начального значения ([4] с. 76)
(СЦтпЩш, п - Ру&Ь - РЩ. + Х1 = 0
. ' (1.7)
- СЕд - ТоР,",3 + ъ = 0
Для выделения единственного решения системы уравнений (1.7) нужно добавить граничные и начальные условия для искомых величин и¡(х, ?), 9(х, например, граничные условия
х еЕ: щ(х, г) = и0(х, г), д(х, г) = (х, г) (1.8)
и начальные условия
г = Щ (х, г) = (х), и.(х, г) = gi(x), д(х, г) = к(х) (1.9)
Задачу (1.7)—(1.9) для неоднородного тела будем называть исходной задачей.
Сопутствующая задача. В случае постоянных коэффициентов термоупругости вместо уравнений (1.7) имеем
С\тпит, п/ - - Р°<Л + Х1 = 0 , - С° у - Т° р ° Ьи + Ъ = 0 (1.10)
где и и у = Т°(х, ?) — Т0 — перемещения и относительная разность температур в однородном теле при тех же входных данных, что и в исходной задаче, т.е. кроме одинаковых объемных сил и источников тепла в задаче для однородного тела совпадают также граничные и начальные условия
х еЕ: и(х, %) = и°(х, %), у(х, %) = &°(х, %) (1.11)
% = и(х, %) = ^(х), щ(х, %) = Е(х), у(х, %) = Н(х) (1.12)
Задачу для тела той же формы и с теми же входными данными, что и в исходной задаче, но из материала с постоянными механическими характеристиками будем называть сопутствующей задачей.
2. Вывод интегральных представлений для перемещений и температуры в неоднородном теле. Введем функции и(х, ?) и ©(х, ?), определенные в области V, занятой телом, непрерывные и обладающие необходимым количеством непрерывных производных по своим переменным. Кроме этого, потребуем, чтобы
и(х, %)|х = у= 0, и(х, %)I = ° = 0, Щх, %)I = ° = 0
®(х, %)\х = уе2 = °, 0(х, %)\% = ° = ° .
Эти функции потребуются для вывода интегральных формул, по которым решение исходной связанной задачи представляется через решение сопутствующей задачи.
Преобразование интегрального тождества, основанного на уравнениях движения. В первом уравнении системы (1.7) произведем замену ? на т, умножим на Щх, ? — т), проинтегрируем по объему тела, а также по т от 0 до Получим % %
Ц[С/тп(х)итг п(х, т—(х, % - т)dVxdт - Ц[Р/(х)д(х, т)]/Щх, % - т)dVxdт -
° V ° V
% %
- Цр(х)й(х, т) Щх, % - т)dVxdт + т) и(х, % - т)dVxdт = ° (2.2)
° V ° V
Преобразуем интегралы в этом тождестве с целью заменить исходные функции иI и 9 на функции ц и у сопутствующей задачи. Опуская подробности несложных, но
довольно громоздких выкладок, приведем результат % %
Ц[(С/тпи,/),п - ри](йт - От)dVxdT + ||Руи,/(& - у)dVxdт =
° V ° V
% %
= Ц[С/тп(х) - С°/тп] и- ищ ^I + Ц[р(х) - р°УЬ^^т -
° V ° V
%
- Ц[Р/(х) - Р°] и,-иdVxdт (2.3)
Преобразование интегрального тождества, основанного на уравнении теплопроводности. Во втором уравнении системы (1.7), точно так же как это было сделано в предыдущем случае, произведем замену t на т, умножим на ©(х, t — т), проинтегрируем по объему тела и по т от 0 до t. Получим после преобразований
г
), - СЕ0](3 - у)dVxdт + То Ц(Р,0),(щ - о,)dVxdт
о V о V
г
= - 4]0,у,^т + То ЦЩх) - Р,]0,^т +
о V о V
г
+
о V
Ц[СЕ(х) - Со]0уdVxdx
о V
Вычтем это равенство из равенства (2.3), получим
г
{{[(С,тпЦ 1 ) „ - рит - То(Рт,0),](ит - ит)dVxdт -
г
г
Ц|(Л,®,-), - Се® - ви, I (3 - у)dVxdт =
о V г
о V
= Ц[С,тп(х) - Со,тп] ииОт, ndVxdТ + Ц[р(х) - ро] и/О^Т ■
о V о V
г г
- 4)0,У^Т - Ц[Се(х) - С°]®ууdVxdт -
о V о V
г
- Ц[Р,(х) - Р,](и,,у + То® Ьи ,)dVxdт (2.4)
о V
Тензоры и функции Грина связанной задачи термоупругости. Функции Ц(х, 0 и ©(х, 0 удовлетворяют однородным начальным и граничным условиям (2.1), в остальном они произвольные. Ограничим теперь произвол этих функций, подчинив их уравнениям, левая часть которых выделена чертой снизу в формуле (2.4). Правую часть нужных в дальнейшем уравнений будем выбирать двумя разными способами. В первом случае в формуле (2.4) положим
(С,тпЧ,)п - Р &т - То(Рт,0) , = -5т*5(х - - Т)
(X,0,), - Се(0 - Р^,, = о .
Решением системы уравнений (2.5) будут функции и(к) (х , г - т) и ©(к)(х, г — т), где Ътк - дельта Кронекера [18], а 8(г — т) и 8(х — — дельта-функции Дирака [19], причем
5(х - £) = 5(х: - ^)5(х2 - ^)5(хз - ^)
В соответствии с условиями (2.1), решение уравнений (2.5) удовлетворяет однородным граничным и начальным условиям
ик(хЛ г)\хе2 = °, и?"(хХ %)\% = ° = °, ик)(хЛ %)\% = ° = ° (2б)
0(к) (хД, %)\х е2 = °, 0(к) (хД, %)= ° = ° .
Во втором случае положим (С/тп и,/),п - Р ит - ЗД/))/ = °
(^■0,)- - СЕ0 - ви- = -5(х - £,)5(% - т) .
Решением этой системы уравнений являются функции Щх, г — т) и ©(х, 1 — т), удовлетворяющие нулевым граничным и начальным условиям
Ц(хЛ %)|хе2 = °, и(х, %)\% = ° = °, и-(хЛ щ% = ° = ° (28)
0(х, Щ х е2 = °, 01% = ° = ° .
Функции ик) (х, г - т), ©(к)(х, г - т) и Щх, 1 - т), ©(х, г - т), по сути дела, являются тензорами Грина связанной краевой задачи термоупругости. Они удовлетворяют связанным уравнениям (2.5) и (2.7) и однородным граничным и начальным условиям. Используя равенство (2.4), с помощью уравнений (2.5), (2.7) можно сразу выписать представляющие интерес интегральные представления
йк(х, %) = и>к(х, %) -%
- Ц[С/тп(^) - С°/тп] и/(х, % - т) ^ п(£, т)dV(dт -
° V %
- Ц[р(^) - Р°] ик(х, % - т) ад, т)dVi.dт +
° V %
+ Ц^-Д) - 4]0 (к) (х, % - т) dVi;dт +
° V %
+ Ц[С(£) - С°]0(к)(х, % - т)у (£, т)dViddт +
° V %
+ ДОР/^) - в°][ и- (х, % - т)у(£, т) + Т° 0(к) (х, % - т) т)] dVd (2.9)
° V
$(х, %) = у(х, %) - ... (2.10)
Многоточием обозначены слагаемые, аналогичные соответствующим слагаемым в представлении (2.9) при замене и^ , ©(к) на Ц, ©.
Заметим, что при выводе интегральных формул (2.9) и (2.10) использовались граничные и начальные условия исходной задачи. Однако в этих формулах входные данные явно не присутствуют, поэтому они остаются справедливы и при других граничных и начальных условиях.
3. Представление решения нестационарной динамической задачи термоупругости в виде рядов по производным от решения сопутствующей задачи. Предположим далее, что функции и и у бесконечно дифференцируемы по всем переменным и раскладываются в ряды Тейлора по времени и координатам:
л (£,т) = £
——
Для дальнейшего прочтения статьи необходимо приобрести полный текст. Статьи высылаются в формате PDF на указанную при оплате почту. Время доставки составляет менее 10 минут. Стоимость одной статьи — 150 рублей.