При освоении новых рудных месторождений и реконструкции действующих, как правило, отсутствуют необходимые экспериментальные данные о напряженном состоянии массива. Такое положение создаёт определённые трудности при разработке проектных решений и приводит к осложнениям для нормального технологического процесса выемки руд. Проведение экспериментальных исследований по оценке напряжённо-деформированного состояния и прогнозу поведения массива в процессе развития горных работ не всегда является целесообразным и экономически оправданным. Кроме того их осуществление требует значительных временных затрат. Это относится в большей степени и к разработке небольших по запасам месторождений, где в создании крупной экспериментальной базы нет необходимости.
В связи с этим возникает необходимость в разработке экспресс-методов прогноза геомеханических процессов, происходящих в массиве горных пород. К числу таких методов относятся инженерные методы расчёта, основанные на математическом решении ряда фундаментальных задач механики горных пород и механики деформируемого твёрдого тела.
При исследованиях напряжённо-деформированного состояния массива горных пород вокруг выработок широко используются численные методы расчёта. Среди них наибольшее распространение получили методы конечных элементов и граничных интегральных уравнений. Достоинство метода конечных элементов заключается в возможности учёта локальных неоднородностей массива, так как каждому из элементов, на которые разбивается массив в области решения, можно присвоить собственные механические характеристики. При решении задач механики горных пород методом конечных элементов, можно осуществлять анализ напряженного состояния каждого из элементов и в соответствии с этим напряжённым состоянием пересчитывать все механические характеристики элементов.
Применение метода конечных элементов позволяет рассматривать произвольную неоднородность среды, анизотропию, учитывать проявление пластичности, вязкости, нелинейность деформационных свойств, учитывать структурные особенности реальных массивов горных пород, такие как трещинноватость, трещиновато-пористость, наличие нарушений сплошности массивов в виде трещин контакта, а также учитывать технологические факторы, например, последовательность очистных и закладочных работ.
Другие направления численных методов решения задач о напряжённо-деформированном состоянии среды являются методы граничных элементов, в которых решение уравнения удовлетворяется на границе исследуемой области.
Наряду с широкими возможностями численных методов, их использование требует специальной подготовки и недоступно большому кругу специалистов.
Это в значительной степени устраняется при использовании предлагаемого метода оценки напряжённого состояния массива горных пород крутопадающих месторождений малой и средней мощности.
Разработанный инженерный метод позволяет прогнозировать нагруженность рудного массива с учётом динамики развития горных работ и тем самым учитывать технологическую последовательность развития горных работ, тип применяемой закладки, прочностные свойства пород. Метод расчёта применим в том случае, когда одно из главных напряжений нетронутого массива нормально к плоскости простирания. Наибольший практический интерес представляют напряжения . Эти напряжения определяют нагруженность рудного и закладочного массивов, которая сопоставляется с несущей способностью этих массивов. Известно, что рассланцевание пород ориентированно вдоль плоскости простирания рудного тела и происходит в условиях высокого давления, действующего перпендикулярно к направлению перемещения пород. Это является основанием для использования метода при оценке напряжённого состояния месторождений жильного типа, приуроченных к сланцеватым текстурам.
Среди разнообразия рудных тел в горнодобывающей практике встречаются такие, которые характеризуются относительно простой формой (без разветвлений) и не имеют резких изменений мощности и условий залегания. Особенность их разработки заключается в том, что очистные камеры имеют вид относительно узких, протяжённых выработок, расположенных примерно в одной плоскости. Поперечный размер камер значительно меньше продольных размеров и (, ).
Выемка ценных руд и руд в сложных горно-геологических условиях сопровождается последующей закладкой выработанного пространства. Определение оптимальных параметров очистных камер, рудных и искусственных целиков, прочности закладочного массива и т.п. является важной практической задачей, решение которой должно основываться на экономической оценке производственных затрат с учётом фактора безопасности горных работ. Фактор безопасности в свою очередь тесно связан с решением проблем горного давления, исследование которых, осуществляемое средствами и методами геомеханики, связанно главным образом с анализом напряжённого состояния массива горных пород. Ниже изложен инженерный метод расчёта изменения нагруженности рудного массива и закладки. Метод достаточно прост и не требует обращения к большим ЭВМ. Например, в приведённых ниже примерах для расчётов использовался карманный калькулятор.
Пусть известно напряжённое состояние нетронутого массива горных пород, которое в простейших случаях может быть определено из анализа геологических условий залегания рудного тела, а в общем случае в результате непосредственных шахтных измерений. Создание выработанного пространства приводит к перераспределению напряжений в массиве горных пород. Суть перераспределения состоит в том, что "вынимаемая" часть рудного массива исключается из силового соотношения структурных элементов массива горных пород. При этом нагрузка, первоначально воспринимаемая "вынимаемой" частью массива, передаётся на оставшуюся часть рудного массива, вмещающие породы или закладочный массив. Нагрузка распределяется также среди специально созданных несущих конструкций, к которым относятся рудные и искусственные целики.
Если руда и вмещающие породы близки по упругим свойствам, а одно из главных напряжений нетронутого массива действует нормально к плоскости простирания рудного тела, приближенная оценка распределения нормальной нагрузки со стороны вмещающего массива между несущими элементами может быть осуществлена на основе известной формулы распределения напряжений в упругой среде при действии в среде сосредоточенных силовых источников. В этой оценке очистные камеры моделируются лежащими в одной плоскости разрезами в однородной упругой среде.
|
Для простоты рассмотрим сначала в условиях плоской деформации одиночную выработку без закладки выработанного пространства.
Рудное тело на рисунке не показано. Напряжения и являются главными.
Распределение напряжений на прямой можно найти, решив задачу теории упругости о разрезе, нагруженном нормальной нагрузкой . Эта нагрузка равна нормальным напряжениям , действовавшим вместе проведения разреза. В частном случае горизонтальной одиночной выработки при однородном исходном поле напряжений на разрезе, моделирующем выработку, распределённая нагрузка постоянна и равна .
Для инженерного расчёта напряжений вдоль оси разобьём разрез и линию его продолжения в среде вдоль оси на линейные элементы (в плоском случае). Выбор размера зависит от степени необходимой точности оценки: чем меньше , тем точнее оценка. Однако следует учесть, что с уменьшением растёт объём вычислений.
Пусть разрез, моделирующий выработку, содержит элементов, которым приписан индекс , а среда элементов с идексом . Последние элементы будем для определённости называть "грузонесущими". Введя линейные элементы будем далее рассматривать средние напряжения в элементах или силы. Задаваемые в центре элементов разреза силы обозначим . Искомые силы, приложенные в центре грузонесущих элементов, обозначим . Очевидно . Таким образом наш расчёт заключается в определении нормальных сил в грузонесущих элементах по известным силам , задаваемым на элементах разреза.
При предположении упругости среды сила Pj в каждом элементе является суммой соответствующих воздействий от совокупности сил
Pj=∑i=1..Nf(Qi)
Функция f отражает закономерность передачи силового возмущения от одиночной силы. Об этой функции можно говорить как о своеобразной функции точечного возмущения для разреза. Используя соответствующее соотношение, полученное в работе "Предельное равновесие хрупких тел с трещинами", Панасюк В.В., Киев, 1968г., напряжение σy(x) в центре грузонесущего элемента j можно записать в виде:
σy(x*j)=∑i=1..N{Qi[a2-x2*i]1/2/(π(x*j-x*i)[x2*j-a2]1/2)}
В целях разумного облегчения расчётов будем искать кусочно-линейное распределение напряжения σy(x). В этом случае в пределах рассматриваемого грузонесущего элемента j напряжения σy(x) можно описать функцией вида:
σy(x)Ij=Aj+Bjx
где Aj и Bj (j=1,2,...,M) - постоянные, которые выбираются таким образом, чтобы было справедливо соотношение:
Aj+Bjx=∑i=1..N{Qi[a2-x2*i]1/2/(π(x*j-x*i)[x2*j-a2]1/2)}
Поскольку в рассматриваемом случае среднее напряжение в грузонесущем элементе j равно σyсрIj=Aj+Bjx*j, то можно записать соотношение для сил:
Pj=∑i=1..NPij
где Pij равно:
Pij=∑i=1..N{LQi[a2-x2*i]1/2/(π(x*j-x*i)[x2*j-a2]1/2)}
Величина Pij представляет собой часть силы Pj, которая обусловлена силой Qi, заданной в i-ом элементе разреза.
Для оценки распределения сил Pj вблизи края разреза, т.е. там, где распределение нагрузки вызывает наибольший интерес, формулу можно упростить и записать в виде:
Pij≈∑i=1..N{LQi[a2-x2*i]1/2/(π[2a]1/2(x*j-x*i)[x*j-a]1/2)}
В другом крайнем случае, т.е. для удалённых точек среды x*j, справедлива другая асимптотика:
Pij≈∑i=1..N{LQi[a2-x2*i]1/2/(πx2*j)}
которая соответствует известному закону убывания напряжения в упругой среде, выносимого парой сил на разрезе (плоский случай). Поскольку на практике наиболее интересно распределение напряжений в целиках и вблизи выработанного пространства, в инженерных расчётах можно использовать только первую формулу. При этом погрешность в определении напряжений σy на достаточном удалении от выработанного пространства не превысит 10%.
Рассмотрим другое силовое соотношение. Сила Qi вносит свой "вклад" во все грузонесущие элементы. Естественно, что сумма этих вкладов по всем грузонесущим элементам должна быть равна Qi. Таким образом можно записать:
∑j=1..MPji=Qi
Приведённые формулы являются основными рабочими формулами расчёта. Целесообразно их записать в ином виде, более наглядно отражающем их структуру:
Pji=Fi/(Rji[rj]1/2)
Соотношение примет вид:
∑j=1..M{Fi/(Rji[rj]1/2)}=Qi
Отсюда можно получить следующее выражение:
Fi=Qi/∑j=1..M{1/(Rji[rj]1/2)}
Методика расчёта действующих сил в грузонесущих элементах (или средних напряжений) заключается в следующем. Линейные величины rj и Rji, входящие в приведённые выше формулы, находятся геометрически в соответствии с условиями разбивки исследуемой области на элементы. Силы Qi определяются по известным нормальным напряжениям σy на контуре будущего разреза, моделирующего новую выработку. Определив по формуле неизвестные величины Fi (i=1,2,...,N), находим затем величины Pji и далее искомые силы Pj в грузонесущих элементах. Используя принцип суперпозиции решений, который справедлив для упругих сред, можно найти нормальные силы в грузонесущих элементах при произвольном числе разрезов среды вдоль оси X.
Отметим, что по указанной методике определяются в грузонесущих элементах силы, а следовательно и средние нормальные напряжения. По средним напряжениям в элементах в общем случае нельзя построить однозначное кусочно-линейное распределение искомых напряжений σy(x). Действительно, для определения 2M неизвестных постоянных Aj, Bj (j=1,2,...,M) имеется M уравнений. Кроме того в однородной среде должно выполняться условие непрерывности напряжений σy на границах элементов, что даёт ещё M-1 уравнение:
Aj+Bj(x*j-L/2)=Aj-1+Bj-1(x*j-1+L/2) j=2,3,..,M
Однако одно уравнение оказывается неопределённым. Для определённости логично принять, что в том элементе, где среднее значение σy минимально, это напряжение постоянно. В этом случае эпюру напряжений σy вдоль оси X легко построить геометрически, проводя соответствующим образом отрезки прямых от элемента, где среднее значение напряжения минимально. Каждый отрезок прямой должен проводиться таким образом, чтобы он, во-первых, в центре элемента соответствовал найденному среднему значению напряжения, во-вторых, совпадал с концом отрезка в соседнем элементе. Нетрудно понять, что при увеличении числа разбиений среды на элементы точность такой кусочно-линейной апроксимации возрастает.
Рассмотрим методику рассчёта напряжений в случае закладки выработанного пространства. Закладочная смесь является более податливым материалом в сравнении с массивом горных пород. При нагружении в режиме заданных смещений горная порода и закладочная смесь естественно деформируются одинаково, однако закладочная смесь нагружается слабее во столько раз, во сколько её модуль деформации Ез меньше модуля деформации горной породы Е. Численными расчётами с помощью МКЭ установленно, что при очередном расширении выработанного пространства зависимость средних напряжений в закладочном массиве вблизи очистной камеры от отношения Ез/Е определяется кривой.
Эту зависимость можно использовать в расчётах средних напряжений в массиве закладочной смеси. Для простоты в инженерных расчётах эту зависимость удобно считать линейной (штриховая линия). Таким образом далее будем считать, что искусственный несущий элемент из закладочной смеси по сравнению с находящимся в тех же условиях рудным целиком воспринимает нагрузку, пропорциональную отношению Ез/Е.
Следовательно, в общем случае, когда часть грузонесущих элементов j, расположенных по оси X, имеет другой модуль деформации, чем массив горных пород, вместо приведённых выше соотношений следует использовать следующие соотношения:
Pji=(Eзj/E)*(Fi/(Rji[ri]1/2)) и Fi=Qi/∑j=1..M{(Eзj/E)*(1/(Rji[ri]1/2))}
Очевидно, при Езj=Е, что эти формулы следуют из предыдущих и наоборот.
Рассмотренный выше метод расчёта перераспределений напряжений в случае выемки с закладкой целесообразно связать с методом последовательных циклов (МПЦ), разработанным С.В.Кузнецовым и В.Н.Одинцевым. Метод последовательных циклов учитывает последрвательность расширения и заложения выработанного пространства. В этом случае исследования перераспределения напряжений можно представить в зависимости от динамики развития горных работ, учитывающей напряжённо-деформированное состояние массивов горных пород и закладки, а также и технологические факторы влияния. При этом имеется возможность учёта изменения деформационных свойств закладочной смеси во времени или под нагрузкой, например, повышение жёсткости и упрочнение закладочного массива. Вычислительные трудности при этом практически не возрастают. Естественно, для учёта геологических свойств закладочной смеси необходимо задать зависимости, отражающие соответствующие свойства, например, Езj(σyср) или Езj(Tj), где Tj - время, прошедшее с момента уложения смеси в очистной камере. Анализ перераспределения напряжений в массиве горных пород и закладочном массиве рассмотренным выше методом в сочетании с МПЦ значительно повышает достоверность исследований в целом.
Рассмотрим теперь пространственный случай. Методически он мало отличается от рассмотренного выше плоского случая. В пространственном случае производится разбивка разреза, моделирующего плоскую очистную выработку, на квадратные элементы с индексом i и оставшуюся часть плоскости y=0 на квадратные грузонесущие элементы с индексом j. Сторона элемента равна L.
Как и в рассмотренном выше случае суть расчёта состоит в определении нормальных сил Pj, действующих на плоскость Y=0. Вводятся в рассмотрение величины Pji, представляющие собой составляющие сил Pj, обусловленные действием сил Qi. В пространственном случае справедливы приведённые выше соотношения, однако выражения для величин Pji и Fi имеют несколько другой вид. Для случая закладки отработанных ранее камер можно записать:
Pji=(Ej/E)*(Fi/([Rji]2[rj]1/2))
Fi=Qi/∑j=1..M{(Ej/E)*(1/([Rji]2[rj]1/2))}
Как можно видеть отличие формул состоит в том, что в пространственном случае величина Rji в квадрате. Зависимость от Rji и rj при Ej=E следует из анализа асимптотического разложения напряжений σy в пространственном случае вблизи границы разреза. Для круговых разрезов приведённые соотношения являются строгими, для разрезов с некруговой границей - приближёнными.
Методика нахождения сил Pj та же самая. Расстояния rj и Rji определяются геометрически, затем по формуле находят величины Fi и далее определяют искомые нормальные силы во всех грузонесущих элементах. Как и в плоском случае, используя принцип суперпозиции решений для упругих сред, можно найти силы в грузонесущих элементах при произвольном числе разрезов, лежащих в одной плоскости.
По известным силам в грузонесущих элементах можно построить кусочно-линейное распределение в элементах напряжений σy. Для этого следует по описанной в плоском случае методике построить кусочно-линейные распределения напряжений σy по осям X и Z. В качестве исходного элемента следует взять элемент, где среднее значение σy минимально.
Изложенную выше методику расчёта распределения нормальных к рудному телу напряжений σy можно связать с развитой в работе "Теория защитных пластов" Петухов И.М., Линьков А.М., Сидоров В.С., Фельдман И.А., М., 1983г. методикой оценки линейного размера αo предельно напряжённой области (пластической зоны) и максимума опорного давления σymax. Оценка величин αo и σymax включает определение коэффициента интенсивности напряжений KI у края разреза. Коэффициент KI, рассматриваемый в теории трещин, определяется из асимптотического представления напряжений у конца разреза (трещины) в упругой среде:
σy(x,0)≈KI/[2π(x-a)]1/2
Этот коэффициент зависит от внешних условий нагружения и размера разреза. Как видно напряжения у края разреза в упругой среде имеют неограниченный рост при x→α, что свидетельствует о реальности перехода материала у края разреза в пластическое состояние.
В упомянутой выше работе считается, что распределение напряжений в пластической зоне, т.е. от края разреза до точки близжайшего максимума напряжений σy, описываются линейным законом:
σy=σ°MAC+2K*ξ/h
За максимумом опорного давления имеет место упругое распределение напряжений. Длина пластической зоны определяется по формуле:
do=(3/4[π/2]1/2)2/3(h/2K*)2/3f(P)
где
P=0,57(σ°MAC)3h/(2KI2K*)
f(P)=(1/4)1/3{([1+P]1/2+1)1/3+([1+P]1/2-1)1/3}
В табличном виде зависимость f(P) приведена в указанной выше работе.
Таким образом, размер пластической зоны вблизи плоской выработки поперечного размера h зависит от величины прочности массива и параметра KI, определяемого из решения соответствующей задачи теории упругости. Величину максимума опорного давления можно найти, положив ξ=do.
Рассмотренный выше метод расчёта распределения нормальных напряжений σy позволяет оценить значение KI. Для этого следует рассмотреть крайний грузонесущий элемент у разреза, для определённости приписав ему индекс j=1. При той же величине среднего напряжения будем считать, что в этом элементе имеет место не линейное распределение напряжений σy, а распределение напряжений у конца разреза σy(x,0). В этом случае на границе рассматриваемого элемента с соседним грузонесущим элементом должен возникнуть скачёк напряжений Δσy. Отношение величины этого скачка к средним напряжениям в элементе невелико. Оно должно стремиться к нулю при L→0. Следовательно, при малом L величина погрешности в крайнем грузонесущем элементе мала. Значение KI по найденной в элементе j=1 силе P1 определяется из соотношения:
KI=P1[π/2L]1/2
Подставив полученное таким образом значение KI в соотношение для длины пластической зоны оценим её размер do. Возможны два случая.
Достоверная информация о состоянии горных массивов, протекающих в них геомеханических процессах, имеет основополагающее значение для решения многих проблем, связанных с созданием безопасных технологий отработки месторождений, рациональным извлечением полезных ископаемых.
Используемые при оценке напряжённого состояния расчётные экспресс-методы, являются инженерными методами нового поколения. В основе относительно простых процедур расчёта лежат математические решения ряда задач механики горных пород и механики твёрдого деформируемого тела.
Оценивая достоверность предложенного метода оценки напряжённого состояния массива и несущих конструктивных элементов, следует остановиться на следующих моментах. Это соответствие расчётной модели реальному состоянию технологических процессов натуры и сама точность расчётного метода. В этой связи необходимо отметить, что достоверность инженерного метода значительно повышается, поскольку при постановке задачи исследований имеется возможность прогнозировать динамику развития горных работ. Использование метода позволяет учесть влияние уже существующих горных работ на изменение геомеханической обстановки при исследовательной отработке горизонтов как в масштабе всего рудного тела или нескольких рудных тел, так и в масштабе очистных камер.
Точность инженерного метода расчёта концентрации напряжений в плоском и пространственном случаях определялась на примерах сравнения решения тестовых задач методом конечных элементов, а также задач по определению нагрузок на ленточные и столбчатые целики, рассчитанные Ж.С.Ержановым и др. методами теории упругости. Расхождения результатов расчёта задачи о распределении напряжения при проведении плоской выработки с прямоугольным уступом для пространственного случая, выполненные инженерным методом и найденное более точным методом не превышали 20%.
Сравнительный анализ показывает, что различия между расчётными данными не превышают 14%. Эти значения допустимы для многих задач и её величину (20%) можно принять в качестве оценки погрешности рассмотренного метода расчёта напряжений.
Опробывание инженерного метода было произведено для оценки размеров зоны опорного давления для случая частичной отработки охранных целиков стволов шахт Джезказганского месторождения и показала хорошие практические результаты.
|
Для оценки опорного давления в предохранительном целике ствола при его частичной отработке параметры ленточных целиков (а) и очистных камер (L) в расчётах варьировались. Установлено, что величина опорного давления в охранном целике в значительной степени зависит от расстояния от ствола до выработанного пространства. Проведённые расчёты позволили определить значение опорного давления, формирующегося вследствии образования камер.
Анализ расчётов показывает, что наиболшее дополнительное давление, достигает 1,5γH, отмечается у границ предохранительного целика. Однако уже на расстоянии 10м от вырабатываемого пространства дополнительная нагрузка на массив не превышает 0,5γH, а на расстоянии 15м - 0,1γH. Таким образом, на основе расчётов можно сделать вывод о том, что частичная отработка охранного целика по фактору горного давления может производиться до тех пор, пока расстояние от границ выработанного пространства ствола не достигнет 20-25м.
На базе проведённых исследований разработаны и согласованы с Госгортехнадзором Респ.Казахстан 29.09.89г. (приложение 29) дополнения к "Временным правилам охраны сооружений от вредного влияния подземных разработок на рудниках Джезказганского горно-металлургического комбината", регламентирующее условия по частичной отработке запасов в охранных целиках стволов шахт и поверхностных сооружений.
безымянный © copyright 2004 |
Опубликовано 12 февраля 2011г.