Гранично-элементное моделирование динамики составных вязкоупругих тел на основе модифицированных методов квадратур сверток и дурбина тема диссертации и автореферата по ВАК РФ 01.02.04, кандидат физико-математических наук Белов, Александр Александрович

• Белов, Александр Александрович
• кандидат физико-математических наук
• 2008

Оглавление диссертации кандидат физико-математических наук Белов, Александр Александрович

ВВЕДЕНИЕ.

Глава I. Математические модели и методы решения.

1.1. Математическая модель.

1.1.1. Постановка упругодинамической краевой задачи.

1.1.2. Определяющие соотношения линейной теории вязкоупругости.

1.2. Фундаментальные и сингулярные решения изотропной теории упругости и вязкоупругости.

1.3. Граничные интегральные уравнения.

1.3.1. Обобщенная формула Сомильяны и гранично-временное интегральное уравнение изотропной теории упругости.

1.3.2. Построение ГИУ для решения задач о колебаниях кусочно-однородных тел.

1.3.3. Граничные интегральные уравнения вязкоупругости.

1.4. Численное обращение преобразования Лапласа.

1.4.1. Интегрирование быстро осциллирующих функций.

1.4.2. Обращение преобразования Лапласа методом Дурбина.

1.4.2.1. Алгоритм метода Дурбина с аппроксимацией трансформанты.

1.4.2.2. Комбинированные формулы метода Дурбина.

1.4.3. Численное обращение преобразования Лапласа на основе метода Дурбина.

1.5. Метод квадратур сверток.

1.5.1. Традиционный метод квадратур сверток.

1.5.2. Модификации метода квадратур сверток.

1.6. Численные результаты на основе метода квадратур сверток.

Глава II. Методика гранично-элементного моделирования.

2.1. Гранично-элементная дискретизация.

2.1.1. МГЭ-схема для однородных тел.

2.1.2. МГЭ-схема для кусочно-однородных задач в изображениях.

2.1.3. Разрешающая система алгебраических уравнений и поэлементное интегрирование.

2.1.4. Учет симметрии задачи.

2.2. Построение дискретных аналогов на основе квадратур сверток.

2.3. Визуализация гранично-элементного моделирования.

2.3.1. Описание формата входных данных.

2.3.2. Пользовательский интерфейс.

2.3.3. Панели инструментов.

2.3.4. Работа мышью.

2.4. Программная реализация.

2.5. Тестовые задачи.

2.5.1. Задача о сферической полости.

2.5.2. Задача о действии скачка давления на торец составного призматического тела.

2.5.2.1. Численное ГЭ-моделирование на основе метода Дурбина.

2.5.2.2. Численное ГЭ-моделирование на основе метода квадратур сверток.

2.5.2.3. Численное ГЭ-исследование кривых деформаций.

2.5.2.4. Сравнение методов квадратур сверток и Дурбина.

Глава III. Решение прикладных задач.

3.1. Задача о действии вертикальной силы на поверхность упругого полупространства.

3.2. Задача о действии вертикальной силы на поверхность вязкоупругого полупространства.

3.3. Задача о действии вертикальной силы на поверхность полупространства с полостью.

3.4. Задача о действии давления внутри сферической полости, расположенной в упругом полупространстве.

3.5. Задача о действии давления внутри кубической полости, расположенной в упругом полупространстве.

3.6. Задача о штампе на полупространстве.

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

Введение диссертации (часть автореферата) на тему «Гранично-элементное моделирование динамики составных вязкоупругих тел на основе модифицированных методов квадратур сверток и дурбина»

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

Метод граничных интегральных уравнений (ГИУ) и метод граничных элементов (МГЭ) являются универсальным численно-аналитическим подходом к решению трехмерных волновых начально-краевых задач теории вязкоупругости. В некоторых случаях (бесконечные и полубесконечные тела и среды) метод ГИУ и МГЭ являются наиболее предпочтительным универсальным подходом к исследованию волн. К настоящему времени сформировались два основных направления развития метода ГИУ и МГЭ для исследования волновых задач: использование интегрального преобразования и построение шаговых процедур. Численные схемы МГЭ на основе интегральных преобразований обладают ограничением на тип решаемых задач (только когда справедлив принцип соответствия). Точность таких схем существенно зависит от параметров метода численного обращения интегрального преобразования. Считается, что для МГЭ наиболее предпочтителен метод Дурбина. Численные схемы МГЭ, опирающиеся на шаговые процедуры (интерполирование по времени), существенно зависят от выбора шага и применимы только к тем задачам, для которых возможно построение тензора Грина. Так как для волновых задач вязкоупругости тензор Грина во времени в общем случае не существует, то такие численные схемы не эффективны. Однако, возможно построение шаговых схем МГЭ на основе сочетания обоих подходов. Таким сочетанием обладает метод квадратур сверток ориентированный на гранично-временные интегральные уравнения вязкоупругости.

Для рассмотрения задачи о распространении волны, необходимы динамические формулировки исходной системы дифференциальных уравнений и соответствующие ей ГИУ. Использование ГИУ в задачах о исследовании волн начинается с работ В. Вольтерра и Ж. Адамара [1], а в динамической теории упругости с работ Ч.Х. Мюнца [135]. Исследование волновых потенциалов начинается с работ С.Г. Михлина и В.Д. Сапожниковой [80]. Построение соответствующих обобщенных потенциалов теории упругости началось с работ A.Y. De Hoop [108], В. Новацкого [82] и др. Существенный результат по исследованию волновых и обобщенных потенциалов теории упругости можно найти у В.Д. Купрадзе с соавторами [88], Т.В. Бурчуладзе и Т.Г. Гегелиа [45], И.Ю.

