Главная » Просмотр файлов » Галлагер - Метод конечных элементов. Основы

Галлагер - Метод конечных элементов. Основы (947497), страница 75

Файл №947497 Галлагер - Метод конечных элементов. Основы (Галлагер - Метод конечных элементов. Основы) 75 страницаГаллагер - Метод конечных элементов. Основы (947497) страница 752013-09-15СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

е. работа этих сил у' на любых кинематически возможных перемещениях зависит только от начальной и конечной конфигураций системы. Тем самым исключаются случаи, когда направление действия силы «отслеживает» направление отклоненного элемента системы, на который эта сила действует. 13. Анализ устойчивости упругих тал Для условий равновесия в невозмущенном состоянии, т, е. когда осевая сила меньше критической нагрузки, применение принг[ипа стас(ионарнослти потпгнциальной энергии в виде равенства нулю первой вариации от П (т. е. 6П„=О) приводит к уравнению жесткости (13.25) При этом (Лу) можно определить пз уравнения (!3.25), используя обычные средства, причем мы придем к результатам, учитывающим взаимодействие осевого и поперечного деформнрования, влияние осевых нагрузок на жесткостные характеристики балки.

Следует отметить, что таким образом можно учесть эффект увеличения изгпбной жесткости при действии растягивающих нагрузок. Чтобы рассмотреть вопрос потери упругой устойчивости, когда интенсивность системы осевых нагрузок, вызывающих выпучивание, еще не известна, инкремептальная матрица жесткости должна быть вначале подсчитана численно при произвольно выбранной интенсивности нагрузки (предполагается, что распределение осевых спл фиксировано). При выпучивании считаем, что интенсивность системы осевых нагрузок в ьт раз больше произвольно выбранной интенсивности сил (Р,), использованных при построении матрицы [[га), поэтому уравнение равновесия принимает вид 6Пг — — [Ку) (Лу )+ от[[та) [Ьу ) — (Р ) =О.

(13.26) В этом месте необходимо рассмотреть условия нейтрального равновесия; онн определяют собственное значение ы и отвечающую ему моду выпучивания (тх„). Необходимую информацию нельзя получить, рассматривая первую вариацию функционала потенпнальной энергии, поэтому следует выписать вторую вариацию 6'П =6(6П ). Йз рис. 6.2, на котором схематически изображена зависимость потенциальной энергии от некоторого представительного параметра перемещения А, следует, что в состоянии устойчивого равновесия бтП )О, а для нейтрального равновесия 6тП =О.

Последнее условие определяет точку бифуркации искомого решения. Теперь, применив эти условия к (13.24), запишем вторую вариацию потенциальной энергии в виде 6'П„= [ 6Ау ) 1 К) (66 у ) =О, и так как в данном случае!К)=[Ку[+!Кх!, то [[К,)+[К [[=О, (!3.27) где символом [ [ обозначен детерминант. Таким образом, имеем условие, согласно которому детерминант матрицы (13.27) равен нулю. Альтернативой этому условию является то, что не существует единственного решения (условие бифур- 13.2, Глобальиая формулироаиа 40! кации) уравнения [13.26), если найдется вектор (Лу) и скаляр оу, такие, что [Ку+оуКа!(Ау)=0. Вычисление детерминанта с целью определения значения оу для системы большого порядка неэффективно. Обычно используют модификацию уравнения [13.27), умножая его на (Ь).

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

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

[ Л ) =) уч ! й ). Дла матРицы [ку), напРимеР, имеем [В действительности, возможно любое разбиение общей совокупности степеней свободы, однако для простоты рассмотрения и обозначения оставим деление на степени свободы (уа) и (8). Согласно равд. 2.8, преобразование от [ уу [ к степеням свободы [ в;й ), основанное на традиционной матрице жесткости [йу), имеет вид ~ (~) = (1'а) (уу) (13.28) [1у! Г 1 '[ й1ое 1о Матрица преобразования [Га! используется теперь для преобразования обычной матрицы жесткости [йу) и геометрической матрицы 14 № 1647 13. Анализ устойчивости упругих тел 402 жесткости Щ. Имеем !й,]-]- ]й,](~)=о, (13.29) где (йг]=(Го]~!(су](Го] ° (]ту]=!Го] !(сх]!Га].

(13.29а, Ь) Собственное значение «т н соответствующая форма моды выпучива- ния выделяются из редуцированной системы (!3.29). Вопросы эф- фективности этой схемы будут описаны в дальнейшем. 13.3. Призматический элемент 13.3.1. Выпуинваниа прн изгиба Целью настоящих рассмотрений является преобразование энергии деформации (13.13) в явную формулировку для матрицы жесткости элемента посредством выбора функционального представления для ш. Точное представление для этого случая можно сформулировать, используя функцию перемещений, которая удовлетворяет соответствующему дифференциальному уравнению !!3.1].

Однако, стремясь, как обычно, при конечно-элементном анализе выбрать простое аппроксимирующее поле перемещений, альтернативно выберем (при $=х/Е) поле перемещений для изгибаемого элемента без осевой нагрузки (см. (5.14а)). Имеем нг= ] (1 — Зсот+2$о) (3$т — 2Р) (1 — 2~+во) х Я вЂ” ~о) х ] т т Использование выписанного выражения в (13.13) приводит к хорошо известному представлению для базисной изгибной матрицы Щ, задаваемой выражением (5,17), а также к следующей явной форме для матрицы (]44]: нг, та 0 0 36: (Симметрично) л„— 36 36 1 41 зи. Зй Зй 4Хт:.

', — Зй ЗŠ— Е,' 4Е' (13.30) В качестве примера применения инкрементальной матрицы жесткости при решении задач о потере устойчивости балок найдем критическую силу для изображенной на рис. !ЗА простой балки, используя один элемент. Здесь нгт=йт=О, Р„= — Р,, (.=тт2.

По- ааз 13.3. Прнаматнчесинй ваемент этому характеристическое уравнение примет вид решая которое получаем Р„=9.9чЕ111а. Зто решение лишь на 0.752ай отличается от точного решения иа(Е11!а1. Йвалалелав Р . 1ЗА. На рис. 13.5 представлена ошибка в процентах при решении этой задачи как функция от числа элементов разбиения. Для сравнения приведены соответствующие результаты для конечно-разностного решения 113.21. Пониженная точность конечно-разностного решения вызвана допущением линейного изменения тв между узлами по сравнению с принимаемым здесь кубическим характером и 1 О 2 4 а З 1О Числа ааеаеапта Рнс.

1З.з. Характериспчки сходиности методов — выпучиваиие колонны Эйлера. 14* 13. Анализ устойчивости упругих тел изменения. Заметим, однако, что конечно-разностные уравнения используют только одну степень свободы в каждом узле, а именно прогиб 2в. На рис. 13.6 изображен рафик зависимости ошибки в процентах от числа элементов в представлении при определении критической нагрузки суживающейся консольной балки. Сплошной линией Ь =и, (-,)" пт 1— Се пенне А-А Рппппнып нпннпнйна пнныйпнпнепне-, в ЯЯ1- - — — -и нне унена пнепгеннтпп Π— 2 и -6 чи -8 -~о — 16 Рис. 13.6. Характеристики сколимости методов — устойчивость суитавииейси балки изображен график для ступенчатого представления 113.31, в котором для задания элемента используются геометрические характеристики в центре элемента.

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

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

Тип файла
DJVU-файл
Размер
3,53 Mb
Тип материала
Учебное заведение
Неизвестно

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

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