К теории блочных и нано структур Текст научной статьи по специальности «Математика»

УДК 539.3

К ТЕОРИИ БЛОЧНЫХ И НАНО СТРУКТУР

© 2007 В.А. Бабешко, О.В. Евдокимова, О.М. Бабешко1

Введение

В работе [1], интегральный и дифференциальный методы факторизации развиваются для совокупности различных областей, граничащих между собой, называемых блочными структурами. При применении этого метода возникают сложности при удовлетворен граничных условий на границе блоков. Из [2] следует, что построенные методом факторизации решения краевой задачи на границе имеют наряду с классическими, также и составляющие в виде обобщенных функций — 6-функций и их производных. В настоящей работе показано, каким образом преодолеваются эти сложности при применении дифференциального метода факторизации в блочных структурах.

1. Ниже под блочными структурами понимаются материалы, занимающие ограниченные, полуограниченные или неограниченные области, называемые контактирующими блоками. Предполагается, что в блочной структуре каждый блок обладает своими специфическими свойствами поведения при воздействии физическими полями различной природы. Считается, что эти поля описываются краевыми задачами для систем связанных дифференциальных уравнений в частных производных с постоянными коэффициентами. Среды такого типа свойственны строению коры Земли [3], конструкционным материалам, находящимся в сложных физико-механических условиях [4], нано материалам, кристаллическим структурам разнотипной компоновки, материалам электроники [5]. Аналогичное строение имеют и различные материалы, в том числе создаваемые на основе комбинации как только нано размерных, так и макро и нано размерных составляющих [6]. В данной работе рассматривается случай структуры, блоки которой являются трехмерными. Отсутствие существенных ограничений на краевые задачи, описывающие свойства отдельных блоков, свидетельствует о том, что рассматриваемые блочные структуры могут иметь большое разнообразие свойств. В общем случае понятие блока включает в себя требова-

1 Бабешко Владимир Андреевич (babeshko@kubsu.ru), Евдокимова Ольга Владимировна, Бабешко Ольга Мефодьевна, Кубанского государственного университета, 350040, Россия, г. Краснодар, ул. Ставропольская, 149.

ние неизменности границы области задания краевой задачи, в том числе многосвязной, и требования ее кусочной гладкости. Каждый блок может быть как ограниченным, так и неограниченным и в нем могут протекать связанные процессы механики деформируемого твердого тела, гидромеханики, электромагнитные, диффузионные, тепловые, акустические и другие процессы. Введенные блочные структуры, являются более общими, чем кусочно-однородные структуры, предполагающие лишь скачкообразное изменение физических параметров среды при переходе от блока к блоку, с сохранением материала среды. Последнее выражается в скачкообразном изменении отдельных коэффициентов дифференциальных уравнений краевых задач, происходящем при переходе из блока в блок с сохранением типа краевой задачи. Естественно, введенные блочные структуры имеют более широкий спектр свойств, чем кусочно-однородные. Это вытекает как из разнообразия свойств блоков, их формы, характера взаимодействия блоков между собой, так и как результат взаимодействия физических полей блоков, ряд из которых ими излучается или трансформируется, проходя через них. Частным случаем блочных структур являются слоистые структуры. Такие структуры с плоскими границами для линейных краевых задач можно считать в настоящее время достаточно полно исследованными. Исследование блочных структур проводится, в основном, численными методами, для которых наличие неограниченных областей всегда создает трудности. Дифференциальный метод факторизации, являясь обобщением метода интегральных преобразований, позволяет, уже на этапе решения краевой задачи отвечать на ряд вопросов относительно свойств физических полей в каждом блоке. Целесообразно заметить, что интегральные преобразования в краевой задаче в области П для дифференциальных уравнений в частных производных оказываются удобным инструментом для исследований в тех случаях, когда дифференциальные уравнения, область П и функции, описывающие интегральное преобразование — согласованы. Под согласованностью понимается возможность трансформации дифференциальных уравнений в частных производных в обыкновенные в результате применения интегрального преобразования, и задание граничных условий на границе, описываемой постоянными значениями геометрических параметров. Такое свойство имеет место, если функции интегральных преобразований являются собственными для дифференциального оператора в рассматриваемой области П. В терминах топологической алгебры это свойство будет иметь место в том случае, если группы преобразований, порождаемые автоморфизмом многообразия П, имеют представления, которые инвариантны относительно дифференцируемого отображения векторного поля, заданного на этом многообразии. Для ряда простых областей, которые будем называть классическими, это преобразования Фурье в областях с плоскими границами, преобразования Бесселя в областях с круговыми границами, преобразования Бесселя-Лежандра в областях со сферическими границами, применяемые, например, к уравнениям с постоянными коэффициентами Гельм-

