Файл: Ониани, Ш. И. Тепловой режим глубоких шахт при гидравлической закладке выработанного пространства и сложном рельефе поверхности.pdf

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

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

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

Добавлен: 17.10.2024

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

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

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

. шичеших

задач,

касающихся

в основном

тел простейших

'форм. До

40-х годов при решении

задач

нестационарной

теплопроводности

преимущественно

пользовались

методом

Фурье. В

горной

теплофизике

классический

метод

успешно

• применил А.

Ф. Воропаев для

описания

процесса

тепло- и

маооообмена

между

горным

массивом

и вентиляционной

струей, протекающей

по выработкам вентиляционной сети

[15].

 

 

 

 

 

Метод операционного .исчисления дли решения дифферен­

циальных уравнений в частных производных второго

порядка

в математической физике широкое распространение

получил

в последние годы- А. В. Лыков в своих работах изложил ос­ новные принципы практического приложения метода преоб­

разования

Лапласа" и успешно применил

его для

решения из­

вестных

и

новых

задач

 

нестационарной

тепло- и масеопро-

водности

[91, 92,

93]. В

теории теплопроводности

этим мето­

дом

широко пользуются

Г. Карелоу и Д. Егер [61]. В обла­

сти

гарной

теплофизики

примером

успешного

применения

метода операционного

исчисления

служат известные работы

А. Н. Щербаня и

О. А.

Крайнева [147,

148, 152].

 

Несмотря на свою распространенность, два точных ана­ литических метода охватывают далеко не все технически интересные задачи нестационарной теплопроводности. Кроме того, точные аналитические решения иногда настолько .слож­ ны и громоздки, что теряют практическое значение. Вследст­ вие этого в последнее время наблюдается стремление к раз­ работке и применению приближенных методов решения крае­ вых задач теплопроводности. В большинстве случаев такие решения при достаточной точности приближения просты и удобны для ведения инженерных расчетов. Среди приближен­ ных методов наибольшее (распространение получил метод ко­ нечных разностей, или так называемый метод сеток, сущность которого заключается в разбивке пространственной перемен­ ной на отдельные отрезки, а временной переменной на от­ дельные периоды. Искомая температура определяется для каждого отрезка и периода времени последовательно. Даль­ нейшим развитием метода конечных разностей является метод элементарных балансов, примененный для решения прост-.

128


ранственных нестационарных -задач теплопроводности яри гра­ ничных условиях первого, второго и третьего рода [9]. Для решения многих технических задач нестационарной теплопро­ водности с граничными условиями всех родов используется и •метод исключения дифференциального уравнения теплового баланса одной из нескольких независимых пространственных переменных. При этом действительная температурная кривая в теле заменяется параболой n-vo порядна [10]. В области гор­ ной теплофизики приближенный метод (метод конечных раз­ ностей) для решения нестационарных задач то прогнозу и регулированию теплового режима глубоких шахт был при­ менен А. Ф. Воропаевым [14].

Классический .метод точного аналитического решения диф­ ференциального уравнения теплопроводности и метод конеч­ ных разностей не дают возможности получить удовлетвори­ тельное решение (поставленной перед нами задачи, а прибли­ женный метод А. И. Вейника [10] приводит к более сложным зависимостям, нежели метод точного аналитического реше­ ния с применением интегрального преобразования Лапласа. Поэтому для точного аналитического решения задачи неста­

ционарного теплообмена между

закладкой и прилегающими

к ней массивами угля и породы

мы пользуемся методам 'Опе­

рационного исчисления, который согласно современной тео­ рии теплопроводности [61, 91], является наиболее удобным методом для решения задач нестационарной теплопроводно­

сти

в

составных

телах-

 

 

 

 

Применим преобразование

Лапласа

к уравнению

(5.1):

 

 

L

 

 

 

L

д \ (х,

х)

(5.13)

 

 

 

дх

 

дх2

 

 

 

 

 

 

 

 

 

 

 

L [t (х,

т)] =

j

t(x, x)e~Pxdx

= t(x,

p),

(5.14)

 

 

 

 

 

0

 

 

 

 

 

 

L

Ma (x,

x)

=

pt3(x,p)-f{0),

(5.15)

 

 

 

дх

 

 

 

 

 

 

 

 

 

 

где

/(0) = fu

 

согласно

уравнению (5.3).

 

 

 

•С учетом (5.15) выражение (5.13) перепишется следую­

щим

образам:

 

 

 

 

 

9. Ш. Ониани

129



 

 

 

д2

 

 

d

 

 

pt3

(x, p)-f(0)

= a3

{L[t3(x, t)]} =

a3

J3(x,

p).

(5.16)

Изображение

функции

t(x, p) не зависит от времени т, по­

этому

частная

производная

от оригинала функции t3(x, х) заме­

няется

обыкновенными производными для изображения

13(х, р).

В результате

дифференциальное уравнение

в частных

производ­

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

для

изображения.

 

 

 

 

 

 

 

 

 

 

 

 

С учетом

условия

(5.3)

можно

написать

 

 

 

ф ,

р )

- ~

1а(х,р)

+

~

= °-

 

 

 

(5Л7)

 

 

 

 

и3

 

 

 

 

и3

 

 

 

 

 

 

 

Теперь применим

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

х

[37,91]:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

р)]-

