Теория катастроф в задачах анализа состояния горнотехнических зданий и сооружений Текст научной статьи по специальности «Строительство и архитектура»

80

Д. И. Назаров

УДК 681.3

Д. И. Назаров ТЕОРИЯ КАТАСТРОФ В ЗАДАЧАХ АНАЛИЗА СОСТОЯНИЯ ГОРНОТЕХНИЧЕСКИХ ЗДАНИЙ И СООРУЖЕНИЙ

Теория катастроф [1], появившаяся во второй половине прошлого века, получила за последние пол века значительное развитие в математике, однако в прикладных науках не нашла должного отражения, только частична отражена в термодинамике [2], квантовой механике, климатологии и экономике. В механических моделях наиболее упоминаемы ферма Мизеса и машина катастроф Зимана.

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

Значительный научный прогресс метода конечных элементов в нелинейном анализе [3] с эффективным использованием в современном программном обеспечении (ANSYS, Algor, MSC/NASTRAN, UAI/NASTRAN, MARC и т.д.) создал мощный инженерно-технический аппарат с достаточно высокой достоверностью результатов вычислений. Однако, нелинейный метод конечных элементов, имея высокую точность для монотонных (квазилинейных) систем, не обеспечен математическим аппаратом для систем, имеющих, при аналитическом расчете или экспериментальном исследовании, бифуркации.

Рис. 1. Пример конструкции с двумя узлами бифуркации

Реализация в методе конечных элементах метода продолжения решения по параметру (в зарубежной литературе - агс-1еп§;Ш), позволяет находить корректное решение в системах имеющих только один узел бифуркации. Реальные инженерные конструкции, в силу изначальной конст-

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

Аналитическое решение представленной конструкции:

У(х) =

(1)

где X - горизонтальная проекция правого элемента, а у(х) - вертикальная проекция левого элемента,

Ry( x) = y(x) ■ EA

пр

x

--------------x

X + X - x

лв пр

f \ 1 1

Vх 2 + Yl 1пр

(2)

где Яу(х) - вертикальная реакция левого элемента (вертикальная проекция усилия в правом элементе).

Несложно убедиться, что все решение фактически зависит от одного параметра - горизонтальной проекции правого элемента (X).

Откладывая по оси абцисс у(х), а по оси ординат Яу(х), строим график зависимости Яу от

Ау=(у-¥лв) (рис. 2.), для удобства анализа совместив с внешней силой (Р).

На графике росту приложенной силы (Р) соответствует кривая Л-Б-Е-Ж-И-І (черная пунктирная линия), напомним, что нами рассматривается статическая внешняя сила, а она (по определению) растет от нуля до указанной расчетной величины.

Таким образом, сила, достигнув значения Рь, не уменьшается (если в бак накапало 95 литров, то меньше не станет).

Реакция левого элемента (Яу), в зависимости от перемещения узла (Ау), представлена кривой Л-Б-Є-Б-Е-Ж-О-И-І (сплошная линия).

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

X

Геотехнология

81

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

I ^ п ,_________,

^піп -^г I

Fтр=^m■g

Рис. 3. Анализ воздействий на тело: процесс деформации при P<Fтр; процесс движения при P>Fтр

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

В физическом представлении - точка бифуркации, что позволяет использовать известный математический метод [3] без учета инерционных характеристик системы, т.е. в нелинейной статической постановке задачи.

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

Для инженерно-технических задач данный подход к подобным задачам может считаться обоснованным.

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

1. АрнольдВ.И. Теория катастроф. М.: Наука, 1990. 128 с.;

2. Басин М.А. // Письма в ЖТФ. 2006. Т. 32. В. 8. С. 30-33;

3. Григолюк Э.И., Шалашилин В.И. Проблемы нелинейного деформирования: Метод продолжения решения по параметру в нелинейных задачах механики твердого деформируемого тела. М.: Наука, 1988. 232 с.;

4. Назаров Д.И. Некоторые особенности геометрически-нелинейных задач// Автоматизация и информатизация машиностроения - 2000: Сб. тр. междунар. науч.-техн. конф. - ТулГУ, Тула, 2000. С. 9699.

□ Автор статьи:

Назаров Дмитрий Иванович - доцент кафедры «Строительство подземных сооружений и шахт» КузГТУ.

E-mail: L01BDV@yandex.ru

Рис. 2. График поведения конструкции

Очевидно, что на кривой С-К-Б и Е-Ь-О (штрих-пунктирная линия), конструкция оказаться не может, хотя бы по причине невозможности (после достижения внешней силой значения Pь) перемещения узла в обратном направлении (создавая отрицательную работу).

Соответственно, раз на участке В-Е отсутствует статическое равновесие, то необходимо учитывать динамические составляющие, которые как раз в существующем методе конечных элементов отсутствуют.

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

Объясняя подобные явления проще провести аналогию с классической задачей теоретической механики - равновесие/движение тела (рис. 3.) лежащего на поверхности.

Ранее [4] указывалось, что для изучения процесса перескока необходимо использовать дина-