1. Выделение конечных элементов (разбиение области на конечные элементы); Разбиение области на элементы обычно начинают от её границы, с целью наиболее точной аппроксимации формы границы. Затем производится разбиение внутренних областей. Часто разбиение области на элементы производят в несколько этапов. Сначала разбивают на крупные части, границы между которыми проходят там, где изменяются свойства материалов, геометрия, приложенная нагрузка (другие физические величины). Затем каждая подобласть разбивается на элементы. Стараются избегать резкого изменения размеров конечных элементов на границах подобластей. После разбиения области на конечные элементы осуществляют нумерацию узлов, причем порядок нумерации имеет существенное значение, так как влияет на эффективность последующих вычислений. Это связано со следующим: матрица коэффициентов системы линейных алгебраических уравнений, к которым приводит метод конечных элементов, является сильно разреженной матрицей ленточной структуры. Ненулевые элементы такой матрицы располагаются параллельно главной диагонали. Обозначим через L число, представляющее наибольшую разность между номерами ненулевых элементов в строке. Число L называется шириной полосы. Чем меньше ширина полосы, тем меньший объем памяти требуется для хранения матрицы при реализации метода конечных элементов в САПР, и тем меньше затраты машинного времени на решение результирующей системы уравнений. Ширина полосы зависит от числа степеней свободы узлов и способа нумерации последних. Если максимальную разность между номерами узлов для конечных узлов обозначить через N, а число степеней свободы – M, то L=(N+1)*M. Информация о способе разбиения на конечные элементы и нумерация узлов является исходной для всех последующих этапов алгоритма метода конечных элементов. При этом требуется указывать не только номер, но и координаты каждого узла, его принадлежность к определенным конечным элементам, информацию о соединении элементов между собой, значения физических параметров объекта в пределах каждого элемента.
2. Определение аппроксимирующей функции для каждого элемента. На этом этапе искомая непрерывная функция аппроксимируется кусочно-непрерывной, определенной на множестве конечных элементов. Эту процедуру можно выполнить один раз для типичного элемента области и затем полученную функцию использовать для остальных элементов области того же вида. В качестве аппроксимирующей функции элементов чаще всего используют полиномы, которые подбираются так, чтобы обеспечить непрерывность искомой функции в узлах и на границах элементов.
3. Объединение конечных элементов в ансамбль. На этом этапе уравнения, относящиеся к отдельным элементам, объединяются в ансамбль, то есть в систему алгебраических уравнений. Полученная система является моделью искомой непрерывной функции. Мы получаем матрицу жесткости.
4. Решение полученной системы алгебраических уравнений. Реальная конструкция аппроксимируется многими сотнями конечных элементов, возникают системы уравнений со многими сотнями и тысячами неизвестных. Решение таких систем уравнений – основная проблема реализации метода конечных элементов. Методы решения зависят от размера разрешающей системы уравнений. В связи с большой размерностью и сильной разреженностью матрицы коэффициентов системы для реализации метода конечных элементов в САПР разработаны специальные способы хранения матрицы жесткости, позволяющее уменьшить необходимый для этого объем оперативной памяти.
Матрицы жесткости используются в каждом методе прочностного расчета, используя конечную элементную сетку. Название матрицы жесткости пришло из строительной механики, где МКЭ начал использоваться раньше, чем в других областях техники.
Для решения систем уравнений применяются методы двух групп: прямые методы (метод Гаусса), косвенные методы, когда решение определяется на основе последовательной аппроксимации (метод Гаусса-Зейделя).
