Векторное преобразование Фурье с разрывными коэффициентами и его применение в теории упругости Текст научной статьи по специальности «Математика»

УДК 539.3+517.44

ВЕКТОРНОЕ ПРЕОБРАЗОВАНИЕ ФУРЬЕ С РАЗРЫВНЫМИ КОЭФФИЦИЕНТАМИ И ЕГО ПРИМЕНЕНИЕ В ТЕОРИИ УПРУГОСТИ

© 2011 А.А. Малышев, О.Э. Яремко1

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

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

1. Предварительные сведения

В предлагаемой работе используется метод векторных интегральных преобразований Фурье. Данный метод широко используется в различных отраслях науки: математическая физика, комбинаторика, теория сигналов и многих других. При использовании метода векторных интегральных преобразований Фурье, дифференциальные уравнения и граничные условия задачи переводятся в уравнения и условия для трансформант, после нахождения которых искомое решение дается формулой обращения для рассматриваемого интегрального преобразования. Метод векторных интегральных преобразований Фурье эквивалентен методу собственных вектор-функций, предложенному А.Ф. Улитко [1], но вместе с тем он может успешно применяться к решению задач теории упругости в кусочно-однородных средах. Более детальное изучение теории интегральных преобразований Фурье с кусочно-постоянными коэффициентами в скалярном случае было осуществлено Я.С. Уфляндом [2; 3], Л.С. Найда [4], В.С. Проценко [5; 6], М.П. Ленюком [7-9]. Векторный вариант метода разработан О.Э. Яремко [10; 11].

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

■^Малышев Алексей Александрович (dovmaa@mail.ru), Яремко Олег Эммануилович (yaremki@yandex.ru), кафедра математического анализа Пензенского государственного педагогического университета, 440026, Российская Федерация, г. Пенза, ул. Лермонтова, 37.

являются динамические задачи теории упругости для неоднородных тел. В этом случае коэффициенты Ламе являются функциями координат, определяющими поле упругих свойств тела. Метод векторных интегральных преобразований Фурье относится к аналитическим методам решения задач теории упругости. Обширный список работ по использованию интегральных преобразований в задачах теории упругости приведен в монографии Я.С. Уфлянда [2].

2. Векторные преобразования Фурье с разрывными коэффициентами

Рассмотрим векторную задачу Штурма - Лиувилля [12] о конструкции ограниченного на множестве 1+ нетривиального решения сепаратной системы линейных обыкновенных дифференциальных уравнений с постоянными матричными коэффициентами

d 2 + ^E + dx2

ym — °j

q2

4.m

A2E + rm, m — 1,n + 1

(1)

по краевым условиям

^(aii + a250J dX + (в01 + a27ii)^) yi

o, Ilyn+iH L=TO <

(2)

x = l 0

и условиям контакта в точках сопряжения интервалов

((j + A2sji) dx + (j + A27ki)^ yk — ^ (j + a2 ^2) dx + (j + A27k2)) yk+i,

(3)

x — lk, k — 1 ,n, j — 1, 2.

где

ym (x, A)

/ yim (x, A)

\ yrm (x, A)

, I ym I — y1m

+... + y2r

m

1, n + 1.

Пусть при некотором A рассматриваемая краевая задача имеет нетривиальное решение

y (x,A) — ^2 0 (x - lk-1) 0 (lk - x) yk (x, A) + 0 (x - ln) yn+i (x, A).

k=i

В этом случае число A называется собственным значением, а соответствующее решение y (x, A) — собственной вектор-функцией. Здесь

_ aii, e0i, Yii, <*?i, aki, j, Yjki, j ak2, j 7k2, j Aj —

(j — 1, 2; m — 1,n +1; k — 1,n) матрицы размера r x r; для матриц

