Введение
Метод конечных элементов обычно используется в тех случаях, когда задачи механики невозможно решить аналитически. При этом часто учитываются нелинейные эффекты, такие как выход из работы при сжатии (геометрическая нелинейность), пластификация (нелинейность материала) и контактные или кинематические степени свободы. Такие эффекты, особенно в случае применения нелинейной модели материала, можно учесть с помощью итерационного метода расчета.
Формулировка МКЭ
Основные этапы формулировки МКЭ (дальнейшая информация находится в [1]):
- Слабая форма равновесия
δu Виртуальное (тестовое) перемещение t0 Коэффициент начальной нагрузки σdV Внутренние силы ρbdV Силы в теле B Интегрированная площадь - Преобразование в нотацию Фойгта с тензором 4-го порядка
C Матрица жесткости δε Разные варианты деформированного состояния B Интегрированная площадь ε Деформация U Удлинение
На данной записи далее основывается приближенное решение МКЭ для нелинейных материалов, - причем поле смещения умножается на функции формы.
u(x,t) Перемещение во времени (приращение нагрузки) H Функция формы û Перемещение узла - Использование производной смещения в слабой форме. Расчет перемещения узлов выполняется с помощью численного интегрирования, а постпроцессинге затем по закону материала определяется напряжения и деформации.
Порядок итерации Ньютона-Рафсона
Из-за нелинейных свойств материала матрица материала C в вышеприведенном равнении 2 меняется с каждым шагом деформации. Поэтому стандартным методом расчета для решения данной задачи является так называемая итерация Ньютона-Рафсона. При ее использовании функция линеаризируется уже в исходной точке. При итерации затем всегда используется матрица жесткости C из предыдущего шага. Только в линеаризованном шаге итерации проходит касательная через нулевую точку функции.
К схеме порядка итерации, которая показана на вышеприведенном рисунке, относятся также следующие уравнения:
- Разделение нагружения на отдельные шаги нагрузки.
fext Внешняя нагрузка t Шаг времени/нагрузки Δf Новое приращение нагрузки - Шаг предсказывания
K Матрица жесткости предыдущего шага времени t+Δtfext Внешние силы, увеличенные последующим шагом нагрузки 0tfint
Внутренние силы предыдущего шага времени ϕ Деформация - В точке 3 - Итерации логической схемы рассчитывается общая деформация с вычетом пластической деформации (шаг коррекции).
Главной целью итерационного расчета всегда является достижение условия, при котором сумма нагрузок равна нулю, но в числовом виде это невозможно. Поэтому задается предел прерывания ε, при достижениее которого будет расчет прерван как достаточно точный.
| R | Предел отказа |
| fext | Внешние силы |
| fint | Внутренние силы |
| ε | Предел отказа эпсилон |
| t | Шаг времени |
В программе предел прерывания настраивается прямо в параметрах расчета.
На следующем рисунке затем показан процесс итерации Ньютона-Рафсона. В первой итерации
предел прерывания R или ε не достигнут. Также во второй итерации (обозначено красным) предел допуска не достигается. Только в третьей итерации становится расстояние касательной жесткости настолько малым, что расчет достигает сходимости.
Как уже упоминалось раньше, в ходе итерации деформации постоянно суммируются.
Заключение
Итерация Ньютона-Рафсона обладает консистентностью, и соответственно, сходимостью второго порядка, причем количество «правильных» точек итерации с каждым шагом удваивается. В результате так итерация Ньютона-Рафсона сходится квадратично, и в случае достижения сходимости ее точность с каждой итерацией увеличивается. Однако, если данным методом нельзя достичь сходимости, то ошибка повторяется бесконечно и вычисление прерывается.
Ошибку расчета чаще всего вызывает либо слишком сильный наклон кривой зависимости деформаций от нагрузки, либо слишком малый наклон кривой в области пластических деформаций. С другой стороны, если кривая зависимости деформации от нагрузки на вышеприведенном рисунке показывает слишком сильный излом на втором шаге итерации, то касательная материала и, следовательно, матрица жесткости не будут правильно отображать наклон упругой области. Потому что в таком случае был бы подъем нулевой точки для пластической области неправильным. Это в принципе и есть одна из причин, по которой увеличение количества шагов нагрузки сопровождается улучшением сходимости.
Dlubal_KohlA_]_LI.jpg?mw=350&hash=ee8d38f1c4853d80307fa156c159b5e78a3fdca9)
.jpg?mw=350&hash=8f312d6c75a747d88bf9d0f5b1038595900b96c1)
.jpg?mw=350&hash=196d91410dd36f58fcc8b0194f4577bb70cf53c3)
.png?mw=600&hash=49b6a289915d28aa461360f7308b092631b1446e)
