Решение задач теории упругости методами Монте-Карло тема автореферата и диссертации по механике, 01.02.04 ВАК РФ

МОСКОВСКИЙ ОРДЕНА ЛЕНИНА, ОРДЕНА ОКТЯБРЬСКОЙ РЕВОЛЮЦИИ И ОРДЕНА ТРУДОВОГО КРАСНОГО ЗНАМЕНИ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ ИМЕНИ М.В. ЛОМОНОСОВА

МЕХАНИКО-МАТЕМАТИЧЕСКИЙ ФАКУЛЬТЕТ

СОРОКИН МИХАИЛ ВЛАДИМИРОВИЧ

УДК 539.3

РЕШЕНИЕ ЗАДАЧ ТЕОРИИ УПРУГОСТИ МЕТОДАМИ МОНТЕ-КАРЛО

01.02.04 — механика деформируемого твердого тела

На правах рукописи

Автореферат диссертации на соискание ученой степени кандидата физико-математических наук

Москва — 2004 г.

Диссертационная работа выполнена на кафедре Механики композитов Механико-математического факультета МГУ им. М. В. Ломоносова

Официальные оппоненты:

Пермский Государственный Технический Университет.

Защита состоится 15 октября 2004 г. в 16:00 час. на заседании диссертационного совета Д501.001.91 при МГУ им. М. В. Ломоносова по адресу: 119899 Москва, ГСП, Воробьевы горы, МГУ, Главное здание, Механико-математическй факультет, ауд. 16-10.

С диссертацией можно ознакомится в библиотеке Механико-Математического факультета МГУ (Главное здание, 14 эт.).

Автореферат разослан 14 сентября 2004 г.

Ученый секретарь диссертационного совета Д501.001.91 при МГУ

д. т. н., профессор Лурье С.А., к. ф.-м. н., с. н. с. Чистяков П.В.

Научный руководитель:

д. ф.-м. н., профессор Победря Б.Е.

Ведущая организация:

Шешенин С. В.

2.006-4 ЮТ

1 ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность теми. Методы статистического моделирования (методы Монте-Карло) являются классическими при решении задач теории переноса. Они также находят широкое применение в таких областях как статистическая физика, теория турбулентности, физико-химическая кинетика, т.е. там, где распространено статистическое описание тех или иных сложных физических процессов. Однако методы Монте-Карло решения краевых задач для уравнений в частных производных применяются значительно реже. Современное состояние исследований по данному вопросу, а также некоторые исторические сведения, можно найти в работах Блепова B.C., Ермакова С.М., Михайлова Г.А., Некруткина В.В., Сабельфельда К.К., Соболя И.М., Победри Б.Е., Чистякова П.В. и других авторов. В этих работах содержится изложение результатов по решению краевых и начально-краевых задач теории потенциала, диффузионных уравнений и уравнения теплопроводности. В них также рассмотрены алгоритмы статистического моделирования для некоторых задач теории упругости как для плоского, так и для трехмерного случаев.

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

Другой важной особенностью этих методов является сильное падение эффективности при росте размерности задачи. Прежде всего, это связано с большим объемом данных, которые необходимо загружать в оперативную память компьютера (или память на жестком диске), скорость работы с которой представляет узкое место практически для всех современных вычислительных систем.

Методы Монте-Карло имеют ряд отличительных особенностей в сравнении с детерминированными разностными и вариационными методами, которые позволяют во многом преодолевать описанные выше проблемы. Во-первых, зависимость требуемой памяти ЭВМ от размерности задачи, как правило, близка к линейной. Это позволяет эффективно решать ^пат» m rrngntt Во-вторых, решение

РОГ. \ Л КМ АЯ

I 'КА

' 1 . ipf

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

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

Наконец, одним из важнейших преимуществ методов Монте-Карло является их автоматическая распараллеленость, которая возникает из самой структуры алгоритмов, представляющих собой многократную реализацию независимых случайных траекторий. Число таких независимых реализаций может быть очень большим (несколько миллионов и больше), что позволяет использовать параллельные вычислительные блоки независимо друг от друга все время решения вплоть до объединения результатов счета в конце работы программы.

