Главная » Просмотр файлов » Брэббия К., Теллес Ж., Вроубел Л. - Методы граничных элементов

Брэббия К., Теллес Ж., Вроубел Л. - Методы граничных элементов (1050609), страница 62

Файл №1050609 Брэббия К., Теллес Ж., Вроубел Л. - Методы граничных элементов (Брэббия К., Теллес Ж., Вроубел Л. - Методы граничных элементов) 62 страницаБрэббия К., Теллес Ж., Вроубел Л. - Методы граничных элементов (1050609) страница 622017-12-27СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 62)

дИдп. Интегрирование по времени в уравнении (1ОА2) можно выполнить аналитически следующим образом: ср ср ~ д(х, Пи»(8, х, (р, !) 81 = — ( 4(х, 1) Л(Г, Ги)д! = — 4(х, !и), 1 (! ОА5) ср ) и (х, 1) Ч' ($, х, Гр, 1) д(=- — 4 [ — ) и (х, !) Л (1 (л) д(— с. ср — — )е и(х, !) дг =- — — 'с — и(х, Ги) + р) дС ) чяе(» -+ — ' "( ' ) ] ~. (10.461 В последнем соотношении интеграл брался по частям. Подставляя соотношения (10.45) и (!ОА6) в уравнение (!0.42), получаем граничное интегральное уравнение для трехмерного случая 1 1 Г1 с(ч) и($, гр) = — ( — ~д (х, Г„)+ п, [ — и(х, 1,) + 1 + — '['"' " ~ '.~дг(.)+ — '[!Л',—, — '(ГМ,)[, ((ОА7) с=си П где Мо и Л', — средние значения соответственно и, и (дигдг 1, на сферической поверхности с центром в точке $ и с переменным радиусом сг 110).

'р(ожно видеть, что коэффициент с (ч) имеет то же значение, что н в сл>чае уравнения Лапласа 18, 9). Уравнение Задачи а р гаришранаиии вони ЗЗЗ Глана га 382 10.8. Двумерные задачи (10.47) можно рассматривать как одно пэ представлений хорошо известно~о интегрального уравнения Кирхгофа (1 — 5).

Процедуры, необходимые для численного решения этого уравнения, подробно обсуждаются в книгах !8, 9). Характерной особенностью исследования трехмерного случая является то, что здесь не требуется интегрирования по времени в отличие от двумерного случаи, который обсуждается в следующем разделе. ге~ вере мав время Рва. 10.18. Прввнвденныи ногевпнвн а точечного игточнииа, помешечного в вшви в форме пврадаевеанаене. Пример 10.5. Гррненбум опубликовал 18! некоторые предварительныс результаты по созданию вычислительной программы, основанной на описанных здесь методах. В качестве специального примера для трехмерного случая было рассмотрено распространение волн от точечного источника, помещенного в ящик в форме параллелепипеда.

Эта задача была ныбраиа потому, что здесь можно получить аналитическое решение, используя метод отображений. Размеры ящика составляли 3 х 4 х 5 (в безразмерных единицах), а его понерхность была разбита на 376 постоянных прямоугольных граничных элементов1 размером 0,5 х 0,5, В качестве функции источника брался набор коротких импульсов с гауссовским распределением Поскольку в уравнении (10.4?) потенциал и выражается через значения величин, найденных на предыдущем шаге по времени, при числевиом решении применялась шаговая по времени процедура. На рнс. 10.16 представлены результаты решения для граничной точки, соответствующие двум шагам по времени (шаг Ы брался равным 0,5 и 2), и дано сопоставление с точным решением.

При болыпем шаге по времени решение для потенциала было не очень точным для столь быстро изменяющегося процесса, поэтому имеются небольшие расхождения с точным решением. Для М = =- 0,5 приближенное решение довольно хорошо описывает картину второго и третьего отражения волны, но затем процесс численного решения становится неустойчивым. Грененбум !8 ! делает вывод, что для столь грубой сетки полученные результаты являются ободряющими. Фундаввентальиое решение двумерного скзлярного яолиового уравнения для нестацпонарных процессов имеет внд (4, 4? ! и* = в,, Н !р (Ги — !) — г), (10.48) 2в [ае (Г . — 1)а — гн! ПВ где Н вЂ” функция Хевисайда.

Нормальная производная функции и' на границе Г равна (9, 10! аи* вг 4*= — =',, „, Н1 (Гн....ф !— аи ~ 2п(ав(г г)а,д)вга — е вав Ь(р((г.- !) г)~пе 110.49) 2п Гае (Г - Г)е .е)врв Поскольку фундаментальное решение для двумерной задачи выражается через функцию Хевисайда (ступенчатую функцию), а не через дельта-функцию, как это имеет место в трехмерной задаче, влияние функции источника, помещенного в точке х, на потенциал в точке (5; 1„) уже не ограничивается временем запаздывания, и интегрирование следует проводить от начального момента времени !о до текущего момента времени !и. Поэтому при численном решении соответствующего граничного интегрального уравнения следует использовать шаговые по времени схемы решения, описанные в гл.

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

Глава [) Колебания 11.1. Введение В то время как в ряде инженерных наук метод граничных элементов быстро завоевывал позиции, исследованию этим методом динамических задач теории упругости было посвящено лишь несколько статей. Впервые формулировка и решение динамической задачи теории упругости для неустановившегося процесса с помощью метода граничных элементов и преобразования Лапласа были представлены в работах Крузе и Риццо [4), а также Крузе [5), в которых была решена задача о распространении волн в полуплоскости. В качестве приложения была рассмотрена задача о распространении одиночной антиплоской волны. В 1978 г. Коле и др.

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

Одной из главных проблем в подобном подходе является способ численного обратного преобразования Лапласа, хотя ряд соображений по этому поводу приводится в работе [8). Динамическая жесткость фундаментов, лежащих на полупространстве прн установившемся движении, недавно исследовалась в работах Домингеца и др. [9, 10). В самое последнее время было опубликовано несколько серьезных исследований динамических задач теории упругости, в том числе работы Хатчинсона [11), Нива и др. [12], Нардини [26), где методом граничных элементов находились собственные значения и собственные векторы.

11.2. Разрешающие уравнения Представленные здесь уравнения соответствуют теории малых перемещений для однородных изотропных материалов с линейно- упругим поведением. Кроме того, здесь используются те же, что и в гл. 5, предложения о регулярности на бесконечности. Здесь рассматривается область Я с границей Г (рис. 1!.1). Введем следующие обозначения для компонент перемещения точки х с координатами х;: и; (х, 1). (11. Ц Заь Конебання При указанных предположениях уравнение движения в пере- мещениях можно записать в виде П вЂ” 3) (С[ — Сз) ип ц+ Сазе, и+ Ь)!р = йи (11.2) где, как и прежде, запятые обозначают пространственные произ- водные, а точки — производные по времени. Скорости распростра- нения продольных волн (воли сжатия) и поперечных (сдвнговых) волн соответственно равны г С, = у'()е-[-2р,)/р, Сн = 1/р!р, й (11.3) где Х и и — постоянные Ламэ: Л = Етl [(1 + т) (1 — 2р) ), я и = б = Е/ [2 (1 -[- т) ), (11.4) р — плотность материала.

Деформации равны х1 вр =* — (и~ + и ) (11.5) (11.8) поведение материала описывается законом Гука он = 2рзц+ )ез„нбц. (11.6) Основными граничными и начальными условиями, которые должны удовлетворяться при решении динамических вадач теории упругости, являются: начальные условия — при 1 = 1, в области 11 -[- Г: и;(х, 1) = и";(х), й~(х, 1) = о~~(х), (11.7) где и) (х) и оер (х) — заданные функции; граничные условия для перемещений — при ! ) 1б на участке Г| границы Г: и;(х, 1) = й~ (х, 1); граничные условия для напряжений — при ! > 1б на участке Г, границы Г: р~(х, 1) = онпз — — р,(х, 1), (1 1.9) где черточкой обозначены заданные величины; Г = Г, -[- Г,.

гя Компоненты внутренних напряжений и напряжений на границе можно представить как функции перемещений: оц = р, (ипе+ иве)+),ин,ьбц. (11.10) Р; = )е(ди;/дп )-из,п!) -[ )илгпп (11.!1) где и, — направляющие косинусы внешней нормали к границе (рис. 11.1). 13 Вреббнн К. н Нр. Гяаяа ~ дояебаааа [2(Р)' — (г/Ся)Р]0(Г, г/С,) ) г,рг,! гр [(Г)2 ( /С )2]п2 (11.19) 1ЗР В отсутствие объемных сил уравнение (11.2) можно записать [3] в следующем виде: '[г е= (1/С1) е, (11.12) 7 га = (1/Се) гзр (11.13) где е — объемное расширение (е =- з„ь), е — вектор вращения, который связан с тензором вращения (см. выражения (5.25)) соотношениями ыг = гаазр газ = гз31р ыз = гз12 (1! .14) Из уравнений (11.12) и (11.13) следует, что скорость С, соответствует волнам сжатия или Р-волнам, а скорость С, соответствует поперечным или сдвиговым волнам, называемым также Я-волнами.

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

Характеристики

Список файлов книги

Свежие статьи
Популярно сейчас
А знаете ли Вы, что из года в год задания практически не меняются? Математика, преподаваемая в учебных заведениях, никак не менялась минимум 30 лет. Найдите нужный учебный материал на СтудИзбе!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6367
Авторов
на СтудИзбе
309
Средний доход
с одного платного файла
Обучение Подробнее