Примечание 1

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

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

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

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

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

Основой математической модели -го КЭ является квадратная матрица жесткости конечного элемента, размер которой равен числу неопределенных коэффициентов в части вектора , относящейся к этому КЭ. В свою очередь, размер этой части вектора есть произведение размерности пространства и числа узлов, выделенных в модели -го КЭ. то следует из того, что каждый неопределенный коэффициент в в задаче анализа механической прочности есть значение перемещения в направлении оси в -м узле. При получении модели КЭ используется матрица координатных функций конечного элемента, у которой число строк равно , а число столбцов × . Например, размер для КЭ в форме параллелепипеда с восемью узлами равен 24×24, так как размеры матриц , и суть 6×3, 3×24 и 6×6. При расчете сначала выбирают вид зависимости перемещений от пространственных координат. Затем преобразуют эти зависимости в выражения и с помощью интегрирования по пространственным координатам определяют матрицу жесткости КЭ.

При наличии библиотеки КЭ применение МКЭ сводится к следующим операциям:

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

2. Выбор библиотечной модели КЭ, задание внешних нагрузок и значений геометрических и физических параметров, формулировка граничных условий. Следующие операции выполняются программой моделирования.

3. Реализация в модели сетки конечных элементов. Тем самым становятся известными координаты узловых точек в модели.

4. Приведение имеющихся объемных сил и поверхностных нагрузок к узловым точкам модели.

5. Объединение моделей КЭ в общую конечно-элементную модель детали (6), в которой матрица жесткости имеет порядок, равный , где — общее число узлов. При объединении элементы матрицы образуются суммированием тех элементов матриц жесткости отдельных КЭ, которые относятся к одному и тому же узлу и направлению перемещения. Если некоторый узел закреплен (его перемещение равно нулю), то соответствующие этому узлу строки в и и столбцы в вычеркиваются.

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

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

К числу известных программ анализа по МКЭ относятся ANSYS, NASTRAN, PATRAN и др.