Двухмерные нестационарные волны в электромагнитоупругой полуограниченной среде Текст научной статьи по специальности «Физика»

Механика деформируемого твердого тела Вестник Нижегородского университета им. Н.И. Лобачевского, 2011, № 4 (4), с. 1423-1424

УДК 539.3

ДВУХМЕРНЫЕ НЕСТАЦИОНАРНЫЕ ВОЛНЫ В ЭЛЕКТРОМАГНИТОУПРУГОЙ ПОЛУОГРАНИЧЕННОЙ СРЕДЕ

© 2011 г. В.А. Вестяк, Д.В. Тарлаковский

Московский авиационный институт (государственный технический университет)

tdv902@mai.ru

Поступила в редакцию 15.06.2011

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

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

Уравнения движения электромагнитоупругой среды

Рассматривается частный случай электромагнитоупругой среды — изотропные проводники без учета температурных полей. Соответствующая замкнутая система линеаризованных уравнений движения может быть записана в виде [1]: ри = (Я + |)grad divu + |Аи + Ее,

гч =

У,и] +У ^

= + 2|ог г

Ее = а[РеоЕ + Ре Ео + у([ ] , Н] + [], Н о])],

Ё + уеЕ = с2(А - grad div)E - Г(11),

(1)

(2)

Г (и) = у у [ и, Но] + у г и,

где и = игег- и Ее = Е'е е { — векторы перемещения и силы Лоренца; Е = Ее{ и Н = И'е^ — векторы на-пряженностей электрического и магнитного полей; ] = /ег- — плотность тока; е1; е2, е3 — ковариан-тный базис некоторой криволинейной системы ко -ординат р — плотность; X и | — упругие

постоянные Ламе; ре — плотность зарядов; а, у, уе, уу, уг и с2 — коэффициенты, выражающиеся через коэффициенты диэлектрической и магнитной проницаемости е и |е, коэффициент электропроводимости а, скорость света с и рео; А — оператор Лапласа; точками обозначены производные по времени ^ нижний индекс «о» соответствует значению величин в начальном состоянии.

При этом компоненты тензоров напряжений и деформаций ац и ец связаны с перемещениями так:

(3)

(4)

(5)

2 - V где Уг- — оператор ковариантного дифференцирования; gij — компоненты метрического тензора.

Напряженности полей, плотность тока и плотность зарядов должны удовлетворять следующим соотношениям электродинамики [2, 3]:

с rotE = —|е Н, с гоШ = 4 л] + еЁЁ, е divE = 4пре, ] = с(Е + | ес _1[и, Н о]) + р

ео11,

Полагается, что в начальный момент времени t = о возмущения в среде отсутствуют:

и t=о = о, и| t=о = о, Е| t=о = о,

Е t=о=о.

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

и времени. При этом на границе ^ = задаются кинематические или силовые и электродинамические условия, и решение в области ^ > считается ограниченным.

Решение методом малого параметра

Даже в частных случаях одномерных волн решение системы уравнений (1), (2) не удается построить аналитически. Поэтому далее исполь-

1424

В.А. Вестяк, Д.В. Тарлаковский

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

да да да

u =Хи* ^, E =ХEm^, H =Х ^am,

m =0

m=0

m=0

j =Z jm a m, Pe = Zpe™-

m =0

m=0

Подстановка этих рядов в (1), (2) приводит к рекуррентной системе уравнений:

р!д0 = (Я + ц)grаd divu0 + цЛи0, (6)

Е m +Т е Е m =

= c2(Л-graddiv)Em - f(iim ), m > 0, (7)

Риm + + m + ЦЛи m + g m'

gm = ре0^-1 + ре( m-1)E0 +

+ У([ЬHm-l] + [^-ЬВД, m > 1.

Связь коэффициентов рядов для величин E, Н, j и ре вытекает из формул (3), (4):

c ГО^ = еН m, rOtHm = 4пjm + ^m,

^^ = 4пPеm,

jm = 0^ +Цеc~1[^i m, 10]) + Р е0 и m, m > 0.

Соответствующие начальные условия, как следует из равенств (5), остаются однородными:

'm\t=0

= 0,

m t=0

= 0,

Em t=0 = 0

Em 11=0 = 0 (m > 0).

(8)

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

Нестационарные поверхностные возмущения в полуплоскости и пространстве со сферической плоскостью

Начально-краевая задача (6)-(8) с соответствующими краевыми условиями при m = 0 явля-

ется задачей о распространении двухмерных нестационарных поверхностных возмущений в упругой среде без учета взаимодействия с электромагнитным полем, а при m > 1 она является упругой задачей при заданных нестационарных объемных возмущениях. Решения всех этих задач существенно зависят от геометрии области. Рассматриваются полуплоскость z > 0 в прямоугольных декартовых координатах = z, = x, = y и пространство r > R со сферической полостью радиуса R. Например, в последнем случае для перемещений имеют место равенства (звездочки означают свертки по координате x и времени):

umk ( xzт) =

ТО

= J Gki( x, z, S, X)**gml( x, S, + 0 ТО

+ J Gk з( x, z, S, t) * *gm3 (x, S, T)dS 0

(k = 1,3).

Здесь um1' gml и ^ gm3 — проекЦии векторов Um,

gm на оси x и z соответственно; Gkl (k, l = 1, 3) — объемные функции Грина соответствующих задач.

Работа выполнена в рамках РФФИ (гранты №0908-00470, 10-08-90412) и ФЦП «Научные и педагогические кадры инновационной России» (проект П2235).

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

1. Вестяк В. А., Тарлаковский Д.В. // Методи роз-в'язування прикладных задач мехашки деформ1вного твердого тша: Зб. наук. праць Дшпропетр. нацюн. ун-та. Дшпропетровськ: 1МА-прес. 2009. Вип. 10. С. 57-62.

2. Гринченко В.Т., Улитко А.Ф., Шульга Н.А. Ме -ханика связанных полей в элементах конструкций. Электроупругость / Под. ред. А.Н. Гузя. Киев: Наук. думка, 1989. Т. 5. 280 с.

3. Ильюшин А. А. Механика сплошной среды. М.: МГУ, 1978. 287 с.

TWO-DIMENSIONAL TIME-DEPENDENT WAVES IN AN ELECTRO-MAGNETO-ELASTIC SEMI-BOUNDED MEDIUM

V.A. Vestyak, D. V Tarlakovsky

A time-dependent problem of two-dimensional waves propagating from the boundary of a semi-infinite medium is described. Linear equations of coherent electromagnetoelasticity for isotropic conductors are used to formulate the problem. To obtain the solution, a small parameter method is used, where the small parameter is the quotient that characterizes the relation between mechanic and electromagnetic fields. The resolving system is formulated as a system of recurrent relations in terms of small parameter series' quotients. The null approximation is actually the solution for a purely mechanic problem. It is shown that in order to obtain higher order approximation one should first construct solutions for mechanical problems with specific time-dependent volumetric disturbances that depend on area's geometry. Solutions for half-plane and infinite medium with spherical cavity are used as examples of this approach.

Keywords: time-dependent waves, coherent electro-magneto-elasticity, two-dimensional problem, small parameter method, recurrent system.