Чудиновича [95, 104], Т. Ha-Duong [122] и др. Первая граничная интегральная формулировка для упругодинамики была опубликована Т.A. Cruse и F.J. Rizzo [106, 107]. Эта формулировка применялась в сочетании с преобразованием Лапласа. Соответствующая формулировка в сочетании с преобразованием Фурье была представлена J. Dominguez [109-111]. Распространение ГИУ с преобразованием Лапласа на численное решение трехмерных динамических задач теории упругости было осуществлено в работе И.З. Ройтфарба и Чу Вьет Кыонга [87]. Конкретных примеров решения трехмерных нестационарных динамических задач авторы не привели. Некоторые результаты были получены в [2]. Первая формулировка граничного элемента непосредственно во временной области была представлена W.J. Mansur и использовалась для скалярного волнового уравнения и упругодинамики с нулевыми начальными условиями [133, 134]. Обобщение этой формулировки для ненулевых начальных условий было представлено Н. Antes [97, 98]. Сравнительное исследование возможностей метода граничных элементов в решении задач упругодинамики приводилось в работе G.D. Manolis [130]. Имеется исследование [136], после которого в сочетании с МГЭ-решением на основе преобразования Лапласа, применяется метод Дурбина, как наиболее подходящий для получения численного МГЭ-обращения. Все формулы численного обращения интегрального преобразования зависят от надлежащего выбора их параметров, поэтому разрабатываются МГЭ-схемы во временной области. Н.М. Хуторянским было осуществлено распространение МГЭ-подхода непосредственно во временной области на решение трехмерных нестационарных динамических задач теории упругости [91] на основе применения сплайна для аппроксимации по временной переменной. При этом в случае конечной границы тела была получена конечно-шаговая численная схема. Дальнейшее развитие методика получила в работах [90, 92, 94], где были рассмотрены различные схемы решения для основных типов краевых задач. Первые трехмерные ГЭ-расчеты на основе постановки для второй основной задачи получены с помощью шаговой схемы Н.М. Хуторянским и В.В. Туриловым [91, 89] и на основе их дискретной модели - в [2]. Детальный обзор по упругодинамическим аспектам граничных элементов и их применению содержится в работах D.E. Beskos [100, 101]. Из последующих работ по построению МГЭ-схем отметим [96, 99, 102, 103, 113, 114, 138, 149, 150, 155]. Все пошаговые процедуры требуют адекватного выбора шага времени. Неверно выбранный шаг по времени приводит к неустойчивости или численному демпфированию. Улучшение устойчивости традиционных пошаговых формулировок граничных элементов возможно на различных основах [21, 97, 102, 105, 123, 132, 139, 141, 143, 146, 147, 151-153].

Ни одна из традиционных пошаговых МГЭ-формулировок не может быть обобщена на вязкоупругий случай. Но можно построить пошаговую МГЭ-формулировку не на основе сплайн-аппроксимации интеграла Вольтерра, а на основе метода квадратур г" сверток. Первая такая формулировка была предложена для динамических задач упругости в [126, 127, 143]. Этот подход использует область Лапласа фундаментального решения, и полученные результаты не только более устойчивы, но и позволяют учитывать действие демпфирования в случае вязкоупругости. Из работ по развитию методов ГИУ и МГЭ в сочетании с методом квадратур сверток отметим следующие: [121, 128, 129].

Принципиально другой подход для динамических задач - это МГЭ с двойным применением теоремы взаимности и метод, опирающийся на линейные регулярные ГИУ. Метод с двойным применением теоремы взаимности предложен D. Nardini и С.А. Brebbia [137, 148]. Построение регулярных линейных ГИУ восходит к работам В.А. Бабешко [8, 13] и развивалось В.А. Бабешко, О.М. Бабешко, О.В. Евдокимовой [7, 16], А.О. Ватульяном [46-51], JI.A. Игумновым [20-22, 64] и др.

Отдельно следует сказать о задачах полупространства. Методы ГИУ и МГЭ в таких задачах являются частью общего направления - интегрального подхода - по их решению [3, 4, 53, 71, 73, 74, 76-79, 119, 140]. Решения интегральных уравнений, их систем являются основой для изучения соответствующих процессов динамического контактного взаимодействия и, в особенности, возникающих при таком взаимодействии резонансных явлений. Вопросы существования и единственности решения интегральных уравнений детально изложены в публикациях и монографиях И.И. Воровича и В.А. Бабешко [5, 6, 9-12, 14, 15, 56-58] и др. К эффективным методам решения таких интегральных уравнений и систем интегральных уравнений относятся метод факторизации [14, 17, 18], метод фиктивного поглощения [11, 17, 59, 84, 85] и др. В работах [10, 12, 15] было предложено обобщение метода фиктивного поглощения для слоисто-неоднородного полупространства и эффективность продемонстрирована в работах [54, 55, 83, 70-72, 83, 86] и др. Основным достоинством методов факторизации и фиктивного поглощения является высокая степень точности учета динамических свойств среды.

