Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Операционное исчисление и обобщенные ряды Лагерра

..pdf
Скачиваний:
0
Добавлен:
12.11.2023
Размер:
9.58 Mб
Скачать

а) т — четное

t mÇfe+i) П—1

 

 

r

t

 

 

д(Ж)

 

 

 

Л m(fe+i)

1

у

1 _______

Г

i

п

e~zdx — с~*1 t

» e dx

j\ tT

 

 

 

e * J x

 

m

G

r|™ (ft+ 1 ) + i j l

o

 

 

 

+

 

 

 

 

 

 

m—2

 

2icr

" - 1

2

 

tcoa—

. 2

V i

V I

 

e

+ n

*= 0

i d

г Г т ( й + 1 ) _

 

r=»l

*

n

 

 

 

 

Г ? !* ± «

 

_ ТС08± -Г

Г2 г.г .

j r

е

C 0 s l_-S- +

 

 

+

(f — 'Osin

 

dx;

 

6) m — нечетное

 

 

 

 

 

 

 

t

 

 

 

 

#

71

(6.3.4)

 

C

 

 

 

 

 

 

|am(+l,t)dt= 2

rfm(fe+l) j.,1 +

 

o *

 

*=°r| ~ — +lJ

 

4 -

,

’t ; 1

 

.

I

m (*+l)

 

i

y _____ e- ____ - \

~

e-*dx +

 

 

 

 

 

 

 

 

 

(6.3.5)

 

nt—l

 

 

 

 

 

 

— 1

А

Jcos-^f-

mCft+1)

_ tcoa?l>

 

 

 

e

с

»

e

