Рассмотрены два алгоритма конечноэлементного моделирования неоднородных свойств, в том числе при связности механических и электрических полей. В первом случае неоднородные механические и пьезоэлектрические свойства тела задаются в виде аналитических зависимостей в некоторых классах функций. Во втором предполагается, что функции, описывающие неоднородные свойства, известны в дискретном наборе точек. Неоднородные свойства учитываются при вычислении элементных матриц. Для этого во втором случае интерполяция неоднородных свойств внутри тела осуществляется двумя способами: с помощью сплайнов, реализующих модель тонкой упругой пластины, и сплайнов Ж. Менге. Предложенные алгоритмы реализованы в виде модулей конечноэлементного комплекса ACELAN. Разработанное программное обеспечение может быть использовано при решении статических и динамических задач для неоднородных упругих и электроупругих тел. Приведены соответствующие примеры.
Two finite element simulation algorithms for heterogeneous properties including the connectivity of mechanical and electrical fields are considered. In the first instance, heterogeneous mechanical and piezoelectric properties of the body are given in the form of analytical relationships in certain classes of functions. The second method assumes that the functions describing the heterogeneous properties are known in a discrete set of points. The changes caused by inhomogeneous properties are taken into account under the calculation of the element matrices. For this, in the second case, the interpolation of heterogeneous properties within the body is implemented in two ways: using thin-plate splines, and J. Meinguet´s splines. The proposed algorithms are implemented as modules of the ACE-LAN finite element complex. Examples of using the developed software for the calculation of the stress strain state of inhomogeneous elastic and electroelastic bodies are presented.
Введение. Функционально-градиентные (ФГМ) материалы широко применяются в современном приборостроении, авиастроении, медицине и т. п. Физическое моделирование таких материалов дорого и не всегда возможно. Поэтому проектирование технических устройств, в которых используются ФГМ, опирается на теоретическое исследование напряжённо-деформированного состояния. В большинстве случаев оно основывается на численном математическом моделировании — например, в рамках метода конечных элементов (МКЭ). Поэтому разработка программного обеспечения (ПО), позволяющего проводить такое моделирование, представляется весьма актуальной. В настоящей статье рассматривается создание такого ПО для конечноэлементного комплекса ACELAN. В современных исследованиях выделяются два класса функционально неоднородных материалов. В первом случае неоднородность возникает при поляризации керамики, когда электродное покрытие имеет сложную геометрию [1]. В работе [2] рассматривается неоднородность другой природы. В неоднородном материале имеется трещина, которая подвергается структурному анализу при помощи метода Галёркина — Петрова. Вблизи трещины все свойства изменяются по одному закону (предполагается увеличение в ё*'5 раз). Исключение — тепловые коэффициенты, которые просто экспоненциально убывают по мере приближения к трещине. В настоящей статье разработан модуль конечноэлементного комплекса ACELAN для моделирования функционально неоднородных материалов, в том числе пьезоэлектрических керамик с неоднородной поляризацией. При этом функция распределения неоднородных свойств материала, зависящая от положения точки в теле, предполагается известной в виде либо функциональной зависимости, либо дискретного набора значений величин соответствующих физических свойств в определённых узлах. Описана методика интерполяции этих свойств на используемую конечноэлементную сетку. Приведены примеры решения задач для тел с неоднородными свойствами.