Адаптация и развитие метода конечных элементов для расчета параметров напряженно-деформированного состояния углепородного массива

Тип работы:
Реферат
Предмет:
Общие и комплексные проблемы естественных и точных наук


Узнать стоимость

Детальная информация о работе

Выдержка из работы

УДК 539. 3
Ю.А. Степанов
АДАПТАЦИЯ И РАЗВИТИЕ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ ДЛЯ РАСЧЕТА ПАРАМЕТРОВ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ УГЛЕПОРОДНОГО МАССИВА
В связи тем, что задачи горной геомеханики и инженерной геологии не могут быть точно сформулированы, то основной целью нахождения напряжений состоит в получении физической картины процесса взаимодействия системы & quot-механизированная крепь — углепородный массив& quot-, а смысл расчетов заключается скорее в качественном моделировании, нежели в количественном анализе. В соответствии с методологией имитационного моделирования одна и та же задача решается несколько раз не для улучшения точности, а для того, чтобы увидеть, как изменяется решение при возможном вероятностном изменении одного или нескольких исходных параметров. В таких условиях требуется гибкий вычислительный инструмент, который позволял бы легко модифицировать задачу, подлежащую решению, а также обеспечивал бы достаточно быстрое нахождение ответов.
Большинство задач геомеханики, решаемых с помощью метода конечных элементов (МКЭ), связаны с анизотропными средами. Так при расчете нагрузок, действующих на механизированную крепь, имеем сложную систему, состоящую из анизотропного углепородного массива, выработанного пространства и секции механизированной крепи.
Принципиальным отличием этой системы является изменение ее геометрических и механических параметров во времени и в пространстве. В этих условиях применение современных численных методов моделирования изменения параметров окружающей среды и изменения положения механизированной крепи требует постоянного изменения формы, размеров и деформационных свойств материала конечных элементов.
Многие исследователи для решения реальных задач используют известные готовые пакеты, характеризующиеся разными возможностями, разными способами постановки задач и представления результатов в зависимости от того, для каких пользователей они предназначены.
Использование любого мощного пакета имеет некоторые недостатки. В первую очередь отметим низкую эффективность его использования. Этот недостаток присущ любому универсальному инструменту, который используют для выполнения немногих из широкого круга функций.
Второе — скрытность алгоритма и невозможность изменения математических моделей, реализующих алгоритм. Существующие пакеты не всегда тривиальным образом могут быть инсталлированы на конкретном ПК. Для использования паке-
та необходимо изучать специальный макроязык, с помощью которого описывается задача.
Третье — желательно иметь полный контроль над проектом и зависеть только от собственных решений.
В качестве альтернативного варианта в ВУЗах России для решения задач горного профиля студентами и аспирантами используются узконаправленные специализированные пакеты, изготовленные собственными силами. Так, например, в Тул-ГУ разработана методика автоматизированного расчета и конструирования крепей базирующаяся на стержневой математической модели [1]. Она позволяет путем численного моделирования с использованием ПЭВМ определять напряженно-деформированное состояние крепей при различных видах воздействий с учетом деформации звеньев, податливости связей, сил инерции и сопротивления среды.
Методика, предлагаемая нами, отличается от имеющихся, прежде всего тем, что нагружение крепи производится не с помощью набора типовых статических или динамических воздействий, а моделированием реальной нагрузки, возникающей в процессе выемки угля механизированным комплексом [2].
Развитием метода конечных элементов для разработки методики прогноза геомеханического взаимодействия механизированных крепей с углепородным массивом в очистных забоях угольных шахт, является разработка модели накопления деструктивных изменений в массиве горных пород.
Для этого, сначала решается упругая задача статического состояния системы & quot-углепородный массив — механизированная крепь& quot- и определяются напряжения в каждом конечном элементе пород и элементов крепи. По паспорту прочности проверяется состояние пород и металла в этих элементах, и запоминаются новые деформационные и механические характеристики.
На следующем этапе производится разгрузка породы от давления гидростоек, т. е. проводится второй этап счета на ЭВМ методом конечных элементов незакрепленной кровли с учетом предыдущего этапа и состояния пород. После этого этапа также определяется состояние боковых пород и металла в конечных элементах и формируется зона разрушения пород кровли и элементов крепи.
На третьем этапе проводится нагружение породы после передвижки секции крепи. Напряжен-
6////^ 24 ^30^ 36 42 48 54 60
у 11^^/ 17 23 29 35 41 47 53 59
4 10^У 16 22 28^/ 34^/ 40^/ 52^/// 58^у
23 У / Л/ /& quot- Л/ ^9^ 45^^ 51^у 57 у/
2 8 У 20/ 32^/ тттттт-гу 38 У 44 У -ГТТТТТТТ-^ 50 у 56 у
66
65
64
63
62
13
19 25 31
37
43
49
55
61
Рис. 1. Схема нумерации узлов конечных элементов
но — деформированное состояние конечных элементов, где порода находится в дискретном состоянии, рассматривается с позиции текущей среды. Если горное давление со стороны кровли превышает распор крепи или впереди крепи преобладает зона сыпучего материала, что определяется путем сравнения по паспорту прочности в каждом конечном элементе, то впереди забоя формируется зона отжима и купол. По результатам третьего этапного расчета проводится оценка режимов нагружения, т. е. величины распора, разгрузки и определяется оптимальный режим передвижки крепи с начальным распором.
Как известно, идея метода конечных элементов основана на аппроксимации непрерывной функции (давления, перемещения и т. д.) дискретной моделью, которая строится на множестве непрерывно-кусочных функций, определенных на конечном числе подобластей, называемых элементами [3]. В качестве функции элемента чаще применяется полином, порядок которого зависит от числа данных о непрерывной функции используемых в каждом узле элемента. Симплекс-элементам соответствуют полиномы, содержащие константу и линейные члены.
Число элементов в таком полиноме на единицу больше размерности координатного пространства. Полином
(р = а1 +а2 х + а3 у (1)
представляет собой симплексную функцию для двумерного треугольного элемента. Этот полином линеен по Х и У и содержит три коэффициента, потому что треугольник содержит три узла. Следовательно, двумерный симплекс-элемент — это
треугольник с прямолинейными сторонами и тремя узлами.
Для того чтобы найти перемещения в узлах конечных элементов, необходимо провести логическую нумерацию вершин треугольника. Правильная нумерация узлов сокращает объем машинной памяти. На рис. 1 показан один из способов нумерации узлов в сечении углепородного массива для исследования процесса взаимодействия механизированной крепи с боковыми породами.
Создание триангуляционной среды может быть выполнено как в автоматическом, так и в интерактивном режиме.
Автором и группой единомышленников разработана программа, позволяющая существенно упростить процесс корректировки формы, размеров и свойств пород любого конечного элемента исследуемой области. Она представляет собой специализированный графический редактор, который используется после машинной процедуры разбиения исследуемой области на конечные элементы, с помощью которого осуществляется корректировка положения вершин конечных элементов относительно выработанного пространства и элементов секции механизированной крепи.
Векторная величина & quot-перемещение"- имеет как величину, так и направление. Поэтому в каждом узле необходимо определять более одной неизвестной (две степени свободы). В этом случае векторную величину представляют ее компонентами, которые рассматриваются как неизвестные скалярные величины. Обозначение компонент вектора изображено на рис. 2.
7
1
Узловые значения скалярной величины ф (рисунок 3) обозначаются через ФІ, Ф], Фк, а координатные пары трех узлов — через (Хі, Уі), (X] ,
Рис. 3. Двумерный симплекс-элемент
У]), (Хк, Ук). Интерполяционный полином элемента имеет вид ф = а1 + а2 х + а3 у.
Решая систему линейных уравнений, получаем интерполяционные функции в узлах элементов. В общей форме интерполяционный полином имеет вид:
Ф,
Ф ,¦
ф
,(е)
= [ N ]*{Ф} = [ ^е), N (-Ы (ке) … N1*
Фк
Фг
(2)
Общепринятая формулировка метода конечных элементов предполагает отыскание поля перемещений путем минимизации потенциальной энергии системы. В результате, уравнения, определяющие элементы, сводятся к системе алгебраических уравнений равновесия, которые можно решить относительно узловых перемещений методом Гаусса.
После того, как будут найдены узловые значения вектора перемещений, производится расчет напряжений в центрах конечных элементов.
Описанная процедура & quot-классического"- метода
конечных элементов была принята в качестве основы алгоритма расчета напряженно-деформированного состояния углепородного массива в окрестности очистного комплексномеханизированного забоя.
Совершенно очевидно, что любая модель может только приближенно отображать свойства моделируемого объекта и процессы, протекающие в нем. Для обеспечения требуемой адекватности модели и повышения точности результатов моделирования модель (алгоритм ее работы) необходимо настраивать.
В данной работе для настройки алгоритма использовали результаты натурных измерений конвергенции основания и перекрытия секций механизированных крепей. Для повышения достоверности результаты натурных измерений подвергали предварительной обработке: осуществляли визуальный анализ данных, отфильтровывали выбросы, производили расчет необходимых статистик. Технология настройки алгоритма состояла в использовании итеративной процедуры корректировки параметров алгоритма на основе отклонений расчетных смещений горных пород в окрестности очистного забоя от предварительно обработанных фактических.
Алгоритм расчета параметров НДС углепородного массива в окрестности очистного забоя может быть реализован на основе классического МКЭ путем добавления адаптирующей части, представляющей собой набор моделей для расчета скорректированных значений модуля упругости и коэффициента Пуассона, обеспечивающих соответствие расчетных данных реальным по формулам:
Е = а о + а I-к ^_-/кр (3)
М= в0 -в1 (4)
где, а о, а I, в о, в I — настраиваемые коэффициенты-
к 8 — коэффициент структурного ослабления-
/ кр — коэффициент крепости породы.
Для настройки моделей использовались результаты замеров смещений пород кровли и почвы в шахтных условиях [4], а также размеры зон разрушенных пород кровли и фактических вывалов при движении очистного забоя в угольных шахтах концерна «Кузнецкуголь».
Упруго-пластичная модель разрушения горных пород реализована путем повторения этапа решения системы линейных алгебраических уравнений для изменененных прочностных свойств конечных элементов, находящихся в углепородном массиве. Изменение прочностных свойств материала конечных элементов производили на основе определения значения коэффициента остаточной прочности (Эо1), значение которого могло изменяться в пределах от 0 до 1. Величину коэффициента остаточной прочности находили по
формулам приведенными в [5]. Корректировка прочностных свойств пород производилась с использованием модели вида:
Енов = Е стар — ехр «'- П (ВЫ), = 1* 5 (5)
где, а , — коэффициенты-
Бо1 — коэффициент остаточной прочности-
, — номер участка (интервала), в котором находится значение Бо1.
Кроме того, адаптация МКЭ к условиям конкретной задачи предполагает:
— выбор формы и размеров конечных элементов, адекватно отражающих конфигурацию исследуемой системы-
— создание алгоритма перестройки триангуляционной сети в соответствии с формой и размерами секции механизированной крепи-
— трансформирование сети конечных элементов в соответствии с периодическими подвижками секций механизированной крепи-
— трансформирование сети конечных элементов в соответствии с периодичностью зависания и обрушения пород кровли-
— изменение сети конечных элементов в соответствии с формой и размерами зон отжимов поверхности забоя и вывалов пород кровли впереди и над секцией механизированной крепи-
— многократное повторение этапа решения системы линейных алгебраических уравнений, для реализации упруго-пластичной модели и изменения прочностных свойств конечных элементов, перешедших в запредельное состояние, для кор-
ректировки положения узлов конечных элементов в соответствии с очертаниями проекций элементов механизированной крепи.
В соответствии с изложенным, можно сделать вывод о том, что большее количество мероприятий, связанных с адаптацией метода конечных элементов, относится к препроцессингу — комплексу мероприятий, проводимых перед вычислением смещений в узлах триангуляционной сети. Это потребовало разработку различного рода сеточных генераторов для автоматического формирования триангуляционной сети и интерактивных графических редакторов [6,7]. При вычислении параметров геомеханического взаимодействия секции механизированной крепи с углепородным массивом использовалась упруго-пластичная модель напряженно-деформированного состояния массива горных пород. Поскольку боковые породы очистного забоя испытывают знакопеременную нагрузку в результате циклического движения секции механизированной крепи, адаптация модуля решения была направлена, на использования алгоритма накопления деструктивных изменений в породах, напряжения в которых превысили порог их прочности.
При разработке программных средств, применялись, преимущественным образом, динамические структуры данных, использование которых, снимает с пользователя дополнительные проблемы, связанные с разбиением и подкачкой данных. Эти функции выполняет сама операционная система.
СПИСОК ЛИТЕРАТУРЫ
1. Применение методики стержневой аппроксимации для расчета крепей различных конструкций. / В. М. Еганов, А. Е. Коряков, С. П. Туляков: Тул. Гос. ун-т. — Тула. 1998. — 12 с. — Рус. — Деп. в ВИНИТИ 27. 11. 98. № 3494-В98.
2. Горский В. Г., Адлер Ю. П. Планирование промышленных экспериментов. — М. & quot-Металлургия"-, 1974. — 264 с.
3. Сегерлинд Л. Применение метода конечных элементов. /Пер. с англ. к.ф.н. Шестакова А. А. -М.: Издательство & quot-МИР"-, 1979. — 392 с.
4. Штумпф Г. Г., Рыжов Ю. А. Шаламанов В.А., Петров А. И. Физико-технические свойствагор-ных пород и углей Кузнецкого бассейна: Справочник. — М.: Недра, 1994 — 447 с.
5. Булычев Н. С. Механика подземных сооружений в примерах и задачах. / Учебное пособие для вузов. — М.: Недра, 1989. — 270 с.
6. Степанов Ю. А. Сеточный генератор. Материалы научно- практической конференции Кузбасса, / Под общ. ред. проф. К. Еафанасьева/. Кемерово: Изд-во & quot-Полиграф"-, часть 2, январь 2001 г. — 220 с. сервер: http: //conference/kemsu. ru/infokuz.
7. Степанов Ю. А., Степанов А. В. Корректировка триангуляционной сети метода конечных элементов с помощью графического редактора ЭВМ. Нетрадиционные и интенсивные технологии разработки месторождений полезных ископаемых: Труды IV Международной конференции / СибГИУ. — Новокузнецк, 1999. — 273 с.
? Автор статьи:
Степанов Юрий Александрович, канд. техн. наук, доцент каф. информационных систем и управления НФИ КемГУ (г. Новокузнецк)
Email: Dambo290@yandex. ru

ПоказатьСвернуть
Заполнить форму текущей работой