Формулировки МГЭ для динамических задач вязкоупругости не так разнообразны: используются представления ГИУ в изображениях по Лапласу и Фурье. Применяется принцип соответствия с упругодинамическими задачами [118, 124, 131]. Ограничения такого подхода: требование численного обращения интегрального преобразования и невозможность задания граничных условий контакта во времени. Формулировки непосредственно во временной области требуют знания вязкоупругих фундаментальных решений, которые не известны для общего случая. Только для асимптотической модели, максвелловой модели и модели Кельвина-Фойгта, модифицированных моделей Кельвина-Фойгта и стандартного вязкоупругого тела могут быть получены аналитически фундаментальные и сингулярные решения [21, 91, 115], а значит, возможно применение МГЭ [21, 67, 91, 145]. Громоздкость и не технологичность таких формулировок исключает развитие метода в этом направлении. Первые ГЭ-расчеты на основе смешанных постановок для осесимметричных задач вязкоупругости с применением шаговой схемы были получены Н.М. Хуторянским и Л.А. Игумновым [67]. Новая шаговая МГЭ-формулировка для задач вязкоупругости опубликована в работах М. Schanz и Н. Antes [142]. Формулировка основана на методе квадратур сверток. Эта формулировка пользуется преимуществом квадратурной формулы, в которой веса интегрирования определяются существующим преобразованным фундаментальным решением по Лапласу.

Таким образом, в развитии МГЭ для применения к трехмерным динамическим задачам упругости и вязкоупругости заметную роль сыграли работы таких ученых, как Ш.М. Айталиев, Л.А. Алексеева, В.А. Бабешко, А.О. Ватульян, И.И. Ворович, Р.В. Гольдштейн, Л.А. Игумнов, М.И. Лазарев, А.Н. Соловьев, А.Г. Угодчиков, Н.М. Хуторянский, И.Ю. Чудинович, J.D. Achenbach, Н. Antes, D.E. Beskos, М. Bonnet, Н.В. Coda, J. Dominguez, L. Gaul, L.J. Gray, S. Kobayashi, M. Kogl, W. Kress, G.D. Manolis, M. Marrero, A.D. Mesquita, Y. Niwa, G.H. Paulino, M. Shanz и др. В развитии методов решения соответствующих интегральных уравнений механики деформируемого твердого тела существенные результаты получили В.М. Александров, В.А. Бабешко, О.М. Бабешко, А.О. Ватульян, И.И. Ворович, Е.В. Глушков, Н.В. Глушкова, Р.В. Гольдштейн, О.В. Евдокимова, Л.А. Игумнов, В.В. Калинчук, A.M. Линьков, А.В. Манжиров, Н.Ф. Морозов, Д.А. Пожарский, О.Д. Пряхина, А.Г. Угодчиков, Н.М. Хуторянский и др.

В работе разрабатывается подход совместного использования МГЭ и интегрального преобразования Лапласа с использованием метода Дурбина для решения проблемы численного обращения интегрального преобразования. Кроме того, разрабатывается МГЭ на основе метода квадратур сверток.

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

Научная новизна работы заключается в следующем. ГЭ-моделирование с использованием, преобразования^ Лапласа впервые реализовано на основе \ комбинированных формул метода Дурбина, учитывающих специфику поведения подынтегрального выражения интеграла Меллина на всем частотном диапазоне. Представлены аппроксимации пошагового (по частоте) кусочно-линейного и кусочно; квадратичного поведения как изображения, так и всего подынтегрального выражения, с возможностью рассмотрения переменного шага.

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

Создано программное обеспечение для компьютерного анализа динамики а ' Ч составных вязкоупругих трехмерных тел. Представлены результаты численного \ моделирования и верификации предложенных методик и схем.

Достоверность исследований основана на эквивалентности исходной краевой/начально-краевой задачи в частных производных математической теории I вязкоупругости системе используемых граничных/гранично-временных интегральных уравнений; на использовании для численных исследований регуляризованных

ГИУ/ГВИУ; на детально проработанных алгоритмах МГЭ-подхода и оттестированном программном обеспечении; на сравнении полученных результатов с решениями других авторов.

I Практическая значимость результатов исследования состоит в развитии методов

Дурбина, квадратур сверток, граничных элементов с целью получения устойчивых f высокоточных численных решений трехмерной теории вязкоупругости; в создании МГЭпрограммного обеспечения для анализа динамики составных трехмерных вязкоупругих тел с использованием интегрального преобразования Лапласа; в создании МГЭ-программного обеспечения для анализа динамики составных трехмерных вязкоупругих i тел на основе шаговой схемы; в создании программного обеспечения граничноэлементной визуализации.

Основные положения, выносимые на защиту

1. Методика численного решения систем ГИУ прямого подхода в сочетании с методом квадратур сверток и соответствующее программное обеспечение для расчета во времени неизвестных волновых полей трехмерных вязкоупругих составных тел;

2. методика численного решения систем ГИУ прямого подхода в сочетании с методом Дурбина и соответствующее программное обеспечение для расчета неизвестных волновых полей трехмерных вязкоупругих составных тел;

3. модифицированные методы квадратур сверток и Дурбина;

4. программное обеспечение по визуализации результатов гранично-элементного моделирования;

5. ГЭ-решение и анализ следующих волновых задач:

- о действии скачка давления на торец составного призматического тела;

- о действии вертикальной силы на поверхность вязкоупругого полупространства;

- о действии вертикальной силы на поверхность полупространства с полостью (сферическая, кубическая);

