Полный текст

В настоящее время существует большое количество работ, посвященных изучению механодиффузионных эффектов в сплошных средах. Экспериментальные исследования этой проблемы были начаты в 30-х годах ХХ века. Было показано, что в результате изгиба пластинки из сплава золота и меди градиент деформации инициировал процесс восходящей диффузии. Это привело к формированию градиента концентрации и, как следствие, перераспределению атомов растворенного вещества. В результате возник перенос вещества из областей сжатия в области растяжения. Результат научных исследований был опубликован в 1936 году [1]. Дальнейшие эксперименты позволили установить, что упругие деформации сжатия лишь немного уменьшают коэффициент самодиффузии в некоторых металлах [2, 3]. Это позволяет сделать вывод о том, что упругая деформация слабо влияет на диффузию и основные механодиффузионные эффекты наблюдаются только в области пластической деформации. Первые модели, описывающие указанные явления, появились в 60-х годах ХХ века. В настоящее время в рамках феноменологического подхода наибольшее развитие получили модели, основанные на теории гомогенных смесей, где компоненты, их составляющие, перемешаны и взаимодействуют на молекулярном или атомарном уровнях (поверхности раздела частей отсутствуют). К ним относят твердые тела с примесями, сплавы, твердые растворы (однофазные кристаллические или аморфные твердые вещества переменного состава из двух или более компонентов). Такой подход обладает рядом преимуществ, основным из которых является то, что физические свойства гомогенных смесей во всех частях одинаковы или меняются непрерывно, без скачков, что позволяет эффективно применять аппарат дифференциального исчисления. Кроме того, с точки зрения термодинамики такие смеси представляют собой однородные термодинамические системы, каждой точке которых в условиях равновесия соответствуют одинаковые значения давления, температуры и концентрации, что тоже существенно упрощает процесс моделирования механодиффузионных явлений. Из наиболее свежих публикаций можно назвать [4-20]. Здесь рассматриваются различные постановки задач термомеханодиффузии, в том числе с конечной скоростью распространения тепловых и диффузионных потоков, а также предлагаются методы решений соответствующих стационарных и нестационарных задач. Обзор публикаций по данной тематике показывает, что проблема анализа взаимодействия механического и диффузионных полей является актуальной. Наибольший практический интерес представляют задачи для балок, пластин и оболочек, которые составляют основу большинства конструкций. Среди немногочисленных публикаций по этой теме можно отметить статьи [21, 22]. В них исследуется влияние диффузионных процессов на несущую способность пологой трансверсально-изотропной оболочки. Контактное взаимодействие стержня с упругим полупространством рассматривается в работах [23, 24]. Публикации [25-27] посвящены исследованию механодиффузионных процессов в пластинах. Расчет сферических оболочек с учетом диффузии рассмотрен в [28]. Следует отметить, что все указанные задачи решаются в стационарной постановке. Постановки задач о нестационарных упругодиффузионных колебаниях балок и пластин и методы их решения в известных на сегодняшний день публикациях отсутствуют. В данной статье исследуются эффекты взаимодействия механического и диффузионного полей в пластине Кирхгофа-Лява. Математическая модель упругодиффузионных колебаний пластины получена на основе вариационных принципов и известных соотношений теории пластин, изложенных в работах [29-32]. 1. Постановка задачи Рассматривается задача о нестационарных упругодиффузионных колебаниях прямоугольной изотропной пластины Кирхгофа-Лява с учетом релаксации диффузионных потоков. Схема приложенных сил и изгибающих моментов, а также ориентация осей прямоугольной декартовой системы координат представлена на рис. 1. Рис. 1. Силы и моменты, действующие на пластину Fig. 1. Forces and moments acting upon the plate Здесь - распределённые по поверхности моменты; - распределённая по поверхности поперечная нагрузка; - распределённые по поверхности плотности объемных источников массопереноса. Для математической постановки задачи используется модель упругодиффузионных процессов в сплошных средах в прямоугольной декартовой системе координат, которая в случае однородной среды имеет вид [4-20, 33-35] (1) где и - компоненты тензора напряжений и вектора диффузионного потока, которые для идеального твердого раствора определяются следующим образом : (2) Здесь точки обозначают производные по времени. Все величины в (1) и (2) являются безразмерными. Для них приняты следующие обозначения: (3) где - время; - прямоугольные декартовы координаты; - компоненты вектора перемещений; - толщина пластины; и - длина и ширина пластины; - приращение концентрации q-й компоненты вещества в составе многокомпонентной среды; и - актуальная и начальная концентрации q-го вещества; и - упругие постоянные Ламе; - плотность; - коэффициенты, характеризующие объёмное изменение среды за счёт диффузии; - коэффициенты самодиффузии; - универсальная газовая постоянная; - температура среды; - молярная масса q-го вещества; - компоненты вектора плотности массовых сил; - плотность источников массопереноса; - время релаксации диффузионных потоков. Полагаем, что поперечные прогибы пластины являются малыми. Считаем также, что прямолинейное нормальное к срединной поверхности волокно после деформации также остается прямолинейным и нормальным к срединной поверхности (пластина Кирхгофа-Лява). Линеаризуя перемещения по толщине пластины, получаем (4) Далее, для построения уравнений изгиба пластины используется принцип Даламбера. В этом случае соотношения (1), (2) записываются в виде [29] (5) Здесь и - виртуальные перемещения и приращения концентраций; , - поверхностные динамические возмущения. В соответствии с физической постановкой задачи полагаем, что и . Подставляя равенства (2) и (4) в вариационное уравнение (5), по аналогии с методикой, изложенной в работе [29], получаем уравнения упругодиффузионных поперечных колебаний пластины Кирхгофа-Лява с учетом релаксации диффузионных потоков ( - оператор Лапласа) (6) Математическая модель задачи о нестационарном изгибе свободно опертой пластины под действием распределенной нагрузки описывается уравнениями (6), которые дополняются нулевыми начальными условиями и граничными условиями, которые также получаются из вариационного уравнения (5) [29]: (7) (8) (9) 2. Метод решения Решения задачи о нестационарном изгибе пластины (6)-(9) под действием распределенных по поверхности возмущений и представляем в виде ( ) (10) где - поверхностные функции Грина, удовлетворяющие уравнениям (11) и однородным граничным условиям (12) Для построения функций Грина применяем к задаче (11) и (12) преобразование Лапласа по времени. Затем домножаем каждое уравнение на , , и интегрируем в прямоугольнике . Получаем (13) Здесь (14) Решение системы (13) имеет вид ( ) (15) где (16) Оригиналы в (15) находим с помощью вычетов и таблиц операционного исчисления (штрих обозначает производную по параметру ) [29, 36]: (17) (18) Здесь , - нули полинома ; - дополнительные нули многочлена . Подставляя теперь (17) и (18) в представления (14), получаем поверхностные функции Грина задачи (6)-(9). 3. Пример Возьмем для расчета трехкомпонентный ( ) материал - дюралюминий, со следующими характеристиками [37]: Нагрузку задаем в виде ( - функция Хевисайда) (19) Подставляя (19) в формулы (10), получаем с учетом равенств (17) и (18) (20) Результаты вычислений представлены на рис. 2-5. Для расчета использовалось 100 членов ряда Фурье. На рис. 3 изображены пространственно-временные распределения приращений концентраций алюминия, инициированные изгибом пластины. Рис. 3, а позволяет также примерно оценить время, когда массоперенос первой компоненты приобретает статический характер ( ). Ниже приводится сравнение полученного решения с решением чисто упругой задачи. В этом случае, полагая в равенствах (16)-(18), (20) и , получаем следующие выражения для поперечных смещений: На рис. 4 продемонстрировано влияние диффузионного поля на поле перемещений. Показано, что начиная с определенного времени упругодиффузионные колебания начинают сдвигаться по фазе относительно упругих колебаний. Примерное начало сдвига ( ) показано на рис. 4, а. а б Рис. 2. Поперечные смещения пластины: а - б - Fig. 2. The plate deflections: a - ; b - а б Рис. 3. Приращение концентрации первой компоненты (алюминий): а - ; б - Fig. 3. The concentration increment of the first component (aluminum): a - ; b - а б Рис. 4. Поперечные смещения с учетом массопереноса (сплошная линия) и без него (пунктирная линия): а - ; б - Fig. 4. The plate deflections with mass transfer (solid line) and without mass transfer (dashed line): a - ; b - а б Рис. 5. Приращение концентрации первой компоненты (алюминий): а - ; б - Fig. 5. The concentration increment of the first component (aluminum): a - ; b - Влияние релаксационных эффектов на кинетику массопереноса продемонстрировано на рис. 5. Хорошо видно, что начиная с некоторого момента времени, который соизмерим с временем релаксации , релаксационные эффекты начинают затухать (рис. 5, б). Заключение В работе построена математическая модель упругодиффузионных нестационарных колебаний прямоугольной изотропной пластины Кирхгофа-Лява с учетом релаксации диффузионных потоков, описывающая взаимосвязь между механическими и диффузионными полями в сплошных средах. Предложен алгоритм построения поверхностных функций Грина, основанный на использовании преобразования Лапласа и разложений в тригонометрические ряды Фурье. На примере изгиба пластины под действием распределенной по поверхности механической нагрузки продемонстрирован эффект взаимодействия механического и диффузионного полей. Показано, что, с одной стороны, нестационарный изгиб инициирует процесс массопереноса. С другой стороны, диффузия влияет на поле перемещений, что проявляется в виде фазового сдвига механодиффузионных колебаний по отношению к чисто механическим. Отмечено также, что релаксационные эффекты, обусловливающие конечную скорость распространения диффузионных возмущений, проявляются только на конечном промежутке времени, соизмеримом с временем релаксации диффузионных потоков. Указанные результаты представлены в аналитической и графической формах.

Об авторах

Юй Гу

Гуандунский университет нефтехимических технологий

А. В Земсков

Московский авиационный институт (национальный исследовательский университет); НИИ механики МГУ имени М.В. Ломоносова

Д. В Тарлаковский

Московский авиационный институт (национальный исследовательский университет); НИИ механики МГУ имени М.В. Ломоносова

Список литературы

Статистика

Просмотры

Аннотация - 200

PDF (Russian) - 128