гольца, Шредингера, Ламе, Навье-Стокса [7]. В [8] показано, что эти и другие интегральные преобразования являются следствием тех или иных отображений на себя многообразий, порождающих группы преобразований пространства, их движений. Представления этих групп достигаются введением специальных функций, о которых говорилось выше. В тех случаях, когда области относятся к классическим, краевые задачи решаются сравнительно просто: сводятся к простым функциональным или обыкновенным дифференциальным уравнениям после интегрального преобразования, затем применяется интегральное обращение.

Если области сложные, то для анализа подобных краевых задач применим, дифференциальный метод факторизации, который позволяет сводить краевые задачи к функциональным уравнениям с понижением размерностей последующих задач.

2. Сформулируем краевую задачу для блочной структуры. Будем считать, что область П блочной структуры состоит из областей Пь, (Ь = 1,2,..., В) с границами дПу. Может оказаться, что часть границы блока является общей с границей другого блока, ее называем контактирующей. Остальная часть, не контактирующая, может быть свободной или подчиненной внешним воздействиям. Предполагается, что в каждой области Пу ставится краевая задача для систем дифференциальных уравнений в частных производных с постоянными, своими в каждой области, коэффициентами.

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

м N к Р

Ку(д х1,дх2,дхэ)фь =

А 5ртпк%Ьр ,х\Л2Хз 0

т=1 п=1 к=1 р=1

s = 1, 2, ..., Ру,

р ,Х1 Х2 хз

1)

Ащшпк = СОШг’ % = {фЬ1, ФЬ2, .... ФЬр }’ Ь = 1 2,..., В,

% = {%Л , % (х) = % (Х1, Х2, Хз) , X = {Х1, Х2, Хз} , X € Пь.

На общей, контактирующей, границе дПЬПдПй задаются следующие граничные условия сопряжения

^Ь(дХ1, дх2, дхз)%Ь + К^(дх1, дх2, дхз)%л = м1 ^ к р

= У У У У [ вь к%т)(п)(к) +вй к%(р)(п)(к) ] =

/ к / к / к / к 1 вршпк^Ьр)Х2хз зршпк^¿р,Х1 Х2хз J Ьа5’

т=1 п=1 к=1 р=1 (2)

5 = 1,2,...^ь0 < Р, х € дПЬ П дПй,

М1 < М, N ,< N’ К1 < К, Ь, й = 1,2,..., В.

Краевая задача исследуется в пространствах медленно растущих обобщенных функций Н5(П), описанных в [1].

Эти граничные условия в общем виде описывают случай контакта блоков, когда на общих границах принимается условие совпадения необходимых компонент физических полей, продиктованные соответствующими физическими законами. Например, в задачах термоэлектроупругости эти соотношения приведены в [9]. В частности, условия (2) могут быть намного проще — не иметь на границе внешних воздействий и представлять лишь требование равенства на общей границе решений и их производных. Но как указано выше, это обстоятельство не позволяет приравнивать производные от решений, записанных в интегральной форме, поскольку их составляющими в методе факторизации являются обобщенные функции [2]. Если речь идет о частях границ блоков, не являющихся общими ни с каким другим блоком, то на них принимаются граничные условия краевой задачи, рассмотренной в работы [1]. Изложим схему применения дифференциального метода факторизации к таким областям.

Следуя дифференциальному методу факторизации [1], сведем краевую задачу к системе функциональных уравнений, рассматривая каждую область Пь в отдельности. В результате приходим к системе функциональных уравнений вида

