Высокоточный тетраэдральный конечный элемент

Объект исследования: тетраэдральный конечный элемент с дополнительными степенями свободы в виде частных производных.

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

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

В данной работе рассматривается трехмерный тетраэдральный конечный элемент. Тетраэдральные конечные элементы широко используются для расчета пространственных тел, имеющих трехосное напряженное состояние.

Недостатки данного конечного элемента

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

2. В случае с распределенной нагрузкой необходимо точное приведение нагрузки к узлам. Это связано с интегрированием, что вызывает значительные сложности.

Достоинства данного конечного элемента

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

В качестве степеней свободы КЭ имеет не только трасляционные перемещения (u,v,w) , но и значения их производных ( Полный кубический полином содержит 20 членов. Каждой из вершин соответсвует по три компоненты трансляционного перемещения и для каждой компоненты перемещения – по три частные производные, всего 12 величин на узел. Это приводит к появлению 48 степеней свободы. Дополнительные 12 степеней свободы можно задать в виде трансляционных перемещений центров каждой из четырех граней. Всего конечный элемент имеет 60 степеней свободы.

Основные формулы.

Полный кубический интерполяционный полином имеет вид (1):

u(x,y,z)=a1+a2x+a3y+a4z+a5x2+a6y2+a7z2+a8xy+a9yz+a10zx+a11x3+a12y3+a13z3+

+a14xy2+a15xz2+a16yx2+a17yz2+a18zx2+a19zy2+a20xyz (1)

Или в тетраэдральных координатах (2):

u=L1a1+L2a2+L3a3+L4a4+L1L2a5+L3L4a6+L4L1a7+L3L2a8+L2L4a9+L3L2a10+(L12L2-L22L1)a11+(L32L4-L42L3)a12+(L42L1-L12L4)a13+(L32L2-L22L3)a14+(L22L4-L42L2)a15+

+(L32L1-L12L3)a16 (2)

где L 1, L2, L3, L4 - тетраэдральные координаты, или объемные L-координаты, находятся по формуле(3):

где – объем, заключенный между линиями, соединяющими точку с вершинами тетраэдра, противолежащими вершине i.

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

Матрица жесткости конечного элемента задается выражением (4):

высокоточный тетраэдральный конечный элемент - student2.ru (4)

где матрица и содержат производные от функций форм.

Матрица имеет размер [6][60] .

Мы использовали формулы численного интегрирования для тетраэдров с 5-ю точками интегрирования. В дальнейшем, для повышения точности решения, необходимо взять большее количество Гауссовых точек интегрирования.

Таким образом, получаем:

высокоточный тетраэдральный конечный элемент - student2.ru (6)

где - весовые коэффициенты (для каждой точки интегрирования), - объем тетраэдра.

В дальнейшем планируется использование данного конечного элемента в программном комплексе DSMFem.

Материал поступил в редколлегию 04.05.2017

УДК 629.1.015

А.С. Кондрашова

Научный руководитель: профессор кафедры «Механика и динамика и прочность машин», д.т.н. В.И. Сакало

[email protected]

Наши рекомендации