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

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

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

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

Добавлен: 17.10.2024

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

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

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

Г Л А В А V I I I

МОД Е Л И РОВАН И Е НЕСТАЦИОНАРНОГО^. ТЕМПЕРАТУРНОГО ПОЛЯ

СИСТЕМЫ УГОЛЬ—ЗАКЛАДКА—ПОРОДА НА ИНТЕГРАТОРЕ ЭИНП-3/66

§ 1. Методика моделирования

Из предыдущих глав (.главы V н VI) следует, что анали­ тическое решение уравнения теплопроводности па.раболичесiKoro типа при неоднородных начальных условиях и .гранич­ ных условиях четвертого рода приводит к очень сложным зависимостям, 'неудобным для практических инженерных расчетов без использования ЭВМ, тем более, что приведен­ ные выше решения были получены после некоторой схема­ тизации задачи, при которой реальная (неоднородная система тел была заменена системой с однородными телами. Если же

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

имеет бо­

лее сложный характер, то поставленную задачу

практически

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

ходящей последовательности выемки слоев. Поэтому

пол­

ное и окончательное решение поставленной задачи было

вы-

190

полнено

методом математического

моделирования (методом:'

электротепловой

аналогии).

 

 

 

Следует отметить, что электрическое-.моделирование уже

широко

используется при

решении

многих технических за­

дач, и

область

применения

этого

удобного,

•сравнительно

простого и наглядного метода быстро расширяется

за счет

унификации существующих

и создания новых,

более

совер­

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

•на электротепловой

аналогии. Решение

краевых

задач

как

стационарной, так' и нестационарной теплопроводности

наи­

более удобно проводить методом электротепловой

аналогии.

Коротко

изложим

некоторые

 

вопросы

моделирования

нестационарных

тепловых

полей

на

интеграторе

ЗИНП-3/66

конструкции Института математаки АН УССР [130,

131].

Для

неоднородной среды и одномерной задачи

дифферен­

циальное

уравнение

 

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

 

можно

переписать

следующим

образом:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д

 

1

 

dt

=

 

 

 

 

ед-тг,

 

 

 

(8-1)

 

 

дх

L Rt(x)

 

дх

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где Rtтермическое

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

среды,

кусочно-постоянная

 

функция

координаты

х;

 

 

 

 

 

 

 

 

 

Ctобъемная теплоемкость

среды,

кусочно-постоянная функ­

 

ция координаты

х;

 

 

 

 

 

 

 

 

 

 

 

1;— время протекания теплового процесса.

 

 

 

 

В такой же среде одномерное

движение

электрического

заряда

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

выражением

 

 

 

 

 

 

 

 

 

 

 

1

ди

=

 

 

ди

ед,-,

 

 

 

(8.2)

 

 

дх

~Rjx) ~д~х~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где Ruудельное

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

электропроводной

среды,

ку­

 

 

сочно-постоянная функция координаты х\

 

 

 

Сиудельная

электрическая

емкость,

кусочно-постоянная

 

функция

координаты

х;

 

 

 

 

 

 

 

 

 

хи~время

протекания электрического

процесса.

 

 

 

Из

анализа уравнений

(8.1)

и

(8.2)

следует, что

электри­

ческое

поле,

которое

описывается

соотношением (8.2), явля-

191.


•етея моделью теплового поля, выраженного зависимостью (8.1), только в том случае, если соблюдено условие подобия

ДО)

[28, 131,

134]

 

 

 

 

 

 

 

 

 

 

xt

 

 

 

CrRrl*t

 

 

 

 

 

 

 

Т ~

=

^

=

C~R

~ ^

= Nc'NR-m.

 

(8.3)

 

Так как

геометрические

размеры

исследуемого

теплового

т о л я

в

большинстве

случаев

значительно

больше

размеров

модели

(/, >

/„),

а

постоянные

 

времени

тепловой

системы

намного

превосходят

произведение

постоянных RC

электри­

ческого

полк

[Rj Ct

> РиСи),

 

 

то

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

электри­

ческие

процессы

в

модели

будут

проходить намного быстрее

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

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

и

громоздких

сеточных интеграторов не всегда

целесообраз­

но

и доступно.

В настоящее время выпущена

опытная пар­

тия принципиально нового, удобного, простого и эффективно­ го прибора ЭИНП-3/66 для электрического моделирования двухмерных нестационарных потенциальных полей на электро­ проводной бумаге с распределенной емкостью [130, 131]. Один экземпляр из этой опытной партии был изготовлен по заказу ИГМ АН ГССР.

Интегратор состоит из двух основных частей: стенда электрической модели и стойки измерительных и питающих устройств. Электрическая модель исследуемой системы

строится

в виде

плоского

