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

Методическое пособие 544

.pdf
Скачиваний:
6
Добавлен:
30.04.2022
Размер:
2.2 Mб
Скачать

мени. В одномерном явлении соответвенные граничные условия ставятся на граничных поверхностях слоя x 0 и x l.

Возможно несколько вариантов постановки граничных услови применительно к адачам теплопроводности.

1. Граничное условие первого типа при условии, что в любой точке поверхности тела известна температура

 

u S P, t , P S, t 0,

 

 

108

где P, t

– заданная функция точки

P

поверхности

S и

времени t.

 

 

 

 

2. Граничное условие второго типа при условии, что на поверхности S тела известен тепловой поток qn q n, здесь

q – вектор плотности теплового потока, а

n

– единичная

внешняя нормаль к поверхности

S.

Согдасно закону Фурье

q

k gradu n k

u

.

 

n

 

n

 

 

 

 

Таким образом, граничное условие второго типа интерпери-

руется на поверхности

S

как нормальная производная темпе-

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

 

 

u

P, t

, P S, t 0,

109

n

 

 

 

 

 

 

 

 

 

 

здесь P, t

q

– заданная зависимость.

 

 

 

n

 

 

 

 

 

 

 

 

 

 

k

 

 

 

P, t 0

 

Если поверхность теплоизолирована, то

и

граничное условие

u

0

задано всюду на поверхности S.

 

n

 

 

 

 

 

 

 

 

3. Граничное условие третьего типа моделирует тепловой режим на поверхности тела, отвечающий конвективному теплообмену согласно закону Ньютона с окружающей внеш-

ней средой с температурой

*

.

Из закона Ньютона следует,

u

что плотность теплового потока на границе тела прямо про-

71

порциональна разности температур тела и окружающей среды, что означает

q

q n

 

u u

*

.

 

 

 

 

 

n

 

T

 

 

 

Коэффициент теплообмена (теплоотдачи)

 

T

 

изменяет-

ся в зависимости от характеристик среды и, вообще говоря,

от приращения температуры

ент

T

является константой,

*

.

На практике коэффици-

u u

свободной от температуры и

имеет одно и тоже значение в любой точке поверхности тела.

 

Следовательно, граничное условие третьего типа выра-

жает зависимость между температурой u

и еѐ нормальной

производной

u

в каждой точке поверхности тела

n

 

 

 

 

 

 

 

 

 

 

 

 

 

u

hu hu

*

P, t , P S, t 0,

110

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

здесь

h

 

T

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

Замечание. Граничные условия первого, второго и третьего типа представляется возможным произвольно записать в форме общего граничного условия

u u P, t , P S, t 0,n

здесь и – произвольные постоянные;

P, t

111

– извест-

ная на поверхности тела зависимость. Задавая в соотношении

111 0,

1, а

P, t P, t ,

будем иметь гранич-

ное условие

первого

типа 108 . В

случае 1,

0,

P, t P, t ,

соотношеиие

111

трансформируется в

граничное условие второго типа

109 .

Задавая 1,

h,

а P, t hu* P, t , получим условие третьего типа 110 .

72

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

где

 

0

 

k

u

u

, P S, t 0,

 

 

 

 

 

4

 

112

 

n

0

0

 

 

 

 

 

 

 

 

 

 

 

 

– мера черноты вещества, вообще говоря,зависящая от

температуры; 0 – константа Стефана-Больцмана.

Правая часть соотношения

112 степенным образом

(нелинейно), зависит от температуры.

5. Математические модели температурных полей в многослойных телах и оболочках на поверхности соприкосновения двух тел имеют граничные условия сопряжения, иначе говоря, граничные условия четвѐртого типа.

В случае идеального теплового соприкосновения данные условия

u

P, t u

, P , t 0;

 

 

1

 

 

 

2

 

 

 

113

 

 

 

 

 

 

 

 

 

 

u

k

 

u

 