« cos[ 5 +

+1-f2i Â2 f [ î ^ + i ] f

 

 

 

 

 

 

+

(f — t)sin

 

dx.

 

 

 

 

if

 

 

 

 

 

2- Для функции J Эт (—1> T)dr справедливы соотноше­

О *

ния:

 

(*_

 

V

(—-D**l)*f

Лл

 

 

 

] э ? ( - 1,т ) * = 2 г ^ (* + 1) + 1т +

 

 

 

О а

 

 

I—

 

 

 

 

 

 

 

 

 

 

Л— 1

 

п m(k+1)

л 7»(&~Н)

 

+*2

(-1)к

е* \| т

п

e~~xdtе

е~~х I х Л

 

+

 

e zd%

т * То r j ” (* + -1 ) ' + l j

 

 

 

 

 

 

т—2

tcos^2:

*

 

 

 

 

71—1

2 (_tyke

 

2лг

 

 

т

(*л m(A+7И(Л-Ь1)

 

 

+ * 2f l

Â2 r [ = ! t t î I + i ]|J Г

— TCO Б------

 

 

* е

т cos Г

-

+

 

 

 

 

 

L »

1

 

+

(? —

x ) s i n

d x ;

 

 

б) л — четное, тл — нечетное

 

 

 

m(H-l)

(6.3.6)

J 3m(~l» x)dx =

Z (~l)*t

п

 

if/»(i+l) +11 +

o »

 

fe=° r [—

+1J

 

я—:

( - l )

V

 

 

 

+ * 2

e~,*

+

■ Й

г[-т|*+ 1>+ Г ] ] т

 

 

 

 

 

(6.3.7)

771—1

П 1

g

.

tC O S —

p£

ш(д4

_ TCOS22L

? n r

i

V?

^

<-Dftg m

г

*t»+H

+ * 2

Z i

г [ ro(A+l)то(А-

 

 

 

cos

771

"+”

2

L

»

ï J T

 

 

A=0

*■->

 

 

 

 

 

 

+

( f - x

) s

i n

2^ J

d x ;

 

 

в) л — нечетное, m — четное

 

 

 

 

 

t

 

 

 

 

 

7П(*+1)

 

 

 

.1

»

 

 

2x ) rd xr m= ( É + l )

. л

 

" Ь

0

 

 

*“° Г[ n

+1J

 

 

 

n—1

» —2

,

2г+1 *

 

 

 

 

 

 

 

-

tcos*--------

t

 

_ rco^Zr+i

 

 

+ * x s

 

( -l)*e

m

f>

 

cos [«¥+

 

T[m(fe+1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.3.8)

 

 

 

 

+

(* — x )sin *^ ± lj <fr;

 

 

г) л — нечетное, m — нечетное

 

 

 

 

 

*

 

 

 

П—1

 

m(ft-fl)

 

 

 

 

 

 

( - 1)kt

n

 

 

 

î Эщ(

1, x)dx =

2

.rj-7»(fe+l)

u j

 

 

 

П n

 

 

i - 01

 

 

 

 

 

 

 

n—1

 

 

/» WUK+I)

 

 

 

 

_

 

j _ y

_

( _ i ) ^

Jr T

»

e7dr-f-

 

 

 

 

m iéorPtf“ii± i!.+ i

 

 

 

 

 

 

 

 

 

L ~ Î

' J]«

 

 

 

(6.3.9)

 

 

 

 

 

 

 

 

 

 

 

B—1

2=5 —1

.

2r+l

,

 

 

 

 

 

2

 

-

<«“ * ------

1

 

- « « s a

, 2r+i

,

+v2

V

 

(-!)»«

~

f s a t s

 

orp < î ± u +1] h

8

« * [ * —

+

*=o

2

 

 

 

 

 

 

 

 

 

 

 

 

+

(f — x) sinit

 

dx.

 

 

 

Таким образом, получены равенства для функции, пред­ ставляющей собой интеграл от функции Работнова дробного индекса, которые можно использовать как рабочие форму? лы при вычислениях значений этой функции.

При проведении расчетов значений функции Работнова и интеграла от нее удобно в формулах (2.4)—(2.9), (3.4)—>

(3.9) с помощью замены г = промежуток интегрирова­

ния [0, i\ свести к отрезку [0, 1 ].

По этим формулам была проведена серия расчетов.на ЭВМ БЭСМ-ЗМ с использованием стандартной программы по квадратурам типа Гаусса (Я. М. Жилейкин, ВЦ, МГУ, вып. 26, 1967). Сравнение..результатов с данными рабоТй [87] показало, что предлагаемые методы расчета функ­ ции Работнова и интеграла от нее весьма, эффективны и но точности не уступают упомянутым таблицам.

§ 4. Построение ядер последействия с заданными функциональными свойствами

Требуется сконструировать ядро последействия K (t—т), обладающее свойствами: а) ô-образности и б) двойственно­ сти. Под 6-образностью понимается выполнение следующих ограничений

Я (f) -з-оо при t -з- + 0,

(6.4.1)

К(f) -> 0 при t -> -f- оо,

00

J К (f) dt < оо.

О

Под двойственностью в широком смысле понимается, что как ядро K (t—т) уравнения Вольтерра II рода, так и его резольвента ô-образны, в узком — резольвента описывается той же функцией (но при других значениях ее параметров), что и ядро K (t—т).

Нетрудно видеть, что ядро К а (X, t) интегрального урав­ нения восстановления удовлетворяет условию двойственно­ сти в узком смысле в том и только том случае, если оно са­ мо является резольвентной некоторого ядра Вольтерра qp«(î)- Пусть Кл{К, t) — резольвента ядра q>«(f) и Ф»(р) — изоб­

ражение Лапласа функции q>«(f), т. е. фа(*)-^-Фа (р). Тогда изображение функции К а (Я, t) имеет вид

К а (К t)

Ф, 0>>

(6.4.2)

1-Х Ф« (р) *

 

t

 

Действительно, если в уравнении 17(f)— ^J

Ÿ«(f—т)17(т)4т=

 

_0

_

=rb(t) перейти к изображениям Лапласа U(p)—ЯФ„(р)С7(р) =

— _ _ Ф (р) — ■чфС*»)» то получим и(р)= 1|>(Р)+ Я j-^ф в (ТрФСр)-

t

Следовательно, H(f)=i|>(f)+A, |лгв (Я, t—х Ш Ф * и таким об-

0

разом формула (2) доказана.

t

Образуем уравнение 17(f)— р jlC« (^» * T:)U{r)dx—f(t).

о

Изображение резольвенты Ra (Я, р, t) ядра К а (Я, t) имеет

ВИД Да (Я, [I, t) "т_ 2

Фа ( р )

В«Да(Я, (.1, Î) —Да(Я"Ьр> t).

1р)Ф

Из этих простых расчетов следует, что ядра последейст­ вия могут задаваться своими изображениями вида Ф„ (р) или (2). В обоих случаях ядра релаксации определяются изображениями типа (2). При этом только во втором случае будет иметь место свойство двойственности в узком смысле.

Вели теперь учесть теоремы тауберова типа для интегра­ ла Лапласа, то для обеспечения свойства ô-образности не­ обходимо потребовать:

1)

р ф“ (р)

J °°’

при р +

00

(6.4.3)

 

1—ХФа (р)

\ 0,

при

+

О

 

2)

1_хф^(р)

= 0 (1)

при

р -> + 0.

 

В частности, если Фа(р) такова, что

рФ«

при р -»■ +

оо

при р-> +

(6.4.4)

 

0,

то условия (3) выполнены. Отсюда заключаем, что свойство двойственности в широком смысле является следствием свойства (а), поскольку ядра релаксации имеют изображе­ ние типа (2).

Таким образом, вопрос о построении ядер последействия, обладающих свойствами (а) и (б), сводится к подбору функ­ ций Фа (р), удовлетворяющих ограничениям (4). Наиболее простым является случай, когда функция Ф« (р) такова, что оригиналы ifW(f), отвечающие изображениям вида [Фа(р)]в,

могут быть вычислены сравнительно просто. Тогда ядро К а(Я, t), определяемое изображением (2), может быть задано рядом

*■«(*.*) = 2 хЧ + 1<*>.

(6.4.5)

*=о

 

Так, функция Работнова отвечает выбору Ф«(р)=1/рв. Дей­

ствительно, ограничения (3) удовлетворяются при 0< а < 1 .

Далее так как [Фа(р)]“

^ал—1

> то ряд (5) определяется

разложением (1.1) функции Эа (Я., t).

Осуществляя аналитическое продолжение по параметру а можно показать, что функция Работнова (уже как обоб­ щенная функция) описывается соотношением (2.3.2),' кото­ рое .подчеркивает тесное родство функции Работнова при малых значениях а с б-функцией.

Замечания. 1. Соотношения (3) можно рассматривать

.как обобщение требований (1).

2. Из (3) следует, что изображение Фв(р) должно быть функцией аналитической в полуплоскости R e p > 0.

3. Если jR(f) — резольвента ядра K {f), то e ^ R ip ) — ре­ зольвента ядра er9* K(t).

Это утверждение вытекает из следующей цепочки опе­ рационных соответствий:

ff.(X, t) -

Ф*(р )

f ) -

Фв (Р+Р)

1-ХФ. 0>) ’

1-ХФв (Р+Р) ‘

§ 5. Частные реализации общего принципа

Ядро Аа

(г, X, t). Положим Фа (р)=(Ур2+ г 2— р)а. Так как

[45] (формула № 11.6)

 

 

« -7 - сГ. (rt)-i-(K Р2+ Г 2 — рУ ,

(а > 0),

то в силу (4.5) изображению (4.2) отвечает ядро

 

во

 

 

A .(r, К t) = - 2 ^ 2 +

l)(X i-)*J.tt+i)(rf),

 

