Файл: Филатов, А. Н. Асимптотические методы в теории дифференциальных и интегро-дифференциальных уравнений.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 21.10.2024

Просмотров: 124

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

причем для

упрощения R будем вновь обозначать через /?. Тог­

да (V.2.1)

примет

вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(*)

-

Е

М О

 

 

 

- х )

tx (x)rfx

]

(V.2.2)

Подставляя (V.2.2)

в уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д2и _ дох

 

 

 

 

 

 

 

 

 

(V.2.3)

 

 

 

 

 

 

 

— д х '

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

находим

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

д 2и

г,

д2и

 

 

 

 

 

 

 

 

 

 

 

 

 

t j p

( i - t )

^

Ф ^

- d t

 

 

(V.2.4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дх 2

 

 

 

 

 

здесь р — плотность стержня; и ( х ,

t)

— перемещение.

 

Зададим граничные и начальные условия:

 

 

 

 

 

и { х ,

t )

 

 

ди

 

= 0, U (Ху t)

 

 

 

ди

 

 

= A (x h

(V.2.5)

 

 

дх

 

= / (* ) .

dt

 

 

 

х=0

 

 

х=1

 

 

t=О

 

 

 

 

*=о

 

где

/ — длина стержня.

 

 

 

 

 

 

 

 

 

 

 

Будем искать решение уравнения (V.2.4), удовлетворяющее

граничным

условиям

задачи, в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

оо

т„(<) sin

 

 

 

 

 

 

 

 

 

 

 

 

и (.к,

t)

= 2

 

 

 

 

 

 

 

(V.2.6)

 

 

 

 

 

 

 

ft=l

 

 

 

 

 

 

 

 

 

 

Подставляя (V.2.6) в (V.2.4), получаем

для

определения

функ­

ции

Tk(t)

уравнение

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

К

+

а \ Тн =

s<4

j’ Ж *

-

t) Tt (т) dx,

k =

 

1, 2........

(V.2.7)

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

CO

( 2 * - l ) , / jE_

 

 

 

 

 

 

 

 

 

 

 

 

 

ft

2/

 

|/

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Опустив для

простоты записи

индексы

в

(V.2.7),

запишем урав­

нение (V.2.7)

в форме

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r ' +

co2r

 

 

 

 

7'(x)flfx.

 

(V.2.

 

 

 

 

 

= e(o2j / ? ( / — т)

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

Заменив

переменные по формулам

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т (/) =

cos ш/

 

с2 sin ct/,

 

 

 

 

 

 

 

 

 

 

 

Т (/) =

( —Ci sin<o/ +

с2cos о)/),

 

 

приведем

уравнение

(V.2.8) к стандартному виду

 

 

181


 

г

 

с \ — £u) sin

J R (t — T) [ (x) cos wx-f- c Q(x) sin (oxjdx

i

^

. (V.2.9)

t

 

C2 = £(1) COS 0)t 'JR (t — x) [ C1 (x) cos on -f-

C2 (x) sin (OxjVx

Усредняя систему (V.2.9) согласно

второй схеме усреднения,

„находим

 

 

(V.2.10)

где

оо

R s =

j R(s) sin usds >

0,

..

о

1

 

Rc =

oo

( ^ •

 

j

(s) cos cosflfs >

0.

Интегрируя систему (V.2.10), решение уравнения (V.2.4) при достаточно малом е представляем в виде

 

 

 

оо

,,, .

 

(2k—-\)Kx

oo

'"ktsk*

 

и (х,

t) =

 

 

 

 

У i k (t) sin '

 

2(

« ■ 2 ‘

-

 

 

 

 

 

k=\

 

 

 

 

ft=l

 

 

 

о

(

1

WkRck \

о

 

1

£0)k^ck\

 

( 2 k - l )кх

1

X cik cos

1

2

 

 

 

 

to

sin

2l

где

 

 

-

• ■

 

 

 

 

 

(V.2.11)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

°

2 Г V ,

 

ч .

(2Л - 1) те* ,

 

 

 

 

 

С\ь — ~ Г ) Л <*) sin

 

----- dx,

 

 

 

C2k ~

 

£MkRsk

 

 

 

 

 

 

 

 

 

2аь

с‘*

°ki

 

 

 

 

 

 

 

 

 

 

 

eco

 

 

 

 

 

 

 

 

а* ="«>* —

2

 

 

 

 

Из (V.2.11) следует, что наличие вязкости в (V.2.4) приводит к тому, что свободные колебания стержня. затухают по экспо­ ненциальному закону и наблюдается сдвиг частот.

2. Нелинейная задача. Свободные колебания вязко-упруго­ го стержня. Продольные колебания вязко-упругого стержня в нелинейной постановке были рассмотрены в работах [85, 145,

182


146]. Пусть связь между напряжением ах и деформацией ех за­ дается нелинейной зависимостью

 

ax ( t ) = £ ( ех - е $ # ( * —

t t t

\

о

 

 

G (t —

t — '2’

 

0 0 0

 

 

где e > 0 — малый параметр;

R (t), G(t, t,

t — время.

 

 

Подставляя

(V.2.12) в (V.2.3), находим

sx* 0 ( x) dx

2di 3 .

(V.2.12)

t) — ядра релаксации;

d2u _ p

d2a

d2u (x , t)

^~dW ~ ^

dx2

 

111

 

 

 

"1> t — x2, t

 

 

d2u (x , Tj)

du (x, т2) du (x, т3)

- е т Ш

 

0 ( *

 

 

dx2

 

дл:

~b

0 0 0

 

 

 

 

 

 

 

 

 

ди (x, tx)

 

d2u (x,

t2) du, (x , t3)

. d a

(x, tx)

du (x, z2)

d2a (x, t3)

d'zld’z2dx3.

dx

 

dx2

dx

dx

'

dx

dx 2

 

 

 

 

 

 

 

 

 

 

 

(V.2.13)

Примем

граничные и начальные

условия

для

уравнения (V.2.13)

в виде (V.2.5). Будем искать решение уравнения (V.2.13), удов­

летворяющее

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

в

виде

 

 

 

 

и (х , t ) =

Т .Х

 

(V.2.14)

 

 

 

 

Т (t) sin 2j".

Подставив (V.2.14)

в (V.2.13), получим для определения функции

T(t)

уравнение

 

 

 

 

 

 

Р

t

 

t tt

 

 

Г ' +

«)2Г = е

 

ш J # ( * — 'О П х) Л + а Щ

G(^— xlt £ — т2, t х3)Х

 

 

L

О

 

0 0 0

 

 

 

 

 

X

т(xt) Т (т2) Г ( х 3) dxxd i2dxz

(V.2.15)

где

 

 

 

 

 

 

 

 

 

 

• = | / т ( * >

“ - Н М

 

Введя новые

переменные по формулам

 

 

 

 

 

 

Т (t ) = Ct COS (at + C2 Sin (at,

 

 

 

 

T

(t ) = a) (— cxsin cot + c2cos (at),

183


приведем

уравнение (V.2.15)

к стандартному виду

 

 

 

е sin a>t

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С0‘ J R (t — х) [ct (х) cos (ОТ -f C2(x) sin cot] X

 

c i =

-

 

 

 

 

 

 

i

t

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*Clt

t — T2t

t ~ X z) [C1 (x1)cO S0)X1

- f

 

 

 

0

0 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

C 2 ( X

i )

Sin c o T j ] •

 

( t 2 ) COS c o t 2

- j - C 2

( x

2 )

sin a n 2 ]

X

 

 

 

X

[^i (x3) cos (ot3 -j- c2(x3) sin C1) X 3 ]

dxxdx2dxz

 

 

 

_

£ COS (lit

 

 

2

l

 

 

 

 

 

 

 

 

 

 

 

(V.2.16)

 

 

 

— x) [<4(x) COS COT

C2(t) sin cot]dx 4-

 

c2 =

a)

 

 

(0‘ J # l *

 

 

/11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

h J J J

G (t t15

t — t2,

* — T3)

(tj) cos cotj -f

 

 

 

0 0 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-j-c2 ( x , ) Sin C D X j ]

. [ c t ( x 2) C 0 S

0 ) X 2 -

f c2(x2) Sin C 0 T 2 ]

( x 3) c o s C 0 X 3

+

 

 

 

 

 

 

+

c 2 (хз) sin an3] dxxdx2dxz

 

 

 

 

 

Усредняя

систему (V.2.16) согласно второй схеме, находим

* ~

со

 

+ ш2# ^

+ а М ?

+

*12) +

a 2'-l(E2+

’l2)

 

, (V.2.17 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£

со

 

 

 

 

ш2/?^

+

a2S (Е2 +

?l2) -

a, ”l(E2 +

Ч2)

 

 

Г1= 2ш~

2/?сЕ -

 

 

где

 

 

 

 

 

h

 

 

 

 

 

 

 

h

 

 

 

 

 

 

0>\—

( 3а

4

"

 

ai)t" f " а2 =

( а 4 "ДЬ6 4 -

Л 8

) ,

 

 

“4“

 

“4~

 

 

 

 

о о

 

 

 

 

 

 

 

 

 

 

ОО

 

 

 

 

 

R s =

J

#

(s) sin cosds >

О,

Rc =

^ R (S) cos coscfs >

0 ,

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 0

0 0

0 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a z =

j

j

j

G (sb

, 2f „ ) sin со

(st 4- s2 — s3)

d stds2ds3,

 

 

 

 

0

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

CO

0 0

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a k =

и

0

I

G ^s"

S2'

cos ш ^S1 ^

s* ~~ S3^dsids2ds3t

 

 

 

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

184


оо оо со

а ь ~III G *S1’

S2’

 

sin 10

~^S3~

 

dslds*ds*'

 

0

0

0

 

 

 

 

 

 

 

 

 

 

 

oo oo oo

 

 

 

 

 

 

 

 

 

 

a &~

J I

j* G (Sl’

S2’

S3 ) C O S (0 ( S 2 +

s 3 sjdsidszds^

 

0

0

0

 

 

 

 

 

 

 

 

 

 

 

oo oo

oo

 

 

 

 

 

 

 

 

 

 

07 =

I11G ^Su

S2’

 

SinW(Sl

S3 _

Sz)

^ 1 ^ 3 ,

 

0

0

0

 

 

 

 

 

 

 

 

 

 

 

oo oo

oo

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos (o (Si -f- s3 s2)ds^ds2d s3.

 

0

0

0

 

 

 

 

 

 

 

 

 

 

Решение системы (V.2.17) имеет вид

 

 

 

 

 

H i)

= a *

1/

 

 

 

------тгл exp

\

 

z

)

 

 

 

у

R s ^ —a { a^e\p[—t(»Rstj

 

J

Xcos |

 

^ In

U)2^

-

a { a\ exp ( — ^ R s t)

j

+

<p0|;

 

 

 

 

 

 

 

 

 

 

 

 

 

(V.2.18>

rl (t) =

 

 

 

 

 

Rs0)2

 

 

 

 

 

t(jiRst \

 

 

 

 

 

 

 

exp |-------o— l X

= « o | / - R s w‘2 —

a\ exp

( — eu>Rs t )

 

(

-

 

{

Eu>Rct

 

a 2

 

“X

- a i a l exP ( - EU)^S 0

j + To } ■

— 2 ~ + W l ]n

где a 0, cp0 — произвольные

постоянные,

определяемые из началь­

ных условий. Решение уравнения (V.2.13) можно аппроксимиро­

вать решением

/??0>2

 

£0>Rs t

U (X, t)

 

exp -

а» | / ‘ R s <*2 -

a ^a\

e x p ( — w

X

 

2R s t )

 

 

£(dRc \

a2

0) 4 “

aia0exP { - ^

R s i )

X COS ^|to------ — 1t 2^ In

 

 

— (Pojsin-^-.

 

(V.2.19>

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

ла стержня.

3. Нелинейные поперечные колебания и динамическая устой­ чивость вязко-упругих стержней. Рассмотрим задачу о попереч­ ных колебаниях прямолинейного вязко-упругого стержня, загру­ женного периодической продольной силой Р (/)[145, 146). Стер-

18S