, P , t 0,

 

 

2

 

 

k

1

 

 

 

 

 

2

 

 

 

 

1

n

 

n

 

 

 

 

 

 

 

 

отражают равенство температур и тепловых потоков на поверхности соприкосновения .

В случае неидеального теплового соприкосновения с

термическим сопротивлением

R

на поверхности соприкосно-

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

 

1

u

 

u

k u1 k

 

u2

, P , t 0,

114

 

 

2

2

 

 

R

1

1

n

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

здесь n

– внешняя нормаль относительно первого тела к по-

верхности соприкосновения

.

 

 

Метод Фурье для уравнения теплопроводности. Однородное уравнение.

73

Предположим, что на сегменте

0, l

оси

Ox

находится

тонкий однородный стержень с теплоизолированной от внешней среды боковой поверхностью. Пусть на концах данного стержня фиксируется одинаковая температура: на левом температура A, на правом – температура B. Начальное рас-

пределение температуры в стержне выражается зависимостью x . Определим температуру в точке x в момент времени t 0, считая, что внутри стержня нет источников и по-

глотителей тепла.

Данная задача заключается в решении дифференциального уравнения

u

 

 

 

 

u

 

 

 

 

 

 

a

 

 

2

 

 

 

,

 

 

 

 

 

2

 

 

 

 

 

 

 

115

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

x

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с начальным условием

 

 

 

 

 

 

 

 

 

116

u

t 0

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и граничными данными

 

 

 

 

 

 

 

 

 

117

u

 

x 0

 

A, u

 

x l

B.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Произведѐм подстановку с условием выполнения для искомой функции однородных граничных данных. Конкрет-

но, зададим зависимость v x, t , выражающуюся через неиз-

вестную зависимость u x, t формулой

 

B A

 

: v x, t u x, t

 

x A .

 

l

 

Зависимость v x, t является решением уравнения

v

 

 

 

v

 

 

 

 

 

 

a

 

2

 

,

 

 

 

 

 

2

 

 

 

 

 

118

 

 

 

 

 

 

 

 

t

 

 

x

2

 

 

 

 

 

 

 

 

 

 

 

 

с начальным условием

 

 

 

 

 

 

 

 

 

 

 

 

 

B A

 

 

v t 0 u

 

t 0

 

 

 

x A

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

74

 

 

 

x

B A

x A

 

 

x

 

 

 

 

l

 

1

 

 

 

 

 

 

 

120 (119)

и однородными граничными данными:

v

x 0

0, v

x l

0.

 

 

 

Для решения однородного уравнения

118

120

с неодно-

родным начальным ными данными 120

условием 119 и однородными гранич-возможно использование метода разде-

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

Конкретно, сперва будем искать нетривиальные решения уравнения 118 , удовлетворяющие граничным данным

120 ,

в форме произведения

 

 

 

 

 

121

 

v

 

x, t

 

X

 

x

 

T

t

0.

 

 

 

 

 

 

 

 

После подстановки 121

в уравнение 118 , будем иметь

 

 

X x T t a

2

X x T t ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

1

T t

 

X x

const.

a

2

T t

X x

 

 

 

 

 

 

 

122

Следовательно, T t a2T t 0;

X x X x 0. Из граничных данных 120 вытекает

X 0 T t 0, X l T t 0,

следовательно, X 0 0, X l 0.

75

123

124

125

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

X x X x 0,

 

 

126

 

X 0 0, X l 0,

 

 

 

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

 

 

 

n

2

 

 

, n N,

 

 

 

n

 

l

 

 

 

 

 

 

127

существуют нетривиальные линейно независимые решения

задачи 124

и 125 :

 

 

 

 

 

 

 

X

 

x sin

nx

.

 

n

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При n

уравнение 123

 

 

 

 

 

 

n

2

 

 

 

 

 

 

 

 

2

t

 

T

t

A e

 

 

 

a

 

 

l

 

 

