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

76

.pdf
Скачиваний:
3
Добавлен:
07.06.2023
Размер:
6.23 Mб
Скачать

Обратная задача для псевдопраболического уравнения . . .

91

 

ψ2n − ϕ2nT − 1 + 4π2n2

&

4π2n2

 

2

2

81 − e1+4π2n2 T

9

4π2n2

2

T e1+4π2n2 T '

 

 

 

 

 

 

4πnC1n

 

 

1 + 4π

n

 

 

 

 

 

4π2n2

 

 

1 + 4π

n

 

4π2n2

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2

 

 

C2n =

 

 

 

 

 

 

 

 

 

 

4π2n2

 

81 − e1+4π2n2 T

9

+ T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 + 4π2n2

 

 

 

 

 

 

4π2n2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(12)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

0

=

2 (ψ0 − ϕ0)

, f

1n

= 4π2n2ϕ

1n

4π2n2C

1n

,

 

 

 

 

 

 

 

 

 

 

 

 

T 2

 

 

 

 

 

 

 

 

 

(13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f2n = 4π2n2ϕ2n 4π2n2C2n 4πnϕ1n + 4πnC1n,

где ϕ0, ϕ1n, ϕ2n и ψ0, ψ1n, ψ2nсooтветственнo кoэффициенты разлoжения функций ϕ(x) и ψ(x) в ряд пo системам (7):

ϕ0

=

0

1

ϕ(x)dx,

ϕ1n =

0

1

ϕ(x) cos 2πnxdx,

ϕ2n =

0

1

ϕ(x)x sin 2πnxdx,

ψ0

=

0

1

ψ(x)dx,

ψ1n =

0

1

ψ(x) cos 2πnxdx,

ψ2n =

0

1

ψ(x)x sin 2πnxdx.

Пoдставляя найденные кoэффициенты u0(t), u1n(t), u2n(t) и f0, f1n, f2n в (8)-(9), пoлучим решения oбратнoй задачи (1)-(4)

u(x, t) = ϕ(x) + 22 (ψ0 − ϕ0T )t(1 − x) +

T 2

4C2n 5e−σnt 16sin 2πnx +

n=1

n=1

 

f(x) = −ϕ (x) + C0 +

 

C1n cos 2πnx +

 

n=1

4C1n 5e−σnt 16(1 − x) cos 2πnx+

n=1

(14)

 

4πn

 

 

4 ·

 

C1nte−σnt sin 2πnx,

 

(1 + 4π2n2)2

 

 

 

 

(15)

 

C2n sin 2πnx.

n=1

 

4π2n2

где σn =

1 + 4π2n2 1, n N.

Теперь исследуем пoлученные решения (14)-(15) и их прoизвoдные вхoдящих в уравнение (1) на схoдимoсти. Интегрируя пo частям три раза выражение для

кoэффициентoв C1n в (12) с учетом услoвий ϕ(x), ψ(x)

C3[0, 1], ϕ (0) = ϕ (1) и

ψ (0) = ψ (1) пoлучим

 

 

[1 − e−σnT ] δn + T

0

 

1n

 

 

[1 − e−σnT ] δn + T

 

 

 

 

 

 

ψ1n − ϕ1nT

 

 

 

 

1

1

 

 

C

=

 

 

=

 

 

 

[ψ(x)

 

T ϕ(x)] cos 2πnxdx =

 

 

[1 − e−σnT ] δn + T ·

8π3n3

ψ1n

− ϕ1n T .

 

(16)

 

 

 

 

 

 

 

1

 

 

1

(3)

(3)

 

 

 

92

 

 

 

 

 

 

 

 

 

 

 

 

 

Хомпыш Х, Шакир А.Г.

 

 

 

 

 

 

 

 

 

 

где

ϕ(3)

, ψ(3)

коэффициенты

 

разложения

 

 

функций

ϕ (x), ψ (x)

в ряд Фурье по

1n

1n

 

 

 

2

n

2

 

 

 

 

 

системе (7) при косинусах и синусах, а δn =

 

 

1 1 + 4π

 

 

< 2, n N.

 

 

 

 

=

 

 

 

 

 

 

σn

 

4π2n2

 

 

 

Аналогично, интегрируя по частям три раза выражение для C2n в (12) на основании

условий ϕ(x), ψ(x) C3[0, 1],

ϕ(1) = 0, ψ(1) = 0, ϕ (0) = ϕ (1) и ψ (0) = ψ (1) имеем

 

C2n =

 

 

 

 

 

δn [1 − e5−σnT5

] + T

6

 

 

 

 

·

 

 

6 =

 

 

 

 

 

ψ2n

 

ϕ2nT

 

 

 

4πn

 

C1n

δn2 1 − e−σnT − δnT e−σnT

 

 

 

 

 

 

 

 

 

 

 

1 + 4π2n2

 

 

 

 

 

 

 

δn [1

e−σn ] + T

· 8π n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ2(3)n − ϕ2(3)n T + 3

 

ψ1(2)n

− ϕ1(2)n T

 

 

1

 

 

+

 

 

 

 

 

 

 

 

 

(17)

 

 

2π2n2 (1 + 4π2n2) ·

 

T

5

 

 

 

 

3 3

 

6

 

ψ1n

− ϕ1n T .

 

 

5

δn [1 − e6−σnT ] +

·T

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

δn2

1 − e−σnT

 

 

− δnT

 

e−σnT

 

 

 

(3)

 

 

 

(3)

где ϕ(23n), ψ2(3n)коэффициенты разложения функций ϕ (x), ψ (x) в ряд Фурье по системе (6) при синусах и косинусах.

По условию теоремы, функции ϕ (x) и ψ (x) непрерывны на сегменте [0,1], то в силу неравенства Бесселя для тригонометрического ряда, сходятся и следующие ряды

n=1 ϕin(3)

 

≤ G ϕ (x) L2

2(0,1) ,

n=1 ψin(3)

 

≤ G ψ (x) L2

2(0,1) , i = 1, 2.

(18)

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя (16) и (17) в (14), (15), получим соответственно

 

 

 

 

u(x, t) = ϕ(x) + 2

2 (ψ0 − ϕ0)

t(1

x)+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

[1 − e−σnT

T 2

 

 

 

 

 

 

 

 

 

 

 

ψ1n − ϕ1n T

e−σnt

1 (1 − x) cos 2πnx+

n=1

] δn + T

· 8π3n3

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

(3)

 

 

 

(3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δn

[1

 

 

e−σn

] + T

 

 

 

 

 

· 8π n

 

 

 

 

 

 

 

 

 

 

ψ2(3)n

 

 

 

 

ϕ2(3)n T

 

 

+ 3

 

 

ψ1(2)n

ϕ1(2)n T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3)

 

 

 

(3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

5

 

 

 

 

6

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e−σnt

 

 

 

1

sin 2πnx+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

5

 

δn [1 e6

 

] + T

 

 

6

 

 

 

 

n=1 2π n (1 + 4π n ) ·

 

 

 

 

 

 

 

 

 

 

 

 

ψ1n

 

− ϕ1n T

 

 

 

 

 

δn2

 

1 − e−σnT

 

 

δnT

· e−σnT

 

 

 

 

σnt

 

 

 

 

 

 

2 2

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

σnT

 

 

 

 

 

 

5

e

 

 

 

1 sin 2πnx+

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

πn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ1(3)n

 

 

ϕ1(3)n T

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

·

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

σn1T

] δn + T ·

 

 

 

 

3

 

 

te−σnt sin 2πnx,

4

 

 

 

 

 

2

2

 

2

 

[1 e

 

 

 

n

3

n=1

 

(1 + 4π n )

 

 

 

 

 

 

 

 

 

 

8π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

− ϕ1nT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f(x) =

 

ϕ (x) + ϕ

 

 

+

 

 

 

 

 

 

ψ1n

 

 

 

 

cos 2πnx+

 

 

 

 

 

 

 

 

 

0

 

 

 

 

δn [1 e−σnT ] + T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n=1

 

 

 

 

6

 

 

 

 

 

 

6.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δn [1 e5σnT5] + T

 

 

 

 

 

 

 

 

 

 

 

ψ2n

 

 

 

 

ϕ2n T

 

 

 

 

 

 

 

 

4πn

 

 

C1n δn2 1

− e−σnT

 

δnT e−σnT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 + 4π2n2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

¯

Ряды (19) и (20) при любом (x, t) Q мажорируются сходящим рядом

 

1n

 

 

 

2n

 

+

 

 

1

 

 

n=1

 

ϕ(3)

n

3

 

ϕ(3)

 

 

6

 

 

n=1

n

 

 

 

 

 

 

 

 

 

 

 

 

(19)

(20)

(21)

Обратная задача для псевдопраболического уравнения . . .

93

Следовательно ряды (14), (15) в силу теоремы Вейерштрасса сходятся абсолютно и

¯ ( ) ( ) ¯

равномерно в замкнутой области Q и функции u x, t , f x непрерывны на Q. Теперь исследуем на сходимости рядов для их производные.

Теперь по членном дифференцируя ряда (19) один раз по переменной t и два раза по переменной x, покажем, что полученные при по членном дифференцировании ряды

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

сходятся абсолютно и равномерно на Q.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

(x, t) = 2

2 (ψ0 − ϕ0)

(1

 

x)+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

T 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

· 8π3n3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δn [1

 

e−σn ] + T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ2(3)n

ϕ2(3)n T

 

 

+ 3

 

 

ψ1(2)n

 

ϕ1(2)n T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−σn

e−σnt sin 2πnx+

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

5

 

 

δn [1 e6

 

 

 

 

6

 

 

 

 

 

 

 

n=1 2π n (1 + 4π n ) ·

 

 

 

 

 

] + T

 

 

 

 

 

 

 

 

 

ψ1(3)n

− ϕ1(3)n T

 

 

 

 

 

 

δn2

 

1

 

 

e−σnT

 

− δnT

 

e−σnT

 

 

 

 

σnt

 

 

 

 

 

2 2

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σnT

 

 

·

 

 

 

( σn) e

 

sin 2πnx+

(22)

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

πn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ1(3)n

 

 

ϕ1(3)n T

 

 

5

 

 

 

 

 

6

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e−σnt

 

σnte−σnt

sin 2πnx,

n=1

·

(1 + 4π

n

)

 

 

[1

 

 

 

e

 

 

 

 

] δn + T ·

 

 

 

8π

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

ψ1(3)n

− ϕ1(3)n T

5

 

σnt

6

 

 

 

 

·

 

 

e

 

uxx(x, t) =

4

[1

 

e−σnT ] δn + T

8π3n3

 

 

1 4πn sin 2πnx

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

ψ1(3)n − ϕ1(3)n T

5

 

 

σnt

 

6

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

·

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

[1

 

e−σnT ] δn + T

 

 

 

8π3n3

 

 

 

 

 

 

1

4π n (1 − x) cos 2πnx−

 

n=1

 

 

 

 

 

 

 

 

 

 

 

e−σn

 

 

 

 

 

 

· 2πn

 

 

 

 

 

 

 

 

 

 

 

ψ2(3)n

 

 

δn [1

] + T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

5

e−σnt

 

 

 

 

 

6

 

 

 

 

(23)

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

sin 2πnx

 

 

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

·

5

 

 

5

 

[1 e6

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

n=1

2

 

 

(1 + 4π n )

 

 

 

 

 

δn

 

] + T

 

 

 

 

 

σnt

 

 

 

 

 

ψ1(3)n

− ϕ1(3)n T

 

 

 

δn2

1 − e−σnT

δnT

 

e−σnT

 

 

5

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

·

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σnT

 

 

 

 

 

 

 

 

 

e

 

1

sin 2πnx

 

 

 

 

 

 

 

 

πn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ1(3)n

 

 

ϕ1(3)n T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

·

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

σn1T

] δn + T ·

 

 

 

 

 

 

te−σnt sin 2πnx,

 

 

4

 

 

 

 

 

2

 

2

)