Исторически, методы Монте-Карло включают в себя три наиболее важных подхода, используемых для решения краевых задач. Первый подход заключается в приближенном вычислении континуальных интегралов, представляющих решение соответствующей краевой задачи. При этом основная сложность состоит в получении такого представления решения в виде интегралов по определенной мере, для вычисления которых может быть применен метод Монте-Карло.

Второй подход основан на использовании теорем о среднем и формул Грина для стандартных областей. Чаще всего в качестве такой области выбирается шар или сфера. Указанный подход приводит к методам блуждания внутри области, в частности, к методу блуждания по сферам.

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

Указанные выше преимущества методов Монте-Карло перед детерминированными методами определили естественный интерес к их применению для решения задач механики деформируемого твердого тела. В работах Ворошко П.П., Елепова B.C., Михайлова Г.А. и других авторов представлены различные алгоритмы решения плоских задач теории упругости методами Монте-Карло, а также некоторых частных

задач для пространственного случая.

Однако наибольший интерес с вычислительной точки зрения представляют пространственные задачи теории упругости. Здесь тоже практически сразу выделились два основных метода Монте-Карло: метод блуждания по сферам и метод блуждания по границе. Описание и различные аспекты применения первого метода можно найти в работах Победри Б.Е., Чистякова П.В., Сабельфельда К.К., при этом круг решаемых задач ограничивался первой и второй внутренними краевыми задачами теории упругости для однородной изотропной среды. Описание второго метода можно найти у Сабельфельда К.К., Симонова H.A., этот метод расширил область применения методов Монте-Карло на внешние краевые задачи тех же типов.

Здесь необходимо отметить ряд особенностей этих методов, которые определяют области их преимущественного применения.

1. Алгоритм блуждания по границе позволяет строить статистические оценки решения одновременно в произвольном множестве точек, используя одни и те же траектории случайных блужданий. Это позволяет применять этот метод и в тех случаях, когда требуется построить поле решения во всей рассматриваемой области.

2.Метод блуждания по границе может быть применен к решению как внутренних краевых задач, так и внешних, причем схема метода практически не изменяется. Кроме того, как будет показано ниже, этим методом решается более широкий класс задач теории упругости.

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

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

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

Карло для решения смешанной задачи теории упругости, когда на одной части границы заданы перемещения, а на остальной — нагрузки.

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

Методика исследований. Исследование свойств методов Монте-Карло и их модификаций, предложенных в диссертационной работе, основано на сочетании аналитических методов (оценивание трудоемкости, построение случайных функционалов) и численных экспериментов с использованием программ, написанных на Microsoft Visual С++ .NET.

Апробация работы. Материалы диссертации докладывались на научном и аспирантском семинарах кафедры механики композитов МГУ; на научных конференциях "Ломоносовские чтения"МГУ (секция механики) в 2003 и 2004 гг.

Публикации. Основное содержание диссертации отражено в пяти статьях и тезисах, указанных в конце данного автореферата.

Структура диссертации. Работа состоит из введения, четырех разделов, разбитых на параграфы, двух приложений и списка литературы.

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

1. Построен и исследован алгоритм блуждания по границе для решения смешанной задачи теории упругости. В трехмерном случае смешанные краевые задачи методами Монте-Карло ранее не решались. Введена новая модификация метода, основанная на сохранении оценок итераций интегрального оператора (см. [3,5]).

2. Впервые построен метод блуждания по границе для случая анизотропной упругости (см. [1,2,4]).

3. Метод Монте-Карло применен к решению задач о напряженно-деформированном состоянии шины. Приведены результаты сравнения этого метода и разностного метода решения интегральных уравнений (см. [5]).

2 ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

Во введении дается краткий обзор литературы по теме диссертационной работы и излагается краткое содержание диссертации.

В первом разделе излагается общая схема метода Монте-Карло, применяемая для решения систем линейных интегральных уравнений

в том виде, в котором она используется в диссертационной работе. В целом, данная схема соответствует работам Сабельфельда К.К. и Симонова H.A.. Приведены различные способы решения задачи как в случае сходящегося ряда Неймана, так и в случае, когда спектральный радиус интегрального оператора больше 1, что приводит к расходящемуся ряду Неймана. Второй случай требует применения методов построения резольвенты на основе ее аналитического продолжения, использования Паде-аппроксимаций или разложения в цепные дроби. Заметим, что эти методы могут применяться и в случае слабо сходящихся рядов Неймана, что практически всегда имеет место для интегральных уравнений, получающихся при решении реальных задач теории упругости. Эффективность методов Монте-Карло может быть существенно повышена при использовании разложения резольвенты в цепную дробь, в этом случае количество итераций интегрального оператора может быть гораздо меньше, чем в классическом случае, но требования к точности их подсчета сильно возрастают. В качестве методов ускорения сходимости рядов в работе применяются и способы обобщенного суммирования.

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