- о действии скачка давления внутри полости (сферическая, кубическая), расположенной в упругом полупространстве; I

- о штампе на полупространстве;

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

Апробация работы

Результаты диссертационной работы докладывались на Всероссийской научной конференции, посвященной памяти профессора А.И. Весницкого (Н.Новгород, 2004); IX-X Нижегородских сессиях молодых ученых (Саров, 2004, 2005) и XI-XIII Нижегородских сессиях молодых ученых (Семенов, 2006, 2007, 2008); IX Всероссийском съезде по теоретической и прикладной механике (Н.Новгород, 2006); Всероссийской научно-технической конференции — фундаментальные проблемы машиноведения: новые технологии и материалы (ИМАШ РАН, 2007); XXII Международной конференции «Математическое моделирование в механике деформируемых тел и конструкций. Методы граничных и конечных элементов» (Санкт-Петербург, 2007); Итоговой научной конференции учебно-научного инновационного комплекса «Модели, методы- и программные средства» (Н.Новгород, 2007); VI Всероссийской научно-практической конференции «Машиностроение: наука, техника, образование» (Саранск, 2007); XI Международной конференции «Современные проблемы механики сплошной среды»

Ростов-на-Дону, 2007); XIV международный симпозиум «Динамические и технологические проблемы механики конструкций и сплошных сред» им. А.Г. Горшкова (Москва, 2008).

Публикации. Основные результаты диссертации опубликованы в работах [24-44, 52, 60-63, 65, 66, 68, 69].

Структура и объем работы. Диссертационная работа состоит из введения, трех глав, заключения, списка литературы из 155 наименований. Общий объем диссертации составляет 180 страниц машинописного текста.

Заключение диссертации по теме «Механика деформируемого твердого тела», Белов, Александр Александрович

ЗАКЛЮЧЕНИЕ

1. Развита методика численного решения систем ГИУ прямого подхода в сочетании с методом квадратур сверток и создано соответствующее программное обеспечение для расчета во времени неизвестных волновых полей трехмерных вязкоупругих составных тел;

2. Развита методика численного решения систем ГИУ прямого подхода в сочетании с методом Дурбина и создано соответствующее программное обеспечение для расчета неизвестных волновых полей трехмерных вязкоупругих составных тел;

3. Разработаны модификации методов квадратур сверток и Дурбина;

4. Создано программное обеспечение по визуализации результатов гранично-элементного моделирования;

5. Получены ГЭ-решения и проведен анализ тестовых и прикладных волновых задач, в том числе:

- о действии скачка давления на торец составного призматического тела;

- о действии вертикальной силы на поверхность вязкоупругого полупространства;

- о действии вертикальной силы на поверхность полупространства с полостью (сферическая, кубическая);

- о действии скачка давления внутри полости (сферическая, кубическая), расположенной в упругом полупространстве;

- о штампе на полупространстве;

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

Список литературы диссертационного исследования кандидат физико-математических наук Белов, Александр Александрович, 2008 год

1. Адамар, Ж. Задача Коши для линейных уравнений с частными производными гиперболического типа / Ж. Адамар. М.: Наука, 1978. - 352 с.

2. Айталиев, Ш.М. Метод граничных интегральных уравнений в задачах динамики упругих многосвязных тел / Ш.М. Айталиев и др.. Алма-Ата: Гылым, 1992. -228 с.

3. Александров, В.М. Аналитические методы в контактных задачах теории упругости / В.М. Александров, М.И. Чебаков. М.: Физматлит, 2004. - 304 с.

4. Александров, В.М. Неклассические пространственные задачи механики контактныхвзаимодействий упругих тел / В.М. Александров, Д.А. Пожарский. М.: Факториал, 1998.-288 с.

5. Бабешко, В.А. Высокочастотный резонанс массивного штампа / В.А. Бабешко //

6. Докл. АН СССР. 1989. - Т. 306, № 6. - С. 1328-1332.

7. Бабешко, В.А. Динамика неоднородных линейно-упругих сред / В.А. Бабешко, Е.В. Глушков, Ж.Ф. Зинченко. М.: Наука, 1989. - 343 с.

8. Бабешко, В.А. Интегральный и дифференциальный методы факторизации в задачахдля сплошных сред / В.А. Бабешко, О.М. Бабешко, О.В. Евдокимова // Тез. докл. IX Всерос. съезда по теоретической и прикладной механике. Н.Новгород, 2006. - С. 12.

9. Бабешко, В.А. К проблеме исследования динамических свойств трещиноподобныхтел / В.А. Бабешко // Доклады АН. 1989. - Т. 304, № 2. - С. 318-321.

10. Бабешко, В.А. К расчету параметров высокочастотного резонанса в трехмерномслучае / В.А. Бабешко // Докл. АН СССР. 1994. - Т. 335, № 1. - С. 55-58.

11. Бабешко, В.А. Метод фиктивного поглощения в задачах теории упругости для неоднородного полупространства / В.А. Бабешко, Т.И. Белянкова, В.В. Калинчук // ПММ. 2002. - Т. 66, вып. 2. - С. 276-284.

12. И.Бабешко, В.А. Метод фиктивного поглощения в плоских динамических задачах /

13. B.А. Бабешко, О.Д. Пряхина // ПММ. 1980. - Т. 44, вып. 3. - С. 477-484.