*1=0

 

 

ограничения (4.4) удовлетворяются при 0 < а < 1 .

 

 

1 _—

Ядро Ва(г, X, t). Положим Фа (р) =^ГС

р . Так как [45]

(формула № 11.42)

 

 

j e

( f ) 2 J.-i(2VTr),

(а > 0,

г > 0),

то в силу (4.5) изображению (4.2) отвечает ядро

В . (Г, х, f) =

Xk

g—1( r(ft+D

 

(* + l)2

Замечание. Функции a y -J -a(rt)е~&; \—^ г Ja-i(^\/~tr)e~^

при 0 < a < l могут рассматриваться как ядра последейст­ вия, при этом соответствующие им ядра релаксации имеют вид

e ^ A .( r , М ); e~ ^ B a(r,l,t).

§ 6. Учет запаздывания в явлениях последействия

Представляет интерес изучение ядер последействия, от­ вечающих случаю, когда максимальный эффект последей­ ствия обнаруживается с некоторым запаздыванием (на ве­ личину т > 0 ) относительно начала отсчета t= 0. Этой мо­ дели отвечает ядро последействия К а (X, t), эффект 6-образ­ ности которого смещен из точки 2=0 в точку 7=т. Она отражает те явления, когда временной (наследственный) характер упругих взаимодействий обнаруживается спустя некоторый промежуток времени. Здесь следует различать два случая:

а) ядро последействия имеет изображения вида

ф (р )

Фа(р)а-'Р или 1_ Хфд-(- )- ег** ;

б) ядро последействия имеет изображение вида

Ф,(У)е tp

(6.6.1)

1—Хфа (р)е~*р'

Во всех указанных случаях ядра релаксации имеют изоб­ ражения вида (1), т. е. характер запаздывающей релаксации не зависит от типа запаздывающего упругого последейст­ вия.

Схема построения ядер последействия, изложенная в § 4, без существенных изменений распространяется на вариант с запаздывающим последействием, при этом требование

Ô-образности накладывается на изображения типа (4.3) или (4.4), но не на (1). Рассмотрим некоторые частные реализа­ ции этой схемы.

Ядро D. (Я, t) (дискретный случай). Положим Фв (р) = = е-чр H-ô(t— а), (а > 0 ). Тогда ядро D a (Я, *), отвечающее изображению вида (4.2), представится рядом

00

=2 Я*5(t — лк — а). *=о

Уравнение Вольтерра с ядром ô(t— а) и резольвентой D„ (Я, t)

вырождается в разностное уравнение U(t)— %U(t— o)=/(t), m(t)

решение которого имеет вид U (t)= /(t)+ 2 **« » — ак — а), А==0

где ira(f)=inax{fe; afe+ aC f}.

