Лекция 15 (1014401), страница 2

Файл №1014401 Лекция 15 (Материалы к лекциям) 2 страницаЛекция 15 (1014401) страница 22017-06-17СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

Покажем теперь решение подобных задач методом конечных элементов (вариационным подходом).

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

Установим зависимость между перемещениями узлов треугольного элемента и действующими на узлы силами (см. рисунок).

Так как максимальная степень (по производным) уравнений (6.38) равна 2 (2m = 2), то минимальная степень интерполирующего полинома для этой задачи будет равна n = 2m –1 = 1. Поэтому примем, что перемещения любой точки внутри треугольного элемента зависят от координат точки линейно, т.е.:

Применяя первую формулу из (6.46) к узлам получим:

Отсюда определяем коэффициенты , выражая их через координаты узлов. Аналогично, записав вторую формулу из (6.46) для перемещений v узлов, выразим через координаты узлов. Подставив найденные значения в формулы (6.46), получим:

Где , а остальные коэффициенты получаются круговой перестановкой индексов i, j, m. Величина f представляет собой площадь треугольника, т.е.

Где означает определитель матрицы.

Имея формулы (6.47) для перемещений, можно по формулам Коши (6.44) найти деформации в точке пластины:

Где , а вектор – вектор перемещений узлов, т.е. .

Напряжения выражаются формулами (6.40) или в матричной записи:

.

Более компактно это запишется так: , где - матрица упругости (жесткости) элемента.

Вектор нагрузок:

и вектор перемещений узлов:

связаны известным уравнением .

Для установления интересующей нас зависимости, воспользуемся обозначением из формулы (6.42), полученной вариационным методом.

Так как , где – толщина пластины, – элемент площади, то считая постоянным и интегрируя, найдем для матрицы жесткости элемента формулу:

Где – матрица геометрии, – толщина пластины.

Матрицу жесткости можно записать в виде:

Если произвести перемножение матриц в формуле (6.48), то для коэффициентов получим:

Напряжения в пластине будут:



Построение матрицы жесткости системы (пластины)



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

Для решения задачи методом конечных элементов разбиваем пластину на треугольные панели – конечные элементы и заменяем распределённые нагрузки сосредоточенными в узлах силами в виде компонентов по осям xи y.

Например, для i- го узла, лежащего на границе, такими компонентами будут силы Pxi и Pyi. Цифры 1,2,3 и т.д – номера узлов, при этом 1,2,3,5 – узлы, лежащие на внешней границы пластины, а узел 4 – внутренний, к которому внешние нагрузки не приложены.

В дальнейшем будем рассматривать узел iи примыкающие к нему конечные элементы Л, П и Н, обозначающие соответственно левый, правый и нижний конечный элемент.

В i – ом узле конечных элементов Л, П и Н действуют силы . Эти же силы будут действовать в обратном направлении и на узел i, вырезанный из пластины, но на него также действуют и внешние силы . Из условия равновесия узла получим:

Такие условия равновесия можно составить для всех узлов пластины.

Условия равновесия в общей матрице будут:

, где - система сил, действующих в узлах конечного элемента I, записанная в столбце, содержащих 2n элементов; - система сил, действующих в узлах конечного элемента II и т.д.

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

Подставляя в условия равновесия в общей матрице получаем:

Где представляет суммарную матрицу жесткости пластины.

Таким образом, образовав матрицу жесткости пластины путём занесения в неё матриц жесткости треугольных элементов, получим систему уравнений (6.50), которая в развернутой записи имеет вид:

Чтобы исключить перемещения пластины как твердого тела, необходимо ввести условия закрепления пластины; если пластина закреплена в каких-то точках, то соответствующие перемещения нужно принять равными нулю и исключить соответствующие этим перемещениям строки и столбцы; если пластина закреплена по всей стороне, то следует приравнять нулю перемещения всех узлов этой стороны. Например, если пластина закреплена по левой стороне, т.е. закреплены узлы 1,3,5 и т.д., то нужно принять и исключить в системе уравнений строки и столбцы 1,2,6,7,10,11, 2n-1, 2n.

Из решения уменьшенной системы уравнений найдутся перемещения всех узлов, а по ним определяются узловые силы .

Напряжения в пластине будут:





МКЭ и МКР



Завершая рассмотрение МКЭ приведем его сравнение с МКР:

Метод конечных разностей

МКЭ

Приближенный метод решения дифференциальных уравнений.

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

Включает решение больших разреженных систем уравнений.

Включает решение больших разрешенных систем уравнений.

Постановка довольно проста, содержит мало выводов и выдач.

Постановка может быть достаточно сложной, с многочисленными вводами и выводами.

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

Метод очень физичен. Структура близка к действительной физической задаче.

Для решения каждой задачи составляется своя программа.

Позволяет создать универсальные программы.

Чтобы получить решение, достаточно описать задачу дифференциальным уравнением.

Требует достаточно высокой квалификации – опыта.

Обычно используются правильные сетки. Локальное сгущение сетки затруднительно. Трудно выбрать экономическую сетку.

Весьма велик выбор типов, форм и размеров элементов. Сгущение сетки не вызывает проблем.

Представление граничных условий связано с трудностями.

Представление граничных условий не вызывает затруднений.

Плохо приспособлен для решения задач с неоднородным и анизотропным распределением параметров.

Позволяет решать задачи с неоднородным и анизотропным распределением параметров.

Каждое решение применимо лишь для конкретной задачи - решения не допускают обобщения.

Каждое решение применимо лишь для конкретной задачи - решения не допускают обобщения.

Приложение. Примеры использования вариационного подхода



Пример 1



Задача с кручением призматического стержня произвольного сечения сводится к решению уравнения в частных производных:

В области Ω при нулевом значении функции на контуре поперечного сечения.

Известно, что решение этой задачи эквивалентно отысканию функции , удовлетворяющей на контуре краевым условиям, и минимизирующей функционал будет:

При использовании МКЭ поперечное сечение стержня разбивается на совокупность конечных элементов, для каждого из которых составляется выражение, аппроксимирующее искомую функцию (x,y):

Подставляя последнее выражение в Э, интегрируя по площади e - го конечного элемента(e ) и дифференцируя полученное выражениеЭ((e )) по каждой из узловых неизвестных, получаем:

Или

Где:



Пример 2



Несколько иной подход – метод Бубнова – Галеркина.

Если записать общее выражение для функционала Э(U) затруднительно, а дифференциальное уравнение рассматриваемой задачи известно, то для построения матриц и можно воспользоваться процедурой метода Бубнова - Галеркина:

От функцииU(x,y,z), входящей в подынтегральное выражение, требуется непрерывность во всей области  вплоть до производных (2m-1) – го порядка.

Напомним, что при использовании вариационного способа для получения уравнений МКЭ достаточно было обеспечить непрерывность функции и ее производных доmго порядка.

18

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

Тип файла
Документ
Размер
198,09 Kb
Тип материала
Высшее учебное заведение

Список файлов лекций

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