14. Бабешко, В.А. Метод фиктивного поглощения в связанных смешанных задачах теории упругости и математической физики для слоисто-неоднородного полупространства / В.А. Бабешко, В.В. Калинчук // ПММ. 2002. - Т. 66, вып. 2.1. C. 285-292.

15. Бабешко, В.А. Новый метод решений краевых задач механики сплошной среды и математической физики для неклассических областей / В.А. Бабешко // Доклады АН.- 1985.-Т. 284, № 1.-С. 73-76.

16. Бабешко, В.А. О методе факторизации в краевых задачах для сплошных сред / В.А. Бабешко, О.М. Бабеппсо // ДАН. 2004. - Т. 399, № 3. - С. 63-68.

17. Бабешко, В.А. О решении одного класса смешанных задач для слоистого полупространства / В.А. Бабешко, Т.И. Белянкова, В.В. Калинчук // ДАН. 2001. -Т. 380, №5. -С. 619-622.

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

19. Бабешко, В.А. Обобщенный метод факторизации в пространственных динамических смешанных задачах теории упругости / В.А. Бабешко. М.: Наука, 1984.-256 с.

20. Бабешко, В.А. Формулы факторизации некоторых мероморфных матриц-функций / В.А. Бабешко, О.М. Бабешко // ДАН. 2004. - Т. 399, № 1. - С. 163-167.

21. Бабешко, В.А. Явление высокочастотного резонанса в полуограниченных телах с неоднородностями / В.А. Бабешко, И.И. Ворович, И.Ф. Образцов // МТТ. 1990. -№ 3. - С. 74-83.

22. Баженов, В.Г. Метод граничных элементов в трехмерной динамической теории упругости и вязкоупругости с сопряженными полями / В.Г. Баженов, JI.A. Игумнов. Учебное пособие. Н.Новгород: Изд-во ННГУ, 2007. - 328 с.

23. Баженов, В.Г. Методы граничных интегральных уравнений и граничных элементов в решении задач трехмерной динамической теории упругости с сопряженными полями / В.Г. Баженов, JI.A. Игумнов. М.: ФИЗМАТЛИТ, 2008. - 352с.

24. Бахвалов, Н.С. Численные методы / Н.С. Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. -М.: БИНОМ. Лаборатория знаний, 2006. 636 с.

25. Белов, А.А. Гранично-элементный расчет динамики составных вязкоупругих тел / А.А. Белов // Проблемы прочности и пластичности: Межвуз. сборник. — Н.Новгород: Изд-во ННГУ. 2008. - Вып. 70. - С. 142-149.

26. Белов, А.А. Расчеты динамических задач теории упругости методом граничных элементов / А.А. Белов, Е.А. Шишкова // 11 Нижегородская сессия молодых ученых. Математические науки./ Н.Новгород: Издание ИП Гладкова О.В., 2006 г. — С. 56.

27. Бурчуладзе, Т.В. Развитие метода потенциала в теории упругости / Т.В. Бурчуладзе, Т.Г. Гегелия. Тбилиси: Мецниереба, 1985.-228 с.

28. Ватульян, А.О. Граничные интегральные уравнения для эллиптических операторов / А.О. Ватульян // Известия высших учебных заведений Северо-Кавказский Регион. 2000. -№3.- С. 34-37.

29. Ватульян, А.О. Новый вариант граничных интегральных уравнений и их применение к динамическим пространственным задачам теории упругости / А.О. Ватульян, В.М. Шамшин // ПММ. 1998. - Т. 62, вып. 3. - С. 112-119.

30. Ватульян, А.О. О граничных интегральных уравнениях 1-го рода в динамических задачах анизотропной теории упругости / А.О. Ватульян // ДАН РАН. 1993. - Т. 333, №3.-С. 312-314.

31. Ватульян, А.О. О новой формулировке граничных интегральных уравнений в задачах о колебаниях анизотропных тел / А.О. Ватульян, Е.В. Садчиков // Механика твердого тела. 1999. - № 2. - С. 78-84.

32. Ватульян, А.О. Обратные задачи в механике деформируемого твердого тела / А.О. Ватульян. М.: ФИЗМАТЛИТ, 2007. - 224 с.

33. Возбуждение упругих волн в слое пьезокерамическими накладками / Е.В. Глушков и др. // Акустический журнал. 2006. - Т. 52, № 4. - С. 470-479.

34. Ворович, Е.И. Некоторые динамические связанные задачи для электроупругого слоя / Е.И. Ворович, И.А. Зайцева, В.В. Калинчук // Современные проблемы механики сплошной среды: Тр. IV Междунар. конф'. Ростов-па-Дону: Изд-во СКНЦВШ, 1999.-Т. 1.-С. 107-110.

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

36. Ворович, И.И. Динамические смешанные задачи теории упругости для неклассических областей / И.И. Ворович, В.А. Бабешко. М.: Наука, 1979. - 320 с.

37. Ворович, И.И. К проблеме низкочастотных резонансов при взаимодействии упругого тела с полуограниченной средой / И.И. Ворович, Т.И. Белянкова, В.В. Калинчук // ДАН. 1998. - Т. 358, № 5. - С. 624-626.

38. Ворович, И.И. Неклассические смешанные задачи теории упругости / И.И. Ворович, В.М. Александров, В.А. Бабешко. М.: Наука, 1974. - 456 с.

39. Игумнов, Л.А Граничные интегральные уравнения трехмерных задач на плоских волнах / Л.А Игумнов // Докл. РАН. 2006. - Т. 409, №5. - С. 1-3.