которые возникают при решении смешанной задачи теории упругости.

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

¥>»(*) = / kv(x>y)Vj{y)dy + /.(*). г" = 1, ■•

G

G

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

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

Кроме того, исследуются вопросы эффективного построения случайных траекторий блуждания в том случае, когда геометрия границы задается по точкам. Тогда значительное время вычислений уходит на поиск точек пересечения прямой, направленной вдоль случайного вектора, и граничной поверхности. В работе предлагается использовать алгоритмы, применяемые в компьютерной графике при расчете освещенности трехмерных объектов методом трассировки лучей (ray tracing). Для произвольного вида геометрии, задаваемой массивом простейших поверхностей, разработана программа эффективного поиска точек пересечения заданной прямой с границей рассматриваемой области, который является важнейшей частью алгоритма блуждания по границе.

Третий параграф раздела посвящен решению методом блуждания по границе двух основных краевых задач для уравнений Ламе:

(1 — 2и)Аи + graddiv« = О задачи Дирихле и задачи Неймана с граничными условиями

и

<т(и) • п|Е = 5°

соответственно. Приводятся результаты решения тестовых задач методом Монте-Карло и при помощи разностного метода решения интегральных уравнений, которые позволяют сделать вывод о высокой эффективности первого при вычислении напряженного состояния в фиксированной точке.

В последнем параграфе второго раздела излагается один из основных результатов диссертационной работы — алгоритм решения смешанной

задачи теории упругости методом блуждания по границе. Граничные условия в этом случае имеют вид:

и(г) = «"(*), * € Е„,

{сг(|.)-п}(*) = в0(*),

Вектор перемещения ищется в виде потенциала простого слоя и(я) = / Т(х,у)<р{у)(И:у

Е

с неизвестной векторной плотностью тогда на части границы Е„

♦ справедливо интегральное уравнение 1-го рода

^ / Г(г, уМ«) <*ЕУ = «»(г), « € Е„,

» Е

а на части границы Е^ — интегральное уравнение 2-го рода *,(*) = - / 1Т«(в>»Му)|Е, + 5°(г), г €

Е

которым удовлетворяет искомая плотность <р(у).

Данные интегральные уравнения можно записать в следующей форме с регуляризирующим параметром И:

( ф) = (I - К£)ф) + Ни0(г), г € £„, \ ф) = -ТГф) + г € Еа.

Эти уравнения представляются в виде одного функционального уравнения на граничной поверхности Е

Ф) = А-Ф) + Ъ(х),

где А — это определенный линейный оператор, действующий из Е в Е, а Ь(х) — известная векторная функция на поверхности Е. Такое представление позволяет использовать описанную в первом разделе схему для нахождения вектора перемещений в произвольной внутренней точке тела методом блуждания по границе.

В данном параграфе также приведен метод подсчета регуляризирующего параметра Л на основе изложенного в первом разделе способа оценки максимального собственного значения оператора.

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

где интеграл берется по дуге единичной окружности, плоскость которой перпендикулярна х, а центр совпадает с началом координат (г это длина вектора х; (ц>) зависит от тензора модулей упругости). Приведены условия, налагаемые на материальные константы, которые позволяют получить вычислительно устойчивый алгоритм. Отдельно рассмотрены варианты метода для трансверсальной изотропии и частного случая ортотропии, при котором задаются шесть материальных констант вместо девяти. В таких предположениях тензор фундаментальных решений может быть построен аналитически.

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

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

На Рис.1, приведена приблизительная форма профиля шины, которая использовалась при решении задач.

Рис.1.

Во втором параграфе раздела приводятся результаты решения данных задач о шине методом Монте-Карло, которые сравниваются с результатами, полученными при использовании разностного метода решения систем линейных интегральных уравнений. При проведении расчетов методом конечных разностей использовались регулярные сетки. Показаны границы эффективного применения алгоритмов статистического моделирования.

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

