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

книги / Оболочки и пластины

..pdf
Скачиваний:
21
Добавлен:
19.11.2023
Размер:
71.66 Mб
Скачать

 

=

 

 

 

 

 

 

А,,

 

 

 

 

 

 

 

 

 

 

 

 

Кп = —~ ~ М0РЛ4,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Kl2=

12

Л4(А,— 1),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J2

v.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/Ci3= — h3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь введены обозначения

 

 

d*w0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vavam = &11

 

 

 

d*w0

 

 

dAw0

 

 

 

 

 

 

дх4 +

(2 b12 + SB) dx2dy2 + b22

~dy4

 

 

 

 

 

 

 

 

 

 

V*<Po =

 

^Фо

Д2фр

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*

дх2

ду*

 

 

 

 

 

 

 

 

 

 

^Фо

 

 

* ч

 

а2ф0

 

 

 

 

 

а2ф0

 

 

 

 

 

 

дх2

 

 

ду2

 

ду2

дх2

 

дхду

 

дхду

 

 

Уравнение (3,22,20) дает связь между

нагрузкой q — vt

и w (t) . Перейдем

к определе­

нию параметров 'решения

(X,

р, ткр, Р)

в случае q = q o = const.

 

 

 

 

 

 

Из уравнения

(3,22,19) получим

 

 

 

 

 

 

 

 

 

 

 

 

' KI ^ /C 2+

^3 +

^4 =

0

 

 

 

 

 

 

 

 

 

при

£ =

0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3,22,21)

K i Ф Къ =

0

 

 

 

 

 

 

 

 

 

 

 

 

 

при

t =

0 0 ,

 

' Ki 4- K ^ aXk+

К3е~2аЧ 4 K«e_ 3 “ x* 4 К ь 4 K*e~aXk+ Кые~аХк 4- К,е~2а4 4-

 

 

 

 

 

 

 

4 - K n e ~ atk4

4 - K u P ~ 2aXkVk = 0

при

t =

tft.

 

 

Таким образом, имеем систему уравнений, которая

вместе

с

условием,

что X—р = 1 ,

полностью определяет значения параметров, входящих в решения.

 

 

 

 

Из

первого уравнения системы

(3,22,21) .получим

упругое

 

решение

 

 

Я=

h3

Г

 

~

 

12

 

-

 

12

.

12

 

9

 

12

 

9 -'I

 

 