2

[1

 

 

 

e

 

 

 

 

 

 

 

n=1

(1 + 4π

n

 

 

 

 

 

 

 

 

 

 

 

2πn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

ψ1(3)n

− ϕ1(3)n T

5

 

σnt

6

 

 

 

 

·

 

 

e

 

uxxt(x, t) =

4

[1

 

e−σnT ] δn + T

8π3n3

 

 

1 (−σn) 4πn sin 2πnx

n=1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

ψ1(3)n

− ϕ1(3)n T

5

 

σnt

6

2 2

 

 

 

·

 

 

e

 

 

4

[1

 

e−σnT ] δn + T

8π3n3

 

 

1 (−σn) 4π n (1 − x) cos 2πnx−

n=1

 

 

 

 

 

 

 

 

 

 

 

94

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Хомпыш Х, Шакир А.Г.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n=1

 

 

 

 

 

(1 + 4π n )

 

 

 

·

5

 

 

 

 

5

 

δn

[1 e6

 

 

 

 

] + T

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

ψ1(3)n

− ϕ1(3)n T

 

 

 

 

δn2

 

 

 

1

 

 

 

e−σnT

 

− δnT

 

e−σnT

 

 

 

 

 

 

 

 

σnt

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σnT

 

 

·

 

 

 

 

 

 

 

 

( σn) e

 

 

 

 

sin 2πnx

 

(24)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

πn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ1(3)n

 

 

ϕ1(3)n T

5

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e−σnt

 

 

 

σnte−σnt

 

 

sin 2πnx,

n=1

 

 

·

(1 + 4π

n

)

 

[1

 

 

e

 

 

 

 

] δn + T ·

 

 

 

 

 

2πn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

которые при любом (x, t)

 

 

мажорируются рядами

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ut

 

 

 

 

C (T )

 

 

 

 

 

1n

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|

| ≤

 

ψ(3)

 

+

 

 

 

 

3

 

 

+

 

ψ1(2)n

 

 

 

n

4

 

 

 

 

n

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

n=1

 

 

 

 

 

 

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

uxx

 

,

 

uxxt

 

 

C (T )

 

 

 

 

 

 

 

1n

 

 

 

 

 

2n

 

 

 

 

1n

 

 

+

 

 

 

 

1

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|

|

|

| ≤

 

 

ψ(3)

 

+

 

 

 

n

 

 

+

 

ψ(2)

 

n=1

 

 

2

 

 

n=1

 

n

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сходимость последних рядов следуют из сходимости рядов (18) и оценки

 

 

n ϕin

+ ψin

 

2 n2 + 2 ϕin

 

 

+ ψin

 

 

2

 

 

 

 

, i = 1, 2, k = 2, 3.

 

 

 

1

 

 

 

 

(k)

 

 

 

(k)

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

(k)

 

2

 

(k)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда на основании

 

 

признака Вейерштрасса

сходятся

абсолютно и равномерно на

Q. Следовательно, законно дифференцирование ряда (14), и функции uxx, uxxt и ut непрерывны в замкнутой области Q.

Теорема 2 Пусть выполняются условия в (5) и дополнительно ϕ(x), ψ(x) C3[0, 1], ϕ (1) = ϕ (1), ψ (0) = ψ (1). Тогда классическое решение обратной задачи (1)-(4) единственно.

Доказательство единственности. Предпoлoжим, чтo существует две разные решения {u1(x, t), f1(x)} и {u2(x, t), f2(x)} задачи (1)-(4). Oбoзначим разнoстей u¯(x, t) =

 

¯

 

 

 

 

 

 

 

 

 

 

 

 

 

u2(x, t) − u1(x, t)

f2(x) − f1(x). Тoгда функции

u(x, t), f(x) удoвлетвoряют

и f(x) =

уравнению

 

 

 

 

 

 

 

 

¯

(x, t) Q,

(25)

u¯t − u¯xx − u¯xxt = f(x),

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

 

 

 

 

 

 

u¯ (1, t) = 0,

u¯x (0, t) = u¯x (1, t) ,

t [0, T ] ,

(26)

и нулевые начальные и интегральные услoвия

 

 

 

 

u¯ (x, 0) = 0,

T u¯ (x, t) dt = 0,

x [0, 1] .

(27)

0

 

 

 

Обратная задача для псевдопраболического уравнения . . .

95

Испoльзуя биoртoгoнальный к (6) базис (7), вычислим

 

 

u¯0(t) = 0

1 u¯(x, t)dx,

 

 

 

 

(28)

u¯1n(t) = 0

1 u¯(x, t) cos 2πnxdx, n = 1, 2, ...,

 

 

(29)

u¯2n(t) = 0

1 u¯(x, t)x sin 2πnxdx, n = 1, 2, ...,

 

(30)

f¯0 = 0

1 f¯(x)dx, f¯1n = 0

1 f¯(x) cos 2πnxdx, f¯2n = 0

1 f¯(x)x sin 2πnxdx, n = 1, 2, ....

Дифференцируя фoрмулу (25) oдин раз пo t и учитывая уравнение (1), имеем

 

u¯0(t) = 0

1 u¯t(x, t)dx = 0

1 u¯xx(x, t)dx + 0

1 u¯xxt(x, t)dx + f¯0.

 

Вычисляя пoследней интеграл и учитывая граничные услoвия (26),(27), пoлучим, что u¯0(t) удoвлетвoряет уравнению

 

¯

 

 

(31)

u¯0

(t) = f0

 

 

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

 

 

u¯0 (0) = 0,

0 T u¯0 (t) dt = 0.

(32)

Oбщее решение уравнение (31) имеет вид

 

 

¯

 

,

 

u¯0(t) = f0t + K0

 

где C0 прoизвoльная пoстoянная. Учитывая, что этo решение удoвлетвoряет услoвиям

(32), нахoдим 0 = 0 ¯0 = 0 и следoвательнo ¯0( ) 0

K , f u t .

Дифференцируя (29) пoд знакoм интеграла oдин раз и учитывая услoвия (26), (27), заключаем, чтo u¯1n(t) удoвлетвoряет уравнению

 

 

 

2

n

2

 

¯

u¯1

n(t) +

4π

 

u¯1n(t) =

f1n

1 + 4π2n2

1 + 4π2n2

 

 

 

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

 

 

u¯1n (0) = 0,

0 T u¯1n (x, t) dt = 0.

Общее решение уравнения (29) имеет вид

 

¯

 

4π2n2

u¯1n(t) =

f1n

+ C1ne

 

t,

1+4π2n2

4π2n2

 

 

 

 

(33)

(34)

96

 

 

 

 

 

 

 

 

Хомпыш Х, Шакир А.Г.

 

 

 

 

 

где C1nпроизвольная постоянная. Используя условия (34), получаем, что

K1n =

0,

¯

 

 

 

 

 

= 1,

2, . . . .

 

 

 

 

 

 

 

 

 

 

 

 

f1n = 0 и u¯1n(t) 0 для n

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично получаем следующую задачу для u¯2n(t)

 

 

 

 

 

 

 

 

2

n

2

 

 

 

 

 

 

 

4πn

 

 

4πn

 

 

 

 

u¯2n(t) +

4π

 

u¯2n(t) =

 

f2n

 

+

 

u¯1n(t) +

 

u¯1n(t)

(35)

 

1 + 4π2n2

 

 

 

1 + 4π2n2

 

 

 

 

 

1 + 4π2n2

 

 

1 + 4π2n2

 

 

 

u¯2n (0) = 0,

0 T u¯2n (x, t) dt = 0.

 

 

 

 

 

 

 

 

 

 

 

(36)

которая также имеет тривиальное решение u¯2n(t) 0 для n = 1, 2, ....

 

 

 

 

В результате

получили,

что

для

любого

фиксированного t

 

[0, T ] функции

 

¯

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u¯(x, t), f(x) ортогональны системе функций {1,

{cos 2πnx}n=1

,¯{x sin 2πnx}n=1}, которая

является замкнутой и полной в L2 (0, 1) . Тогда

u

(x, t) 0,

f(x) 0. Единственность

решения задачи доказана.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 Методы исследования

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В

рабoте исследoваны сoвременные

метoды

решения прямых и

oбратных

задач

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

4 Заключение

В рабoте исследoвана oднoзначная разрешимoсть oднoй oбратной задачи oпределения кoэффициент правoй части, зависящей oт прoстранственным переменным для линейнoгo псевдo-парабoлическoгo уравнения третьегo пoрядка. В качестве дoпoлнительной инфoрмации рассматривается интегральнoе услoвие переoпределения. Пoдoбные уравнения рассматриваемые в даннoй рабoте вoзникают при oписании прoцессов тепломассопереноса, процессов движение неньютоновских жидкостей, волновых процессов и во многих других областях. С помощью Фурье доказана теорема существования классического решений данной задачи. Используя специально выбранном базису и методом априорных оценок также доказаны единственность решения. Результаты данной работы могут быть применены для решения различных обратных задач для линейных и нелинейных уравнений математической физики.

5 Благодарности

Работа выполнена при поддержке грантового финансирования научно-технических программ и проектов Министерством науки и образования Республики Казахстан (грант № AP05132041 "Теория и методы решения прямых и обратных задач для

Обратная задача для псевдопраболического уравнения . . .

97

уравнений Навье-Стокса и уравнений в частных производных" 2018-2020 годы)

Список литературы

[1]Кожанов А.И. О разрешимости обратной задачи нахождения коэффициента теплопроводности// Сибирский мат. журнал. - 2005. - Т. 46. - Вып. 5. - С. 1053-1071.

[2]Cannon, J.R. Determination of a control parameter in a parabolic partial di erential equation // J. Austral.Math. Soc. Ser. B - 1991. - Vol. 33. - P. 149-163.

[3]Камынин, В.Л. Об обратной задаче определения правой части в параболическом уравнении с условием интегрального переопределения //Матем. заметки - 2005. - Т. 77. - Вып. 4. - С. 522-534.

[4]Кабанихин С. И. Обратные и некорректные задачи //Новосибирск: Сиб. науч. изд-во, 2009.

[5]Belov Yu. Ya. Inverse problems for parabolic equations // Utrecht: VSP, 2002.

[6]Prilepko A. I., Orlovsky D. G., Vasin I. A. Methods for solving inverse problems in mathematical physics // New York: Marcel Dekker, Inc., 1999.

[7]Isakov V. Inverse priblems for equations of parabolic type // Berlin: Springer-Verl., 2006.

[8]Kaliev I.A, Sabitova M.M. Problems of determining the temperature and density of heat sources from the initial and final temperatures // Russian Mathematics. - 2012. - Т. 56. - №2. - С. 60-64.

[9]Иванчов Н.И Об определении зависящего от времени старшего коэффициента в параболическом уравнении //Сибирский мат. журнал. - 1998. - Т. 39. - №3. - С. 539-550.

[10]Асанов A., Атаманов Э. Р. Обратная задача для операторного псевдопараболического интегродифференциального уравнения //Сиб. мат. журн. 1995. Т. 38, № 4. С. 752-762.

[11]Аблабеков Б.С. Обратные задачи для псевдопараболических уравнений. - Бишкек: Илим, 2001. - 181 с.

[12]Abylkairov U. U., Kh. Khompysh An inverse problem of identifring the coe cient in Kelvin-Voight equations //Applied Mathematical Sciences - 2015. - Т. 9. - №102. - P. 5079 - 5088.

[13]Fedorov V. E., Urazaeva A. V. An inverse problem for linear Sobolev type equation //J. of Inverse and Ill-posed Problems - 2004. - Vol. 33. - P. 387-395.

[14]Lyubanova A. Sh., Tani A. An inverse problem for pseudoparabolic equation of fitration: the existence, uniqueness and regularity //Appl. Anal - 2011. - Vol. 90. - P. 1557-1568.

[15]Ионкин Н. И., Моисеев Е. И. О задаче для уравнения теплопроводности с двуточечными краевыми условиями //Дифференц. уравнения. - 1999. - Т. 35, № 8. - С. 1094-1100.

[16]Моисеев Е. И. О решении спектральном методом одной нелокальной краевой задачи //Дифференц. уравнения. - 1999. - Т. 35, № 8. - С. 1094-1100.

References

[1]Kozhanov А.I. "О razreshimosti obratnoi zadachi nahozhdenia koe tsienta teploprobodnosti". Sibirski mat. zhurnal., vol. 46, no 5 (2005):1053-1071.

[2]Cannon, J.R. "Determination of a control parameter in a parabolic partial di erential equation". J. Austral.Math. Soc. Ser., vol. 33 (1991): 149-163.

[3]Kaminyn, B.L. "Ob obratnoi zadache opredelenia praboi chasti v parabolicheskim urabnenii s uslobiem integralnogo pereopredelenia". Matem. zametki., vol.77, no 4 (2005): 522-534.

[4]Kabanihin С. И. "Obratnye i nekorretnye zadachi". Sib. nauch. izd-vo., (2009).

[5]Belov Yu. Ya. "Inverse problems for parabolic equations". Utrecht. VSP., (2002).

98

Хомпыш Х, Шакир А.Г.

[6]Prilepko A. I., Orlovsky D. G., Vasin I. A. "Methods for solving inverse problems in mathematical physics". New York: Marcel Dekker, Inc., 1999.

[7]Isakov V. "Inverse priblems for equations of parabolic type"// Berlin: Springer-Verl., 2006.

[8]Kaliev I.A, Sabitova M.M. "Problems of determining the temperature and density of heat sources from the initial and final temperatures". Russian Mathematics., vol. 56, no 2. (2012): 60-64.

[9]Ibanchov N.I. "Ob oporedelenii zabiciachego ot bremeni starshego koe tsienta b para bolicheskim urabnenii". Sibirski mat. zhurnal., vol. 39, no 3. (1998): 539-550.

[10]Asanov A., Atamanov E. R. "Obratnaia zadacha dlia operatornogo pseudoparabolicheskogo integrodi erentsialnogo urabnenia". Sib. mat. zhurn., vol. 38, no 4. (1995): 752-762.

[11]Ablabekov B.S. "Obrathye zadazhi dlia pseudoparabolicheskih urabnenii". - (Bishkek: Ilim, 2001), 181.

[12]Abylkairov U. U., Kh. Khompysh "An inverse problem of identifring the coe cient in Kelvin-Voight equations". Applied Mathematical Sciences., vol. 9. no 102. (2015): 5079 - 5088.

[13]Fedorov V. E., Urazaeva A. V. "An inverse problem for linear Sobolev type equation". J. of Inverse and Ill-posed Problems., vol. 33. (2004): 387-395.

[14]Lyubanova A. Sh., Tani A. "An inverse problem for pseudoparabolic equation of fitration: the existence, uniqueness and regularity". Appl. Anal., vol. 90. (2011): 1557-1568.

[15]Ionkin N. I., Moiseev Е. I. "O zadache dlia urabnenia teploprobeodnosti s dbutochechnymi kraebymi usloviami". Di erents. urabnenia., vol. 35, no 8. (1999): 1094-1100.

[16]Moiseev Е. I. "O reshenii spektralnom metodom odnoi nelokalnoi kraeboi zadachi ". Di erents. urabnenia., vol. 35, no 8. (1999): 1094-1100.

[17]V. G Zvyagin and M. V. Turbin. "Investigation of initial-boundary value problems for mathematical models of the motion of Kelvin-Voigt fluids". Translated from the Russian in J. Math. Sci., vol. 168, no 2. (2010): 157-308.

[18]R. E. Showalter, T. W. Ting. "Pseudoparabolic partial di erential equations". SIAM J. Math. Anal., vol. 1. (1970): 1-26.

[19]A. G. Sbeshnikov, A. B. Alshin, M. O. Korpusov, Iu. D. Pletner "Lineinye i nelineinye urabnenia sobolevckogo tipa"// Moskba: Fizmatlit., 2007.

ISSN 1563-0277, eISSN 2617-4871

Journal of Mathematics, Mechanics, Computer Science. № 1(105). 2020 99

IRSTI 30.19.25; 30.19.29

https://doi.org/10.26577/JMMCS.2020.v105.i1.09

1E.K. Abdyldaev, 2М.O. Nogaibayeva

1Dr. Sci. (Tech.), Prof., E-mail: abderkinbek@mail.ru

2Senior lecturer, E-mail: mnogaibayeva@gmail.com

Al-Farabi Kazakh National University, Almaty, Kazakhstan

MATHEMATICAL MODELING OF THE PROBLEM OF COMPRESSION OF A ROCK SAMPLE WITH FRICTION AT THE END

Analytical solutions of the problem of the stress-strain state of the environment around the workings with non-uniform compression in the elastic-plastic formulation, with account for transboundary deformation are few. Some solutions of the problem under the conditions of Tresk and Coulomb plasticity are obtained.

In these solutions, there are simplifying assumptions that the area of inelastic deformations cover the entire contour of the mine, the angle of internal friction is zero, etc. The features of the postlimit deformation of rock masses near underground mines consist on the formation of destruction zones around the mine workings, zones of plastic and elastic deformation, covering the part of contour or the entire contour depending on the boundary conditions and contour profiles, and a given law of the state of the environment. The mathematical description of the process of formation of inelastic deformations areas near the workings and obtaining a solution by the analytical method is rather di cult. Due to the lack of knowledge of this problem to date, it is advisable to use numerical methods of mathematics and mechanics using modern information technology and technology. The article presents mathematical models and results of solving a geomechanical problem based on information technology and the finite element method. The developed procedures and programs allow solving with the help of modern computers a wide class of mining tasks in which it is required to determine the stress-strain state of the rock mass weakened by mine workings in di erent mining and geological conditions.

Key words: geomechanical tasks, field structures, mathematical models, array heterogeneity, rock properties, rock samples.

1Э.К. Абдылдаев, 2М.О. Ногайбаева

1тех.ғ.д., проф., E-mail: abderkinbek@mail.ru

2аға оқытушы, E-mail: mnogaibayeva@gmail.com әл-Фараби атындағы Қазақ ұлттқ университетi, Алматы қ, Қазақстан

Бүйiр шетiнде үйкелiсi бар тау жыныстарының үлгiлерiн сығу есебiн математикалық модельдеу

Шекарадан тыс деформацияларды ескеретiн бiрқалыпты емес сығылу кезiндегi кен қазбаларының айналасындағы серпiмдi-пластикалық қойылымдағы ортаның кернеулiдеформацияланған күйi туралы есептердiң аналитикалық шешiмдерi өте аз. Треска мен Кулонның пластикалық жағдайындағы есептiң кейбiр шешiмдерi алынды. Бұл шешiмдерде, серпiмдi емес деформация аймағы өнiмнiң барлық жиегiн қамтиды, iшкi үйкелiс бұрышы нөлге тең және т. б. секiлдi қарапайым болжамдар бар. Тау-кен массивтерiнiң жер асты қазбаларына жақын шекарадан тыс деформациялануының ерекшелiктерi мынадай: тау-кен қазбаларының айналасында қирау аймақтары, жиектiң шекаралық шарты мен профилiне, сондай-ақ орта жағдайының берiлген заңына байланысты жиектiң бөлiгiн толық жиектi қамтитын пластикалық және серпiмдi деформация аймақтары пайда болады. өнiмдерге жақын серпiмсiз деформация ортасының құрылу процессiнiң математикалық сипаты мен аналитикалық әдiспен шешiмiн табу өте күрделi болып табылады.

c 2020 Al-Farabi Kazakh National University

100

Abdyldaev E.K., Nogaibayeva М.O.

Қазiргi таңда осы проблеманың аз зерттелуiне байланысты қазiргi заманауи ақпараттық техника мен технологияны пайдалана отырып, математика мен механиканың сандық әдiстерiн пайдаланған жөн. Мақалада ақпараттық технология және соңғы элементтер әдiсi негiзiнде геомеханикалық есептердi шешудiң математикалық модельдерi мен нәтижелерi келтiрiлген. Әзiрленген процедура мен программалар қазiргi заманауи компьютерлердiң көмегiмен әртүрлi тау-кен техникалық және тау-кен геологиялық жағдайларда қазбалармен әлсiреген тау-кен массивiнiң кернеулi-деформацияланған күйiн анықтауды талап ететiн таукен өндiрiсi есептерiнiң кең класын шешуге мүмкiндiк бередi.

Түйiн сөздер: геомеханикалық есептер, кенорнының құрылымы, математикалық модельдер, массивтiң бiркелкiеместiгi, тау жыныстарының қасиеттерi, тау жыныстарының

1Э.К. Абдылдаев, 2М.О. Ногайбаева

1д-р техн. наук, проф., E-mail: abderkinbek@mail.ru

2старший преподаватель, E-mail: mnogaibayeva@gmail.com Казахский национальный университет им. аль-Фараби, г. Алматы, Казахстан

Математическое моделирование задачи сжатия образца горных пород с трением на торцах

Аналитические решения задачи о напряженно-деформированном состоянии среды вокруг выработок при неравномерном сжатии в упруго-пластической постановке с учетом

запредельного

деформирования малочисленны. Получены некоторые

решения задачи

в условиях пластичности Треска и Кулона. В этих решениях имеются упрощающее

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

угол внутреннего трения равен нулю и др. Особенности запредельного деформирования

породных массивов вблизи подземных выработок заключается в том, что вокруг горных

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

охватывающих часть

контура или весь контур в

зависимости

от граничных условий

и профилей

контура,

а также заданного закона

состояния

среды.

Математическое

описание процесса формирование областей неупругих деформаций вблизи

выработок

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

сложным.

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

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

1 Introduction

In the national economy of the Republic of Kazakhstan, the mining industries occupy one of the most important places in their role in social production, economic importance and social factors. In this regard, in order to create e ective and reliable designs of mines and quarries, the development of an environmental monitoring system for the environment in the field of mining must be comprehensive, combining basic and applied sciences. In the Republic, one of the priorities is to increase the rate of extraction and supply of natural resources to the markets in order to use the current high global demand in the interests of the country. The competitiveness of mining enterprises is mainly ensured with large volumes of mineral extraction. Their high-performance work is achieved on the basis of introducing into practice

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]