Full Text

Введение Благодаря способности электромагнитоупругих (MEE) материалов преобразовывать один вид энергии в другой (между механической, энергией электрического и магнитного полей) они находят широкое применение в различных областях науки и техники. Они используются во множестве элементов интеллектуальных устройств, таких как преобразователи, сенсоры, актуаторы. Кроме того, они применяются при разработке акустических устройств, гидрофонов, лазеров и т.д. Повышенное внимание к электромагнитоупругим материалам объясняется тем, что из-за взаимосвязанности электрического, механического и магнитного полей они обладают улучшенными свойствами по сравнению с традиционными интеллектуальными материалами. Однако ярко выраженная связанность между полями различной физической природы значительно усложняет решение статических и динамических задач линейной электромагнитоупругости. Численное моделирование играет очень важную роль при разработке, производстве и исследовании ответственных элементов инженерных конструкций, деталей машин и оборудования, особенно при использовании композитных материалов со связанными полями. Подобные материалы, как правило, являются хрупкими и существенно анизотропными, что ведет к увеличению трудоемкости математического моделирования их поведения. Таким образом, разработка высокоточных численных подходов для решения статических и нестационарных динамических трехмерных задач линейной теории электромагнитоупругости является актуальной научной задачей. За последние два десятилетия были предложены некоторые аналитические и полуаналитические подходы и разнообразные численные методы решения статических и динамических двумерных и трехмерных задач линейной электромагнитоупругости [1-10]. Тем не менее относительно малое количество работ посвящено использованию метода граничных элементов (МГЭ) для численного моделирования динамики электромагнитоупругих тел. Ключевыми особенностями МГЭ являются: уменьшение размерности задачи (требуется дискретизация только границы рассматриваемой области); более высокая точность и простое применение к задачам с бесконечными и полубесконечными границами. МГЭ обладает некоторыми преимуществами перед широко используемыми методами конечных элементов и конечных разностей, однако при применении для решения задач анизотропной теории упругости со связанными полями возникают известные трудности при выводе и реализации фундаментальных решений. Насколько известно автору, в научной литературе отсутствуют примеры использования традиционной гранично-элементной схемы для решения нестационарных динамических задач электромагнитоупругости. В этой работе, в развитие предыдущих работ [11, 12], прямой подход трехмерного метода граничных элементов в пространстве Лапласа распространён с анизотропной упругости, вязкоупругости и электроупругости на анизотропную линейную электромагнитоупругость. Трехмерные статические фундаментальные решения (функции Грина) выражены в неявном виде как интеграл по единичной окружности [13, 14], а динамические функции Грина - как поверхностный интеграл по единичной полусфере [15]. Для пространственной дискретизации интегрального уравнения используется классический узловой метод коллокаций совместно со смешанными граничными элементами [16]. Для уменьшения времени вычислений динамические части фундаментальных решений интерполируются по граничным элементам [17]. Для получения решения во временной области используется модифицированный метод Дурбина [18, 19] численного обращения интегрального преобразования Лапласа. Для демонстрации надежности и точности предложенной гранично-элементной формулировки приведены примеры решения статической и нестационарной динамической задач. 1. Постановка краевой задачи Рассматривается однородное, полностью анизотропное трехмерное электромагнитоупругое твердое тело, занимающее объем , ограниченный поверхностью . При использовании квазистатического приближения для электрического и магнитного полей, подразумевая нулевые начальные условия и отсутствие массовых сил, свободных зарядов и тока проводимости, уравнения движения сплошной электромагнитоупругой среды в изображениях по Лапласу имеют следующий вид: (1) где s - это параметр интегрального преобразования Лапласа; - плотность материала. Здесь использована стандартная система сокращенных обозначений для обобщенных перемещений и обобщенного тензора модулей упругости : (2) (3) где , и - упругие перемещения, электрический и магнитный потенциалы, соответственно; , и обозначают тензоры упругих модулей, пьезоэлектрических и пьезомагнитных свойств соответственно; , и - магнитоэлектрические коэффициенты, тензоры диэлектрической и магнитной проницаемостей соответственно. Справедливы следующие соотношения: (4) Обобщенные граничные условия на границе выражаются следующим образом: (5) (6) где (7) Здесь через , и обозначены тензор напряжений Коши, вектор электрической индукции и вектор магнитной индукции соответственно; - компоненты единичного вектора внешней нормали к границе. 2. Гранично-элементная схема Гранично-интегральное представление обобщенных перемещений в точке коллокации x следует из расширенной формулы Сомильяны [9, 21, 22] и имеет вид: (8) где - матрица коэффициентов, равная (где - дельта Кронекера) на гладких участках границы а в угловых точках определяемая геометрией поверхности в окрестности точки коллокации x [16]; через и обозначены обобщенные динамические фундаментальные и сингулярные решения в изображениях по Лапласу в точке наблюдения y, вызванные действием сосредоточенного единичного источника, приложенного в точке коллокации x. Обобщенные сингулярные решения определяются следующим соотношением: (9) где - компоненты вектора внешней единичной нормали к границе в точке наблюдения y. Обобщенные фундаментальные решения в пространстве Лапласа могут быть записаны как сумма сингулярной и регулярной частей следующим образом: (10) где сингулярная часть является фундаментальным решением статической линейной электромагнитоупругой задачи. В соответствии с [13-15] сингулярная и регулярная части определяются как (11) (12) где (13) (14) (15) Здесь - собственные значения матрицы ; Q - число различных собственных значений, - фазовые скорости упругих волн; - волновые числа; n - единичный вектор направления распространения волн. Для пространственной дискретизации интегрального уравнения (8) используются смешанные граничные элементы [16]. Граница аппроксимируется совокупностью квадратичных четырехугольных элементов. Для того чтобы одновременно корректно моделировать непрерывность обобщенных перемещений и разрывный характер обобщенных поверхностных усилий, на каждом граничном элементе используются линейные и постоянные функции формы соответственно. Используется специальная схема интерполяции динамических частей фундаментальных решений по граничным элементам [17] на основе пятиузловых элементов. Это позволяет сократить время вычислений почти в два раза. Учет граничных условий и применение классической коллокационной схемы на узлах аппроксимации неизвестных граничных функций позволяет сформировать разрешающую систему линейных алгебраических уравнений, параметризованную комплексным параметром интегрального преобразования Лапласа s. Для получения решения во временной области используется схема численного обращения преобразования Лапласа. Если изображение определено на некотором наборе частот , то оригинал изображения может быть аппроксимирован следующим образом [18, 19]: (16) (17) (18) (19) определяются из следующей трехдиагональной системы: (20) определяются из системы (21) где , , , , , , и - некоторое вещественное число, большее чем действительные части всех особенностей . 3. Численные примеры Пример 1. Рассматривается электромагнитоупругий прямоугольный параллелепипед с размерами Центр параллелепипеда совпадает с началом координат. Используется материал со следующими параметрами [9]: Заданы следующие граничные условия: В таблице приведено сравнение полученных на различных сетках результатов с аналитическими решениями [9]. Следует отметить, что продемонстрировано очень хорошее соответствие аналитических и гранично-элементных решений, наибольшая относительная погрешность не превышает 0,45 %. Упругие перемещения , , электрический потенциал φ и магнитный потенциал ψ в точке Elastic displacements , , electric potential and magnetic potential in the point Сетка u1·10-10, м u3·10-9, м φ, V ψ·10-2, А 4×12×4 -6,351148 1,137324 1,896969 4,294634 6×18×6 -6,344309 1,136904 1,896607 4,295536 8×24×8 -6,341515 1,136717 1,896484 4,295896 10×30×10 -6,340156 1,136617 1,896430 4,296073 12×36×12 -6,339358 1,136556 1,896404 4,296168 Аналитические решения [9] -6,333160 1,136676 1,899100 4,278120 Пример 2. Рассматривается электромагнитоупругий куб с длиной ребра 1 м. На нижней грани заданы следующие граничные условия: ; на верхней грани х3 = 1 м приложена равномерно распределенная одноосная нагрузка , где - функция Хевисайда. Остальные грани куба свободны от обобщенных поверхностных усилий. В качестве материала взят композит состава феррит кобальта-титанат бария BaTiO3-CoFe2O4 (с объемной долей титаната бария BaTiO3 50 %) плотностью и со следующими параметрами [20]: Рис. 1. Упругие перемещения u3 Fig. 1. Elastic displacements u3 Рис. 2. Электрический потенциал φ Fig. 2. Electric potential φ Рис. 3. Магнитный потенциал ψ Fig. 3. Magnetic potential ψ На рис. 1-3 представлены гранично-элементные решения в серединной точке нагруженной грани для четырех сеток с 96 (сетка 1), 216 (сетка 2), 384 (сетка 3) и 600 (сетка 4) элементами. Полученные результаты отчетливо демонстрируют наличие быстрой сеточной сходимости, даже для магнитного потенциала. Заключение В данной работе представлена формулировка прямого подхода метода граничных элементов в пространстве Лапласа для решения статических и нестационарных динамических задач трехмерной линейной теории электромагнитоупругости. Обобщенные фундаментальные решения в изображениях по Лапласу представлены в виде суммы сингулярной и регулярной частей. Динамическая часть выражена как интеграл по единичной полусфере, сингулярная статическая часть - как интеграл по единичной окружности. Граничные элементы выбраны таким образом, чтобы можно было одновременно моделировать непрерывные обобщенные перемещения и разрывные обобщенные поверхностные усилия. Пространственная дискретизация основана на классической коллокационной схеме. Для получения решений во временной области используется метод численного обращения интегрального преобразования Лапласа. Представлены два численных примера: статическая задача о прямоугольном параллелепипеде под действием заданной нагрузки и задача о нестационарном отклике единичного куба, под действием равномерно распределенной одноосной нагрузки в виде функции Хевисайда по времени. Полученные результаты очень хорошо соответствуют точным решениям. Для каждой задачи продемонстрировано наличие сеточной сходимости.

About the authors

I P Markov

Research Institute for Mechanics, National Research Lobachevsky State University of Nizhni Novgorod

Email: teanku@gmail.com

References

Statistics

Views

Abstract - 115

PDF (Russian) - 42