Вид и описание всех, участвующих в этой системе обозначений, заимствованы из [1], с добавлением индексов Ь. В частности, Юь — вектор внешних форм краевой задачи в области Пь.

Сопоставляя настоящий случай с рассмотренным в [1], заметим, что граничные условия (2) в общем случае содержат значения решений и их производных на границе, по крайней мере, из двух соседних областей. Этим блочные структуры существенно отличаются от изучавшихся в [1].

3. В соответствии с правилами дифференциального метода факторизации дальнейшее исследование задачи предполагает факторизацию матрицы-функции Кь(аз), даваемой формулой (3). Для этого выберем матрицу-функцию КЬ(аз, т) порядка Р — 1, получающуюся вычеркиванием строки и столбца под номером т у сопряженной матрицы-функции КЬ(аз), такую, что нули ^п ее определителя Qь(a3) = ¿е! КЬ^, т) не совпадают с нулями

дПь

ь (а) = ^-Ь ( гаі, іа2, газ) \\кьпт(а)\1,

(3)

Ь = 1,2,..., В.

[ВД’ т)]-1 = \\Q~b Qpsь\\.

Тогда элементы матрицы-функции К-1 (а?, -), имеющей вид

К-1 (а?.- =

1 0

1

$ т1 $ т2 с тт С тИ

0 1

допускают интегральное представление в форме

1 (2рзь(цз)М

$ГП

(Мз)й?Мз 11 ^тр(^з)

= аь(из)К(и3)(и3-а;)~{2 + Г К(о%) ’

т Ф р,

^тр(^з)

Кь(«р

^тр(аз)

+

¿тр(Ц)

&(а3)К„(ар ' ^ - а^)’

(5)

$тт(а?) = К-Х(а3), а3 е

2тр(а?) = £ '2р.ъ(а3)Мт(а3).

5=1

Здесь замкнутый контур Г+ занимает положение, при котором область Л+ содержит только нули ¿+, ¿5-, а область Х_ —только нули ^п. Замкнутый контур Г- охватывает область, содержащую все нули ¿5+, ¿1_, ^ Из этого представления следует, что элементы матрицы-функции К-1(а?, —) являются рациональными функциями, единственными особенностями которых являются нули ¿+, ¿—, причем член, содержащий их, К—1(а?, —), явно выделен.

Граничные условия в соответствии с алгоритмом применения дифференциального метода факторизации в случаях не контактирующих границ применяются в соответствии с правилами, оговоренными в [1] Удовлетворение граничных условий осуществляется по следующей схеме. Вначале граничные условия на не контактирующей границе для каждого отдельного блока вносятся в соответствующие векторы внешних форм функциональных уравнений (3). При контакте блоков, на общих границах соседних блоков выполняются условия сопряжения (2), которые могут, в зависимости от свойств описываемых полей, включать некоторые соотношения для решений и их производных. В простейшем случае — это равенство на общей границе решений и их производных при переходе из одного блока в соседний. Эти соотношения должны быть внесены в соответствующие векторы внешних форм функциональных уравнений (3), предварительно разрешенных относительно неизвестных на границе производных по нормали. Последнее обеспечит выполнение контактных граничных условий (2) при решении псевдодифференциальных уравнений, что доказывается по схеме, изложенной в [2].

Опуская дальнейшие преобразования, изложенные в [1], получаем в случае блоков, занимающих выпуклые области, представление решения в каж-

дом блоке в виде

ОО

%(xv) = JJJ К^М-И [[ ^ъе~1{а^da^daldal, xv е Пъ.

д&

Решение можно сделать более наглядным, вычислив интеграл по пара-

метру а3 по теории форм-вычетов Лере. В результате имеем ,^\ _ ^ || \ ' „-¡(а^+а^х*)

dalda2,

Л Хе"'<“:У;“;'=)[кл'(‘^)т«Н'02.С)«-г

— ТО

~к»(‘ц) Т-ь(«'„а'2,г1У''-4

р

1т+ъ{а\,а1Х+) = -У (Г , ь=1,2,...,В,

[ 1 2 г±) ^ Л <2ъЮ+)К№+)’ ’ ’ ’ ’

р=1 дП±ь

Т±ь = {0,0,..., 0, т±ь, 0,..., 0}.

В этой формуле граница дПь для выбранного х^ < 0, ху е Пь, разбита по следующему правилу:

И Ы = IЫь + ff Ыь •

ddb дО.+ь дО.-ь

JJ' Ыь exp(-mlx3) ^ 0, Imal

Ш+ь

JJ" Ыь exp(-ialx3) ^ 0, Imal

дП—ь

В случае, если блок вырождается в полупространство или слоистую среду, псевдодифференциальные уравнения, появляющиеся в процессе решения краевой задачи, вырождаются в алгебраические, после обращения которых решение строится в конечном виде [1].

В том случае, если рассматриваемый блок не является выпуклым телом, для исследования краевой задачи применяется метод обобщенной факторизации [10].

Замечание. Вопрос ’’сшивания” решений, получаемых в каждом блоке, в методе факторизации осуществляется автоматически при удовлетворении граничных условий. При наличии трещин, разломов или включений меньших размерностей, последние надо рассматривать как границы блоков. В результате получается однотипный алгоритм исследования блочных структур с указанными неоднородностями.

В том случае, когда блоки вырождаются в слои с плоско-параллельными границами, изложенный алгоритм исследования блочных структур приводит к решениям, полностью совпадающим с получаемыми традици-

онными методами интегральных преобразований по координатам, лежащим в граничной плоскости [11].

Отдельные фрагменты работы выполнены при поддержке грантов РФФИ, (06-01-00295), (06-01-08017), (06-08-00671), программа Юг России, проекты — (06-01-96802)-(06-01-96805), (06-05-96806),

(06-01-96634)-(06-01-96638), проект НШ-4839.2006.1, программ отделения ЭММПУ и Президиума РАН, выполняемых Южным научным центром РАН.

Литература

[1] Бабешко, В.А. Об интегральном и дифференциальном методах факторизации / В.А. Бабешко, О.М. Бабешко, О.В. Евдокимова // ДАН. -2006. - Т. 410. - №2. - С. 168-172.

[2] Бабешко, В.А. Выполнение граничных условий в дифференциальном методе факторизации / В.А. Бабешко, О.М. Бабешко, О.В. Евдокимова. - ДАН. - 2007. (в печати)

[3] Садовский, М.А. Деформирование геофизической среды и сейсмический процесс / М.А. Садовский, Л.Г. Болховитинов, В.Ф. Писаренко. -М.: Наука, 1987. - 104 с.

[4] Новацкий, В. Электромагнитные эффекты в твердых телах / В.Новацкий. - М.: Мир, 1986. - 160 с.

[5] Нефедов, Е.И. Широкополосные микрополосковые управляющие устройства СВЧ / Е.И. Нефедов, А.С. Саидов, А.Р. Тагилаев. - М.: Радио и связь, 1994. - 168 с.

[6] Гуткин, М.Ю. Физическая механика деформируемых наноструктур / М.Ю.Гуткин, И.А. Овидько. СПб.: Изд-во ”Янус”, 2003 - Т. 1. - С. 192; Т. 2. - С. 352.

[7] Тихонов, А.Н. Уравнения математической физики / А.Н. Тихонов, А.А. Самарский. - М.: Наука, 1966. - 724 с.

[8] Виленкин, Н.Я. Специальные функции и теория представления групп / Н.Я. Виленкин. - М.: Наука, 1991. - 576 с.

[9] Ворович, И.И. Динамика массивных тел и резонансные явления в деформируемых средах / И.И. Ворович, В.А. Бабешко, О.Д. Пряхина. -М.: Научный мир, 1999. - 248 с.

[10] Бабешко, В.А. Обобщенная факторизация в краевых задачах в многосвязных областях / В.А. Бабешко, О.М. Бабешко // Докл. РАН. -2003. - Т. 392. - №2. - С. 185-189.

[11] Бреховских, Л.М. Акустика слоистых сред / Л.М. Бреховских, О.А. Годин. - М.: Наука, 1989. - 412 с.

Поступила в редакцию 15/У/2007; в окончательном варианте — 15/У/2007.