40. Калинчук, В.В. Динамические контактные задачи для предварительно напряженных тел / В.В. Калинчук, Т.И. Белянкова. М.: Физматлит, 2002. - 240 с.

41. Калинчук, В.В. Динамические контактные задачи для предварительно напряженных электроупругих сред / В.В. Калинчук, Т.И. Белянкова. М.: Физматлит, 2006. - 272 с.

42. Контактные задачи теории упругости для неоднородных сред / С.М. Айзикович и др.. М.: Физматлит, 2006. - 237 с.

43. Линьков, A.M. Комплексный метод граничных интегральных уравнений теории упругости / A.M. Линьков. СПб.: Наука, 1999. - 382 с.

44. Литвинчук, С.Ю. Гранично-элементное моделирование нестационарных трехмерных динамических задач теории упругости и вязкоупругости: автореф. дис.канд. физ.-мат. наук: 01.02.04 / Литвинчук Светлана Юрьевна. Н. Новгород, 2006. - 24 с.

45. Манжиров, А.В. Контактные задачи для неоднородных стареющих вязкоупругих тел / А.В. Манжиров // Механика контактных взаимодействий. М.: Физматлит, 2001.-С. 549-565.

46. Механика контактных взаимодействий. М.: Физматлит, 2001. - 672 с.

47. Михлин, С.Г. Интегральные уравнения в теории упругости / С.Г.Михлин, Н.Ф.Морозов, Н.В.Паукшто. СПб., 1994. - 272 с.

48. Михлин, С.Г. Метод потенциалов в смешанной задаче для волнового уравнения / С.Г. Михлин. В.Д. Сапожникова // Изв. ВУЗов. Математика. 1977. - № 10. - С. 100-112.

49. Николаев, О.П. Разработка и применение модифицированной методики граничных элементов для трехмерных смешанных задач упругого равновесия: автореф. дис. канд. ф.-м. наук / Николаев Олег Петрович. — Горький, 1983.

50. Новацкий, В. Теория упругости / В. Новацкий. М.: Мир, 1975. - 872 с.

51. Нигул У.К. Волны в слоистых линейных средах / У.К. Нигул, А.С. Стулов // АН Эстонской ССР. Припринт. Таллин, 1985. - 66 с.

52. Пряхина, О. Д. Метод фиктивного поглощения в динамических задачах электроупругости / О.Д. Пряхина, А.В. Смирнова, О.М. Тукодова // ПММ. 1998. -Т. 62. Вып. 5. - С. 834-839.

53. Пряхина, О.Д. Связанная нестационарная задача о возбуждении электроупругого слоя массивным электродом / О.Д. Пряхина, М.Р. Фрейгейм // Изв. РАН, МТТ. -1998.-№2.-С. 111-118.

54. Ройтфарб, И.З. Численный метод решения пространственных динамических задач теории упругости на основе метода потенциала / И.З. Ройтфарб, Чу Вьет Кыонг // Сопротивление материалов и теория сооружений. Киев: Будивельник. - 1976. - С. 32-38.

55. Трехмерные задачи математической теории упругости и термоупругости / В.Д. Купрадзе и др.; ред. В.Д. Купрадзе. Изд. 2-е. М.: Наука, 1976. - 664 с.

56. Турилов, В.В. Расчет нестационарного динамического деформирования трехмерных упругих элементов конструкций методом гранично-временных элементов: автореф. дис. канд. тех. наук: 01.02.04 / Турилов Валерий Вячеславович. Горький, 1986. - 20 с.

57. Угодчиков, А.Г. Граничные интегро-дифференциальные уравнения нестационарных динамических задач теории упругости / А.Г. Угодчиков, Н.М. Хуторянский // Актуальные проблемы механики деформируемых сред. -Днепропетровск: ДГУ, 1979. С. 197-204.

58. Угодчиков, А.Г. Метод граничных элементов в механике деформируемого твердого тела / А.Г. Угодчиков, Н.М. Хуторянский. Казань: Изд-во Казанск. ун - та, 1986. — 295 с.

59. Хуторянский, Н.М. Метод гранично-временных элементов в пространственных задачах нестационарной динамики упругих и вязкоупругих тел: автореф. дисс.доктора техн. наук: 01.02.04 / Хуторянский Наум Маркович. Рига, 1988. -32 с.

60. Чудинович, И.Ю. Метод граничных уравнений в динамических задачах теории упругости / И.Ю. Чудинович. Харьков, 1990. - 121 с.

61. Aimi, A new space-time energetic formulation for wave propagation analysis in layered madia by BEMs / A. Aimi, M. Diligenti // Int. J. Numer. Meth. Engng. 2008.

62. Antes, H. The boundary integral Approach to static and dynamic contact problems / H. Antes, P.D. Panagiotopoulos // Int. Series of Numerical Mathematics 108. Birkhauser, Basel, 1992.-313 p.

63. Beer, G. The Boundary Element Method with Programming: For Engineers and Scientists / G. Beer, I. Smith, C. Duenser. Springer New York, 2008.

64. Beskos, D.E. Boundary Element Methods in Dynamic Analysis / D.E. Beskos // Applied Mechanics Review. 1987. - Vol. 40, № 1. - P. 1-23.

65. Beskos, D.E. Boundary element methods in dynamic analysis: Part II 1986-1996 / D.E. Beskos // Appl. Mech. Reviews. 1997. - Vol. 50. - P. 149 - 197.

