Численные методы расчета - Алгоритмы компьютерного моделирования

Решить задачу для математической модели - значит указать алгоритм для получения требуемого результата из исходных данных.

Алгоритмы решения условно делятся на:

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

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

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

Решение задач строительной механики связано с большими математическими трудностями, которые преодолеваются с помощью численных методов, позволяющих с применением ЭВМ получать приближенные, но удовлетворяющие практическим целям решения [2].

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

В зависимости от способов выполнения дискретизации и алгебраизации выделяют различные методы.

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

Алгоритм МКР состоит из трех этапов:

    1. Построение сетки в заданной области. В узлах сетки определяются приближенные значения функции (узловые значения). Набор узловых значений - сеточная функция. 2. Частные производные заменяются разностными выражениями. При этом непрерывная функция аппроксимируется сеточной функцией. В результате получают систему алгебраических уравнений. 3. Решение полученной системы алгебраических уравнений.

Еще одним численным методом является метод граничных элементов (МГЭ). Он основывается на рассматривании системы уравнений, включающей только значения переменных на границах области. Схема дискретизации требует разбиения лишь поверхности. Граница области делится на ряд элементов и считается, что нужно найти приближенное решение, которое аппроксимирует исходную краевую задачу. Эти элементы называются граничными. Дискретизация только границы ведет к меньшей системе уравнений задачи, чем дискретизация всего тела. МГЭ уменьшает размерность исходной задачи на единицу.

При проектировании различных технических объектов широко используется метод конечных элементов (МКЭ). Возникновение метода конечных элементов связано с решением задач космических исследований в 1950-х годах [3]. В настоящее время область применения метода конечных элементов очень обширна и схватывает все физические задачи, которые могут быть описаны дифференциальными уравнениями. Наиболее важными преимуществами метода конечных элементов являются следующие:

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

Решение задач по МКЭ содержит следующие этапы:

    1. Разбиение заданной области на конечные элементы. Нумерация узлов и элементов. 2. Построение матриц жесткости конечных элементов. 3. Сведение нагрузок и воздействий, приложенных к конечным элементам, к узловым силам. 4. Формирование общей системы уравнений; учет в ней граничных условий. Решение полученной системы уравнений. 5. Определение напряжений и деформаций в конечных элементах.

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

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

Похожие статьи




Численные методы расчета - Алгоритмы компьютерного моделирования

Предыдущая | Следующая