Full Text

Введение В настоящее время при проектировании различных объектов новой техники актуальными являются вопросы взаимодействия полей различной физической природы. Во многих случаях, особенно для изделий аэрокосмической отрасли, возникает необходимость учета взаимовлияния электромагнитных и механических полей. Основные результаты в этом направлении получены для стационарных процессов (см., например, [1-4]). Численно-аналитические подходы к решению стационарных задач размерностью больше единицы, в том числе для тел со сферическими границами приведены, например, в [5]. Точное решение задачи о статическом воздействии магнитного поля на тело со сферической полостью получено в работе [6]. Соответствующие нестационарные задачи исследованы в основном в одномерном случае при частичном учете влияния электромагнитного поля. Нахождение точных решений нестационарных задач приобретает особую ценность в связи с тем, что численное обращение преобразования Лапласа по времени, к которому часто прибегают при решении, является некорректной задачей. В связи с этим, как показано, например, в [7], представление решения в виде ряда Лорана по параметру преобразования Лапласа в окрестности бесконечно удаленной точки позволяет получить фундаментальное решение для пространства на начальном промежутке времени. Основные принципы построения таких решений применительно к задачам нестационарной линейной электроупругости изложены в работе [8]. Ими можно эффективно воспользоваться, например, при использовании метода граничных элементов. В [9] приведено решение одномерной, но уже связанной задачи для толстостенной сферы в случае воздействия температуры на поверхности. Явные решения для пространства со сферической полостью для случая задания магнитного поля на границе сферы, причём в случае зависимости свойств материала от температуры, приводятся в работе [10]. Вопросы нахождения аналитических решений нестационарных задач для тел сферической формы, взаимодействующих с различными средами, в том числе с учётом связанности полей и наличием пьезоэффектов, исследованы в работах [11-15]. В работах [16-18] получено решение связанной нестационарной задачи электромагнитоупругости о распространении радиальных возмущений в толстостенной сферической оболочке и пространстве со сферической полостью. Таким образом, двумерные нестационарные задачи, в которых бы учитывалась связь электромагнитных и механических полей, в настоящее время изучены недостаточно, не говоря уже о получении точных решений подобного рода проблем. Количество публикаций, посвященных исследованию подобных двумерных задач, сравнительно мало. Например, в статье [19] дано решение лишь одной составляющей этой проблемы - определение нестационарного электромагнитного поля по заданному полю перемещений в пространстве со сферической полостью. А в работе [20] построено решение второй части этого вопроса - определение напряженно-деформированного состояния той же области, заполненной упругой средой и находящейся под действием объемных сил, под которыми можно понимать возмущения, создаваемые внешним электромагнитным полем. В данной статье для этого геометрического объекта предлагается решение полной связанной нестационарной двумерной задачи электромагнитоупругости. В настоящее время авторам не известны работы, посвящённые решению двумерных связанных задач электромагнитоупругости в нестационарной постановке применительно к пространству со сферической полостью. Предложенный метод позволяет находить точные решения подобного рода задач для любого момента времени, а не только на начальном этапе и в одномерном случае, как это предлагается некоторыми авторами в предложенном выше обзоре. 1. Постановка задачи В сферической системе координат рассматривается заполненное изотропным проводником пространство со сферической полостью радиусом . Осесимметричное движение среды описывается линеаризованной моделью [21]. Она включает в себя уравнения движения (1) уравнения Максвелла (2) и линеаризованные относительно начального состояния (его компоненты обозначаются дополнительным нижним индексом «0») формулы для радиальной и тангенциальной координат силы Лоренца и обобщенный закон Ома [22]: (3) (4) Здесь и , и , и - радиальные и тангенциальные координаты векторов перемещения, напряженности электрического поля, плотности электрического тока; - ненулевая координата напряженности магнитного поля; - плотность поверхностных зарядов; точками обозначены производные по времени. В формулах (1)-(4) и далее использованы следующие безразмерные параметры (при одинаковом начертании их размерные аналоги обозначены штрихом): где - размерное время; и - некоторые характерные линейный размер и напряженность электрического поля; и - плотность и упругие постоянные Ламе среды; и - скорости распространения волн расширения-сжатия и сдвига; и - коэффициенты электропроводимости, диэлектрической и магнитной проницаемости; - скорость света. В начальный момент времени среда является невозмущенной: (5) Искомые функции предполагаются ограниченными. Поскольку методы решения задач при всех возможных граничных условиях на поверхности полости идентичны, то далее ограничимся вариантом задания кинематических возмущений и напряженности электрического поля: . (6) Отметим, что из уравнений электромагнитодинамики (2) и (4) вытекают следующие соотношения для компонент векторов напряженностей магнитного и электрического полей, а также для плотности зарядов: (7) (8) . (9) Далее будем полагать, что начальное электромагнитное поле является стационарным, радиальным и удовлетворяет условиям (10) 2. Представление решения в виде рядов Решение начально-краевой задачи (1)-(6) представляем в виде рядов по полиномам Лежандра и Гегенбауэра [23]: . (11) Тогда с учетом (10) соотношения (1), (3), (4) и (7)-(9) переходят в следующие равенства для коэффициентов этих рядов: (12) ; (13) ; (14) (15) (16) . (17) Соответствующие начальные условия вытекают из (5): (18) При этом граничные условия (6) с учетом разложений и первого соотношения в (16) переходят в следующие равенства: (19) К ним добавляются условия ограниченности искомых функций. Для решения начально-краевых задач (12)-(19) используем преобразование Лапласа по времени ( - его параметр, верхний индекс « » соответствует изображению) [24]. При этом разрешающие уравнения (12) и (15) с учетом (13) трансформируются так: (20) . (21) Изображения коэффициентов других компонент электромагнитного поля и граничные условия согласно (16), (17) и (19) определяются следующим образом: (22) ; (23) (24) К ним опять же добавляются условия ограниченности изображений. Как показано в [16], аналитически найти оригиналы решения краевых задач (20)-(24) даже при невозможно. Поэтому будем использовать разложения искомых функций в степенные ряды по малому параметру : (25) Подставляя изображения по Лапласу этих рядов в (20)-(24), приходим к следующим соотношениям: при (26) (27) (28) ; (29) при (30) (31) ; (32) (33) . (34) Соответствующие граничные условия вытекают из (24) и разложений (25): ; (35) ; (36) ; (37) . (38) Соотношения (26)-(38) являются рекуррентной по индексу последовательностью краевых задач относительно ограниченных функций. 3. Интегральные представления решения Задача (26), (30), (35) является чисто упругой. Поскольку эти вопросы подробно исследованы в работе [25], далее в граничных условиях (6) положим, что . При этом эта задача становится однородной. Следовательно, ее решение тривиальное: . (39) Решение же задачи (27), (31), (36) записываем в интегральном виде ( ): ; (40) (41) где Ядра этих представлений - функции Грина краевых задач, соответствующих уравнениям (27), (31) и граничным условиям (36), а именно ограниченные решения следующих задач ( - дельта-функция Дирака): ; Эти функции найдены в [20] и имеют следующий вид ( - функция Хевисайда): (42) (43) Здесь (44) где Отметим, что числитель и знаменатель дробей в (44) являются экспоненциальным многочленом и просто многочленом аргумента соответственно. Подробный анализ этих формул показывает, что степени числителей и меньше степени знаменателя соответственно на единицу и на двойку, что позволяет точно находить их оригиналы с помощью методов компьютерной алгебры и теорем операционного исчисления. Также в интегральном виде при и записывается решение задачи (32), (38): (45) Здесь - соответствующая объемная функция Грина, а именно ограниченное решение следующей краевой задачи: . Решение задачи (21) при , (37) с учетом (39) имеет следующий вид: , (46) где - поверхностная функция Грина, т.е. ограниченное решение следующей краевой задачи: . Функции и , а также их оригиналы в квазистатическом приближении при найдены в [19] и имеют следующий вид: (47) . Формулы для коэффициентов рядов изображений координат напряженности электрического поля следуют из (22), (45) и (46): ; (48) (49) где Соответствующие коэффициенты для изображения плотности зарядов при определяются равенствами (34), а при в силу (39) имеет место равенство . 4. Разрешающая рекуррентная система уравнений Как следует из п. 3, коэффициенты рядов (25) при каждом определяются независимыми рекуррентными системами интегральных соотношений. При эта система включает в себя соотношения (34), (40) и первое равенство в (49) с начальными условиями (39). Она является однородной и имеет тривиальное решение: . При каждом рекуррентные системы образовываются соотношениями (41), (45), (49) и (34) при , которые в пространстве оригиналов с учетом свойств преобразования Лапласа записываются так: (50) ; (51) (52) . (53) Здесь и далее звездочка обозначает свертку по времени, а дополнительный нижний индекс « » у функции соответствует результату применения к ней следующего оператора: . Начальные условия к системе (50) - (53) следуют из (39), (46), (48) и (53): (54) В соотношения (50)-(53) входят производные по времени, а также, как следует из (15) и (17), производные по радиусу. Для того чтобы избежать численного дифференцирования, необходимо их модифицировать. Прежде всего, с помощью интегрирования по частям преобразовываем формулы (51) и (52). При этом полагаем, что начальная плотность поверхностных зарядов удовлетворяет условию , а также учитываем вытекающие из формул (47) и (49) равенства В результате приходим к следующим равенствам: ; (55) (56) Здесь где Для устранения производной в формуле (53) замечаем, что согласно (17) имеет место следующее равенство: . Веденная здесь функция имеет смысл коэффициента разложения в ряды по полиномам Лежандра, коэффициента объемного расширения для поля перемещений с компонентами и . При этом равенство (53) можно преобразовать так: . (57) Тогда необходимо дополнительно построить интегральное представление для . Его получаем из (50): , (58) где . Используя результаты работы [20] и формулы (42)-(44) для ядер в (58), получаем следующий результат: Здесь где Эти формулы аналогичны (42)-(44). Отличие состоит в том, что, как показывает подробный анализ, у функций степени числителя и знаменателя совпадают. Поэтому необходимо учитывать, что их оригиналы могут содержать слагаемые , (59) которые находятся методами компьютерной алгебры и должны быть в свертках в (58) в соответствии со свойствами дельта-функции. Далее, дифференцируя равенства (50), получаем следующие интегральные представления для производных по времени: (60) где Явный вид этих ядер следует из (42)-(44): где Изображения последних функций имеют структуру (44). При этом в знаменателе степень аргумента уменьшается на единицу. Поэтому функции и могут содержать аналогичные (59) слагаемые, которые также должны быть учтены при вычислении сверток в (60). Таким образом, разрешающая рекуррентная система уравнений состоит из соотношений (50), (55), (56), (57), (58), (60) при и начальных условий (54). Она позволяет находить коэффициенты рядов (11) и (25) для перемещений, напряженностей магнитного и электрического полей, а также плотности зарядов. Соответствующие коэффициенты для плотности тока могут быть найдены по формулам (14). 5. Пример расчета Полагаем, что материал пространства - алюминий, что соответствует следующим параметрам [26] (характерная напряженность электрического поля ): . Радиус полости единичный: , начальные параметры электрического поля следующие: , а на границе полости напряженность электрического поля имеет вид , где , что соответствует таким коэффициентам: . Интегралы в рекуррентных соотношениях находились численно. Распределение по радиусу нетривиальных коэффициентов рядов (11) при для перемещений и компонентов электромагнитного поля представлены на рис. 1-8: сплошные кривые соответствуют , пунктирные - , а штрихпунктирные . Расчеты проводились с учетом первых трех членов рядов (25). Учет последующего члена практически не приводит к изменению результатов. На рис. 1-8 видно характерное убывание механических характеристик и характеристик электромагнитного поля с ростом времени, при этом для тангенциальных перемещений этот процесс менее выражен. Для тангенциальной составляющей вектора напряжённости электрического поля характерно «опрокидывание» в отрицательную область. Для плотности электрических зарядов характерна концентрация вблизи границы полости. Рис. 1. Радиальные перемещения Fig. 1. Radial displacements Рис. 2. Тангенциальные перемещения Fig. 2. Tangential displacements Рис. 3. Напряженность магнитного поля Fig. 3. Magnetic field strength Рис. 4. Радиальная координата напряженности электрического поля Fig. 4. Radial coordinate of electric field strength Рис. 5. Тангенциальная координата напряженности электрического поля Fig. 5. Tangential coordinate of electric field strength Рис. 6. Плотность зарядов Fig. 6. Density of charges Рис. 7. Радиальная координата плотности тока Fig. 7. Radial coordinate of current density Рис. 8. Тангенциальная координата плотности тока Fig. 8. Tangential coordinate of current density Таким образом, в отличие от публикаций других исследователей, предложенный в работе метод связи механических и электромагнитных полей и впервые полученная с помощью этого подхода реккурентная система уравнений позволяют находить перемещения и компоненты электромагнитного поля новой нестационарной осесимметричной связанной задачи электромагнитоупругости для проводящего пространства со сферической полостью.

About the authors

V A Vestyak

Moscow Aviation Institute (National Research University)

Email: v.a.vestyak@mail.ru

E L Kuznetsova

Moscow Aviation Institute (National Research University)

Email: vida_ku@mail.ru

D V Tarlakovski

Moscow Aviation Institute (National Research University); Lomonosov Moscow State University

Email: tdvhome@mail.ru

References

Statistics

Views

Abstract - 356

PDF (Russian) - 41