Автор благодарен профессору Б.Е.Победре за постановку задач и постоянное внимание к работе.

СПИСОК ПУБЛИКАЦИЙ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Победря Б.Е., Сорокин М.В. Применение методов Монте-Карло к решению задач теории упругости. // Тезисы докладов Ломоносовских чтений.—М.: Изд-во МГУ, 2003.-С.110.

2. Сорокин М.В. Решение второй краевой задачи методом Монте-Карло в случае анизотропной среды. // Вест. Моск. ун-та. Сер. 1, Математика. Механика. 2004. N 2. с.74—76.

3. Сорокин М.В. Решение смешанной задачи теории упругости методом Монте-Карло. // Вест. Моск. ун-та. Сер. 1, Математика. Механика. 2004. N 5. (в печати).

4. Сорокин М.В. Применение метода Монте-Карло к решению задач 1 линейной теории упругости для анизотропной среды. // Изв. РАН. Механика твердого тела. 2004. N 3.— С.38—44.

5. Сорокин М.В. Метод Монте-Карло в задаче об упругом • деформировании шины. // Тезисы докладов Ломоносовских чтений,—

М.: Изд-во МГУ, 2004— (в печати).

I

t

\

У

]

Подписано в печать 3Û, £>£ pif Формат 60*84/16. Усл.печ.л. 0,7f Тираж Л)0 экз. Заказ У3£ Отпечатано в Отделе печати МГУ

«

1

1 »

о/. ¿?S. ' ¿PS. .

РНБ Русский фонд

2006-4 1377

11 'гл 'к

Введение 1

1 Общая схема метода Монте-Карло для решения линейных интегральных уравнений. 9

1.1 Схема Неймана-Улама в случае сходящегося ряда Неймана. 9

1.2 Аналитическое продолжение резольвенты. Обобщенное суммирование ряда Неймана. . . 13

1.3 Паде-аппроксимации резольвенты и разложение ее в цепную дробь. 15

1.4 Определение максимального собственного значения интегрального оператора и регуляризация аналитического продолжения для интегральных уравнений 1-го рода.17

2 Решение основных краевых задач линейной теории упругости для однородной изотропной среды. 19

2.1 Интегральные уравнения теории упругости.19

2.2 Метод блуждания по границе.23

2.3 Задачи Дирихле и Неймана.39

2.4 Смешанная задача теории упругости.42

3 Решение задач линейной теории упругости для однородной анизотропной среды. 51

3.1 Частный случай трансверсальной изотропии.51

3.2 Численное построение фундаментального решения и оператора напряжений.

Произвольный тип анизотропии.55

4 Применение метода блуждания по границе для решения задач о напряженно-деформированном состоянии шины. 62

4.1 Постановка задач. Задание геометрии шины.62

4.2 Решение задач о шине методом блуждания по границе. Анализ результатов.66

4.3 Замечание о возможности применения метода блуждания по границе в случае неоднородной среды.69

Актуальность темы. Методы статистического моделирования (методы Монте-Карло) являются классическими при решении задач теории переноса [7, 15]. Они также находят широкое применение в таких областях как статистическая физика, теория турбулентности, физико-химическая кинетика, т.е. там, где распространено статистическое описание тех или иных сложных физических процессов. Однако методы Монте-Карло решения краевых задач для уравнений в частных производных применяются значительно реже. Первая книга [5], посвященная данному вопросу, вышла в 1980 г., вторая [8] - в 1984 г. Монография [27], вышедшая в 1989 г., содержит наиболее полное изложение результатов по решению краевых и начально-краевых задач теории потенциала, диффузионных уравнений и уравнения теплопроводности. В ней также рассмотрены алгоритмы статистического моделирования для некоторых задач теории упругости как для плоского, так и для трехмерного случаев.

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

Другой важной особенностью этих методов является сильное падение эффективности при росте размерности задачи. Прежде всего, это связано с большим объемом данных, которые необходимо загружать в оперативную память компьютера (или память на жестком диске), скорость работы с которой представляет узкое место практически для всех современных вычислительных систем.

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

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

