Геометрическая дискретизация на конечные элементы, конгруэнтные по форме и размерам структурным блокам разрушаемого углепородного массива

Павлова Лариса Дмитриевна

УДК 622.831.232 Ё.Д. Павлова

ГЕОМЕТРИЧЕСКАЯ ДИСКРЕТИЗАЦИЯ НА КОНЕЧНЫЕ ЭЛЕМЕНТЫ, КОНГРУЭНТНЫЕ ПО ФОРМЕ И РАЗМЕРАМ СТРУКТУРНЫМ БЛОКАМ РАЗРУШАЕМОГО УГЛЕПОРОДНОГО МАССИВА

Семинар № 10

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

Кроме того, размеры области исследования следует выбирать исходя из ожидаемого характера НДС таким образом, чтобы задаваемые граничные условия мало влияли на результаты решения задачи. Разумеется, что увеличение рассматриваемой области повысит точность решения, однако это будет сопровождаться значительным ростом затрат вычислительных ресурсов, особенно при решении пространственных задач, поэтому необходимо использование эффективных ресурсосберегающих алгоритмов.

Метод вложенных областей. Первый этап численного решения задачи методом конечных элементов (МКЭ) заключается в дискретизации модели массива горных пород на конечные элементы.

Дискретизацию области исследования

необходимо проводить с учетом того, что решение МКЭ дает в пределах конечного элемента постоянные значения напряжений. Поэтому, с одной стороны, в местах ожидаемых высоких градиентов напряжений сеть элементов следует сгущать. Однако при построении нерегулярных сеток необходимо избегать использования узких и длинных элементов. С другой стороны, применение равномерных сеток или построенных по какому-либо закону позволяет использовать близкие к правильным конечные элементы, автоматизировать расчет узловых координат и сократить объем вводимой информации. Использование же нерегулярных сеток как правило, экономичнее в отношении затрат машинного времени.

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

Количество узлов сетки определяется имеющимися вычислительными ресурсами, и в первую очередь - объемом оперативной памяти. Если используемое число узлов не обеспечивает необходимой густоты сетки в подобласти высоких градиентов, то для достижения заданной точности можно использовать метод вложенных областей [1].

Первоначально область исследования разбивается на крупные элементы и проводится расчет. Затем выделяется подобласть, которая разбивается на более мел-

Рис. 1. Схема вылеления вложенных областей

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

В промежуточных узлах на контуре подобласти, появившихся в результате ее разбиения на более мелкие элементы, перемещения вычисляются по методу источника:

Хєхі Рі Хєуі Рі Хєі Рі

єхі* =-----; єуҐ =-------; є%* =-----, (1)

X Рі ХРі X Рі

где Зш = (1 . Рд/ж рш = (Чо - Чи) + (Но -Нщ)2 +(яо - яш)2б ш^ож чоб наб яоб (о = 1бь) - исходные узловые координаты х, уі, г, (і = 1,к) - промежуточные узловые координаты; єхі, єуі, єгі - узловые перемещения исходной области.

Описанная процедура может повторяться многократно до получения необходимой степени детализации объекта исследования.

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

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

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

Рис. 2. Схема выделения призматических элементов: а - левосторонний элемент; б - правосторонний элемент; в, г - верхние половины призм; д, е - нижние половины призм

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

Локальная нумерация выполняется последовательно против часовой стрелки, начиная с некоторого узла элемента, который выбирается произвольно. Критерием соблюдения выбранного обхода является положительный знак объема тетраэдра.

Для глобальной нумерации узлов возможны различные варианты обхода элементов, но для эффективного решения необходим поиск оптимальной нумерации, которая позволит получить матрицу жесткости элемента в виде ленточной матрицы с минимальной шириной полосы ее ленты. В данной работе предлагается обход элементов в последовательности ОХ

^ оу ^ ог.

Для моделирования угла падения пласта определяется положение плоскости, проходящей через три точки с заданными координатами М0 (х0, у0, ?0), Мх (хх, ух, z1), М2 (х2, у2, г2) уравнением вида [2]

А X + В У + С г + Б = 0. (2)

Если три точки не лежат на одной прямой, то проходящая через них плоскость представляется уравнением

х - х0 у - у0 г - г0

х - х1 у - у1 г - г1 = 0. (3)

Рис. 3. Схема разбиения призматических элементов на тетраэдры: а, б - левосторонний элемент; в, г - правосторонний элемент

х - Х2 У - У2 z - Z2

Вычислением определителя находятся значения А, О, С и О.

Для произвольной точки М(х, у, z), используя уравнение плоскости (2), координата zi с учетом угла падения пласта вычисляется по формуле

^ = - (О + А х,■ + В у) / С. (4)

Дискретизация области исследования в декартовой системе координат применяется при построении модели блочного разрушения горных пород по природным контактам с использованием метода вложенных областей, который позволяет уменьшить размеры объекта исследований без потери точности вычислений.

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

Базовый вариант объекта исследований. В качестве базового варианта объекта исследований принят слоистый массив горных пород размером 100x100x500 м, включающий выработку прямоугольной формы размером 40x40x2 м. Исследуемая область дискретизируется на призматические конечные элементы размером 5x5 м. Вертикальный размер элемента совпадает с мощностью соответствующего породного слоя.

На рис. 4 представлены параметры НДС пород непосредственной кровли. Компоненты напряжений в тетраэдрах усреднены и отнесены к центру призматического элемента. В силу симметричности вертикальных напряжений приводится 1/4 часть изображения изолиний полных и дополнительных вертикальных напряже-

Рис. 4. Напряжения в породах непосредственной кровли: а - полные вертикальные напряжения; б - дополнительные вертикальные напряжения; в - полные горизонтальные напряжения по оси ОХ

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

Разработанный алгоритм пространственной дискретизации модели массива горных пород в декартовой системе коор-

1. Фадеев А.Б. Метод конечных элементов в геомеханике / А.Б. Фадеев. - М.: Недра, 1987. -221 с.

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

Алгоритм пространственной дискретизации модели массива горных пород используется автором для моделирования геомеханических процессов блочного обрушения горных пород, подрабатываемых высокопроизводительными очистными забоями на угольных шахтах, при сложной форме выработанного пространства.

-------------- СПИСОК ЛИТЕРАТУРЫ

2. Корн Г Справочник по математике для научных работников и инженеров / Г. Корн, Т. Корн. - м.: Наука, 1973. - 832 с.

— Коротко об авторах--------------------------------------------------------------

Павлова Лариса Дмитриевна - кандидат технических наук, доцент, Сибирский государственный индустриальный университет, г. Новокузнецк.

УДК 622.232.72:622.26 А.Г. Скуров

ТЕХНОЛОГИЧЕСКИЕ ТРЕБОВАНИЯ К ПРОХОДЧЕСКИМ КОМПЛЕКСАМ НОВОГО УРОВНЯ

Семинар №10

■ Я рименение современных механи-Л л зированных комплексов для выемки пологих угольных пластов обеспечивает высокий уровень добычи и, как следствие, высокие скорости подвигания очистных забоев. Для своевременного воспроизводства очистного фронта требуется значительное увели-

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

Похожие темы научных работ по наукам о Земле и смежным экологическим наукам , автор научной работы — Павлова Лариса Дмитриевна