Файл: Ониани, Ш. И. Тепловой режим глубоких шахт при гидравлической закладке выработанного пространства и сложном рельефе поверхности.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^ , |
С = — • |
|
|
|
|
|
|
|
|
2р |
2р |
|
|
|
|
|
|
Подставляя |
разложение (5.24) с полученными |
|
значениями |
||||||
коэффициентов |
А, |
В, |
С в выражение (5. 22), |
получаем |
|
|||||
|
|
Ф - ' Р ' ^ ^ _ ( ] |
/ Т Г + ^ + |
|
|
Н |
7 = ^ |
(5.25) |
Для каждого члена последнего выражения обратное изо бражение берется непосредственно из соответствующей таб лицы. Следовательно обратное изображение функции Ф3 (Р, р) по оператору Р будет
131