[ vV »* -

 

 

 

■м»в -

 

М*А-

 

V№A -

 

-jjT ^

фоВ J • (3-22-22)

Упругое решение в первом приближении имеет вид

 

 

 

 

 

 

 

 

 

Я=

2я4

(ВцВ22В^2)

 

 

Зяа

(xi Ф х г) (В ц В 22 В *2)

 

 

 

В ц

(#п Ф С ф В 22)

 

 

В11

 

ц ф С ф £ 22)

 

 

 

 

 

 

 

 

 

 

 

 

 

л4(5Ц +

25

85 -ф- 5 la)

(xi - f x-J2 (5ц 5 аа— 5f2)

]«■

(3,22,23)

 

 

[

 

 

125,:

 

 

 

^11 (^11 Ф С ф B22)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где q= q

а4

 

. ___ f_

 

 

kia*

 

k2b2

 

 

 

 

 

 

 

 

 

X 2 =

 

 

 

 

 

 

 

 

 

5u fc4

 

5 _

h

Xl

h

 

 

 

 

 

 

 

 

 

 

 

Подставляя значения £ г:> ,в .(3,22,23), получим

У = 22,345£3— 3,352 (хх ф Ха) £а ф [0,115 (хх ф х2)а 4* 22,075] g .

dq

Из условия —— = 0 найдем границы устойчивости

«Е

= 67,035ga— 6,704 (х,, -ф- х2) I ф 0,115 -ф- х2)г ^ 22,075 = 0.

<*£

Таким образом,

 

3,352 (х, -f х 2) ± У" 3,6769 (х, + х2)2— 1479^6476

 

Sl- 2 —

 

67,035

 

 

Например, для сферической

оболочки

с кривизной XI = K2=|20

имеем £ц=1,04,

£ в =2,96,

им соответствуют

верхняя

и нижняя

критическая «нагрузка

к = 94,34, q ni<='14,82. Из

упругого решения

найдем величину хлопка р, которую будем

считать равной

величине

хлопка в упруго-вязкой задаче, поскольку хлопок, как правило, происходит за столь малый промежуток «времени, что ползучесть за время хлопка не проявляется. Таким образом,

Р = Рву — Ру •

Из второго уравнения системы (3,22,21) находим

X3 4 - № D ± 4- X D 2 + D 3 == 0,

(3,22,24)

где

 

Ву|фо + Л40(А + ЗЬр)

Di =

М 0В

у 2Фо (А + Щ ) + М 0 (2Ар + Щ * ) - V

М 0Ь

V jfto (А Р 4 - W

+ М0(Ар* +

Bp*) -

v 2V ^ ftp +

(1 — у)

 

D3 = -

 

 

 

 

 

 

 

М 0В

 

 

 

 

 

 

 

 

 

Подставляя цифровые

значения

в

(3,22,24), получим (значение Р снимается

с

упругой диаграммы)

 

 

 

 

 

 

А,3 — 1 ,276А* - f 0 ,406Я. — 0,274 = 0.

 

 

Решая по формуле «Кардана,

найдем

для

X одно действительное значение X ='1,136

и

два комплексных.

 

 

 

 

 

 

Таким образом, имеем

Х=«1,'136;

ц=А,— 1=0,136. Из третьего

уравнения системы

(3,22,21) найдем ткр в зависимости от

q0.

 

 

 

 

 

При ^о=17 кг имеем уравнение

 

 

 

13

— 2,9е

If5Tfe — 54,8е

°'5т<!4>44,4е Tft_

 

 

— 34,9Tfte_tfc 4 -1 9 ,ЗтАе~0,5т* = 0,

 

 

решая которое, графически найдем ткр= 0,1.

 

2. Цилиндрическая оболочка

Рассмотрим устойчивость [154] сжатой в осевом направлении замкнутой круговой ци­ линдрической оболочки, изготовленной из ортотропного стеклопластика. (рис. 3.75). Кром­ ки оболочки закреплены шарнирным образом, так, что должны удовлетворяться следующие условия:

О, y,t) = w(x = L, у, t) = /о,

d*w

d*w

дх2 х = 0 , у, t

~д&

а*ф

а2ф

д:=0, у , t

ду2

= о,

x = L , у, t

— Р>

x = L , у, t

д ф = 0

при х = О, х = L,

дхду

 

W (х, у, 0 ) = w0 (х, у),

ф(*. у, 0) = ф0 (х,у).

Будем считать, что кривые релаксации напряжений в материале обо­ лочки могут быть аппроксимированы функциями вида

<*ц( 0 = М 1 Уu— ( 1 —

Опыты на релаксацию растянутых образцов стеклопластиков показы­ вают, что некоторые параметры Yij и atj для образцов, по-разному ориентированных к основе ткани, могут быть с достаточной для инже­ нерной практики точностью приняты одинаковыми:

Yu = Y22*» а1Х= а12 — а22 — а -

Так, для полиэфирного стеклопластика (ПН-1, Т-1) имеем Оц = а12 = а22 = а = 0,5,

Yu = Y*a=-g-’ Y12 = M'l’Yii> Y= 3Yu.

(3,22,25)

*ц = 1,49-10»-^-, 622 = 1,14-Ю5- ^ - ; bn — bn = 0,18-10®

6 = 0,13-105

Примем

Rn = R22 = y№~at

Ri2= H i*Я н

R = 3Ru,

(3,22,26)

а в преобразовании Лапласа

 

 

 

=

=

Зуа

(3,22,27)

где \i\ есть неизвестный параметр. Представим прогиб w(x, у, /) в виде

w (х, у, 0 = ш0 (л:, у) о»!(0,

(3,22,28)

где w0(x, у) есть упругое решение, а W\(t) выбирается в форме

w1(t) = 'K— \ie-at.

(3,22,29)

Параметры X и р, подлежат определению. Очевидно, что X—ц=1. Функ­ цию напряжений возьмем в виде

ф(*, У, 0 = Фо(*, y)-q>i(0,

где ф0 (JC, у) есть упругое решение задачи, а <pi(0 при заданной a>i(f) найдем из уравнений (2,24,16) и (2,24,12). Таким образом, задача сво­ дится к построению упругого решения w0, <р0 и определению параметров jxi, X и ji. Полагая в уравнениях (2,24,12) и (2,24,16) &i= 0 и <7=0, полу­

чим уравнения для цилиндрической оболочки, к которым необходимо добавить условие замкнутости

2яг

(3,22,30)

до

ду

1

д2фо

+

lhwo-

^11^22 — ^12^2L

дх2

 

 

Для построения упругого решения опустим в уравнениях (2,24,12) и (2,24,16) аддитивные члены, зависящие от времени (что сведется к опу­ сканию звездочек), и добавим к ним условие (3,22,30). Примем, что прогиб оболочки можно аппроксимировать функцией [94]

w0= /о + fi sin

sin - f - + /. sin2

,

(3,22,31)

где f0— постоянная, учитывающая радиальное обжатие концевых диа­ фрагм; fi — стрела выпучины при х = Ц 2; /2—параметр, отражающий симметричную (/2 = 0) или несимметричную (/г¥=0) форму выпучивания цилиндра; m — число полуволн вдоль образующей; п — число волн вдоль направляющей. Заметим, что при таком выборе w граничные условия удовлетворяются в «среднем». Внося (3,22,31) в (2,24,16) и проинтегрировав полученное линейное уравнение относительно <ро, по­ лучим частное решение для функции напряжений

<р0 = Dxcos2

-f D2C O S

-^ L -)_ £)3sin

R

sln-^- -f

 

1 2

 

 

R

R

 

+

Di sin

R

sin — ------ ,

2

(3,22,32)

 

 

R

\ • > r

•где

D2 = —

(Ьг1Ь22 — b12b21) Л202

,

 

 

 

 

 

 

t?22

 

 

 

 

 

 

u Z

 

 

 

D3 =

 

 

(Р1Ф22 ----612621) Cl? 2,

 

 

 

 

 

 

 

g l ( 0)

 

 

 

 

 

 

 

 

 

 

A

=

 

fi

( * i A z —

^12^21) h2Q2t,i f ------------ ^

 

 

 

 

£ г(°)

 

 

 

 

 

 

V

T)

/

 

Здесь введены безразмерные величины

 

£2, £о, Л

ЧJ

и обозначения

§ и С\I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5-

 

/ l

 

У

f 2

 

V

/о .

 

 

 

 

 

 

t l

 

. h » Ь2

, h > гО —

, h у

 

 

 

о

? r

,

у

_

 

я2/*

у

л

p

т я г

 

т я г

»

*

,

Л

 

г

0

п

~ У Р

 

"

 

ОцП

 

 

 

 

 

 

nL

 

L

 

£ = “гт“ (611^12 b12b21 46612)| СА= (рхФхъ b12b2^)'r

а функции gk (6) даны формулами

£ 1 (6) =

8161104 +

9«е2 +

&22,

Ы е) = М

4+ * 6

2 + &22,

ё з (0) =

М

4 + 2М 2 +

Ь22 + 4626.

Нетрудно видеть, что ср0 удовлетворяет граничным условиям в «сред­ нем».

Для определения параметров /о, /г воспользуемся методом Ритца. Найдем выражение полной энергии системы Э, состоящей из по­

тенциальной энергии деформации срединной поверхности

"

-

тО о -

£

П

[

Ч

З

)

, +

Ч

З

)

-

*

К-£кУ\дхду'

потенциальной энергии изгиба

“ - - - я -

L2nR

* \ дх2 )

12 ду2

дх2

+Ч^у-

оШо

1 L

Г Г

Гh..

+ 2b..—

 

+ib(£k)]ixiy-

работы внешних сил

L2nR

0 0

После соответствующих преобразований получим

5 = nL.fi

4Сгг 3

[p46u (32Z)2 + 81D2 + DD+ b22n*(32D2+ D\ -±D\ + 4P2+;i6PD2) -

- 2B12p2na (9D§ + D\) + -g - p2n2 (9D2 + Щ j +

+ ^ ( ^ f P4/? + 46uP^ + b ^ f \ +

n*f\ + 2np25/2) _

nPLh |^2622n2P + - ^ p 2(/M-2/2)]

Cxr

Составляя уравнения метода Ритца

 

дЭ = о,

-2 L = О

dfi

[df.

и учитывая условие замкнутости (3,22,30), получим три уравнения для определения трех параметров fu /2 и f0:

- J -$ h

{b22+

м

4) C

? + e r , H

- ^

+

^

S

J

-

160* *tl622

 

 

 

 

 

*11 V g l (0)

 

g 2 (0)

У

 

Cl

f 1

_|_

864fcU Л г

-L-

rCl9!'

+

(9) 11

 

= P , (3,22,33)

462ie2

V1

+

ы

е ) / £

'

6uft(0)n

 

12bn&

 

 

 

v

 

 

 

р _

 

еал

 

 

 

 

 

 

 

 

40аТ1*п

 

1

1

к

(3,22,34)

 

 

C l

 

Схб»

\

еатА

 

 

 

 

 

[ - ( • i66?,e2

2^11^2 (в) /

 

gi (в)+

g2(в)

 

 

 

_ М

й р

+ Л .С2 __ 1 _ Сг==Со>

 

 

(3,22,35)

 

 

 

Ci

 

8

2

 

 

 

 

Из

уравнений

(3,22,33) —(3,22,35)

можно построить

зависимость Р

для

различных

значений

параметров

волнообразования т) и 0.

 

Результаты вычислений для материала (3,22,25), (3,22,26) пред­ ставлены на графике (рис. 3.76) в случае 0= 0,5 и -л=0,22. Определяя критические нагрузки, получим для верхнего значения

Рв = 0,558,

а для нижнего

РИ= 0,171.

Из анализа значений вторых производных от да0 и <р0 по х и у в точке (x=L/2, у) легко видеть, что для m-нечетного решение задачи о влия­ нии времени на прогиб получается без затруднений. Рассмотрим случай /n-четного. Преобразуя по Лапласу уравнение (3,22,28) с учетом форму­ лы (3,22,29), имеем

w'(p) = w0( - --------

(3,22,36)

Ч р

р + о /

Внесем (3,22,36) в уравнение (2,24,12). Решая полученное уравнение относительно <pi, получаем после инверсии

<Pi = Ах-|- A1e~at + A3e~'>ai + A4e~2al + Abte~at,

(3,22,37)

где

Л = А ( * 1 » . “'о, Фо. К ц. Hi).

Внося ср* и

w\ в уравнение

(2,24,16)

и

учитывая (3,22,27), получаем

после инверсии преобразования Лапласа

 

 

Я 1 + К

2е - “ +' K3e~‘/,at +

Кле-'Ч^ +

Kse~2ai + K6e~™ +

K7te~^ +

 

 

+ K e t e - * *

=

Q ,

(3,22,38)

где Ki = Ki(bn, w0, фо, К ц, |xi), причем в К\, /С2, Кз, Кб и К7 входит ли­ нейным образом нагрузка Р. Используя метод коллокации и полагая

в уравнении

(3,22,38)

Р = 0,171, £i = 5,396, £2=2, получаем для

материа­

ла

(3,22,25),

(3,22,26)

значения Ц1=0,12, К= 1,242, (х=0,242.

 

на

Таким образом, величина рм имеет порядок коэффициента Пуассо­

[155] для

рассматриваемого материала (для полиэфирного стекло­

пластика ПН-1, Т-1,

v = 0,11; vi2 = 0,12 [153]). В дальнейшем

значение

jii = 0,12 сохраним для всех Р. Рассматривая прогиб для различных зна­

чений ^=const на участке 0<Р<Рь (£i = 0, £2= 0), получаем методом коллокации (полагая /= 0, 1, 2, 3, 4, с») систему уравнений, из которой

с помощью формулы Крамера найдем Х=1,2. На участке 0<Р<Рь из­ менение прогиба во времени при постоянных нагрузках будет выражать­ ся формулой

w (x,y t) = w0(x,y)(l,2 — 0,2е-°.5'),

т. е. прогиб может изменяться по сравнению с упругим на 20%, если /—>-оо. Предполагая нагрузку изменяющейся во времени

 

P = PoPi(t),

 

 

 

 

где P\{t) есть безразмерная

функция времени,

такая,

 

что

P j(0 )= l,

имеем

 

 

 

 

 

 

 

д2<Ро(*. У)

= Фх(0 = — P*Pi(t),

 

 

д у* x = 0 ,L

ду*

 

 

 

 

 

 

откуда 9i(f) —Pi(t), или, учитывая (3,22,37),

 

 

 

 

- ± t

-

±

t

-

L

t

(3,22,39)

Р — Ро {А -}- Л2е 2

А$е

12

-)-

1-f- A$te

2

).

Время критического состояния можем найти из условия dP/dw = 0, учтя, что прогиб задан в форме

 

W (х, у, t) =

w0 {х, у) (К — р,е

4 * )•

 

 

Полагая

оо, получаем уравнение для tкр:

 

 

 

 

А4е~‘ + Аье 2 *

1

Y

t

2

 

2

= 0 .

12

 

 

 

 

 

Это уравнение можно решить

приближенно,

представив

eat в виде

ряда

 

 

 

 

 

 

 

l + a t +

Ш

+ H W

 

 

 

 

 

2!

3!

 

 

 

Для Х= 1,242, ц = 0,242, Ро = 0,171, ш0=1,76 (получим кубическое уравне­

ние, из которого по формулам Кардана найдем

fKp = 7,058. Далее най­

дем

 

 

7,058

рн = 0,1374, wH= w0(1,242 — 0,2426

2 ) = 2,17.

Таким образом, учет линейной наследственности снижает значение кри­ тических нагрузок для оболочек из стеклопластика.

Упругое решение РВу = 0,558; с наследственностью Рвн=0,470. Упру­ гое решение Рну=0,171; с наследственностью Лш=0,137. На рис. 3.76 дан график зависимости нагрузка—прогиб: а) упругое решение; б) из­ менение прогиба во времени при t->оо и Р = const; в) с учетом наслед­ ственности при нагрузке, изменяющейся во времени. Заметим, что кри­ тические нагрузки оболочек, материал которых обладает линейной на­ следственностью, существенно зависят от режимов нагружения, увели­ чиваясь с увеличением скорости нагружения.

3. Применение табулированных сингулярных ядер в задачах

устойчивости упруго-вязких оболочек

В предыдущих пунктах при решении задач использовалась простей­ шая функция влияния — ядро релаксации в виде экспоненты. Однако, как уже отмечалось в § 23 гл. II, в начале релаксационных процессов наблюдаются весьма большие скорости, которые не могут быть описа­ ны экспоненциальной функцией или их суммой и довольно хорошо ап­ проксимируются сингулярными функциями влияния типа (2,23,17) или (2,23,18). Однако их применение в расчетах ползучести элементов кон­ струкций ограничивалось трудностью вычисления интегралов от них, появляющихся в решениях. Эти трудности в значительной мере снижа­ ются при использовании таблиц сингулярных функций влияния и их интегралов [156]. Покажем, что решение задачи об устойчивости гиб­ ких упруго-вязких оболочек может быть выражено через табулирован-

t

t

ные ядра R, К, ^ R(u)du,

j*K(u)du, где R и К даны формулами

6

о

(2.23.12) и (2,23,17). Задача сводится к решению системы уравнений (2.24.12) и (2,14,16) при определенных граничных условиях. Примем их такими же, как и в п. 1 для шарнирно закрепленных кромок обо­ лочки.

Приближенное представление прогиба во времени примем в виде

w(x, у, t) = w0(x,

w±(t) = X — \ie-<x”t — yQ(t — tk),

(3,22,40)

где

 

 

 

9 < < - у = { °

 

Считая возможным

представить функцию напряжений

в виде

ф(*>У, t) =фо(*>#)<Pi(0

и внося w*(p) в ур.авнение неразрывности

(2.24.12) , получим уравнение, разрешая которое относительно ф1*, по­

лучим зависимость,

аналогичную

(3,22,13), которую здесь

запишем

в виде

 

 

 

ф; = Qi (х, У) Ю * -

Qi (х, у) R* (w*y -

Q2 (*, У) vf[ + Q2 (х, у) R‘w\,

(3,22,41)

где

Q2 у) —

 

(Ь цЬ п - Ь М | \

Hr h ^

(3,22,42)

д*Фо

1

 

^Фо

 

 

^Фо

 

bil ~д*

* 2Ь

{bllh* ~ 612621“

466ls) д*д? +

622~ду*

Q l ( x , у ) =

 

(6ll622~

6 12 6 2 0

[ ( l a y )

~ ( дх2

ду2

)]

 

 

 

 

 

 

 

d2w0

*d2w0

 

,

^Фо

,

1

,

t ,

ч ^фо ' , ,

^Фо

611

д х *

+

2b

( и

22~

12 21 “ 6

l2) дх* ду> ф

22 ^

 

величины, определяемые по упругим решениям. Используя теорему умножения операционного исчисления [157], согласно которой, если интегралы

Ш = ] м )е п * < и и r 2 ( p ) =

] f 2 ( t ) e - p ‘ d i

6

о

сходятся абсолютно в области Rep>yi, то функция

П (Р ) = П ( Р ) П ( Р )

есть преобразование Лапласа функции

 

 

/ з

( 0 = j f i ( i

— s)ft (s)ds,

 

 

 

 

 

 

О

 

 

 

 

 

 

а сходимость интеграла

 

 

 

 

 

 

 

 

 

 

n(p) = ] fs(t)e -p‘dt

 

 

 

в области Rep>yi абсолютная, — запишем (3,22,42)

в виде

 

где

ф; ср) = Qi к г -

Q i^; -

Q2«>;+

Q2^

;

(3.22.43)

t

 

 

 

 

<

 

 

 

 

s) ш2,(s) ds,

J£2 =

 

 

 

 

= j R (i _

j R (/ — s) tj»j.(s) ds.

(3.22.44)

Так как

c j k (f)j* =

£

c„/‘ (p),

то, обозначая

 

 

 

 

A=1

A=1

 

 

 

 

 

 

 

ф*(p) = ф* (p) — Qi (cei?)* H- Qi-%1 +

Q2W1 — QH$2 = 0,

(3,22,45)

и пользуясь теоремой [157] о том, что «если интеграл

 

 

 

 

ф* (р) =

j e-pt?i(0 dt

 

 

 

сходится в точке ро и во всех точках pk=Po + kl

(k—0, 1, 2,...),

функция

Тр, (рА) =0, то <pi(0 равна

нулю в интервале (0,

оо)

всюду, за

исключе­

нием точек разрыва функции ф!* (р)"», — получим из (3,22,45) уравнение

Фх (0 = <Pi(0 — Qi (*. у) wi (0 +

Qi (*. у)

(0 +

Qz (•*. у) щ (0 —

откуда

- - -0.2 (,Х >

у )

( 0 =

О ,

 

 

 

 

 

 

 

 

 

 

Фх (0 = Qx(*. lJ)w\(t) — Q2 (х , у) wl {t) Qi (л:, у)

(0 + Q2

(х, у)

(t).

 

 

 

 

 

 

 

(3,22,46)

Внося сюда выражения

для,#]

из

(3,22,44)

и учитывая (3,22,41),.

получим решение для <pi (t) в виде

 

 

 

 

 

 

Фх (0 = Qi (Ь -

- Y0 (t -

У )2 -

(Q^2 -'Л ) j Я (i -

s) ds

f

 

 

 

 

 

о

 

 

+ (2pX — p) J

 

 

f

 

 

— 2[ry I e~a\sQ(s — sk) R(t — s) ds.

(3,22,47)

о

 

 

В это выражение входят интегралы

от произведения ядра

R на неко-

 

t

 

торые функции. Поскольку функции

R и . | Rds протабулированы, по*

 

о

 

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

t

^ R (t — s) Ф (s) ds,

о

где в качестве* R (t)

выбрана сингулярная функция типа

 

R(t) = Ае-Щ«-'

(0 <

а < 1).

В самом деле, встречающиеся в

(3,22,47) интегралы, если учесть, чтск

Q2(s—sk) = 0(5—5ь), будут равны

 

 

 

t

 

t

R(t — s)ds— j* R (t — s) ds,

I*02(s — sk) R (t—*s)ds=

j*

6

 

0

 

0

t

t

 

 

*k

I*g-ais 0 (s — sk)R(t s)ds = f e~a'sR (/,— s)ds — j e~aiSR (t — s) ds.

6

o

 

o

t

 

 

 

 

Последние интегралы выражаются через/? и

J /?^следующим образом:

 

 

 

 

о

= j e-a'sR (t s)ds

 

<)

 

 

о

 

 

 

 

0