,

 

n

 

 

n

 

 

 

 

 

128

имеет решение в форме

129

здесь

An – коэффициенты, которые нужно найти. Следова-

тельно, в силу 121 , любая зависимость

 

 

 

n 2

2

 

nx

 

 

 

 

a

t

130

vn x, t X n x Tn t Ane

l

 

sin

l

 

 

 

 

 

 

является решением уравнения

118 , для которого выполня-

ются граничные данные

120 .

Решением того же уравнения

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

 

 

 

 

 

n

2

 

 

 

 

 

2

t

 

n

 

 

 

 

a

v x, t

x, t

n

l

 

 

v

 

A e

 

 

 

 

 

 

 

 

 

 

n 1

 

n 1

 

 

 

 

sin

nx

.

l

 

 

131

Найдѐм коэффициент An при условии удовлетворения начального условия 119 . Считая t 0, будем иметь

76

 

 

 

nx

 

1

x

n

 

 

 

A sin

l

.

 

 

n 1

 

 

 

 

 

132

Представленный

 

ряд

есть

 

 

разложение

 

этой

зависимости

1 x

в ряд Фурье по синусам кратных аргументов. Коэффи-

циенты

An вычисляются по формуле Эйлера-Фурье элемен-

тарным интегрированием:

 

 

 

 

 

 

B A

 

 

 

nx

 

 

2

l

x sin

nx

 

 

 

2

l

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

dx

 

 

 

 

 

A

 

 

x

 

sin

 

dx

n

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

0

 

 

l

 

 

 

 

l

0

 

 

 

 

 

 

l

 

 

 

l

 

 

 

2

l

 

nx

 

 

 

 

2

 

A 1

 

B .

 

 

 

 

 

 

 

x sin

dx

 

 

n 1

 

 

 

133

 

l

l

 

n

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для ряда 131

 

с коэффициентами

133

выполняются все

условия задачи 118 – 120 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Элементы

vn x, t

ряда

131

являютя моделью темпе-

ратурных волн. Амплитуда температурных волн сокращаается с течением времени, так как она содержит множитель

 

 

n

2

 

 

2

t

 

 

 

a

e

 

l

 

.

 

 

 

 

 

 

Итоговое решение заданной краевой задачи

117 :

 

115

u x, t A

B A

l

 

здесь коэффициенты

x An

 

 

n

2

 

nx

 

2

 

 

 

 

 

a

t

 

n

l

 

sin

 

,

 

A e

 

l

n 1

 

 

 

 

 

 

 

 

 

 

 

определяется соотношением

134133 .

Пример 1. Решить однородное уравнение теплопроводности на отрезке:

u u

, x 0, 1 , t 0, ,

t

xx

 

 

u

t 0

sin3

2 x, x 0, 1 ,

 

 

 

 

77

u

x 0

u

x 1

0, t 0, .

 

 

 

Решение.

1. Находим вспомогательные решения исходного урав-

нения в форме v x, t X x T t

при условии v x 0 v x 1 0,

что означает

X 0 X 1 0.

 

С этой целью подставим зави-

симость v x, t X x T t

в заданное уравнение и разделим

аргументы:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

,

X

 

x

T

t

 

 

 

 

 

X

x

T

t

 

 

 

,

 

 

 

 

 

v

xx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

x

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T t X x ,

 

 

 

 

 

X x T t T t X x x

 

 

 

 

 

 

 

 

 

 

 

 

 

T t

 

X

x

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T t

X x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

зависимости

X x

и

 

T t

представля-

ют собой решения пары задач:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а)

X

 

x X x

0, X

0 X 1 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б)

 

t T t 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2. Находим решение задачи а).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

X x X x 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

0;

k

 

 

i.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1, 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Общее решение:

 

X x C1 cos

 

 

x C2 sin

x. Исполь-

зуя граничные данные, будем иметь:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X 0 C 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X 1 C sin

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

n

 

