Конечные элементы для деформационной задачи
При моделировании деформации деталей применяют конечные элементы с тремя и более узлами, показанные на рис. 13.3. В плоском трехузловом конечном элементе перемещения и внутренних
точек выражаются через перемещения Uj узлов согласно линейной интерполяционной формуле
(13.15) |
где |
N;=4-
1 A
- функции формы, зависящие от расположения внутренней точки М (рис. 13.12); А[ - площадь одного из треугольников, на которые делят треугольный конечный элемент площадью А линии, соеди - I няющие точку М с узлами элемента. На
Рис. 13.12. Схема определения значений функции формы для точки в трехузловом конечном элементе |
рис. 13.12 показан конечный элемент и заштрихован треугольник, площадь которого А подставляется в (13.15) для определения значения функции формы N в точке М. Из формулы (13.15) и рис. 13.12 следуют некоторые свойства функций формы. Когда точка М совпадает с узлом
/, значение функции N равно единице, а значения двух других функций формы равны нулю. При этом, независимо от перемещений узлов 2 и 3, и = U. Сумма трех функций формы для любой точки М равна единице. Поэтому если перемещения всех трех узлов одинаковы, то и = U = U2 = t/3.
С помощью функций формы легко определить компоненты деформации в точке Мпо формулам (11.6), например:
8N2 dN3 х2~^~ + ихЪ' |
дх |
дх |
дх |
дих ехх ~ Л — их1 _ + U |
дх |
Перемещения узлов не зависят от координат, поэтому дифференцировать необходимо только функции формы. Если эти функции зависят от координат линейно, то их производные являются константами и не зависят от координат точки М. Поэтому в рассмотренном трехузловом конечном элементе компоненты деформации во всех точках одинаковы.
Применяя закон Гука, можно по деформациям рассчитать напряжения, используя формулы (11.18), (11.19) и (11.24). Напряжения также одинаковы во всех точках конечного элемента.
Упругая энергия в единице объема материала пропорциональна сумме произведений компонент деформации и напряжения. Чтобы найти упругую энергию для всей модели, необходимо сложить энергии, рассчитанные для всех конечных элементов. Она оказывается квадратичной функцией перемещений всех узлов модели.