Уравнение Вольтерра с ядром Da(Я, t) и резольвентой

 

m(t)

Dü(k+ n , t) вырождается в уравнение 17(0—

U(t— ak

a)= f(t), решением которого является сумма

m ( f)

04*) — f(f) + | * 2 +{*■)*/ (t — ak — a). fe=0

Замечание. Необходимо иметь в виду, что U{t) и f(t) равны 0 для отрицательных значений аргумента.

Ядро Н« (Я, т, t) (непрерывный случай). Положим Фа(р)=

Ч ^ )Ч -г)%-?т£г-=*■<'>•

Тогда ядро Не, (Я, т, t), отвечающее изображению вида (4.2),.

имеет вид

 

 

в*+в—1

 

_о*

 

 

 

-

Н а O ', i , t ) = —

e

2

r(ajfe+a)

 

 

А=0

И Л И

 

 

 

H .(X ,T ,f) =

- 2- e

 

г Эь ( * , ? ) .

Эффект б-образности в точке f = т > 0 как функции q>«, (f), так и функции Я .(Я , т, t) обнаруживается тем острее, чем больше величина параметра а. В пределе при а-*-оо функ­

ция

вырождается в функцию ô(£—а), а На

т, t)

В Da (Àj t)»

 

 

Ядро я« (Л,, т, t). Положим Ф«(р) —

Tj(t—х),.

0< а < 1 , т > 0,

где т\(t) — единичная функция Хевисайда. Тогда ядро За (Л, х, t), отвечающее изображению вида <4.2), представит­ ся рядом

Я. х, г) = 2 Х*(* rfcfe+Ii------~ ftx — х)-

*=0

Г(аА+а)

 

Отметим следующие важные свойства функции Я« (X, х, t).

1). Для каждого конечного значения

t функция:

За (À, х, t) представляется конечной суммой.

 

Пусть n x^ f^ '(H + l)x, тогда

 

 

д

 

 

Яв (^» х» t) — 2 х*

Г(аА+а)

 

А=0

 

 

2) Будем говорить, что функция t*~1 имеет в точке t O'

степенную особенность порядка л, если

^ а <

При любом фиксированном

а(0< а < 1 ) функция Яо (^»

х, г) имеет конечное число степенных особенностей, точнее,

есл и ^ г^ ^ а<

, то

Я« (Л,

х, t) имеет ровно п

степенных

особенностей в

точках

£=х,

2т, . . . , тгх, порядки

которых

равны соответственно п, п— 1 , . . . , 1.

3) Интеграл от Я« (Л, т, t) как функция верхнего предела

непрерывен по t

для любого 0< а < 1 и определяется выра­

жением

 

 

 

 

J я« (*, х, ï) =

2

"X~£V fe+a+ij—

’ если

nt< t< (n_r 1)t-

Замечание. Ясно,

что перечисленные

свойства могут

быть без труда распространены на ядра вида

a-î*D a (A, t),

е-^Н а (А, т, £),

е~р<Я» (Л, х, £).

Обширный обзор литературы по функциям Миттаг — Лефлера, род­ ственным функции Эа (Я, 0» дан в работах [22, 38].

Впервые функция Эа (Я, t ) была введена Ю. Н. Работяовым [85]. На

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

Метод интегральных операторов Э£ , содержащих в качестве ядра функции Эв (Я, £), получил дальнейшее развитие в другой работе

Ю. Н. Работнова [86].

Асимптотические формулы для функции Ю. Н. Работнова были изу­ чены Б. Д. Анниным [18]. Расшифровка получаемых операторных выра­ жений, а также всевозможные приложения метода интегральных опера­ торов Э* к решению задач наследственной ползучести анизотропных и

неоднородных сред проделаны М. И. Розовским [88, 89], Ж. С. Ержанов [51] приложил метод интегральных операторов к вопросам механики горных пород. В дальнейшем аппарат этой теории был применен для ма­ тематического описания явлении складкообразования в толще горных по­ род [52]. Интересное исследование свойств рассматриваемой функции с привлечением операционного исчисления, поиски путей ее дальнейших приложений принадлежат Я. В. Быкову й А. И. Боташеву [24]. Отметим также, что ими указывается на возможность табулирования этих функ­ ций.

Вопросам по1М1роения интегральных соотношений для функций Эа (А,

О и интеграла от нее, а также конструирования ядер последействия с за­ данными функциональными свойствами, в частности со свойствами, близкими к свойствам функции Э а (Я, t ) , были посвящены работы [10,

11]. Ю. Н. Работновым с сотрудниками [87] приводятся таблицы функ­ ции Эа (Я, t ) и интеграла от нее.

Обширную библиографию но затронутым вопросам можно найти в работе [86].

Соседние файлы в папке книги