Файл: Ониани, Ш. И. Тепловой режим глубоких шахт при гидравлической закладке выработанного пространства и сложном рельефе поверхности.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