ЬЛШ

Р)\ + L x

t3

,

0,

(5.18)

 

 

 

 

а з

 

 

 

 

 

 

 

 

 

 

 

' о

 

 

 

 

 

 

 

 

С1-,

 

 

 

 

 

 

 

 

 

 

 

Ьх[ф,

р)] = рзф3(Р,

 

р)-Р<р(0)-

ф'(0),

 

 

(5Л9)

где

ср(0) —значение

искомой потенциальной функции на нейтраль­

 

 

 

ной плоскости в данный момент времени.

 

 

 

Согласно

условию

 

(5.6),

ср'(О) =

0 и <р(0) = const =

tx (где

^—некоторая

постоянная,

равная

температуре

закладки

иа ней­

тральной плоскости), поэтому уравнение (5.18)

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

в

виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р'Ф3(Р,

p)-Ptx-

 

 

" ^ Ф з ( Л

Р)+

=

0

(5.20)

или

 

 

 

 

 

 

 

 

 

 

 

 

 

Решая уравнение (5.21) относительно Ф3 (Р, р), получаем окончательное изображение функции t(x, х) для обеих переменных:

ФДР, Р) = —

V " ^

 

I

 

 

Р N '

С 5 ' 2 2 *

Р 2 _

Р

\

Р 2

-

 

 

а3

 

 

 

а3

 

Обратное изображение по оператору Р первого члена правой части последнего выражения берется непосредственно из таблицы изображений [37, 61, 911

130


I

-

 

 

 

 

^ch у

±-

Xm

(5.23)

 

Для нахождения обратного изображения по Р

второго

чле­

на

того же выражения

воспользуемся

теоремой разложения,

сог­

ласно

которой

 

 

 

 

 

 

 

 

 

tlla3

А

 

В

С

 

 

 

'

(

- - i f

т

 

 

 

 

 

( 5 ' 2 4 '

та,к как корни

характеристического

уравнении

 

 

 

 

 

 

 

 

 

 

г—

 

 

 

 

 

 

 

 

 

 

р

 

 

 

 

Определение значений постоянных коэффициентов при­

водит к следующим выражениям:

 

 

 

 

 

 

 

 

 

/3

^3

^3

 

 

 

 

 

 

А =

 

р

, В = 7^ ,

С = — •

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя

разложение (5.24) с полученными

 

значениями

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

А,

В,

С в выражение (5. 22),

получаем

 

 

 

Ф - ' Р ' ^ ^ _ ( ]

/ Т Г + ^ +

 

 

Н

7 = ^

(5.25)

Для каждого члена последнего выражения обратное изо­ бражение берется непосредственно из соответствующей таб­ лицы. Следовательно обратное изображение функции Ф3 (Р, р) по оператору Р будет

131