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

•в силе. Но (в решениях для угольного массива с невозмущен­ ным температурным полем и слоя закладки толщиной, допу­

(6.27)

При более сложном виде начального распределения тем­ пературы в прилегающих к закладке массивах (в практике' встречаются случаи, когда даже для приближенного описания начального температурного поля требуется сложная матема­ тическая зависимость) приведенными решениями не может быть описана температурное поле системы трех тел.

169>

Г Л А В А V I I

РЕЗУЛЬТАТЫ РАСЧЕТА ПО ПОСТРОЕНИЮ ТЕМПЕРАТУРНОГО ПОЛЯ СИСТЕМЫ У Г О Л Ь - ЗАКЛАДКА — ПОРОДА

§ 1. Температурное поле системы при экспоненциальном

начальном распределении температуры в среде

Для решений (5.60) — (5.63) Вычислительным центром АН ГССР составлена специальная программа. По этой про­ грамме на ЭВМ БЭСМ-2 работниками Вычислительного центра были проведены расчеты с целью построения темпе­ ратурного поля системы уголь—закладка—порода для крае­ вых условий, приведенных в главе V.

верхние песчаники, а с другой — липтобиолнты и 'угли пласта I I I . Так как полученные нами результаты не учитывают неод­ нородности материала, то для тепловой характеристики иолуограниченных тел взяты теплофизич'еокие .свойства аргил­ литов и липтобиолитов шл'аста «0,90». При этом приняты сле­

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

170

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

4цм

7

1

JO

s

i

• 3

i

10

Of.n -u -to -a - 6 - 4 - s o 9 4 в в ar.«

Рис. 49. Результаты расчета на ЭВМ БЭСМ-2: 1 - х = 500 час; 2—1=1000 час; 3—1=3000 час; 4—т=6000 час; 5—т=12000 час; 6—начальное распре­ деление температуры в охлажденной зоне; 7—естгствзиноз распределение температуры; X—точка минимальной температуры

вильно 'Описывают температурное поле данной системы тел. При малых значениях времени (т = 500—1000 час) на границе раздела тел закладка—липтойиолиты происходит резкое изме­ нение градиента температурыЗатем с увеличением времени кривые постепенно оглаживаются и изменение температуры по .координате происходит совершенно плавно. Это свидетель­ ствует о том, что в 'начале процесса на характер теплообме­ на существенное влияние 'Оказывает неоднородность системы. С уменьшением интенсивности теплообмена это влияние умень­ шается и при больших значениях времени (т> 10000 час) ста- н овитоя несущественным.

171

Высказанное нами выше соображения по поводу величин ны температурных возмущений .и глубины проникновения этих возмущений полностью подтверждаются приведенными на рис. 49 данными. В этом случае нейтральная плоскость, проходящая через точку минимальной температуры, переме­ щается довольно медленно. В начале процесса она располо­ жена почти в середине закладки. Затем медленно переме­ щается в сторону липтобколитов и при т^г3000 час занимает практически стабильное положение. Этим подтверждается вывод о том, что влиянием перемещаемости системы отсче­ та можно пренебречь.

Сравнение результатов расчета на поверхности раздела системы (x — li,2) и при нулевом и больших значениях абсо­ лютной величины координаты показывает однозначность ре­

не соответствует реально возможной величине (заведомо не­

бок и погрешность расчета. По мнению програмистов ЭВМ, причиной этого является при рода самих решений. В решениях (5.60) — (5.63) при av2— 1 получаются произведения очень больших величин на очень малые, и поэтому незначительная погрешность в определении малых величин может вызывать большую погрешность в конечном результате расчета. С та­ ким объяснением поведения машины БЭСМ-2 трудно согла­

ся, а малые — уменьшаются) погрешности вовсе исчезают и. если х— / 1 ! 2 >7, 5 м, машина выдает результаты расчета, име­ ющие вполне приемлемую точность. В случае справедливости замечания программистов ЭВМ картина, полученная при рас­ смотрении результатов расчета, должна быть обратной (по­ грешности расчета с увеличением координаты должны были

172

повышаться, а не уменьшаться). Кроме того, при т-»-оо как это было показано в главе V, тепловое состояние системы стремится к восстановлению естественного распределения тем­ пературы (произведения очень больших величин на очень ма­ лые при t-yoo имеют нулевой предел) и при достаточной точности расчета температура любой точки системы для лю­ бых значений времени имеет вполне определенную конечную величину.

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

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

§ 2. Метод приближенного расчета

Разработка более рациональных методов окончательного решения задачи производилась в двух направлениях — изыс­ кания путей упрощения полученных громоздких точных ана­ литических решений и применения аналоговой техники.

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

173

[69, 97]. На первой стадии температурное поле системы оп­ ределяется начальным распределением температуры и тепло­ вой активностью тел. Эта стадия не зависит от граничных ус­ ловий теплообмена. Стадия регулярного теплового режима устанавливается после истечения довольно большого проме­ жутка времени, и при этом влияние начального распределе­ ния температуры незначительно, а зависимость температуры любой точки системы от времени имеет простой экспоненци­ альный характер. На процесс теплообмена определяющее вли­ яние оказывают теплофизические .свойства тел. Третья стадия

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

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

174