66. Birgisson,B. Elastodynamic direct boundary element methods with Enhanced numerical stability properties / B. Birgisson, E. Siebrits, A.P. Peirce // International journal for numerical methods in engineering. 1999. - Vol. 46. - P. 871-888.

67. ЮЗ.Саггег, J.A.M. Alternative time-marching schemes for elastodynamics analysis with the domain boundary element method formulation / J.A.M. Carrer, W.J. Mansur // Comput. Mech. Springer-Verlag. - 2004. - Vol. 34. - P. 387-399.

68. Chudinovich, I. Boundary equations in dynamic problems of the theory of elasticity /1. Chudinovich // Acta Applicandae Mathematicae. 2001. - Vol. 65. - P. 169-183.

69. Coda, H.B. Three-Dimensional Transient BEM Analysis / H.B. Coda, W.S. Venturini // Computers & Structures. 1995. - Vol. 56(5). - P. 751-768.

70. Cruse, T.A. A direct formulation and numerical solution of the general transient elastodynamic problem. Part I. / T.A. Cruse, F.J. Rizzo // J. Math. Anal, Applic. 1968 -№ 22. - P. 244-259.

71. Cruse, T.A. A direct formulation and numerical solution of the general transient elastodynamic problem. Part II / T.A. Cruse, F.J. Rizzo // J. Math. Anal, Applic. 1968 -№22, 2.-P. 341-355.

72. De Hoop, A.Y. Representation theorems for the displacement in an elastic olid and their application to elastodynamic diffraction theory / A.Y. De Hoop // Delft: Tech. Hogeschoof, 1958, Dr. Sci. Thesis.

73. Dominguez, J. Boundary elements in Dynamics / J. Dominguez // Computational mechanics publication, Southampton, 1993.

74. Dominguez, J. Dynamic Stiffness of Rectangular Foundations / J. Dominguez // Report no. R78-20, Department of Civil Engineering, MIT, Cambridge MA, 1978.

75. Dominguez, J. Response of Embedded Foundations to Traveling Waves / J. Dominguez // Report no. R78-24, Department of Civil Engineering, MIT, Cambridge MA, 1978.

76. Durbin, F. Numerical inversion of Laplace transforms: an efficient improvement to Dubner and Abate's method / F. Durbin // The Computer Journal. 1974. - Vol.17, 4. -P. 371-376.

77. Frangi, A. "Causal" shape functions in the time domain boundary element method / A. Frangi // Comput. Mechanics / Springer-Verlag. 2000. - Vol. 25. - P. 533-541.

78. Freguency- and Time-domain BEM Analysis of Rigid Track on a Half-Space with Vibration Barriers / W. Hubert et al. // Mechanica. 2001. - Vol. 36. - P. 421-436.

79. Gaul, L. Boundary Element Methods for Engineers and Scientists / L. Gaul, M. Kogl, M. Wagner. Berlin Springer, 2003. - 488 p.

80. Gaul, L. Dynamic boundary element analysis of foundation slabs on layered soil / L. Gaul, P. Klein, M. Plenge // In Boundary elements X. (Brebbia C.A., ed.). Computational mechanics publications. Southampton, Boston, 1988. Vol. 4. - P. 29-44.

81. Guiggiani, M. A General Algorithm for Multidimensional Cauchy Principal Value Integrals in the Boundary Element method / M. Guiggiani, A. Gigante // Journal of Applied Mechanics, ASME. 1990. -№ 57. - P. 906-915.

82. Ha-Duong, T. On Retarded Potential Boundary Integral Equations and their Discretisation / T. Ha-Duong // In: Topics in Computational Wave propagation ( Eds. M. Ainsworth, P. Davies et at. Springer-Verlag, Berlin, 2003. - P. 301-336.

83. Junior, D.S. A time-marching scheme based on implicit Green's functions for elastodynamic analysis with the domain boundary element method / D.S. Junior // Comput. Mech. 2007. - № 40. - P. 827-835.

84. Lachat, J.C. Effective numerical treatment of boundary integral equations: a formulation fot three-dimensional elastostatics / J.C. Lachat, J.O. Watson // Int. J. Numer. Mech. Eng. 1976. -№ 10.-P. 991-1005.

85. Lubich, C. Convolution Quadrature and Discretized Operational Calculus. I. / C. Lubich //Numerische Mathematik. 1988. -№ 52. - P. 129-145.

86. Lubich, C. Convolution Quadrature and Discretized Operational Calculus. II. / C. Lubich // Numerische Mathematik. 1988. -№ 52. - P. 413-142.

87. Lubich, C. On the multistep time discretization of linear initial-boundary value problems and their boundary integral equation / C. Lubich // Numer. Math. 1994. - № 67. - P. 365-389.

88. Lubich, C. Time discretization of parabolic boundary integral equations / C. Lubich, R. Schneider // Numer. Math. 1992. - № 63. - P. 455-481.

89. Manolis, G.D. A comparative study on three boundary element method approaches to problem in elastodynamics / G.D. Manolis // Int. J. Num. Meth. Eng. 1983. -Vol. 19, 1. -P. 73-91.

90. Manolis, G.D. Dynamic stress concentration studies by the boundary integrals and Laplace transform / G.D. Manolis, D.E. Beskos // Int. J. Num. Meth. Eng. 1981.- Vol. 17, 4.-P. 573-599.

