Содержание
- 2. Метод конечных элементов (МКЭ) - основной метод современной строительной механики, лежащий в основе подавляющего большинства современных
- 3. Суть метода заключается в том, что область (одно- , двух- или трехмерная), занимаемая конструкцией, разбивается на
- 4. Виды МКЭ По способу получения основных, т. е. разрешающих, уравнений различают четыре основных вида метода конечных
- 5. Вариационный метод основан на принципах стационарности некоторой переменной, зависящей от одной или нескольких функций (такая переменная
- 6. Метод невязок представляет собой наиболее общий подход к построению основных соотношений МКЭ. Этот метод целесообразно применять
- 7. Метод энергетического баланса (метод Одена) основан на балансе различных видов энергии, записанном в интегральной форме. Этот
- 8. Формы МКЭ В МКЭ, аналогично классическим методам строительной механики, за основные неизвестные могут приниматься величины разного
- 9. Метод перемещений – в настоящее время наиболее распространенная форма МКЭ. Это объясняется тем, что для заданной
- 10. Принцип минимума дополнительной энергии и связанные с ним схемы МКЭ в форме метода сил, а также
- 11. Аппроксимация МКЭ относится к методам дискретного анализа. Однако в отличие от численных методов, основывающихся на математической
- 12. Аппроксимирующие функции Аппроксимация, как правило, дает приближенное, а не точное, описание действительного распределения искомых величин в
- 13. Выделенный элемент конструкции имеет те же свойства, что и конструкция в месте расположения элемента. Благодаря такому
- 14. Для каждого элемента можно выбрать оси, в которых можно записывать геометрические или статические условия. Эти оси
- 15. Для системы в целом задаются произвольно выбранные оси, которые называются глобальными осями. На рисунке это оси
- 16. Задачи статики, решаемые с помощью МКЭ, сводятся к следующим основным этапам: 1. Разбиение системы на конечные
- 17. 4. Построение матрицы жесткости конечного элемента и вектора узловых нагрузок. С помощью принципа Лагранжа на основе
- 18. 7. Решение системы алгебраических уравнений. Для решения системы линейных алгебраических уравнений (СЛАУ) используются как точные, так
- 19. РАЗБИЕНИЕ СИСТЕМЫ НА КОНЕЧНЫЕ ЭЛЕМЕНТЫ Данный этап включает в себя разбиение системы на конечные элементы, выбор
- 20. РАЗБИЕНИЕ СИСТЕМЫ НА КОНЕЧНЫЕ ЭЛЕМЕНТЫ Важным моментом в процессе решения задачи по МКЭ является нумерация узлов
- 21. ПОСТРОЕНИЕ ИНТЕРПОЛИРУЮЩИХ ФУНКЦИЙ МКЭ основан на аппроксимации непрерывной функции, определенной на всей области, дискретной моделью с
- 22. ПОСТРОЕНИЕ ИНТЕРПОЛИРУЮЩИХ ФУНКЦИЙ Критерий совместности. Интерполирующая функция должна быть непрерывна вместе со своими производными до (p-l)-гo
- 23. ВЫВОД ОСНОВНЫХ ГЕОМЕТРИЧЕСКИХ И ФИЗИЧЕСКИХ СООТНОШЕНИЙ На основе выбранной системы интерполирующих функций выводятся зависимости между деформациями
- 24. ВЫВОД ОСНОВНЫХ ГЕОМЕТРИЧЕСКИХ И ФИЗИЧЕСКИХ СООТНОШЕНИЙ Физические соотношения, определяющие зависимость между напряжениями и деформациями, имеют вид
- 25. ПОСТРОЕНИЕ МАТРИЦЫ ЖЕСТКОСТИ КОНЕЧНОГО ЭЛЕМЕНТА И ВЕКТОРА УЗЛОВЫХ НАГРУЗОК Полная потенциальная энергия упругой системы определяется по
- 26. ПОСТРОЕНИЕ МАТРИЦЫ ЖЕСТКОСТИ КОНЕЧНОГО ЭЛЕМЕНТА И ВЕКТОРА УЗЛОВЫХ НАГРУЗОК Это условие записывается в виде: Значение полной
- 27. ПОСТРОЕНИЕ МАТРИЦЫ ЖЕСТКОСТИ КОНЕЧНОГО ЭЛЕМЕНТА И ВЕКТОРА УЗЛОВЫХ НАГРУЗОК где матрица есть вектор узловых нагрузок. представляет
- 28. ПЕРЕХОД ОТ ЛОКАЛЬНЫХ КООРДИНАТ К ГЛОБАЛЬНЫМ КООРДИНАТАМ СИСТЕМЫ В процессе решения задачи по МКЭ в большинстве
- 29. ПЕРЕХОД ОТ ЛОКАЛЬНЫХ КООРДИНАТ К ГЛОБАЛЬНЫМ КООРДИНАТАМ СИСТЕМЫ В матричной форме зависимость между перемещениями узлов отдельного
- 30. ПЕРЕХОД ОТ ЛОКАЛЬНЫХ КООРДИНАТ К ГЛОБАЛЬНЫМ КООРДИНАТАМ СИСТЕМЫ Значения sinα и cosα вычисляются через координаты узлов
- 31. ПЕРЕХОД ОТ ЛОКАЛЬНЫХ КООРДИНАТ К ГЛОБАЛЬНЫМ КООРДИНАТАМ СИСТЕМЫ Для получения зависимости между матрицами жесткости в системах
- 32. ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ Для выполнения операции суммирования необходимо преобразовать векторы
- 33. ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ Полученная система представляет собой систему линейных алгебраических
- 34. ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ Перед началом формирования глобальной матрицы всем ее
- 35. ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ Система уравнений МКЭ для конструкции может быть
- 36. ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ Выражение для полной потенциальной энергии может быть
- 37. ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ или в матричной форме: Поскольку матрица жесткости
- 38. ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ Учет граничных условий может быть осуществлен различными
- 39. РЕШЕНИЕ СИСТЕМЫ АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ Все методы решения систем линейных алгебраических уравнений делятся на две группы -
- 40. ОПРЕДЕЛЕНИЕ УСИЛИЙ В результате решения системы уравнений МКЭ определяется вектор узловых перемещений всей конструкции. Для вычисления
- 41. ОПРЕДЕЛЕНИЕ УСИЛИЙ Ранее было принято правило знаков, согласно которому положительные направления усилий совпадают с положительными направлениями
- 42. КОНЕЧНЫЙ ЭЛЕМЕНТ СТЕРЖНЯ ПРИ ЦЕНТРАЛЬНОМ РАСТЯЖЕНИИ И СЖАТИИ. ПРИМЕР РАСЧЕТА. Этот конечный элемент представляет собой прямолинейный
- 44. Скачать презентацию
Метод конечных элементов (МКЭ) - основной метод современной строительной механики, лежащий в
Метод конечных элементов (МКЭ) - основной метод современной строительной механики, лежащий в
Метод конечных элементов впервые был применен в инженерной практике в начале 50-х гг. XX в. Первоначально он развивался по двум независимым один от другого направлениям – инженерному и математическому. На раннем этапе формулировки МКЭ основывались на принципах строительной механики, что ограничивало сферу его применения. И только когда были сформулированы основы метода в вариационной форме, стало возможным распространение его на многие другие задачи. Быстрое развитие МКЭ шло параллельно с прогрессом современной компьютерной техники и ее применением в различных областях науки и инженерной практики.
Метод конечных элементов, как и многие другие численные методы, основан на представлении реальной континуальной конструкции ее дискретной моделью и замене дифференциальных уравнений, описывающих НДС сплошных тел, системой алгебраических уравнений. Вместе с тем МКЭ допускает ясную геометрическую, конструктивную и физическую интерпретацию.
Суть метода заключается в том, что область (одно- , двух- или
Суть метода заключается в том, что область (одно- , двух- или
В зависимости от типа конструкции и характера ее деформации КЭ могут иметь различную форму. Так, при расчете стержневых систем (фермы, балки, рамы) КЭ представляют собой участки стержней; для двумерных континуальных конструкций (пластины, плиты, оболочки) чаще всего применяют треугольные и прямоугольные (плоские или изогнутые) КЭ; а для трехмерных областей (толстые плиты, массивы) – КЭ в форме тетраэдра или параллелепипеда. В отличие от реального сооружения в дискретной модели конечные элементы связываются между собой только в отдельных точках (узлах) определенным конечным числом узловых параметров.
Виды МКЭ
По способу получения основных, т. е. разрешающих, уравнений различают четыре
Виды МКЭ
По способу получения основных, т. е. разрешающих, уравнений различают четыре
Прямой метод аналогичен матричному методу перемещений для стержневых систем, в основе его лежат положения, которые использовались на ранней стадии развития МКЭ. Этот метод удобен своей простотой и очевидным геометрическо-физическим значением отдельных шагов аппроксимации. Соотношения для КЭ здесь строятся непосредственно на основе трех групп уравнений (трех сторон задачи): статической, геометрической и физической.
Однако область применения прямого метода весьма ограничена: его можно использовать лишь для конечных элементов простой геометрии с малым числом степеней свободы в узле.
Вариационный метод основан на принципах стационарности некоторой переменной, зависящей от одной или
Вариационный метод основан на принципах стационарности некоторой переменной, зависящей от одной или
Метод невязок представляет собой наиболее общий подход к построению основных соотношений МКЭ. Этот метод
Метод невязок представляет собой наиболее общий подход к построению основных соотношений МКЭ. Этот метод
Метод энергетического баланса (метод Одена) основан на балансе различных видов энергии, записанном
Метод энергетического баланса (метод Одена) основан на балансе различных видов энергии, записанном
Из приведенных видов МКЭ в строительной механике особенно актуальны вариационный метод и метод взвешенных невязок Галеркина, которые для рассматриваемой задачи представляют собой два взаимно дополняющих метода одинаковой точности. Широкое применение этих методов обусловлено тем, что выражения в функционале или во взвешенном интеграле, как правило, имеют низший порядок производных по сравнению с производными в соответствующем дифференциальном уравнении для данной задачи. Это позволяет выбирать аппроксимирующие функции из более широкого семейства простых функций. Можно сказать, что вариационный вид МКЭ вышел из классического метода Ритца, а метод Галеркина – из обобщенного метода Бубнова-Галеркина. В принципе, из других методов также выводятся соответствующие виды МКЭ, однако их применяют значительно реже.
Формы МКЭ
В МКЭ, аналогично классическим методам строительной механики, за основные неизвестные
Формы МКЭ
В МКЭ, аналогично классическим методам строительной механики, за основные неизвестные
Метод перемещений – в настоящее время наиболее распространенная форма МКЭ. Это объясняется
Метод перемещений – в настоящее время наиболее распространенная форма МКЭ. Это объясняется
В основе математической формулировки МКЭ в форме метода перемещений лежит вариационный принцип Лагранжа, т. е. принцип минимума потенциальной энергии системы. Основными неизвестными здесь являются перемещения узловых точек дискретной схемы, напряжения же вторичны и определяются путем численного дифференцирования перемещений.
К достоинствам метода относятся: простота реализации; удовлетворительные точность и устойчивость решения с гарантированной сходимостью к нижней границе. Минусы: точность определения напряжений намного ниже, чем перемещений, хотя именно значения напряжений важны при прочностных расчетах, к тому же поскольку приближенное решение отвечает нижней границе, то значения и перемещений, и напряжений оказываются заниженными.
Принцип минимума дополнительной энергии и связанные с ним схемы МКЭ в форме метода сил, а
Принцип минимума дополнительной энергии и связанные с ним схемы МКЭ в форме метода сил, а
Главным плюсом МКЭ в форме метода сил является то, что основные неизвестные здесь – напряжения. И если бы в реализации метода сил не было определенных сложностей, значения напряжений можно было получать той же степени точности, что и перемещения в методе перемещений. Кроме того, использование принципа Кастилиано дает верхнюю границу приближенного решения (т. е. напряжения завышены), что в принципе лучше при расчетах на прочность, нежели заниженная оценка. Тем не менее, пока нет алгоритмов, в той же степени простых и устойчивых, имеющих гарантированную сходимость в обширном классе задач, подобно МКЭ в форме метода перемещений.
Аппроксимация
МКЭ относится к методам дискретного анализа. Однако в отличие от численных методов,
Аппроксимация
МКЭ относится к методам дискретного анализа. Однако в отличие от численных методов,
Сущность аппроксимации сплошной среды по МКЭ состоит в следующем:
1) рассматриваемая область разбивается на определенное число КЭ, семейство элементов по всей области называется системой или сеткой конечных элементов;
2) предполагается, что КЭ соединяются между собой в конечном числе точек – узлов, расположенных по контуру каждого из элементов;
3) искомые функции в пределах каждого КЭ (например, распределение перемещений, деформаций, напряжений и т. д.) с помощью аппроксимирующих функций выражаются через узловые значения, представляющие собой основные неизвестные МКЭ;
4) для анализа и расчета полученной системы конечных элементов действительны все принципы и методы, применяемые для любых дискретных систем.
Аппроксимирующие функции
Аппроксимация, как правило, дает приближенное, а не точное, описание действительного
Аппроксимирующие функции
Аппроксимация, как правило, дает приближенное, а не точное, описание действительного
Под точностью понимается отклонение приближенного решения от точного или истинного решения.
Устойчивость, прежде всего, определяется ростом ошибок при выполнении отдельных вычислительных операций. Неустойчивое решение является результатом неудачного выбора аппроксимирующих функций, «плохой» разбивки области на КЭ, некорректного представления граничных условий ит.п.
Под сходимостью подразумевается постепенное приближение последовательных решений к предельному, по мере того как уточняются параметры дискретной модели, такие как размеры элементов, степень аппроксимирующих функций и т. п. В этом смысле понятие сходимости аналогично тому значению, которое оно имеет в обычных итерационных процессах.
Выделенный элемент конструкции имеет те же свойства, что и конструкция в
Выделенный элемент конструкции имеет те же свойства, что и конструкция в
Неизвестные величины перемещений и усилий приводятся к узлам конечно-элементного разбиения. Например, в задачах, где в элементах конструкции реализуется растяжение (сжатие), стержневой элемент имеет две степени свободы при выборе неизвестных в виде узловых усилий (рис. а) или перемещений (рис. б). В зависимости от постановки задачи в каждом узле n-го элемента (1 или 2) имеется неизвестная величина с номером, соответствующим номеру узла.
Для каждого элемента можно выбрать оси, в которых можно записывать геометрические
Для каждого элемента можно выбрать оси, в которых можно записывать геометрические
При разбивке системы на конечные элементы необходимо учитывать возможные формы деформирования. Например, в задаче изгиба (рис. а) минимально возмож-ное количество конечных элементов 3. В задаче устойчивости (рис. б) минималь-но возможное количество элементов составляет 4, так как необходимо описать возможные отклонения стойки вследствие потери устойчивости.
Теоретически, чем чаще разбиение, тем точнее решение. Но в этом случае возникают трудности при решении системы уравнений, когда матрица жесткости системы становится плохо обусловленной.
Для системы в целом задаются произвольно выбранные оси, которые называются глобальными
Для системы в целом задаются произвольно выбранные оси, которые называются глобальными
Благодаря унификации и типизации элементов удается записать уравнения состояния конечного элемента и конструкции единообразно и кратко. Метод конечных элементов позволяет формализовать задачу в матричном виде, удобном для расчетов на ЭВМ. Большим достоинством МКЭ является его алгоритмичность. Наибольшее распространение в практике расчета и проектирования конструкций получил вариант МКЭ в перемещениях (в качестве искомых неизвестных принимаются перемещения рассчитываемой системы), который будет рассмотрен ниже.
Задачи статики, решаемые с помощью МКЭ, сводятся к следующим основным этапам:
Задачи статики, решаемые с помощью МКЭ, сводятся к следующим основным этапам:
2. Построение интерполирующих функций. Выбирается система функций (чаще всего кусочно-полиномиальная), однозначно определяющая перемещения внутри каждого конечного элемента через перемещения узловых точек. Интерполирующие функции подбираются таким образом, чтобы обеспечить непрерывность искомых величин (перемещений и их производных) в узлах элементов.
3. Вывод основных геометрических и физических соотношений. На основе выбранной системы интерполирующих функций выводятся зависимости между деформациями и перемещениями (геометрические соотношения), а также между напряжениями и деформациями (физические соотношения).
4. Построение матрицы жесткости конечного элемента и вектора узловых нагрузок.
4. Построение матрицы жесткости конечного элемента и вектора узловых нагрузок.
5. Переход от локальных координат к глобальным координатам системы. Перед тем, как включить матрицу· жесткости КЭ в глобальную матрицу жесткости всей конструкции, ее необходимо преобразовать из местной системы координат в общую для всей конструкции систему координат.
6. Получение системы уравнений метода конечных элементов. Каждая матрица жесткости отдельного конечного элемента включается в глобальную матрицу жесткости в цикле по элементам. Таким образом, формируется система алгебраических уравнений всей конструкции (уравнения равновесия), которая имеет вид:
где R - матрица жесткости системы (ансамбля) конечных элементов; z — вектор неизвестных узловых перемещений, Р - вектор узловых нагрузок.
В матрице жесткости R системы уравнений необходимо учесть граничные условия, так как в противном случае эта матрица будет вырожденной.
7. Решение системы алгебраических уравнений. Для решения системы линейных алгебраических уравнений
7. Решение системы алгебраических уравнений. Для решения системы линейных алгебраических уравнений
8. Определение деформаций и напряжений. Деформации, напряжения и усилия в конструкции определяются с помощью найденных узловых перемещений на основе геометрических и физических соотношений.
РАЗБИЕНИЕ СИСТЕМЫ НА КОНЕЧНЫЕ ЭЛЕМЕНТЫ
Данный этап включает в себя разбиение системы
РАЗБИЕНИЕ СИСТЕМЫ НА КОНЕЧНЫЕ ЭЛЕМЕНТЫ
Данный этап включает в себя разбиение системы
Разбиение стержневой системы на конечные элементы представляет собой первый и весьма важный шаг в процедуре решения задачи по МКЭ, поскольку от него во многом зависит точность получаемого решения. Успех на этом этапе зависит, в первую очередь, от имеющихся инженерных навыков. Неудачно выполненное разбиение конструкции на конечные элементы может привести к ошибочным результатам. При назначении размеров КЭ возникает задача оптимального разбиения конструкции на подобласти. При этом следует иметь в виду, что размеры элементов должны быть достаточно малыми для того, чтобы обеспечить приемлемую точность решения, с другой стороны, использование мелких КЭ приводит к большим системам алгебраических уравнений, решение которых связано со значительным объемом вычислительной работы.
В процессе разбиения конструкции на конечные элементы необходимо учитывать некоторые общие представления об окончательных результатах расчета для того, чтобы уменьшить размеры конечных элементов в зонах концентрации напряжений, где искомые величины быстро меняются, и увеличить размеры КЭ там, где искомые величины меняются медленно.
РАЗБИЕНИЕ СИСТЕМЫ НА КОНЕЧНЫЕ ЭЛЕМЕНТЫ
Важным моментом в процессе решения задачи по
РАЗБИЕНИЕ СИСТЕМЫ НА КОНЕЧНЫЕ ЭЛЕМЕНТЫ
Важным моментом в процессе решения задачи по
В настоящее время разработаны сервисные программы автоматизированной разбивки конструкции на конечные элементы и рациональной нумерации узлов.
ПОСТРОЕНИЕ ИНТЕРПОЛИРУЮЩИХ ФУНКЦИЙ
МКЭ основан на аппроксимации непрерывной функции, определенной на всей
ПОСТРОЕНИЕ ИНТЕРПОЛИРУЮЩИХ ФУНКЦИЙ
МКЭ основан на аппроксимации непрерывной функции, определенной на всей
Перемещения, являющиеся функциями координат произвольной точки конечного элемента, записываются через компоненты вектора узловых перемещений с помощью интерполирующей функции (функции формы или базисной функции):
где - матрица функции формы, - вектор узловых перемещений КЭ; s - количество степеней свободы КЭ.
Данные функции должны удовлетворять критериям полноты и совместности:
Критерий полноты. Интерполирующая функция должна обеспечивать постоянные значения рассматриваемых величин при уменьшении размеров элемента.
Для выполнения этого условия интерполирующая функция должна представлять собой полный полином как минимум степени р, где р - наивысший порядок производной в подынтегральном выражении функционала энергии.
Условие полноты удовлетворяется в том случае, когда
ПОСТРОЕНИЕ ИНТЕРПОЛИРУЮЩИХ ФУНКЦИЙ
Критерий совместности. Интерполирующая функция должна быть непрерывна вместе со
ПОСТРОЕНИЕ ИНТЕРПОЛИРУЮЩИХ ФУНКЦИЙ
Критерий совместности. Интерполирующая функция должна быть непрерывна вместе со
Критерии полноты и совместности представляют собой достаточные условия сходимости метода конечных элементов. При их выполнении с уменьшением размеров конечного элемента приближенные решения МКЭ монотонно сходятся к точному решению. Сказанное отнюдь не означает, что нарушение данных критериев приводит к невозможности получения достоверного решения. Существуют несовместные, и даже неполные элементы, которые дают высокую точность и быструю сходимость.
ВЫВОД ОСНОВНЫХ ГЕОМЕТРИЧЕСКИХ И ФИЗИЧЕСКИХ СООТНОШЕНИЙ
На основе выбранной системы интерполирующих функций
ВЫВОД ОСНОВНЫХ ГЕОМЕТРИЧЕСКИХ И ФИЗИЧЕСКИХ СООТНОШЕНИЙ
На основе выбранной системы интерполирующих функций
В общем виде зависимость между деформациями и перемещениями (геометрические соотношения) записывается следующим образом:
где ε - вектор, содержащий компоненты тензора деформаций; z - вектор узловых перемещений; В - матрица, связывающая вектор узловых перемещений с вектором, содержащим компоненты тензора деформаций.
Таким образом, полагается, что зависимость между деформациями и узловыми перемещениями носит линейный характер. Линейной зависимости отвечают такие условия работы конструкции, когда деформации и углы поворота малы по сравнению с единицей, а квадраты углов поворота малы по сравнению с соответствующими компонентами деформации.
ВЫВОД ОСНОВНЫХ ГЕОМЕТРИЧЕСКИХ И ФИЗИЧЕСКИХ СООТНОШЕНИЙ
Физические соотношения, определяющие зависимость между напряжениями
ВЫВОД ОСНОВНЫХ ГЕОМЕТРИЧЕСКИХ И ФИЗИЧЕСКИХ СООТНОШЕНИЙ
Физические соотношения, определяющие зависимость между напряжениями
где σ - вектор, содержащий компоненты тензора напряжений; D - матрица упругости.
Уравнения физических соотношений (уравнения состояния) представляют собой обобщенный закон Гука, устанавливающий прямо пропорциональную зависимость между напряжениями и деформациями и справедливый для определенного класса материалов на определенном участке графика σ-ε.
ПОСТРОЕНИЕ МАТРИЦЫ ЖЕСТКОСТИ КОНЕЧНОГО ЭЛЕМЕНТА И ВЕКТОРА УЗЛОВЫХ НАГРУЗОК
Полная потенциальная энергия
ПОСТРОЕНИЕ МАТРИЦЫ ЖЕСТКОСТИ КОНЕЧНОГО ЭЛЕМЕНТА И ВЕКТОРА УЗЛОВЫХ НАГРУЗОК
Полная потенциальная энергия
где W - потенциальная энергия деформации; А - потенциал внешних сил.
Потенциальная энергия деформации упругой системы определяется соотношением:
где р - вектор внешних распределенных нагрузок; S - площадь, на которой приложена нагрузка.
При этом деформации и напряжения, входящие в выражение для потенциальной энергии, так же выражаются через узловые перемещения.
Получение уравнений МКЭ в перемещениях основано на одном из фундаментальных энергетических принципов механики - принципе Лагранжа, согласно которому для системы, находящейся в состоянии равновесия, полная потенциальная энергия принимает стационарное значение.
где V — объем, занимаемый телом, а потенциал внешних распределенных нагрузок определяется по формуле:
ПОСТРОЕНИЕ МАТРИЦЫ ЖЕСТКОСТИ КОНЕЧНОГО ЭЛЕМЕНТА И ВЕКТОРА УЗЛОВЫХ НАГРУЗОК
Это условие записывается
ПОСТРОЕНИЕ МАТРИЦЫ ЖЕСТКОСТИ КОНЕЧНОГО ЭЛЕМЕНТА И ВЕКТОРА УЗЛОВЫХ НАГРУЗОК
Это условие записывается
Значение полной потенциальной энергии для всей области V равно сумме энергий отдельных конечных элементов:
Рассмотрим отдельный конечный элемент, индекс k при этом опускаем. Полная потенциальная энергия КЭ с учетом полученных формул запишется в виде:
где m - количество конечных элементов.
Тогда
ПОСТРОЕНИЕ МАТРИЦЫ ЖЕСТКОСТИ КОНЕЧНОГО ЭЛЕМЕНТА И ВЕКТОРА УЗЛОВЫХ НАГРУЗОК
где матрица
есть вектор
ПОСТРОЕНИЕ МАТРИЦЫ ЖЕСТКОСТИ КОНЕЧНОГО ЭЛЕМЕНТА И ВЕКТОРА УЗЛОВЫХ НАГРУЗОК
где матрица
есть вектор
представляет собой матрицу жесткости конечного элемента, а
ПЕРЕХОД ОТ ЛОКАЛЬНЫХ КООРДИНАТ К ГЛОБАЛЬНЫМ КООРДИНАТАМ СИСТЕМЫ
В процессе решения задачи
ПЕРЕХОД ОТ ЛОКАЛЬНЫХ КООРДИНАТ К ГЛОБАЛЬНЫМ КООРДИНАТАМ СИСТЕМЫ
В процессе решения задачи
в двух системах координат - местной (локальной), связанной с осью стержня, и общей (глобальной). Рассмотрим стержень, входящий в состав плоской стержневой системы, который имеет три перемещения (три степени свободы) в каждом узле: u, ν, φ. На рис. показан конечный элемент до и после деформации и отмечены компоненты вектора перемещений. Начальный узел конечного элемента обозначен через «i», конечный узел - через «j».
Учитывая, что углы поворота в начальной и конечной точках стержня φ′i и φ′j при преобразовании координат не меняются, получим:
ПЕРЕХОД ОТ ЛОКАЛЬНЫХ КООРДИНАТ К ГЛОБАЛЬНЫМ КООРДИНАТАМ СИСТЕМЫ
В матричной форме зависимость
ПЕРЕХОД ОТ ЛОКАЛЬНЫХ КООРДИНАТ К ГЛОБАЛЬНЫМ КООРДИНАТАМ СИСТЕМЫ
В матричной форме зависимость
где
ПЕРЕХОД ОТ ЛОКАЛЬНЫХ КООРДИНАТ К ГЛОБАЛЬНЫМ КООРДИНАТАМ СИСТЕМЫ
Значения sinα и cosα
ПЕРЕХОД ОТ ЛОКАЛЬНЫХ КООРДИНАТ К ГЛОБАЛЬНЫМ КООРДИНАТАМ СИСТЕМЫ
Значения sinα и cosα
Матрица преобразования координат V является ортогональной матрицей и для нее V-1 = VT . Тогда
Первоначально матрица жесткости для стержня строится в местной системе координат х'у’, связанной с осью стержня. Перед тем, как включить эту матрицу в единую матрицу жесткости всей конструкции, ее необходимо преобразовать из местной системы
координат в общую для всей конструкции систему координат ху
ПЕРЕХОД ОТ ЛОКАЛЬНЫХ КООРДИНАТ К ГЛОБАЛЬНЫМ КООРДИНАТАМ СИСТЕМЫ
Для получения зависимости между
ПЕРЕХОД ОТ ЛОКАЛЬНЫХ КООРДИНАТ К ГЛОБАЛЬНЫМ КООРДИНАТАМ СИСТЕМЫ
Для получения зависимости между
Сравнивая левую и правую части данного соотношения, находим :
Учитывая получим
ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ
Для выполнения операции
ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ
Для выполнения операции
Дифференцирование данного выражения по z, в соответствии с условием
приводит к системе уравнений:
где использовано правило дифференцирования матричных соотношений:
Подстановка данных формул в выражение для полной потенциальной энергии конечного элемента дает:
ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ
Полученная система представляет
ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ
Полученная система представляет
Матрица жесткости отдельного элемента , фигурирующая в полученном уравнении, представляет собой расширенную матрицу, размерность которой равна размерности глобальной матрицы. Поэтому использование процедуры суммирования в форме
при численной реализации МКЭ неэффективно. В практических расчетах выполняется прямое построение глобальной матрицы жесткости. В этом случае строится матрица R для отдельного конечного элемента по формуле
имеющая размерность sxs. Затем строкам и столбцам этой матрицы приписываются номера глобальных степеней свободы, что позволяет определить местоположение элементов матрицы жесткости КЭ в глобальной матрице жесткости.
ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ
Перед началом формирования
ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ
Перед началом формирования
Пусть, например, система состоит из двух КЭ, содержащих по два узла, в каждом из которых имеется одно неизвестное (одна степень свободы). Общее количество узлов системы - 3, размерность глобальной матрицы 3x3, элементы соединены между собой во 2-м узле.
Матрицы жесткости 1-го и 2-го КЭ с соответствующей нумерацией элементов и глобальная матрица имеют вид:
где верхний индекс указывает номер КЭ, нижние индексы соответствуют номерам строки и столбца в глобальной матрице.
ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ
Система уравнений МКЭ
ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ
Система уравнений МКЭ
Воспользовавшись соотношением
запишем выражение для полной потенциальной энергии в виде:
Пусть, как и в предыдущем примере, система состоит из двух КЭ, содержащих по два узла, в каждом из которых имеется одно неизвестное (одна степень свободы). Общее количество узлов - 3, общее количество узловых неизвестных - 3, размерность глобальной матрицы 3x3, элементы соединены между собой во 2-м узле.
ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ
Выражение для полной
ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ
Выражение для полной
Дифференцирование представленной выше квадратичной функции, позволяет получить систему линейных алгебраических уравнений. Учитываем при этом, что
Тогда
ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ
или в матричной
ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ
или в матричной
Поскольку матрица жесткости конечного элемента симметрична, то есть
то система уравнений (без учета граничных условий) будет иметь вид:
В матрице жесткости R системы уравнений необходимо учесть граничные условия, в противном случае она будет вырожденной, т. е. ее определитель будет равен нулю.
ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ
Учет граничных условий
ПОЛУЧЕНИЕ СИСТЕМЫ УРАВНЕНИЙ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ. УЧЕТ ГРАНИЧНЫХ УСЛОВИЙ
Учет граничных условий
1. Из матрицы R удаляются i-я строка и i-й столбец, соответствующие перемещению zi=0. После этого строки и столбцы матрицы перенумеровываются. Соответственно уменьшается размерность вектора узловых перемещений.
2.Уравнение zi=0, соответствующее граничному условию, формируется в составе матрицы R. Для получения zi=0 в матрице R i-я строка и i-й столбец, а также соответствующий элемент в векторе внешних нагрузок Р заполняются нулями. На место диагонального элемента rii в матрице R ставится единица. В результате порядок матрицы не изменяется, а заданные перемещения получают нулевые значения.
3. Для получения уравнения zi=0 диагональный элемент rii умножается на большое число. Порядок матрицы при этом не меняется.
РЕШЕНИЕ СИСТЕМЫ АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ
Все методы решения систем линейных алгебраических уравнений делятся
РЕШЕНИЕ СИСТЕМЫ АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ
Все методы решения систем линейных алгебраических уравнений делятся
Большинство известных алгоритмов прямого решения системы алгебраических уравнений МКЭ основано на методе Гаусса. При этом учитываются симметрия, положительная определенность, ленточная структура глобальной матрицы жесткости системы. В процессе исключения по Гауссу матрица коэффициентов с помощью прямого хода сводится к верхней треугольной форме, из которой обратным ходом находятся неизвестные перемещения.
ОПРЕДЕЛЕНИЕ УСИЛИЙ
В результате решения системы уравнений МКЭ определяется вектор узловых перемещений
ОПРЕДЕЛЕНИЕ УСИЛИЙ
В результате решения системы уравнений МКЭ определяется вектор узловых перемещений
Далее, с помощью матрицы жесткости конечного элемента R‘k` в местной системе координат определяются усилия на концах стержня по формуле:
где Ni, Qi, Mi - соответственно продольное и поперечное усилия и изгибающий момент в начальной точке КЭ; Nj, Qj, Mj - то же в конечной точке КЭ.
ОПРЕДЕЛЕНИЕ УСИЛИЙ
Ранее было принято правило знаков, согласно которому положительные направления усилий
ОПРЕДЕЛЕНИЕ УСИЛИЙ
Ранее было принято правило знаков, согласно которому положительные направления усилий
КОНЕЧНЫЙ ЭЛЕМЕНТ СТЕРЖНЯ ПРИ ЦЕНТРАЛЬНОМ РАСТЯЖЕНИИ И СЖАТИИ. ПРИМЕР РАСЧЕТА.
Этот конечный
КОНЕЧНЫЙ ЭЛЕМЕНТ СТЕРЖНЯ ПРИ ЦЕНТРАЛЬНОМ РАСТЯЖЕНИИ И СЖАТИИ. ПРИМЕР РАСЧЕТА.
Этот конечный
Здесь и далее апостроф, указывающий на то, что матрица жесткости построена в местной системе координат опускается.
Интерполирующая функция для продольного перемещения принимается в виде:
Линейная функция будет точно описывать продольные перемещения на участках стержня между точками приложения сосредоточенных сил.
Коэффициенты α0 и α1, определяются с помощью граничных условий:
при x=xi u=ui;
при x=xj u=uj.