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

Материалы и методы. Для построения матриц жесткости объемных КЭ применен математический аппарат метода двойной аппроксимации, суть которого состоит в раздельном представлении функций распределения перемещений и деформаций внутри элемента. Хранение и обработка результирующей системы уравнений реализованы в алгоритмических терминах разреженных матриц. Разработка программного обеспечения и проведение вычислительных экспериментов осуществлены с использованием 64-х разрядной вычислительной платформы Microsoft Visual Studio 2013 и компилятора Intel® Parallel Studio XE 2019 со встроенным текстовым редактором Intel® Visual Fortran Composer XE 2019. Визуализация результатов расчетов выполнена с помощью дескрипторной графики пакета компьютерной математики Matlab. В качестве тестового образца использован объемный восьмиузловой КЭ SOLID185 программного комплекса ANSYS Mechanical.

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

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