книги / Операционное исчисление и обобщенные ряды Лагерра
..pdfа) т — четное
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].