конденсатора, обкладкой которого

с одной

стороны

является

основная модель, построенная из

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

другой модель начальных

условий, изготовленная

из низко-

- омной электропроводной .бумаги. В качестве диэлектрика

при­

меняется

конденсаторная

пленка разной толщины

(от

1

до

20 мкм).

Средняя удельная емкость модели определяется

за-

192


.ранее по специальной методике [132]. Изготовленная модель помещается на стенде и прижимается к нему специальными планками. Измерительное устройство интегратора имеет пять фиксированных скоростей отклонения луча, которые соответ­ ствуют длительности процесса, моделирующего исследуемое

явление в натуре (1, 2, 10, 50 и 500 мсек),

причем

наиболь­

шая точность достигается

на пятой шкале

(500 мсек).

С уменьшением продолжительности горизонтальной раз­

вертки электроннолучевого

индикатора снижается

точность

моделирования. На первой и второй скоростях отклонения лу­ ча (1 и 2 мсек) оказывается инерционность измерительного устройства и качественная картина процесса существенно искажается. Поэтому при работе на этих скоростях не всег­ да обеспечивается требуемая точность. Ощутимые погреш­ ности получаются и на третьей шкале. Таким образом, при моделировании с помощью интегратора ЭИНП-3/66 целесо­ образно пользоваться пятой или, по крайней мере, четвертой шкалой. А вообще, чем медленнее протекает процесс в мо­ дели, тем точнее будут результаты моделирования. В данном случае процесс теплообмена между закладкой и средой про­

должается бесконечно долго. Но для нас наибольший

интерес

представляет подстадия собственно

начальных

условий

(т < 5000—6000 час), а с учетом запаса

по времени

можно

брать т = 20000 час. Следовательно, масштаб времени выби­

рать не приходится его значение ограничивается

условиями

моделирования

и возможностью

моделирующего устройства.

В данном

случае при работе на пятой шкале

масштаб

времени

 

 

 

 

 

Nx =

20000-3600

 

 

 

^

= 144-100,

 

т. е. исследуемое явление в натуре протекает в 144 миллиона раз медленнее, чем в модели.

 

Размер

стенда составляет 0,75 X 1,20 м. Температурные

возмущения,

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

закладки,

в

течение 20000 час могут распространяться в глубь массива

на

расстоянии 30—40 м, т. е. общая минимальная

толщина

горного массива, вовлеченного в одномерный нестационарный

13. Ш. Ониани

-

193


процесс теплообмена с закладкой, 'Составляет 60 м. Следо­ вательно, линейный масштаб моделирования не может быть меньше 50, т. е. = 50. Максимальное удельное сопротив­ ление электропроводных бумаг, применяемых в настоящее время для моделирования, не превышает 1000 ком на квад­ рат. Допустим, модель угольного массива пласта IV изготов­ лена из такой бумаги, тогда масштаб сопротивления составит

R R

149,25

^ = X" =

T6ooooo-= 1 4 'S 2 5 -1 { ^B '

где Ry = 149,25 см-грах/вт —удельное термическое сопротивле­ ние угольного массива (см. табл. 13);

/?б =1000 ком—электрическое сопротивление моделирующей бу­ маги на квадрат.

Следовательно, масштаб емкостей, согласно выражению

(8.3),

будет

 

 

 

 

 

 

 

 

 

 

NT

_

14,4-10-

 

 

 

 

 

-

лг

\W

~

l i r o - i n - ; о с л п

 

=" 3 3 , 5 9 5 - 1 0 '

 

 

N^NY~

14,925- 10"5-2500

 

 

откуда

удельная

емкость модели

составит

 

 

С"=

Щ =

 

38,5С5 8 10^ -

4663 . 8 . 10 - ^/см - ,

 

где

Ct = 1,8

дж/см3 град—объемная

теплоемкость

угольного

 

 

массива (см. табл. 13).

 

 

 

 

 

Существующие

конденсаторные

пленки толщиной 1—20

М'км

при помещении между

высокоомной (RQ= 1,0 Мом) и

низкоомной

 

 

0,6 ком)

электропроводными

бумагами

на стенде электрической модели интегратора создают элект­ рическую емкость не более 30,0 пер/ом2. Требуемую по усло­ виям моделирования удельную электрическую емкость не мо­ гут создать самые высококачественные конденсаторные плен­ ки. Таким образом, задачи, подобные рассматриваемым, не могут быть решены на стенде электрической модели .интегра­

тора. Поэтому аналоги .исследуемой

системы строились нами

с помощью

моделирующей цепи

(#С-.контура)

с сосредо­

точенными

параметрами. Это тем

более целесообразно, что

поставленную задачу можно рассматривать как

одномерную

194