91. Mansur, W.J, Time discontinuous linear traction approximation in time-domain BEM scalar wave propagation / W.J. Mansur, J.A.M. Carrer, E.F.N. Siqueira // International journal for numerical methods in engineering. 1998. Vol.42, № 4. - P. 667-683.

92. Mansur, W.J. A Time-Stepping Technique to Solve Wave Propagation Problems Using the Boundary Element Method / W.J. Mansur // Phd thesis, University of Southampton, 1983.

93. Mansur, W.J. Transient Elastodynamics Using a Time-Stepping Technique / W.J. Mansur, C.A. Brebbia // In Boundary Elements. Berlin, Springer-Verla, 1983. - P. 677698.

94. Muntz, Ch.H. Sur la resolution du probleme dynamique de l'elasticite / Ch.H. Muntz // C. r. Acad. Sci.- 1932.-Vol. 194, № 17.-P. 1456-1459.

95. Narayanan, G.V. Numerical operational methods for time-dependent linear problems / G.V. Narayanan, D.E. Beskos // Int. J. Num. Meth. Eng. 1982. - Vol.18 (12). - P. 1829-1854.

96. Nardini, D. A New Approach to Free Vibration Analysis Using Boundary Elements / D. Nardini, C.A. Brebbia // In Boundary Element Methods. (Brebbia C.A., Ed.). Springer-Verlag: Berlin, 1982. - P.312-326.

97. Numerical computation of internal stress and velocity in time-domain BEM formulation for elastodynamics / Jr. D Soares et al. // Computational Mechanics / Springer-Verlag. -2002.-Vol. 30.-P. 38-47.

98. Peirce, A. Stability Analysis and Design of Time-Stepping Schemes for General Elastodynamic Boundary Element Models / A. Peirce, E. Siebrits //. International Journal for Numerical Methods in Engineering. 1997. - Vol. 40(2). - P. 319-342.

99. Polyanin, A.D. Handbook of Integral Equations, Second Edition / A.D. Polyanin, A.V. Manzhirov // Boca Raton: Chapman&Hall / CRC. Precss, Boca Raton. 2008. - 1144 p.

100. Rizos, D.C. An Advanced Direct Time Domain BEM Formulation for General 3-D Elastodynamic Problems / D.C. Rizos, D.L. Karabalis // Computational Mechanics. -1994. -№ 15.-P. 249-269.

101. Schanz, M. A new visco- and elastodynamics time domain boundary element formulation / M. Schanz, H. Antes // Computational mechanics. 1997. - Vol. 20, № 5. -P. 452-459.

102. Schanz, M. Application of 'Operational quadrature methods' in time domain boundary element methods / M. Schanz, H. Antes // Mechanica. 1997. - Vol. 32, № 3. - P. 179186.

103. Schanz, M. Boundary element Analysis / M. Shanz, O. Steinbach // Lecture Notes in Applied and Computational Mechanics. Springer-Berlin. - Vol.29. - 2007. - 352 p.

104. Schanz, M. Implementation of viscoelastic behaviour in a time domain boundary element formulation / M. Schanz, L. Gaul // Applied mechanics review. 1993. - Vol.46, № 11 (part 2).-P. S41-S46.

105. Schanz, M. Numerical damping and instability of a 3-D BEM time-stepping algorithm / M. Schanz, L. Gaul, H. Antes // Extended abstracts of IABEM 93, Braunschweig, 1993.

106. Schanz, M. Wave Propogation in Viscoelastic and Poroelastic Continua / M. Schanz // Berlin Springer, 2001.- 170 p.

107. Schclar, N.A. Anisotropic analysis using boundary elements / N.A. Schclar. Southampton, Boston: Computational Mechanics Publications. 1994. - 152 p.

108. Soares, Jr. D. An efficient stabilized boundary element formulation for 2D time-domain acoustics and elastodynamics / Jr. D. Soares, W.J. Mansur // Comput. Mech. Springer-Verlag. - 2007. - Vol. 40. - P. 355-365.

109. Wave Propagation in Infinite Domains // Springer-Verlag (D. Haank, M.Mos, P. Hendriks eds) 2007. - Vol. 31.

110. Yu, G. A Linear 9 for 2-D Elastodynamic BE Analysis / G. Yu, W.J. Mansur, J.A.M. Carrer // Computational Mechanics. 1999. - Vol. 2A. - P. 82-89.

111. Yu, G. Stability of Galerkin and collocation time domain boundary element methods as applied to the scalar wave equation / G. Yu, W.J. Mansur, J.A.M. Carrer, L. Gong // Computers and Structures. 2000. Vol. 74, № 4. - P. 495-506.

112. Yu, G. Time weighting in time domain BEM / G. Yu, W.J. Mansur, J.A.M. Carrer, L. Gong // Engineering analysis with boundary elements. 1998. Vol. 22, № 3. - P. 175181.

113. Zhao X. An efficient approach for the numerical inversion of Laplace transform and its application in dynamic fracture analysis of a piezoelectric laminate. // Int.J. of Solids and Structures. 2004. V. 41. P. 3653-3674.

114. Zhong, M. The analysis of dynamic stress intensity factor for semi-circular surface crack using time-domain BEM formulation/ M. Zhong, Y. Zhang // Appl. Mathem. and Mechan. Engl. Edition. - 2001. - Vol. 22, № 11. - P. 1344-1351.