Наконец, одним из важнейших преимуществ методов Монте-Карло является их автоматическая распараллеленость, которая возникает из самой структуры алгоритмов, представляющих собой многократную реализацию независимых случайных траекторий. Число таких независимых реализаций может быть очень большим (несколько миллионов и больше), что позволяет использовать параллельные вычислительные блоки независимо друг от друга все время решения вплоть до объединения результатов счета в конце работы программы.

Исторически, методы Монте-Карло включают в себя три наиболее важных подхода, используемых для решения краевых задач. Первый подход заключается в приближенном вычислении континуальных интегралов, представляющих решение соответствующей краевой задачи. При этом основная сложность состоит в получении такого представления решения в виде интегралов по определенной мере, для вычисления которых может быть применен метод Монте-Карло.

Второй подход основан на использовании теорем о среднем и формул Грина для стандартных областей. Чаще всего в качестве такой области выбирается шар или сфера. Указанный подход приводит к методам блуждания внутри области, в частности, к методу блуждания по сферам.

Третий подход связан с использованием глобальных интегральных уравнений на границе рассматриваемой области. Чаще всего, интегральные уравнения записываются не для самого решения, а для некоторой вспомогательной функции. Примером применения данного подхода может служить метод блуждания по границе на основе граничных интегральных уравнений теории потенциала. Также существуют различные модификации методов Монте-Карло, которые можно интерпретировать как частные случаи трех основных подходов. К ним относятся методы блуждания по решетке и методы решения систем линейных алгебраических уравнений [23, 31].

Данная классификация дана в монографии [27], там же можно найти подробный обзор работ по этому вопросу. Более подробно второй и третий подходы описаны в книгах [46] и [47]. Таким образом, в области общей теории решения краевых задач методами статистического моделирования выделились основные направления исследований, были получены эффективные алгоритмы для решения классических краевых задач уравнений математической физики.

Указанные выше преимущества методов Монте-Карло перед детерминированными методами определили естественный интерес к их применению для решения задач механики деформируемого твердого тела. В работе [4] был предложен метод решения краевой задачи для уравнения Ламе в плоском случае. В [5] и [17] методы Монте-Карло были применены к решению задачи об изгибе пластин. Более подробное изложение этого вопроса можно найти в [9]. Работа [28] посвящена алгоритмам блуждания по малым сферам для решения метагармонических уравнений, а в [30] рассматриваются неоднородные задачи теории пластин. В работе [48] были приведены методы решения частных задач для трехмерного случая (плоский и антиплоский сдвиги), которые могут быть сведены к классическим задачам теории потенциала.

Однако наибольший интерес с вычислительной точки зрения представляют пространственные задачи теории упругости. Здесь тоже практически сразу выделились два основных метода Монте-Карло: метод блуждания по сферам и метод блуждания по границе. Описание и различные аспекты применения первого метода можно найти в [23, 25, 27, 41, 42, 46], при этом круг решаемых задач ограничивался первой и второй внутренними краевыми задачами теории упругости для однородной изотропной среды. Описание второго метода дано в [26, 27, 29, 47], он расширил область применения методов Монте-Карло на внешние краевые задачи тех же типов.

Здесь необходимо отметить ряд особенностей этих методов, которые определяют области их преимущественного применения.

1.Алгоритм блуждания по границе позволяет строить статистические оценки решения одновременно в произвольном множестве точек, используя одни и те же траектории случайных блужданий. Это позволяет применять этот метод и в тех случаях, когда требуется построить поле решения во всей рассматриваемой области.

2.Метод блуждания по границе может быть применен к решению как внутренних краевых задач, так и внешних, причем схема метода практически не изменяется. Кроме того, как будет показано ниже, этим методом решается более широкий класс задач теории упругости.

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

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

К настоящему времени в рамках обоих подходов построен целый ряд алгоритмов для решения пространственных задач линейной теории упругости. Однако круг задач, решаемых методами статистического моделирования, еще недостаточно широк для их внедрения в вычислительную практику. Пространственные задачи теории упругости рассмотрены только для однородной изотропной среды. Также остается открытым вопрос о возможности применения методов Монте-Карло для решения смешанной задачи теории упругости, когда на одной части границы заданы перемещения, а на остальной — нагрузки.

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

