х
I
X, |
0,3 |
0303 |
\ |
У] |
0,9 |
0399 |
j |
У yz |
IJ |
0,915 \ |
|
х2 |
1,3 |
0,995 |
|
|
2,0 |
1,000 |
|
vz =aм ( |
-Kxi) (куг ~Kyj) j |
Рис. 112.
-Схема к расчету оседания точки земной поверхности по таблице значений функции К21 со ставленной по аргументу i = r/R [24, 414, 457]
работ, а затем суммированием этих площадей наклона от нуля до соответству ющего положения забоя находится объем отделенной части области влияния для всех положений забоя, подробно описан в литературе [370].
Оседание можно определить как произведение табличных значений раз ностей, если определить значения i = r/R при помощи листа кальки с нане сенными на нем окружностью площади полной подработки, центральной точ кой Р, для которой производится расчет, и двумя взаимно перпендикулярными осями, разбитыми на десятые доля радиуса полной подработки i?, и взять из таблицы соответствующие значения К г, Вычисление производится по фор муле
Vz — vz п (КХ 2 &xl) (Ку2 |
K yi)0 |
|
(169) |
|
Для показанного на рис. 112 примера эта формула примет вид |
|
|||
vz = 120 (0,995 - 0,303) (0,915 - 0,399) = 21,1 |
см. |
|
||
Необходимость поиска табличных значений |
[290] отпадает, если на круге, |
|||
изображающем |
площадь полной подработки, |
нанесена сетка К г |
[457]. |
|
Процесс расчета оседаний можно упростить и ускорить при помощи а н а |
||||
л о г о в о г о |
м е т о д а , |
выполняя подсчет |
площадей влияния, |
например, |
фотоэлектрическим способом |
[207, 383]. Для этой цели используется специаль |
|||
ный прибор (рис. ИЗ), проецирующий световую плоскость при помощи объек тива на оптический фильтр, расположенный перед фотоэлементом. При помощи кольцеобразных поляроидных пленок оптического фильтра интенсивность светового потока регулируется таким образом, что степень ослабления освещен ности в разных точках изменяется в соответствии с изменениями коэффициента влияния е, так что проходящий через различные по площади, но равноценные по степени влияния участки фпльтра световой поток возбуждает в фотоэлементе одинаковый по силе фототок. Если выполненный на прозрачной основе план горных работ, на котором отработанные участки пласта закрыты черной бума гой, передвигать по светопульту так, чтобы нанесенная на плане расчетная
Рис. ИЗ.
Принцип действия оптического аналогового при бора для расчета оседаний и горизонтал!>иых сдвижений [205, 207]:
1 — фотоэлемент; 2 — фильтр; з — объектив; |
бу |
мага; 5 — освещаемая поверхность |
|
5
точка земной поверхности совместилась с нанесенным на светопульте цецхром проекции, то присоединенный к фотоэлементу гальванометр покажет отСдет, соответствующий оседанию, а если применен дифференциальный элемент^ то и соответствующий горизонтальному смещению. На шкале гальванопара предварительно должно быть установлено значение полного оседания.
Этот метод может быть применен и для расчета сдвижений нри наклоном залегании пласта. В этом случае применяется способ центрального проециро вания, при котором прибор для измерения светового потока, снабжений широкоугольным объективом, наклоняется на угол ц, а изображающая по^ерХность очистной выработки непрозрачная бумага устанавливается в плоск^сХИ наклоненной на угол, равный углу падения пласта [207]. В ПНР также разра ботан фотоэлектрический аналоговый прибор, при помощи которого могут §Ь1ТЬ вычислены и построены кривые профиля мульды оседания [109]. В кач^стве оптического фильтра в этом приборе используется вращающийся диск ^ 30_ нально расположенными щелями.
Известен также сконструированный в виде опыта механический и т е р а тор, по устройству сходный с планиметром и действующий по принципу KoJieCa трения [167]. Счетчик этого прибора показывает величину коэффициента в^иЯ_ ния е, если обводный штифт прибора провести по контуру очистной выраб^тКИ на плане. Прибор оборудован тремя передачами изменения масштаба и функциональными передачами, при помощи которых он может быть устано^яеи „ на выбранную функцию распределения.
8.3.5.
Расчет оседаний с помощью ЭВМ для горизонта^, залегания пласта
Процесс подсчета элементов площади очистной выработки, оказывадццих влияние на земную поверхность, может быть запрограммирован для Рас^ета оседаний с помощью ЭВМ, если интеграционная сетка выполнена в виде ^ ткя
Ь |
Ь |
У |
Ряс. 115.
Интегрирование по площади очистной выработки с контуром в виде много угольника, разбитой на трапеции
Рис. 114.
Полосы равного влияния для интеграционной сетки квадратной формы, построенной для расчета оседаний [75]:
1 — разрез тела влияния
из 200 элементарных площадок равного влияния, каждая из которых соответ ствует 0,5% е, и введена в память вычислительной машины. Координаты вер
шин треугольников, на которые разбивается площадь |
выемочного участка, |
фиксируются в закодированном виде на перфокартах. |
После поиска интегра |
ционной сетки, соответствующей заданным значениям |
угла падения пласта |
и глубины разработки, а также после преобразования координат ЭВМ сопо ставляет отдельные участки площади очистной выработки с центрами соответ ствующих элементарных площадок сетки и суммирует число центров тяжести, оказавшихся в пределах контура очистной выработки. Полученное число, умноженное на 0,5% и на значение коэффициента влияния, даст искомую вели чину оседания в данной точке земной поверхности [193].
Была также разработана программа вычислений при помощи ЭВМ для метода расчета оседаний, основанного на представлении о симметричном рас пределении влияния элементарных площадок очистной выработки на оседание вышележащих элементарных площадок породного массива (см. рис. 105) [75]. Если применяемую при графическом способе расчета квадратную сетку, соот ветствующую площади полной подработки, разбить на 100 прямоугольных площадок равного влияния, каждая из которых соответствует 1 % оседания, то мы тем самым разобьем на такое же число частей вертикальными секущими плоскостями колоколообразное «тело влияния» (рис. 114), после чего подсчет числа влияющих элементарных площадок может быть выполнен способом интегрирования но площади при помощи ЭВМ. Этот метод расчета оседаний, принятый на шахтах акционерного общества «Рурколе» в качестве стандарт ного [261], основывается на уравнении (159), дополненном параметрами с и ку
в следующем виде: |
(170) |
dvz= с ехр ^ — ■% -] . / ' ) dx du |
7 заказ 74 Л
для бесконечно малой составляющей оседания dvz, обусловленной влиянием элементарной площадки очистной выработки на точку Р земной поверхности с координатами х = 0 и у = 0. Если ввести полярные координаты
£ = pcos(p и y = psincp, |
(171) |
то при интегрировании бесконечно малых составляющих dvz для очистной выработки, имеющей в плане форму круга с радиусом г, получим
2л г г
(г) = с j |
j |
р exp ( — |
dp d(p = 2яc J p exp ( |
- |
dp. |
(172) |
||
0 |
0 |
|
|
|
0 |
|
|
|
Если ввести |
обозначения |
• |
Z/C* |
dp==~r-dp |
и dp = |
— |
то оседание |
|
в данной точке |
|
' i t * |
|
p |
|
|||
|
|
|
|
г |
|
|
||
|
|
|
|
|
|
|
|
|
у2 (г) = 2лс&2 ^ exp ( — р) dp = |
— 2ясА:2 exp ( — -т ^ ) | = |
|
|
|||||
|
|
|
|
|
|
О |
|
|
= 2яс/с2[ 1 - е х р ( - - £ • ) ] • |
|
|
|
|
(173) |
|||
Теперь можно приближенно определить параметры с и А, входящие в вы ражение (170), через граничный угол. Подставив в выражение (173) г = получим
Уг (оо) = уг п = 2лс/г2, откуда с = , (174)
а подставив вместо г величину «практического» радиуса площади полной под работки i?np, измеренного при площади очистной выработки, равной 0,99у2п (т. е. для упоминавшегося выше «практического» граничного угла), получим571
|
|
|
|
R 2 |
|
vz (7?пр) = 0,99 иг n = аМ |
1 |
|
пр |
(175) |
|
|
2к* |
||||
|
|
|
|
|
|
и если принять ехр ( — |
^= 0,01, то получим |
|
|||
|
\ 2к2 ) |
|
|
|
|
откуда |
1 0 0 = 2 ’ |
|
|
||
|
|
|
|
|
|
к = 4 -V Tie Д = 0,3295 |
tgY |
. |
|
|
|
2 |
|
|
|
|
|
Oh6 щ е е |
и п т е г р и р о в а н и е экспоненциальной функции в |
выра |
|||
жении (173) |
даст в прямоугольных координатах |
|
|||
Уг=сИехр(~ ~ m * ) dxdy |
(176) |
|
или, поскольку при / / = const можно считать, что к = const,
пг -=с J exp ( - JL -} dx j exp ( - J O } dy. |
(177> |
Для показанной на рис. 115 очистной выработки с контуром в виде много угольника, разбитой на трапеции с отсчетом координат как в положительном, гак и в отрицательном направлении (на рисунке — заштрихованные площади за пределами контура выработки), интегрирование по площади для отдельных
трапеций выполняется н о |
ф о р м у л е |
С и м п с о н а |
|
|
А — У^ ’2т (Xl ^Х2 |
~'r ^г4 |
+ |
/iX2m4* x2m+l), |
(178)* |
причем трапеции разбиваются через интервалы, равные 2 т, с шагом шири
ной |
2т |
При этом |
|
||
|
|
* |
|
||
где |
х {у )^ а у + Ъ, |
(179) |
|||
|
|
|
|
|
|
|
а = |
Гг~ Д1 |
и b ^ Xl- Х'2~ Х1 ух. |
(180> |
|
|
|
Уг— У1 |
Уг~У1 |
|
|
|
Интеграл |
составляющей оседания над площадью трапеции, |
отнесенный |
||
к точке земной поверхности Р (0,0), в соответствии с формулой (177) имеет вид
Уг |
* (У) |
|
Vz(трап) = J |
j ехР ( — - £ - ) dz] exp ( - -J^-) dy. |
(181) |
Ух |
О |
|
При использовании таблиц значений вероятностной функции, имеющей вид
х
= у = - j*exp( — t2)dt,
О
заключенный в квадратные скобки интеграл, входящий в выражение (181), записывается в виде
(182>
о |
|
|
|
если подставить значения |
|
||
Х = х(у), |
t —— |
и dt |
—4 — dx. |
У’ |
1 |
2 А* |
1 2 к |