Об одном способе аппроксимации производных в гексаэдрических 8-узловых конечных элементах

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

On the Approximation of Derivatives in Hexahedral 8-Node Finite Elements

In numerical solutions of elasticity and plasticity problems, finite elements with a reduced integration are often used, especially in solving dynamic problems. In this case, for 8-node 3D elements, one point of the numerical integration is used instead of 8. In this case, it is actually assumed that the strains and stresses are constant within an element. In this case, the traditional technique of constructing the stiffness matrix for an element of a standard shape in the form of a cube with a subsequent mapping of the actual finite elements of an arbitrary shape and size onto the standard one is not necessary. Instead, the stiffness matrix can be constructed directly for a finite element of an arbitrary shape. In this case, it is expressed through the coefficients of grid operators approximating the first partial derivatives of the displacement field in the finite element. The paper considers a new approach to approximating derivatives when constructing the stiffness matrix for a 3D 8-node finite element with one integration point. The theoretical basis for this approach is the further development of the class of rare mesh FEM schemes. The obtained formulas allow one to construct incompatible FEM schemes with improved properties. The paper discusses the problems of hourglass instability, shear and volume locking. A new effective approach to solving the hourglass problem is proposed. The possibility of applying new derivative approximation formulas to finite elements of a degenerate form with a number of nodes less than eight is also discussed. It is shown that they remain applicable in a standard way in this case too. The results of the study are confirmed by the presented numerical solution results of the model static elasticity problems.