Методика исследований. Исследование свойств методов Монте-Карло и их модификаций, предложенных в диссертационной работе, основано на сочетании аналитических методов (оценивание трудоемкости, построение случайных функционалов) и численных экспериментов с использованием программ, написанных на Microsoft Visual С++ .NET.

Апробация работы. Материалы диссертации докладывались на научном и аспирантском семинарах кафедры механики композитов МГУ; на научных конференциях "Ломоносовские чтения "МГУ (секция механики) в 2003 и 2004 гг.

Публикации. Основное содержание диссертации отражено в пяти статьях и тезисах, указанных в конце данного введения.

Структура диссертации. Работа состоит из введения, четырех разделов, разбитых на параграфы, двух приложений и списка литературы.

Основные результаты диссертационной работы состоят в следующем:

1. Построен и исследован алгоритм блуждания по границе для решения смешанной задачи теории упругости. Предложена новая модификация метода, основанная на сохранении оценок итераций интегрального оператора (см. [33, 35]).

2. Построен метод блуждания по границе для случая анизотропной упругости (см. [24, 32, 34]).

3. Метод Монте-Карло применен к решению задач о напряженно-деформированном состоянии шины. Приведены результаты сравнения этого метода и разностного метода решения интегральных уравнений (см. [35]).

Заключение.

1. Александров А.Я. Решение основных трехмерных задач теории упругости для тел произвольной формы путем численной реализации метода интегральных уравнений. ДАН СССР, 208, N2, - 1973.

2. Бурчуладзе Т.В., Гегелгш Т.Г. Развитие метода потенциала в теории упругости. Тбилиси: Мецниереба, 1985.—226с.

3. Бухин Б.Л. Введение в механику пневматических шин. Москва: Химия, 1988, 224с.

4. Ворошко П.П., Квитка А.Л., Цыбенко А. С. Применение метода случайных блужданий для решения задач теории упругости.— Проблемы прочности, 1973.—N4.—С.53—57.

5. Елепов Б.С., Кронберг А.А., Михайлов Г.А., Сабельфельд К.К. Решение краевых задач методом Монте-Карло. Новосибирск: Наука, 1980.—173с.

6. Ермаков С.М. Метод Монте-Карло и смежные вопросы. М.: Наука, 1975.—471с.

7. Ермаков С.М., Михайлов Г.А. Статистическое моделирование. М.: Наука, 1982—295с.

8. Ермаков С.М., Некруткин В.В., Сипин А.С. Случайные процессы для решения классических уравнений математической физики. М.: Наука, 1984—205с.

9. Купрадзе В.Д. Методы потенциала в теории упругости. М.: Физматгиз, 1963.—472с.

10. Лаврентьев М.М., Романов В.Г., Шишатский С. П. Некорректные задачи математической физики и анализа.— М.:Наука, 1980.—269с.

11. Михайлов Г.А. Некоторые вопросы теории методов Монте-Карло. Новосибирск: Наука, 1974.—142с.

12. Михайлов Г.А. Оптимизация весовых методов Монте-Карло. М.: Наука, 1987.—239с.

13. Михайлов Г.А., Сабельфельд К.К., Ченцов Н.Н. Векторно-статистические модели некоторых задач математической физики.—В кн.: Актуальные проблемы прикладной математики и математического моделирования.—Новосибирск: Наука, 1982.— С.69—82.

14. Михлин С. Г. Лекции по линейным интегральным уравнениям. М.: Физматгиз, 1959—232с.

15. Михлин С. Г. Многомерные сингулярные интегралы и интегральные уравнения. М. Физматгиз, 1962.

16. Партон В.З., Перлин П. И. Интегральные уравнения теории упругости. М.: Наука, 1978.

17. Порее В.Н. Компьютерная графика.— СПб.: БХВ-Петербург, 2004.— 432с.

18. Победря Б.Е. Задача в напряжениях для анизотропной среды // ПММ. 1994. Т.58. Вып.1. С.77-85.

19. Победря Б.Е. Численные методы в теории упругости и пластичности. М.: Изд-во МГУ, 1981.—343с.

20. Победря Б.Е., Сорокин М.В. Применение методов Монте-Карло к решению задач теории упругости. // Тезисы докладов Ломоносовских чтений.—М.: Изд-во МГУ, 2003.—С.110.