, n N;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

n

x C sin nx C sin nx, C

 

C

, n N.

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

n

 

 

 

2

 

 

 

3. Найдѐм решение задачи б). Для

 

2n2

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

78

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T t

 

n T t 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dT

 

2

2

 

 

 

dT

 

 

 

2

 

2

 

 

dT

 

 

 

 

 

2

 

2

 

 

 

dt

n T ;

 

T

 

 

n dt;

T

 

n

 

 

dt;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln T

 

2

n

 

 

 

 

~

;

 

 

ln T

2

n

2

 

 

 

 

 

 

 

~

 

;

t ln An

 

 

t ln An

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

t

 

~

 

 

 

 

 

 

 

 

~

 

 

 

 

2

 

2

t

 

 

ln T

ln e

 

n

ln An

;

 

ln T

 

ln An e

 

 

n

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

2

 

2

t

 

 

 

 

 

 

 

 

 

Общее решение данного уравнения:

Tn

An e

 

n

, n N.

 

 

 

 

 

 

 

4. Таким образом, вспомогательные решения заданного уравнения представляются в форме

 

 

 

 

~

2

2

 

 

2

2

 

 

 

 

 

 

x, t C An e

t

 

t

 

 

 

 

v

 

n

sin nx A e

 

n

sin nx,

 

 

 

n

n

 

 

 

n

 

 

 

 

 

 

Cn

~

 

 

 

 

 

 

 

 

 

где

An

An .

 

 

 

 

 

 

 

 

5. Решение исходного уравнения теплопроводности ищем в виде

 

 

 

 

 

 

 

u x, t

n

x, t

n

2

2

t

v

A e

 

n

 

 

 

 

n 1

 

n 1

 

 

 

sin nx.

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

при которых ряд для функции u x, t сходится.

6. Находим коэффициенты

A

,

n

 

такие, что u x, t удо-

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

 

u

 

sin

3

2 x sin

2

2 x sin 2 x

1 cos 4 x

sin

2 x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t 0

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin 2 x

1

sin 6 x sin 2 x

 

 

sin 2 x sin 2 x cos 4 x

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

79

 

 

 

 

 

 

С другой стороны,

3

sin 2 x

1

sin 6 x.

 

 

 

 

4

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

1

 

u

 

 

A sin nx

sin 2 x

sin 6 x.

t 0

 

 

 

 

 

n

 

4

 

4

 

 

 

 

 

n 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда

A

3

,

 

2

4

 

 

 

A

1

,

 

6

4

 

 

 

A A

A

A 0,

1

3

4

5

An 0, n 7.

Подставляя эти коэффициенты в формулу для решения исходного уравнения теплопроводности из пункта 5, получаем:

u x, t

3

 

2

t

 

1

 

2

t

 

e

4

sin 2 x

 

e

36

4

 

 

4

 

 

 

 

 

 

 

 

 

 

sin 6 x.

Пример 2. Решить однородное уравнение теплопроводности на отрезке:

 

u

1

u , x 0, 3 , t 0, ,

 

 

 

 

t

9

 

 

xx

 

 

u

 

 

 

 

3 x 4sin12 x, x 0, 3 ,

 

8sin

3

 

t 0

 

 

 

 

 

 

 

 

 

u

x 0

u

x 3

0, t 0, .

 

 

 

 

 

Решение.

1. Находим вспомогательные решения исходного урав-

нения в форме

v x, t X x T t

при условии v

x 0

v

 

x 3 0,

что означает

X 0 X 3 0.

 

С этой целью подставим зави-

симость v x, t X x T t

в заданное уравнение и разделим

аргументы:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

xx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

9

v

 

,

 

X x T t

 

 

t

 

9

 

 

X x T t

 

x x

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X x T t

 

1

T t

 

 

 

 

 

 

 

 

 

 

1

T t X x ,

 

9

X x

x

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

t

 

X

x

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

T t

X x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

80