M = ( ekm + A27lm alm + ^lm mk ^2m + A272m a2m + A2^2m

— 1,2; k —1,

потребуем их невырожденность, т. е.

n

det Mmk — 0, A e [0, to ).

Матрицы Am и rm,m = 1 , n +1 — положительно-определенные [13]. Обозначим Ф„+1 (x) = e9n+ixi; Ф„+1 (x) = e-qn+ixi; q^+i = A-+i (А2Е + Г2). Индукци-

онными соотношениями определим остальные n пар матричнозначных функций ($fc, ^k), k =1 ,n :

,,k i \2xk ) d I (ak I \2„,k

[(ajfi + АДУ dX + (j + A27|i)] (Фk, *k) =

vk j_ \2xM A _i_ f »k \ 2

dx

[(aj2 + Aj dx + (j + A27k2)] (^k+1, *k+i), k = 1,n, j = 1, 2.

Введем также обозначения

„О , \2r0

О

Ф (А) =

О

Ф(А) =

(a0i + — + (e0i + a270i)

(a°i + А2^0^ + (в°1 + a270i)

Ф1 (x, A) Ф1 (x, A)

X = lo

X=lo

Qk —

Фk ^k

ф/k ^'k

i = 1, n +1.

Теорема 1. Спектр задачи (1)—(3) непрерывен и заполняет всю полуось (0, ж). Задача Штурма - Лиувилля r раз вырождена, т. е. каждому собственному значению А соответствует ровно r линейно независимых собственных вектор-функций, в качестве последних можно взять r столбцов матричнозначной функции

П

u (x, А) = ^2 в (x - 1k-i) в (Ik - x) Uk (x, A) + в (x - ln) Un+i (x, A),

k=1

00

Uj (x, A) = Ф^' (x, А) Ф-1 (A) - Ф2 (x, А) Ф-1 (A).

^ u1m (x, A)

ym (x, A) =

V urm (x A)

Двойственная задача Штурма - Лиувилля состоит в нахождении нетривиального решения сепаратной системы линейных обыкновенных дифференциальных уравнений с постоянными матричными коэффициентами

(a™ dx2 + A2E + rm) ym = 0, qm = ^ + 1^ m =1,n +1

по краевым условиям

(dxy* Ai + А^) +y* (a?i + A%) ^

yn+iN <

(4)

(5)

и условиям контакта в точках сопряжения интервалов:

\2х k

— y* уЛ ( в11 + А2^11 ak1 + А^1

dxyk,ykJ I eki + A272ki a2ki + A252k

d

в?2 + A27k2 ak2 + A25k

12

" dxyk+1, yk+1 11 dk _l \2^,k T i \ 2 xk } , x lk, k 1,n.

Решение рассматриваемой краевой задачи будем записывать в виде

(6)

У* (£, А) = £ о (£ - lk-i) в (lk - О у* (£, А) + в (li - О у? (£, А) + в (£ - ln) уП+i (£, А),

k=2

т. е.

1

1

у*(£,а) = ( у*т1 (е,л) ••• y*mr(е,А)),

\\У*т II = у/(Ут)2+--+1УГт)2,та = 1,П + 1

Теорема 2. Спектр задачи (4)—(6) непрерывен и заполняет полуось (0, ж). Задача Штурма - Лиувилля r раз вырождена, т. е. каждому собственному значению А соответствует ровно r линейно независимых собственных строк-функций, в качестве последних можно взять r строк матричнозначной функции

П

u* (x, А) = £ О (х - lk—1) О (Ik - x) и\ (x, А) + О (x - In) <+1 (x, A),

k= 1

u* (x, в) = (ф (в), f (в)) ^—1 (x, в ) ( E ) А—2,

т. е.

у** х а) = (и*1 (е, а) ••• u*r (е, а) ) ,j = ir

Наличие спектральной функции и (x, А) и сопряженной спектральной функции и* (x, А) позволяет написать на множестве 1+ векторную теорему разложения.

Теорема 3. Пусть вектор-функция f(x) определена на 1+, непрерывна, абсолютно интегрируема и имеет ограниченную вариацию. Тогда для каждого x G 1+ справедлива формула разложения

СЮ СЮ

f (x) = - nj уu (x, а) (Ju*(е а) f (е) de+

0 lo

n / о \

+ (7О1Л (lo) + °0f (lo)) + £ f (A) ^ (A) ^—1 (Ik,A) M—11 (A) x

k=1 A 1 /

Y2k1

Y2k2

°k

°21

Xk

°22

fk + 1 (lk )\_ / Ylk1 °kM f fk (lk)

fk+1 (ik w v Yik2 °kw v fk (ik)

Установленные теоремы разложения позволяют ввести прямое Fn+

) AdA.

и обратное

F—+ матричные интегральные преобразования Фурье на действительной полуоси с n точками сопряжения:

Fn+ [f] (а) = J u* (е, а) f (е) de+ + (y0±^1 (io) + xOf (io)) +

lo

где

+ Ё Ё (A) ,£ (A)) ^—1 (lk,A) M—11 (A) x k=1 1 1

x

Y2k1

Y2k2

°21 \ ( fk+1 (lk) A A Yl1 °П \ A fk (lk)

°22 / V fk+1 (lk) ) V Y12 °k2 J V fk (lk)

F

-1

(x)

1

nj

СЮ

J Au (x, A) f (A) dA

0

f (x),

f (A),

f (x) = О (lk - x) О (x - lk—1) fk (x) + О (x - ln) fn+1 (x).

k=1

(7)

(8)

С целью применения полученных интегральных формул для решения рассматриваемой задачи теории упругости приведем основное тождество интегрального преобразования дифференциального оператора

B

^ в (x - lj-1) в (j - x) (a?dX + j + в (x - In) (

A”+1 dx? + гП+1

Теорема 4. Если вектор-функция

n

f (x) = ^ в (x - lk-i) в (lk - x) fk (x) + в (x - ln) fn+1 (x), k = 1

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

fkm) (lk-1)

lim

х —1+0

f(m) fk

(x),

m

0,1, 2, 3;

k = 1, n + 1,

удовлетворяет краевому условию на бесконечности

lim

X——

u* (x, A)

dxf (x) - dxu (x,A) f (x)

0,

то имеет место основное тождество интегрального преобразования дифференциального оператора B

Fn+ [B (f)] (A) = -A2/ (A) - {(++ (lo) + «Of (lo)) -- (7O1 A1f// (lo) + ^01A1f/// (lo)) } - £ + ,ф) О-1 (lk, A) M-11 (A) x

k = +

21

+2 a22

+ik1 +1

+2 +2

fk + 1 (lk ) ffc + 1 (lk )

721 S21

722 S22

fk (lk )

/ k

fl (lk)

7i1 7ik2 +

sk

°11

12

Afc+1ffc+1 (lk )

A2k+1fl+l1 (lk)

A2kfk/l (lk)

A!fk//l (lk)

(9)

Доказательство теорем 1-4 проводится аналогично представленному в работе [11].

k

3. Динамические задачи теории упругости

Рассмотрим задачу о распределении напряжений в полубесконечном n +

n+1

+ 1-слойном упругом теле 1+ x R = {(x,y) : x G Z+,y G R}, где 1+ = U (li-1,li).

j=1

В случае плоской деформации вектор смещения у» имеет компоненты uj,vj, 0. Если согласно [14] ввести две функции напряжения у» (x,y,t) и ^ (x,y,t), определяемые соотношениями

дуч d^j d^j d^j

ui “Д I о , vj "о о ,

dx dy dy dx

d+ dy

то выражения для компонент напряжения принимают вид [15]

^jX

Aj^^j + 2^j

d2 у» + d+j dx2 dxdy

ajy

Aj^yj + 2Д»

d+j

dy2

d+j

dxdy

Tixy

}^i

d2Pi

2

oxoy

dVi + д2Ц dx2 dy2

(И)

где Ai, pi — упругие постоянные Ламе [16]. При этом уравнения движения будут удовлетворяться, если функции напряжений Pi и Ц выбрать в виде решений системы волновых уравнений:

dVi _ 2 Л d2^i

dt2 = CliЛ^i, Of2

c^iA^i t > 0, —ж < y < ж, /i_i < x < li

(12)

с нулевыми начальными условиями

Pi (x, y, 0) = 0, Ц (x, y, 0)

0,

d<Pi (x, y, 0) dt

0,

дЦ (x,y, 0) dt

0.

На границе тела прилагается изменяющееся со временем давление p (y,t). Если считать, что касательное напряжение равно нулю, то граничные условия принимают вид

<7ix = —Р (y, t), rixy = 0 при x = 0.

Считая непрерывными компоненты вектора смещения щ и компоненты тензора напряжений aix, Tixy, приходим к следующим внутренним граничным условиям, т. н. условиям сопряжения [16; 17]

ui ui+1, vi vi+1, ®ix Tixy Ti+1xy, x li.

Для решения задачи (10)—(12) применим по переменной y преобразование Фурье [18], а по переменной x векторное интегральные преобразование Фурье с разрывными коэффициентами из п.1. В образах Фурье по переменной y задача (10)—(12) примет вид системы уравнений

dVi _ 2 d2¥i 2 С2- д2ц

dt2 Cli dx2 Clit Pi dt2

с начальными условиями

2 d2i2i C2i dx2

с1Л2 i’i

t > 0, li_1 < x < li (13)

Pi (x,y, 0) = 0, Ц (x,y, 0)

0,

dpi (x, y, 0) dt

0,

d'tpi (x,y, 0) dt

0,

(14)

где Pi, Ц — изображения Фурье по переменной y функций напряжения

Pi

Pi (x,y,t) e>iydy,

Ц

Ц (x,y,t) ejydy

с граничными условиями

0"ix = Ai — Alt2 ‘Pi + 2pi

92yi

dx2

+

d'tpi

dx

-p (t, t) при x = 0,

I 0-cdPi d2ipi 2 i | n n

Tixy = Mi ( 23t~ox HP — t Ц I =0,x = 0,

dx2

с внутренними условиями сопряжения

Px + jt^i = + jt'^i+ ^ jtPi — lip = jtPi+i —

при x = li

= li

Ai^ — A’P'P’ + ЦX + =

\ d2Pi+i

Ai+i dx2

Ai+it2Pi+i + 2Mi+i

d 2Pi+1 + ^дЦг+Л dx2 j dx )

при x = li

.-df д2ф,

М 2Я—---

i - e2ii) = Mi+н т

1 х J s dx dx2 Обозначим c = max {cii, C2i}.

^d<p.

i+1

d2ip

i+1

dx

dx2

- eVi+O x = li. (16)

Применим к задаче (13)—(16) векторное интегральное преобразование Фурье с разрывными коэффициентами, определенное формулами (7)—(8). Положим в системе уравнений (1)

2 _ ( c21 0

r 2 A l 0 c2

i2

Г2 = ( (c2 - c?o e2

(c2 - 4) e2

в краевых условиях (2) будем считать

«11 =

0

2j>1e

2j>1e о

e?1 = -

Y11 =

Л1 + 2^1 0

0

"М1

А-2Г^ / Л4 0 Г I 0 М1£2

Л1 + 2М1 0 -2 хО = ( 0 0

0 -М1 J А- , д11 =10 0

в условиях сопряжения (3) положим

«и =

0

2jMfce 0

2j>fc e

0k

( Лк + 2Mk 0 \ л-2р2 ( Лкe2 0

1 0 -aJ k k v 0 Mke2

Y11 =

Лк + 2Mk 0 0

Mk

A-2 Xk =

Ak , °11 =

00

00

к _ I 0 2jMk + 1e

“12 i 2jMk+1e 0

ek2 = -

Yl2

Лk + 1 + 2Mk + 1 0

0 —Mk+1

Л1+1 + 2Mk+1 0

д—2 “p2 l Л1+^ 0

Ak+lik+1 — t 0 Mk+1e2

0

"Mk+1

A-2 Xk

Ak+1, °12

00

00

«2ki

1 0

0 -1

eki

Y2ki

00

00

xk

°2i

00

= 1, 2.

0 me

me 0 )1 ,2i ^ 0 0 у ’ ”2i ^00

Переходя к образам Фурье по переменной x и учитывая тождество (9), запишем задачу (13)—(16) в виде

d2 ( f \ 2е2 ( f ^ — f

dt2 у 1

—c e

1

+

jp (e, t)

0

l)(e,n0) = 0, 1(1 )(e,Y,0) = 0,

(17)

(18)

здесь принято обозначение

(n,e) = F

n+

(П):

0

(j) = ^6 (lk x) 6 (x 1k-l) ( Jk ) +6 (x 1n) ( Jk )

Приведем решение задачи (15),(16)

* wo= i'sin-T’) r p«-»I*.

n2

sin | \/c2£2 + n2

J J 4 ^ Jo л/c2^2 + n2 ^ 0

Возвращаясь к оригиналам по переменной у, получим

г' /• у+с('-т)

* J'n^W2

0 Jy-c('-r)

J0 I П\ (t - T) -

2 (У - s) W P (S,T)

0

dsdr,

где Jo — функция Бесселя [19]. Применив в заключение обратное интегральное преобразование Фурье F-+ из (8), найдем функции напряжений

U: j (x,y,t)=

1 /*' /• y+c('-T)

р/2к

o Jy-c('-T)

H I x,\ (t - T) -

2 (У - s)2

P (s,T)

dsdT, (19)

0

где

H: (x, z)= nu: (x, n) Jo (nz) dn•

o

Наличие выражений (17) для функций напряжения позволяет по формулам

(10),(11) найти компоненты вектора перемещений u:,v:, 0 и компоненты тензора напряжений <7:х,<х:у,пХу.

Литература

[1] Гринченко В.Т., Улитко А.Ф., Шульга Н.А. Динамика связных полей в элементах конструкций. Электроупругость. Киев: Наукова думка, 1989. 279 с.

[2] Уфлянд Я.С. Интегральные преобразования в задачах теории упругости. Л.: Наука, 1967. 402 с.

[3] Уфлянд Я.С. О некоторых новых интегральных преобразованиях и их приложениях к задачам математической физики // Вопросы математической физики. Л., 1976. С. 93-106.

[4] Найда Л.С. Гибридные интегральные преобразования типа Ханкеля-Лежанд-ра // Мат. методы анализа динам. систем. 1984. T. 8. С. 132-135.

[5] Проценко В.С., Соловьев А.И. Некоторые гибридные интегральные преобразования и их приложения в теории упругости неоднородных сред // Прикладная механика. 1982. T. 13. № 1. C. 62-67.

[6] Проценко B.C., Головченко А.В. Обобщенное интегральное преобразование типа Фурье-Лежандра // Мат. методы анализа. Харьков, 1982. № 6. С. 26-28.

[7] Ленюк М.П. Гибридные интегральные преобразования (Бесселя, Лежандра, Бесселя) // Укр. матем. журнал. 1991. Т. 43. Вып. 6. С. 770-779.

[8] Ленюк М.П. Гибридные интегральные преобразования (Бесселя, Фурье, Бесселя) // Матем. физика и нелинейная механика. 1989. Вып. 12(46). С. 68-74.

[9] Ленюк М.П. Интегральное преобразование Фурье на кусочно-однородной полупрямой // Изв. вузов. Сер. Математика. 1989. T. 4. С. 14-18.

[10] Яремко О.Э. Матричные интегральные преобразования Фурье для задач с разрывными коэффициентами и операторы преобразования // Доклады РАН. 2007. T. 417. № 3. С. 323-325

[11] Баврин И.И., Матросов В.Л., Яремко О.Э. Операторы преобразования в анализе, математической физике и теории распознавания образов. М.: Прометей, 2006. 292 с.

[12] Судаков Р.С. Простые методы прикладной теории матриц. М.: РХД, 2005. 450 с.

[13] Брейсуэлл Р. Преобразование Хартли. M.: Мир, 1990. 584 с.

[14] Снеддон И. Преобразование Фурье. М.: Иностр. лит., 1955, 668 с.

[15] Ландау Л.Д., Лифшиц Е.М. Теоретическая физика: Теория упругости. М.: Наука, 1987. T. 7. 247 с.

[16] Тихонов А.Н., Самарский А.А. Уравнения математической физики. 7-е изд. М.: Наука, 2004. 743 с.

[17] Ахтямов А.М., Садовничий В.А., Султанаев Я.Т. Обратные задачи Штурма-Лиувилля с нераспадающимися краевыми условиями. М.: Изд-во Московского университета, 2009.

[18] Оболашвили Е.И. Преобразования Фурье и его применение в теории упругости. Тбилиси: Мецниереба, 1979. 230 с.

[19] Снеддон И.Н., Бери Д.С. Классическая теория упругости. М.: Вузовская книга, 2008. 215 с.

Поступила в редакцию 6/V/2011; в окончательном варианте — 6/V/2011.

FOURIER VECTOR TRANSFORMATION WITH DISCONTINUOUS COEFFICIENTS IN THE THEORY OF

ELASTICITY

© 2011 A.A. Malyshev, O.E. Jaremko2

The article deals with the description of the method of Fourier vector transformation with discontinuous coefficients. The dynamics of the theory of elasticity illustrates the technique of method in solving problems of mathematical physics of heterogeneous media.

Key words: Fourier transformation, theory of elasticity, the Sturm — Liouville vector.

Paper received 6/V/2011.

Paper accepted 6/V/2011.

2Malyshev Alexey Alexandrovich (dovmaaamail.ru), Jaremko Oleg Emanuilovich

(yaremkiayandex.ru), the Dept. of Mathematical Analysis, Penza State Pedagogical University,

Penza, 440026, Russian Federation.