21. Победря Б.Е., Чистяков П. В. Решение пространственных задач теории упругости методом Монте-Карло // ПММ. 1988. Т2. С.341— 345.

22. Сабельфельд К. К. Векторные алгоритмы метода Монте-Карло для решения систем эллиптических уравнений 2-го порядка и уравнений Ламе // Докл.АН СССР.-1982.-Т.262, N5.-C. 1076-1080.

23. Сабельфельд К.К. Методы Монте-Карло в краевых задачах. Новосибирск: Наука, 1989.—280с.

24. Сабельфельд К.К., Макаров С.Е. Алгоритм блуждания по малым сферам для решения метагармонических уравнений // Методы и алгоритмы статистического моделирования.—Новосибирск: Наука. Сиб. отд-ние, 1983.—С.40-47.

25. Сабельфельд К.К., Симонов Н.А. Решение пространственных задач теории упругости в детерминированной и стохастической постановках методом Монте-Карло // Докл.АН СССР.—1984.— Т.285, N4.—С.802—805.

26. Сабельфельд К.К., Ультан А.Е. Решение некоторых неоднородных задач теории пластин методом Монте-Карло // Численные методы механики сплошной среды.—Новосибирск: Наука. Сиб. отд-ние, 1981.-Т.12, N 2.—С.117-124.

27. Соболь И.М. Численные методы Монте-Карло. М.: Наука, 1973.— 131с.

28. Сорокин М.В. Решение второй краевой задачи методом Монте-Карло в случае анизотропной среды. // Вест. Моск. ун-та. Сер. 1, Математика. Механика. 2004. N 2. с.74—76.

29. Сорокин М.В. Решение смешанной задачи теории упругости методом Монте-Карло. // Вест. Моск. ун-та. Сер. 1, Математика. Механика. 2004. N 5. (в печати).

30. Сорокин М.В. Применение метода Монте-Карло к решению задач линейной теории упругости для анизотропной среды. // Изв. РАН. Механика твердого тела. 2004. N 3.— С.38—44.

31. Сорокин М.В. Метод Монте-Карло в задаче об упругом деформировании шины. // Тезисы докладов Ломоносовских чтений.—М.: Изд-во МГУ, 2004.—(в печати).

32. Трехмерные задачи математической теории упругости. // Купрадзе В.Д., Гегелиа Т.Г., Башелейшвили М.О., Бурчуладзе Т.В. — М.: Наука, 1976.—790с.

33. Трусов П.В. Механика сплошной среды. // Пермь: Изд-во ПГТУ, 1996.—ч.З.—141с.

34. Трусов П.В., Дударь О.И., Келлер Н.Э. Тензорные алгебра и анализ. // Пермь: Изд-во ПГТУ, 1998.—132с.

35. Фоли Док., ван Дэм А. Основы интерактивной машинной графики. В 2-х книгах.—М.: Мир, 1985.

36. Фридман В.М. Метод последовательных приближений для интегрального уравнения Фредгольма 1-го рода // Успехи мат. наук.—1956.—Т. 11, вып. 1(67).— С.233-234.

37. Чистяков П. В. О смещенной оценке решения пространственной задачи теории упругости методом Монте-Карло. // Теория вероятн. и ее применен., 1984, т.29., N 1, — С.174—175.

38. Чистяков П. В. О решении пространственной задачи теории упругости методом Монте-Карло. // Вест. Моск. ун-та. Сер. 1, Математика. Механика. 1984. N 6, — С.88—90.

39. Fair W. Continued fraction solution to Fredholm integral equations // Rocky Mountains J.Math—1974—Vol.4.—P.357—360.

40. Fichera G. Existence theorems in elasticity. Springer-Verlag. Berlin-Heidelberg, N.Y., 1972. Перевод: Фикера Г., Теоремы существования в теории упругости. "Мир", 1974.]

41. Sabelfeld K.K. Shalimova I.A. Spherical means for PDEs. VSP, Utrecht, The Netherlands, 1997.

42. Sabelfeld K.K., Simonov N.A. Random walks on boundary for solving PDEs. Utrecht: VSP, 1994. 137 p.

43. Shia D., Hui G. Y. A Monte Carlo solution method for linear elasticity // Int. Journ. of Solids and Structures. 2000. N37. P.6085—6105.