CC BY
8
Юй Гу, Земсков А.В., Тарлаковский Д.В. Упругодиффузионные колебания изотропной пластины Кирхгофа-Лява под действием нестационарной распределенной поперечной нагрузки // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2021. - № 3. С. 48-57. DOI: 10.15593/perm.mech/2021.3.05
Yu Gu, Zemskov A.V., Tarlakovskii D.V. Elastic Diffusion Vibrations of an Isotropic Kirchhoff-Love Plate under an Unsteady Distributed Transverse Load. PNRPUMechanics Bulletin, 2021, no. 3, pp. 48-57. DOI: 10.15593/perm.mech/2021.3.05
ВЕСТНИК ПНИПУ. МЕХАНИКА
№ 3,2021 PNRPU MECHANICS BULLETIN
https://ered.pstu.ru/index.php/mechanics/index
Б01: 10.15593/регш.шееЬ/2021.3.05 УДК 539.3, 539.8
УПРУГОДИФФУЗИОННЫЕ КОЛЕБАНИЯ ИЗОТРОПНОЙ ПЛАСТИНЫ КИРХГОФА-ЛЯВА ПОД ДЕЙСТВИЕМ НЕСТАЦИОНАРНОЙ РАСПРЕДЕЛЕННОЙ ПОПЕРЕЧНОЙ НАГРУЗКИ
Гу Юй1, А.В. Земсков2'3, Д.В. Тарлаковский3'2
Гуандунский университет нефтехимических технологий, МаоМинг, Китай
2
Московский авиационный институт (национальный исследовательский университет), Москва, Россия
3
НИИ механики МГУ имени М.В. Ломоносова, Москва, Россия
О СТАТЬЕ
АННОТАЦИЯ
Получена: 20 октября 2020 г. Принята: 13 июля 2021 г. Опубликована: 22 октября 2021 г.
Ключевые слова:
упругая диффузия, преобразование Лапласа, ряды Фурье, функции Грина, пластина Кирхгофа-Лява, нестационарные задачи.
Исследуются нестационарные упругодиффузионные колебания свободно опертой прямоугольной изотропной пластины Кирхгофа-Лява, находящейся под действием распределенной поперечной нагрузки. Для математической постановки задачи используется модель, описывающая связанные упругодиффузионные процессы в сплошных многокомпонентных средах с учетом релаксации диффузионных потоков. Из нее с помощью вариационного принципа Даламбера получены уравнения поперечных колебаний прямоугольной изотропной пластины Кирхгофа-Лява с учетом диффузии. На основе полученных уравнений сформулирована постановка начально-краевой задачи об изгибе свободно опертой изотропной прямоугольной пластины под действием распределенных по поверхности упругодиффузионных возмущений.
Решение задачи о нестационарных упругодиффузионных колебаниях пластины ищется в интегральной форме. Ядрами интегральных представлений являются поверхностные функции Грина, для нахождения которых используется преобразование Лапласа по времени и разложение в двойные тригонометрические ряды Фурье по пространственным координатам. Трансформанты Лапласа функций Грина представлены через рациональные функции параметра преобразования Лапласа. Переход в пространство оригиналов осуществляется аналитически с помощью вычетов и таблиц операционного исчисления. Получены аналитические выражения для поверхностных функций Грина рассматриваемой задачи.
В качестве расчетного примера рассмотрен изгиб свободно опертой механодиффузионной пластины, находящейся под действием внезапно приложенных, распределенных по поверхности нестационарных изгибающих моментов. На примере трехкомпонентного материала выполнено численное исследование взаимодействия нестационарных механического и диффузионного полей в изотропной пластине. Исследовано влияние релаксационных эффектов на кинетику массопере-носа. Решение представлено в аналитической форме и в виде графиков зависимости искомых полей перемещения и приращений концентрации компонент среды от времени и координат.
В заключение приведены основные выводы о влиянии связанности полей и релаксационных эффектов на напряженно-деформированное состояние и массоперенос в пластине.
©ПНИПУ
© Гу Юй - PhD, профессор, e-mail: guyu@mail.buct.edu.cn.
Земсков Андрей Владимирович - д.ф.-м.н., доц., e-mail azemskov1975@mail.ru, ¡D: 0000-0001 -6725-4361. Тарлаковский Дмитрий Валентинович - д.ф.-м.н., проф., завлаб.; e-mail: tdv902@mai.ru, : 0000-0002-5694-9253.
Yu Gu - PhD, Professor, e-mail: guyu@mail.buct.edu.cn.
Andrei V. Zemskov - Doctor of Physical and Mathematical Sciences, Associate Professor, e-mail: azemskov1975@mail.ru, : 0000-0001 -6725-4361.
Dmitry V. Tarlakovskii - Doctor of Physical and Mathematical Sciences, Professor, Head of Laboratories, e-mail: tdv902@mai.ru, : 0000-0002-5694-9253.
Эта статья доступна в соответствии с условиями лицензии Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
ELASTIC DIFFUSION VIBRATIONS OF AN ISOTROPIC KIRCHHOFF-LOVE PLATE UNDER AN UNSTEADY DISTRIBUTED TRANSVERSE LOAD
Yu Gu1, A.V. Zemskov2,3, D.V. Tarlakovskii3'2
Guangdong University of Petrochemical Technology, MaoMing, China
Moscow Aviation Institute (National Research University, Moscow, Russian Federation
Research Institute of Mechanics MSU, Moscow, Russian Federation
ARTICLE INFO
ABSTRACT
Received: 20 October 2020 Accepted: 13 July 2021 Published: 22 October 2021
Keywords:
elastic diffusion, Laplace transform, Fourier series, Green's functions, Kirchhoff-Love plate, unsteady problems.
We investigated an unsteady elastic diffusion vibration of a simply supported rectangular isotropic Kirchhoff-Love plate. The plate is under the action of a distributed transverse load. A model that describes coupled elastic diffusion processes in a multicomponent continuum is used for the mathematical problem formulation. The model is taking into account the diffusion fluxes relaxation. The transverse vibration equations of a rectangular isotropic Kirchhoff-Love plate with diffusion were obtained from the model using the d'Alembert variational principle. The initial-boundary value problem of a freely supported isotropic rectangular plate bending is formulated on the basis of the obtained equations. The plate is under the action of elastic diffusion perturbations distributed over the surface.
The problem solution of an unsteady elastic diffusion plate vibration is sought in an integral form. The surface Green's functions are the kernels of the integral representations. To find the Green's functions, we used the Laplace transform in time and the expansion into double trigonometric Fourier series in spatial coordinates. Green's functions in the image domain are represented in the form of rational functions and depend on the Laplace transform parameter. The transition to the original domain is done analytically through residues and tables of operational calculus. The surface Green's function analytical expressions are obtained.
As a calculation example, we considered a freely supported elastodiffusive plate under the action of suddenly applied unsteady bending moments distributed over the plate surface. By using a three-component continuum, a numerical study of interactions between unsteady mechanical and diffusion fieldsis done for an isotropic plate. The influence of relaxation effects on the kinetics of mass transfer is investigated. The solution is presented in the analytical form, as well as in the graphs of the displacement fields and concentration increments on time and coordinates.
At the end of the publication, the main conclusions are given about the fields coupling effect and the relaxation of diffusion fluxes on the stress-strain state and mass transfer in the plate.
©PNRPU
В настоящее время существует большое количество работ, посвященных изучению механодиффузионных эффектов в сплошных средах. Экспериментальные исследования этой проблемы были начаты в 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
Здесь m = {ml, m2} - распределённые по поверхности моменты; q - распределённая по поверхности поперечная нагрузка; z - распределённые по поверхности плотности объемных источников массопереноса.
Для математической постановки задачи используется модель упругодиффузионных процессов в сплошных средах в прямоугольной декартовой системе координат, которая в случае однородной среды имеет вид [4-20, 33-35]
д.xv dxt
N „ _ч
л(N+')=_£ n(q) ( q = 1, N),
q=1
где а^ и - компоненты тензора напряжений и вектора диффузионного потока, которые для идеального твердого раствора определяются следующим образом
(ч = ЦУ) :
[X8A +ц(5л Sji +8 а Sk j aq^(q);
CTj =
cX,
l q=1
' q • u 4 dXj q ,J " dxjdx,
(2)
Здесь точки обозначают производные по времени. Все величины в (1) и (2) являются безразмерными. Для них приняты следующие обозначения:
T(q)
x,. и ,. ct
X = . , U = , , т = ~' h* h* h*
X* ц* с
X = --7, Ц = --7, \ = —~
X* + 2ц* X* + 2ц* q h*
f п*Ы
l = D = , m h q ch
m(q)D'(q)a*(q)n(q) a*(q)
Aa =-, aq =—--
q pRT0ch q X*+ 2ц*
(3)
F =■
0
h* F *
X* + 2ц* 2 _ X* + 2ц*
Y (q) =
*Y*(q)
h Y
P
h=hL=i, h*
, *
где I - время; х* - прямоугольные декартовы коорди-
* ^ 1 *
наты; и* - компоненты вектора перемещений; п -
* *
толщина пластины; I* и I* - длина и ширина пластины; ^(ч) = п(ч) - п(ч) - приращение концентрации q-й компоненты вещества в составе многокомпонентной
(ч) (ч)
среды; П ' и щ' - актуальная и начальная концентрации q-го вещества; X* и ц* - упругие постоянные Ламе; р - плотность; а*(ч) - коэффициенты, характеризующие объёмное изменение среды за счёт диффузии;
- коэффициенты самодиффузии; Я - универсальная газовая постоянная; Т0 - температура среды;
т(ч) - молярная масса q-го вещества; ^ * - компоненты вектора плотности массовых сил; У ' - плотность источников массопереноса; т( ч) - время релаксации диффузионных потоков.
c
Полагаем, что поперечные прогибы пластины являются малыми. Считаем также, что прямолинейное нормальное к срединной поверхности волокно после деформации также остается прямолинейным и нормальным к срединной поверхности (пластина Кирхгофа-Лява). Линеаризуя перемещения по толщине пластины, получаем
dw (x, x,т)
dx,
dw (x, X,т)
dx„
, (x, x, x,т)= w( x, x,x);
"П(г)(х,х2,Х3,т)= N (х„х2,т) + хНч (х,х2,т). (4)
Далее, для построения уравнений изгиба пластины используется принцип Даламбера. В этом случае соотношения (1), (2) записываются в виде [29]
i
(
«,--—Ft
dx.
v j
\
5udG +
q= g
Sil1+'. dd
(
8j\q) dx.
-Y
(q)
5^(q)dG -
+
+ ii( jj -P)8".dS +
n0
(5)
ц=1 п..
Здесь Ъщ и - виртуальные перемещения и приращения концентраций; р , 1(^ - поверхностные динамические возмущения. В соответствии с физической постановкой задачи полагаем, что р = 0 и = 0 .
Подставляя равенства (2) и (4) в вариационное уравнение (5), по аналогии с методикой, изложенной в работе [29], получаем уравнения упругодиффузион-ных поперечных колебаний пластины Кирхгофа-Лява с учетом релаксации диффузионных потоков (А - оператор Лапласа)
N
Дм>-12м> = ДДи' + ^а^АЯ^ -12(сЦут + д);
q=1
1/2
m = i x3Fdx3 (i = 1,2);
-1/2
Hq+^qHq=DqMiq+hq^w + Uzq,
(6)
12
= | 7хъёхъ.
-12
Математическая модель задачи о нестационарном изгибе свободно опертой пластины под действием распределенной нагрузки описывается уравнениями (6), которые дополняются нулевыми начальными условия-
ми и граничными условиями, которые также получаются из вариационного уравнения (5) [29]:
lx1 =0
d2 w „ —w +x 52 w dx? N + Sa q=1
+ X dx2 d2 w dx? N + Va — q=1
x^.+ dx? d2 w N + Va q=1
„ d2 w x^+ dx d2 w dx2 N + Vaq q=1
= 0, w| x1 ='1 = 0, w| x
= 0;
x =0
= 0;
= 0;
(7)
= 0;
x2 =l2
x2 =0
HA = 0, H\ = 0, Hl = 0, Hl
qlx1 =0 qlx1 =4 q\x2 =0 qlx2 =2
= 0; (8) = 0. (9)
2. Метод решения
Решения задачи о нестационарном изгибе пластины (6)-(9) под действием распределенных по поверхности возмущений ^ (X, X, т) = -12 (ц + div т) и
(х, X, т) = 12г(?) представляем в виде (ц = 1, N +1)
V (х, X, т) =
N+1 т 11 12
= ЕЛ/ (х1, х,С, т- /) ^ с, /) аа
к=! 0 0 0
Нц (^ Х2 , т) =
N+2 т 11 12
= Е///(х1,х2,С,т-/)^с,/)ада/, (10)
*=1 0 0 0
где Отк - поверхностные функции Грина, удовлетворяющие уравнениям
AGli-12Gli=AAGli +
+ Zaq AGq+1,i -51k S( x1 -^)8( x2 -C)S(x) ;
q=1
^q+l,k + Tq^q+l,k
= DqAG— +AqAAG1k + 5,++5(x -i;)5(x -£)5(T) (11)
q
и однородным граничным условиям
N Л
fd2 Glt+k d2G\k
v
dx;
dx:
rd2Glk+x d2 G,k
2 q=1
N
q+1, k
= 0;
x =0
dx
dx2 q=1
SaqGq
q+1,k
^ d2G,k , d2G,k
v dxf
dx.
/
N
+ SaqGq
= 0;
q+1,k
q=1
/
= 0;
x2 =0
M
G
+
f. д2 Glk д2 Glk N ^ X-^ +-^ + ZaqGq+l, k
дх1 дх2 q=l
= G;
_4_ lll2
Sqp ^ЛapПp (Vnm , s)
A+l (Vnm , s) Qq (Vnm , s)
sin Ц-mС sin Xn1
G
= G, Gj , = G, G,| = G, G,| , = G; (12)
' lklxl =4 ' lklx2 =g ' lklx2 =l2 ' v y
G
q+l,k
X = G
= G, G„
G
q+l,k
1x2 = G
= G, G„
= G; = G.
Для построения функций Грина применяем к задаче (11) и (12) преобразование Лапласа по времени. Затем домножаем каждое уравнение на sin Xn x sin |ira x2, Xk = лп / Zj, !m = лот / /2 и интегрируем в прямоугольнике D = [0, Z ] х [0,4 ]. Получаем
К (Vnm , *) GLS (Xn , , С, s)"
N
vlG%u k (Xn , Mm , I, С, s ) =
q=1
4
= 778iksin im Сsin Xn
-AqvlGfr (Xn, im, I,С, s) +
+ kq+1 (vnm , s) GL:i,i (Xn , im , I, С, s) = 4
= T^8q+1,k sin lmC sin Xn|-
lll2
(13)
Здесь
/ ,
nm
kl (vnm , s)= s2 (V^m + 12)+ V kq+l (Vnm , s)= s + V' + Dq VL , Vnm = X2 + MÍ ,
GLss (Xn, ^, I, С, s) =
4 ll l2
= — jjGft (X,X,I,С,s)sinXnX sinцтx2dx2dx;
lll2 G G
GL (X1, X2, |, С, s) = = ZZGLss (Xn,Ц-m,I,С,s)sinXnXsmMmx2. (14)
Решение системы (13) имеет вид (q, p = 1,N)
Gif (Xn,im,I,С,s) = 4-sinImCsmXnI;
lll2 P (v„, s )
gls+l (xn, цга, i, c, s ) =
4 v2 a П (v , s )
4 nm q q \ nm5 /
lll2 P (v„ , s)
G^l (Xn, Mm, I,С, s) =
4 Vnm^(Vnm , s) .
sin ЦтС sin X„|:
sin Mm С sin Xn1
lll2 Qq (vnm , s)
GL+1,p+1 (Xn , Mm , I, С, s) =
где
N
^^ (Vnm , s) = П kJ+1 (Vnm , s); J=1
N
ПJ (Vnm , s)= П kr+1 ( Vnm , s);
r = 1, r
(16)
P ( Vnm , s) = kl ( Vnm , s) П (Vnm , s) _ V™ Z aJЛJПJ (Vnm , s) ;
J=1
(2ч (У пт , 5) = кч+1 (У пт , 5) Р (^ пт , 5) •
Оригиналы в (15) находим с помощью вычетов и таблиц операционного исчисления (штрих обозначает производную по параметру 5 ) [29, 36]:
Gk (Xn , Ц-m , I, С, *) =
_4 .
lll2
2 N +2
sin MmC sin Xn1 Z Ak (Vnm , sl (Vnm )) eXP (sl (Vnm Ь)
l=1
Gq"+U (Xn, Mm, I, C, t) =
^ 2 N +4
= T7sin MmC sin Xn1 Z 4+1,1 ( Vnm , sl (Vnm )) eXP ( sl (Vnm Ь) lll2 l=l
Gq"+l,p+l (Xn , Mm , I, C, T) =
4 . . Л exP (Xl (v nm ^
= JJ- sin Mm C sin Xn -/ +
lll2 l=1 kq+l ( Vnm , Xl (Vnm ))
¿J 2N+4
^T7sin M„C sin X„| Z (V»m, sl (V«m )) eXP ( sl (Vnm (17)
lll2
A(l)(V s (V ))= "(V°m , sl (Vnm )) ,
Aiz (Vnm , sl (Vnm ))= P,(vnm , SI (vnm )) ;
Al,q+l (vnm , sl (vnm )) =
V2maqП, (Vnm , sl (Vnm )) , P'(Vnm , sl (Vnm )) ;
A)u (Vnm , sl (Vnm = ; (18)
Xl (Vnm )
Q'q (Vnm , sl (Vnm ))
-l + 7l-4x/g+1(vBm,0)
2x„
(15) 4+1,p+1 (Vnm , sl (Vnm )) =
VnmЛдapПp (Vnm , sl (Vnm )) Qq (Vnm , sl (Vnm ))
Здесь (\пт ), ] = 1,2И + 2 - нули полинома Р (Улт, я); %1 (У„т )= Я2М+2+1 (Vпт ) - дополнительные нули многочлена 0,ц (vлm , 5) .
Подставляя теперь (17) и (18) в представления (14), получаем поверхностные функции Грина задачи (6)-(9).
i=i
3. Пример
Возьмем для расчета трехкомпонентный (N = 2) материал - дюралюминий, со следующими характеристиками [37]:
Н
Н
X* = 6,93-1010—, ц* = 2,56-10 —;
T = 800 К, р = 2780 , l* = l* = 0,01 м;
м3
D*(1) = 7,74 -10-14 —, D*(2) = 6,67 -10-14 —;
с с
n® = 0,935, n(2) = 0,045, И = 0,0005 м;
а*«= 1,55-107 ДЖ, а*(2) = 6,14-107 ^;
кг кг
ш(1) = 0,027 —, ш(2) = 0,064 —.
моль моль
Нагрузку задаем в виде (H (т) - функция Хевисайда)
F (x, X, т) = -12 (q + div m) = H (т) sin nx2
Fq+1 (^ X2, т)= 0.
(19)
Подставляя (19) в формулы (10), получаем с учетом равенств (17) и (18)
2 nx ■» 2 N+ 2
w (^ x2, т) = -sin ~XL Z Z (Vn1, Sl (Vn1 ))X
l1 l2 n=1 l=1
[exp(s (vm)т)-1][1 -(-1)n] _
-:—7—ч-sin Xnxl,
XnS (Vn1 )
2 TCx ™ 2 N+4
Hq (X , x2 , T) = 2 Sin ~xL Z Z 4+U (Vn1 , S (Vn1 ))>
/1
l2 n=1 /=1
[exp(s, (уя1 )x)-l][l-(-1)"] , x---——[-] sin(20)
XnS/ (Уя1 )
Результаты вычислений представлены на рис. 2-5. Для расчета использовалось 100 членов ряда Фурье. На рис. 3 изображены пространственно-временные распределения приращений концентраций алюминия, инициированные изгибом пластины. Рис. 3, а позволяет также примерно оценить время, когда массоперенос первой компоненты приобретает статический характер
(х ~3 -1015).
Ниже приводится сравнение полученного решения с решением чисто упругой задачи. В этом случае, полагая в равенствах (16)-(18), (20) D = 0 и а = 0, получаем следующие выражения для поперечных смещений:
4 sin |_ц x ^ sin \nxy 1 -(-1)" 1 - cos (s01x)
w =-—7-г--;
/1 n=1 k1(y n 1, % ) Xn %
2/2 \-1/2 S01 =ynm (У nm + 12) ■
На рис. 4 продемонстрировано влияние диффузионного поля на поле перемещений.
Показано, что начиная с определенного времени упругодиффузионные колебания начинают сдвигаться по фазе относительно упругих колебаний. Примерное
начало сдвига (х ~ 1010) показано на рис. 4, а.
Рис. 2. Поперечные смещения пластины: а - w(x, 0.5/2, т); б - x ,X, 2.25 -102) Fig. 2. The plate deflections: a - w(x, 0.5/2, т) ; b - x ,X, 2.25 -102)
x
б
а
а б
Рис. 3. Приращение концентрации первой компоненты (алюминий): а - Их (x, 0.5/2, т); б - H (x ,X, 5 -1015)
Fig. 3. The concentration increment of the first component (aluminum): a - H (x, 0.5/2, т); b - H (x ,X, 5-1015)
а б
Рис. 4. Поперечные смещения w(0,5^, 0,5/2, т) с учетом массопереноса (сплошная линия) и без него (пунктирная линия): а - те[1010 - 2•102,1010 + 2,5-102 ] ; б - те[ - 2 •102,1011 + 2,5-102 ] Fig. 4. The plate deflections w(0.5^, 0.5/2, т) with mass transfer (solid line) and without mass transfer (dashed line): a - те[1010 - 2-102,1010 + 2.5-102]; b - те[1011 - 2 •102,1011 + 2.5 402]
2.0
1.5
10
0.5
10 / / / / /
/ / /
/ / / / /
/ / / / /
/ / / / z—-
2.0 4.0 6.0 8.0 T.ltf
1.0 0.8 0.6 0.4 0.2 0
10 8 /■ /./ / /
/./ У
// /./ // / /
// // /./
/./ /■/
/ /., /у/ /
/./ /У Л/
0 1. 0 2.0 3.0 4.0 T,10'°
а б
Рис. 5. Приращение концентрации H (0,Ц, 0,5/2, т) первой компоненты (алюминий): а - те[0,109 ] ; б - т е [0, 5 • 101(
Fig. 5. The concentration increment H (0Ц, 0.5/2, т) of the first component (aluminum): a - те[0,109 ] ; b - те[ 0, 5-1010
Влияние релаксационных эффектов на кинетику массопереноса продемонстрировано на рис. 5. Хорошо видно, что начиная с некоторого момента времени, который соизмерим с временем релаксации т , релаксационные эффекты начинают затухать (рис. 5, б).
Заключение
В работе построена математическая модель упруго-диффузионных нестационарных колебаний прямоугольной изотропной пластины Кирхгофа-Лява с учетом релаксации диффузионных потоков, описывающая взаимосвязь между механическими и диффузионными полями в сплошных средах. Предложен алгоритм построения поверхностных функций Грина, основанный
Библиографический список
1. Горский В.С. Исследование упругого последействия в сплаве Cu-Au с упорядоченной решеткой // Журнал экспериментальной и теоретической физики. - 1936. - Т. 6, № 3. - С. 272-276.
2. Nachtrieb N.H., Handler G.S A relaxed vacancy model for diffusion incrystalline metals // Acta Metallurgica. - 1954. -Vol. 2, no. 6. - P. 797-802.
3. Petit J., Nachtrieb N. H. Self-diffusion in liquid gallium // Journal of Chemical Physics. - 1956. - Vol. 24. - P. 1027.
4. Еремеев B.C. Диффузия и напряжения. - М.: Энергоатомиздат, 1984. - 182 с.
5. Князева А.Г. Введение в термодинамику необратимых процессов. - Томск: Изд-во «Иван Федоров», 2014. - 172 с.
6. Afram A.Y., Khader S.E. 2D problem for a half-space under the theory of fractional thermoelastic diffusion // American Journal of Scientific and Industrial Research. - 2014. - Vol. 6, no. 3. - P. 47-57.
7. Aouadi M., Boulehmi K. Partial exact controllability for inhomogeneous multidimensional thermoelastic diffusion problem // Evolution Equations and Control Theory. - 2016. - Vol. 5, no. 2. -P. 201-224.
8. Atwa S.Y., Egypt Z. Generalized thermoelastic diffusion with effect of fractional parameter on plane waves temperature-dependent elastic medium // Journal of Materials and Chemical Engineering. - 2013. - Vol. 1, iss. 2. - P. 55-74.
9. Boulehmi K., Aouadi M. Decay of solutions in nonhomogeneous thermoelastic diffusion bars // Applicable Analysis. - 2014. - Vol. 93, no. 2. - P. 281-304.
10. Choudhary S., Deswal S. Mechanical loads on a generalized thermoelastic medium with diffusion // Meccanica. -2010. - Vol. 45. - P. 401-413.
11. Deswal S., Kalkal K.K., Sheoran S.S. Axi-symmetric generalized thermoelastic diffusion problem with two-temperature and initial stress under fractional order heat conduction // Physica B: Condensed Matter. - 2016. - Vol. 496. - P. 57-68.
12. Elhagary M.A. Generalized thermoelastic diffusion problem for an infinitely long hollow cylinder for short times // Acta Mech. - 2011. - Vol. 218. - P. 205-215.
13. El-Sayed A.M. A two-dimensional generalized thermoelastic diffusion problem for a half-space // Mathematics and Mechanics of Solids. - 2016. - Vol. 21, no. 9. - P. 1045-106.
на использовании преобразования Лапласа и разложений в тригонометрические ряды Фурье.
На примере изгиба пластины под действием распределенной по поверхности механической нагрузки продемонстрирован эффект взаимодействия механического и диффузионного полей. Показано, что, с одной стороны, нестационарный изгиб инициирует процесс массопереноса. С другой стороны, диффузия влияет на поле перемещений, что проявляется в виде фазового сдвига механодиффузи-онных колебаний по отношению к чисто механическим. Отмечено также, что релаксационные эффекты, обусловливающие конечную скорость распространения диффузионных возмущений, проявляются только на конечном промежутке времени, соизмеримом с временем релаксации диффузионных потоков. Указанные результаты представлены в аналитической и графической формах.
14. Hwang C.C., Huang I.B. Diffusion in hollow cylinders with mathematical treatment // International Journal of Engineering Research and Development. - 2012. - Vol. 3, iss. 8. -P. 57-75.
15. Knyazeva A.G. Model of medium with diffusion and internal surfaces and some applied problems // Mater. Phys. Mech.
- 2004. - Vol. 7, no. 1. - P. 29-36.
16. Kumar R., Chawla V. A study of Green's functions for three-dimensional problem in thermoelastic diffusion media // African Journal of Mathematics and Computer Science Research.
- 2014. - Vol. 7, no. 7. - P. 68-78.
17. Kumar R., Devi S., Sharma V. Plane waves and fundamental solution in a modified couple stress generalized thermoelastic with mass diffusion // Materials Physics and Mechanics. - 2015. - Vol. 24. - P. 72-85.
18. Othman M.I.A., Elmaklizi Y.D. 2-D problem of generalized magneto-thermoelastic diffusion, with temperature-dependent elastic moduli // Journal of Physics. - 2013. - Vol. 2, no. 3. - P. 4-11.
19. Sharma J.N., Sharma N.K., Sharma K.K. Transient waves due to mechanical loads in elasto-thermo-diffusive solids // Advances in Applied Mathematics and Mechanics. - 2011. -Vol. 3, no. 1. - P. 87-108.
20. Sherief H.H., El-Maghraby N.M. A Thick plate problem in the theory of generalized thermoelastic dffusion // Int. J. Thermophys. - 2009. - Vol. 30. - P. 2044-2057.
21. Швец Р.Н., Флячок В.М. Уравнения механодиф-фузии анизотропных оболочек с учетом поперечных деформаций // Математические методы и физико-механические поля. - 1984. - Вып. 20. - С. 54-61.
22. Швец Р.Н., Флячок В.М. Вариационный подход к решению динамических задач механотермодиффузии анизотропных оболочек // Математическая физика и нелинейная механика. - 1991. - № 16. - С. 39-43.
23. Aouadi M. Copetti M.I.M. Analytical and numerical results for a dynamic contact problem with two stops in thermoelastic diffusion theory // ZAMM Z. Angew. Math. Mech. -2015. - P. 1-24. DOI 10.1002/zamm.201400285.
24. Copetti M.I.M., Aouadi M. A quasi-static contact problem in thermoviscoelastic diffusion theory // Applied Numerical Mathematics. - 2016. - Vol. 109. - P. 157-183.
25. Aouadi M. Miranville A. Smooth attractor for a nonlinear thermoelastic diffusion thin plate based on Gurtin-Pipkin's model // Asymptotic Analysis. - 2015. - Vol. 95. -P. 129-160.
26. Aouadi M. On thermoelastic diffusion thin plate theory // Appl. Math. Mech. -Engl. Ed. - 2015. - Vol. 36, no. 5. -P. 619-632.
27. Aouadi M., Miranville A. Quasi-stability and global attractor in nonlinear thermoelastic diffusion plate with memory // Evolution Equations and Control Theory. - 2015. - Vol. 4, no. 3. -P. 241-263.
28. Bhattacharya D., Kanoria M. The influence of two temperature generalized thermoelastic diffusion inside a spherical shell // International Journal of Engineering and Technical Research (IJETR). - 2014. - Vol. 2, iss. 5. - P. 151-159.
29. Zemskov A.V., Tarlakovskii D.V. Modelling of rectangular Kirchhoff plate oscillations under unsteady elastodiffusive perturbations // Acta Mechanica. -2021. DOI: 10.1007/s00707-020-02879-1
30. Afanasieva O.A., Zemskov A.V. Unsteady elastic-diffusion oscillations of a simply supported kirchhoff plate under the distributed transverse load action // Gdoutos E., Konsta-Gdoutos M. (eds) Proceedings of the Third International Conference on Theoretical, Applied and Experimental Mechanics. ICTAEM 2020. Structural Integrity. - Vol 16. - Springer, Cham, 2020. - P. 181-186.
References
1. Gorskij V.S. Issledovanie uprugogo posledejstviya v splave Cu-Au s uporyadochennoj reshetkoj [Study of elastic aftereffect in Cu-Au alloy with an ordered lattice]. Zhurnal e'ksperimental'noj i teoreticheskoj fiziki, 1936, vol. 6, no 3, pp. 272-276 [in Russian].
2. Nachtrieb N.H., Handler G.S. A relaxed vacancy model for diffusion incrystalline metals. Acta Metallurgica, 1954, vol. 2, no 6, pp. 797-802.
3. Petit J., Nachtrieb N. H. Self-Diffusion in Liquid Gallium. Journal of Chemical Physics, 1956, vol. 24, p. 1027.
4. Eremeev V.S. Diffuziya i napryazheniya [Diffusion and stress] //Moscow, Energoatomizdat, 1985, 182 p. [In Russian].
5. Knyazeva A.G. Vvedenie v termodinamiku neobratimy'x processov [Introduction to the thermodynamics of irreversible processes], Tomsk, Ivan Fedorov Publishing House, 2014, 172 p. [in Russian].
6. Afram A.Y., Khader S.E. 2D Problem for a Half-Space under the Theory of Fractional Thermoelastic Diffusion. American journal of scientific and industrial research, 2014, vol. 6, no. 3, pp. 47-57.
7. Aouadi M., Boulehmi K. Partial exact controllability for inhomogeneous multidimensional thermoelastic diffusion problem. Evolution equations and control theory, 2016, vol. 5, no. 2, pp. 201-224.
8. Atwa S.Y., Egypt Z. Generalized Thermoelastic Diffusion With Effect of Fractional Parameter on Plane Waves Temperature-Dependent Elastic Medium. Journal of Materials and Chemical Engineering, 2013, vol. 1, iss. 2, pp. 55-74.
9. Boulehmi K., Aouadi M. Decay of solutions in nonhomogeneous thermoelastic diffusion bars. Applicable Analysis, 2014, vol. 93, no. 2, pp. 281-304.
10. Choudhary S., Deswal S. Mechanical loads on a generalized thermoelastic medium with diffusion. Meccanica, 2010, vol. 45, pp. 401-413.
31. Le K.C. Vibrations of shells and rods. - Berlin, Springer Verlag, 1999. - 425 p.
32. Михайлова Е.Ю., Тарлаковский Д.В., Федотенков Г.В. Общая теория упругих оболочек. - М.: МАИ, 2018. - 112 с.
33. Igumnov L.A., Tarlakovskii D.V., Zemskov A.V. A two-dimensional nonstationary problem of elastic diffusion for an orthotopic one-component layer // Lobachevskii Journal of Mathematics. - 2017. - Vol. 38, no. 5. - P. 808-817.
34. Земсков А.В., Тарлаковский Д.В. Двумерная нестационарная задача упругой диффузии для изотропного однокомпо-нентного слоя // Прикладная механика и техническая физика. -2015. - Т. 56, № 6. - С. 102-110 = Zemskov A.V., Tarlakovskiy D.V. Two-dimensional nonstationary problem elastic for diffusion an isotropic one-component layer // Journal of Applied Mechanics and Technical Physics. - 2015. - Vol. 56, no. 6. - P. 1023-1030.
35. Tarlakovskii D.V., Vestyak V.A., Zemskov A.V. Dynamic processes in thermoelectromagnetoelastic and thermoelastodiffusive media // Encyclopedia of Thermal Stress. -Vol. 2. - Springer Dordrecht Heidelberg; New York, London, Springer reference, 2014. - P. 1064-1071.
36. Диткин В.А., Прудников А.П. Справочник по операционному исчислению. - М.: Высшая школа, 1965. - 586 с.
37. Физические величины: справочник / А.П. Бабичев, Н.А. Бабушкина, А.М. Братковский [и др.] / под общ. ред. И.С. Григорьева, И.З. Мелихова. - М.: Энергоатомиздат, 1991. - 1232 с.
11. Deswal S., Kalkal K.K., Sheoran S.S. Axi-symmetric generalized thermoelastic diffusion problem with two-temperature and initial stress under fractional order heat conduction. Physica B: Condensed Matter, 2016, vol. 496, pp. 57-68.
12. Elhagary M.A. Generalized thermoelastic diffusion problem for an infinitely long hollow cylinder for short times. Acta Mech, 2011, vol. 218, pp. 205-215.
13. El-Sayed A.M. A two-dimensional generalized thermoelastic diffusion problem for a half-space. Mathematics and Mechanics of Solids, 2016, vol. 21, no. 9, pp. 1045-106.
14. Hwang C.C., Huang I.B. Diffusion in hollow cylinders with mathematical treatment. International Journal of Engineering Research and Development, 2012, V. 3, Is. 8, pp. 57-75.
15. Knyazeva A.G. Model of medium with diffusion and internal surfaces and some applied problems. Mater. Phys. Mech., 2004, vol. 7, no. 1, pp. 29-36.
16. Kumar R., Chawla V. A study of Green's functions for three-dimensional problem in thermoelastic diffusion media. African journal of mathematics and computer science research, 2014, vol. 7, no 7, pp. 68-78.
17. Kumar R., Devi S., Sharma V. Plane waves and fundamental solution in a modified couple stress generalized thermoelastic with mass diffusion. Materials Physics and Mechanics, 2015, vol. 24, pp. 72-85.
18. Othman M.I.A., Elmaklizi Y.D. 2-D Problem of Generalized Magneto- Thermoelastic Diffusion, with Temperature-Dependent Elastic Moduli. Journal of physics, 2013, vol. 2, no 3, pp. 4-11.
19. Sharma J. N., Sharma N. K., Sharma K. K. Transient Waves Due to Mechanical Loads in Elasto-Thermo-Diffusive Solids. Advances in Applied Mathematics and Mechanics, 2011, vol. 3, no. 1, pp. 87-108.
20. Sherief H.H., El-Maghraby N.M. A Thick Plate Problem in the Theory of Generalized Thermoelastic Diffusion. Int. J. Thermophys., 2009, vol. 30, pp. 2044-2057.
21. Shvets R.N., Flyachok V.M. Uravneniya mexanodiffuzii anizotropny'x obolochek s uchetom poperechny'x deformacij [The equations of mechanothermodiffusion of anisotropic shells taking account of transverse strains]. Mat. Met. Fiz.-Mekh. Polya, 1984, no 20, pp. 54-61. [In Russian].
22. Shvets R.N., Flyachok V.M. Variacionny'j podxod k resheniyu dinamicheskix zadach mexanotermodiffuzii anizo-tropny'x obolochek [A variational approach to solving dynamic problems of mechanothermal diffusion of anisotropic shells]. Mat. fiz. i nelinejnaya mexanika, 1991, no 16, pp. 39-43. [In Russian].
23. Aouadi M., Copetti M.I.M. Analytical and numerical results for a dynamic contact problem with two stops in thermoelastic diffusion theory. ZAMM Z. Angew. Math. Mech., 2015, pp. 1-24. DOI 10.1002/zamm.201400285.
24. Copetti M.I.M., Aouadi M. A quasi-static contact problem in thermoviscoelastic diffusion theory. Applied Numerical Mathematics, 2016, vol. 109, pp. 157-183.
25. Aouadi M. Miranville A. Smooth attractor for a nonlinear thermoelastic diffusion thin plate based on Gurtin-Pipkin's model. Asymptotic Analysis, 2015, vol. 95, pp. 129-160.
26. Aouadi M. On thermoelastic diffusion thin plate theory. Appl. Math. Mech. Engl. Ed., 2015, vol. 36, no. 5, pp. 619-632.
27. Aouadi M., Miranville A. Quasi-stability and global attractor in nonlinear thermoelastic diffusion plate with memory. Evolution equations and control theory, 2015, vol. 4, no 3, pp. 241-263.
28. Bhattacharya D., Kanoria M. The influence of two temperature generalized thermoelastic diffusion inside a spherical shell. International Journal of Engineering and Technical Research (IJETR), 2014, vol. 2, iss. 5, pp. 151-159.
29. Zemskov A.V., Tarlakovskii D.V. Modelling of rectangular Kirchhoff plate oscillations under unsteady
elastodiffusive perturbations. Acta Mechanica., 2021, DOI: 10.1007/s00707-020-02879-1.
30. Afanasieva O.A., Zemskov A.V. Unsteady Elastic-Diffusion Oscillations of a Simply Supported Kirchhoff Plate Under the Distributed Transverse Load Action. In: Gdoutos E., Konsta-Gdoutos M. (eds) Proceedings of the Third International Conference on Theoretical, Applied and Experimental Mechanics. ICTAEM 2020. Structural Integrity, vol 16. Springer, Cham, 2020, pp. 181-186.
31. Le K.C. Vibrations of shells and rods. Berlin, Springer Verlag, 1999, 425 P.
32. Mikhailova E.Yu., Tarlakovskii D.V., Fedotenkov G.V. Obshchaya teoriya uprugikh obolochek [General theory of anisotropic shells]. Moscow, MAI, 2018, 112 p. [In Russian].
33. Igumnov L.A., Tarlakovskii D.V., Zemskov A.V. A two-dimensional nonstationary problem of elastic diffusion for an orthotropic one-component layer. Lobachevskii Journal of Mathematics, 2017, vol. 38, no. 5, pp. 808-817.
34. Zemskov A.V., Tarlakovskiy D.V. Two-dimensional nonstationary problem elastic for diffusion an isotropic one-component layer. Journal of Applied Mechanics and Technical Physics, 2015, vol. 56, no. 6, pp. 1023-1030.
35. Tarlakovskii D.V., Vestyak V.A., Zemskov A.V. Dynamic Processes in Thermoelectromagnetoelastic and Thermoelastodiffusive Media. Encyclopedia of thermal stress, volume 2, Springer Dordrecht Heidelberg New York London, Springer reference, 2014, pp. 1064-1071.
36. Ditkin V.A., Prudnikov A.P. Spravochnik po operacionnomu ischisleniyu [Operational Calculus Handbook]. Moscow, «Vy'sshaya shkola», 1965, 586 P.
37. Babichev A.P., Babushkina N.A., Bratkovskij A.M. ets. Edited by Grigorev I.S., Melixov I.Z. Fizicheskie velichiny: Spravochnik [Physical quantities: Handbook] // Moscow, Energoatomizdat, 1991, 1232 P.
