Файл: Лекция Приближенные методы расчета в строительной механике.doc
Добавлен: 10.02.2019
Просмотров: 517
Скачиваний: 4
Лекция 17
Приближенные методы расчета в строительной механике
Учебные вопросы:
17.1. Введение
17.2. Идея метода конечных элементов (МКЭ)
17.3. Метод конечных элементов в расчетах плоских стержневых систем
17.4 Матрица жесткости КЭ в местной системе осей координат
17.5 Матрица преобразований (направляющих косинусов)
17.6 Матрица жесткости для сооружения в целом
17.6 Порядок расчета стержневых систем методом конечных элементов
17.7 Формирование матрицы жесткости плоского треугольного конечного элемента в локальной системе осей координат
17.1. Введение
Большинство задач строительной механики, связанных с исследованием напряженно-деформированного состояния конструкций и их элементов (стержней, пластин, оболочек) сводится, как правило, к решению одного или нескольких дифференциальных уравнений равновесия элемента, соответственно с одним или несколькими неизвестными.
Точное решение таких уравнений (решение в замкнутом виде) не представляет затруднений лишь в некоторых элементарных случаях. При решении реальных задач часто приходится сталкиваться с таким объемом вычислительных работ, что от точного решения отказываются, а во многих случаях точное решение задачи вообще невозможно, т.к. граничные условия или условия на контуре просто не выражаются в аналитической форме. Поэтому, как правило, при решении практических задач приходится прибегать к приближенным методам решения.
Приближенные методы решения задач могут быть разбиты на две основные группы:
1. Вариационные методы, которые дают приближенные аналитические выражения искомой функции (функции перемещений или функции внутренних усилий).
2. Численные методы, которые дают значения искомой функции при тех или иных значениях аргумента.
К первой группе относятся вариационные методы Ритца, Бубнова-Галеркина, метод Треффца и другие.
Ко второй группе относятся метод сеток и его более совершенная модификация - метод конечных элементов, а также ряд графических и полуграфических методов таких, как, например, метод прямых, метод коллокаций и другие.
Преимущество вариационных методов заключается в том, что задача сводится обычно к решению системы двух, трех, редко четырех уравнений, которые дают хорошее приближение к действительному состоянию сооружения. К их недостаткам следует отнести то, что возможности вариационных методов ограничены сложными контурами и сложными законами распределения внешней нагрузки, т.к. применение вариационных методов требует, чтобы было, хотя бы в приближенной форме, определено аналитическое выражение внешней нагрузки, деформированной упругой поверхности элемента и др. условий задачи.
Численные методы, в сравнении с вариационными, имеют более универсальный характер, т.к. не требуют аналитических выражений условий задачи. Однако численные методы обладают рядом недостатков. Так, для получения удовлетворительного решения они требуют нанесения на исследуемую область густой сетки или разбиения на достаточно большое число элементов, что неизбежно влечет за собой решение систем алгебраических уравнений с большим числом неизвестных, что становится возможным только при наличии ЭВМ. Кроме того, численные методы часто приводят к неточности решений, особенно в местах приложения сосредоточенных сил, при наличии острых углов, подкреплений и т.д., т.е. там, где нарушается гладкость полей переменных.
Весьма существенным недостатком численных методов является то, что они не дают аналитического выражения искомой функции, а, следовательно, для определения параметров напряженно-деформированного состояния в данной области приходится вычислять эти величины во всех узлах стыковки элементов, т.е. получать массу ненужной информации для тех областей, которые нас совершенно не интересуют.
Таким образом, выбор метода расчета при рассмотрении конкретной упругой системы зависит от постановки задачи и исходных условий, а также от вооруженности расчетчика вычислительной техникой.
17.2. Идея метода конечных элементов (МКЭ)
Исторически возникновение МКЭ связано с идеей применения хорошо разработанных процедур для расчета статически неопределимых стержневых систем к решению континуальных задач.
Первоначально эта идея была высказана еще в 1933 году И.М. Рабиновичем, но развитие получила только в 70-х годах, с появлением ЭВМ.
Метод конечных элементов основан на мысленном представлении сплошного тела в виде совокупности отдельных конечных элементов, взаимодействующих между собой в конечном числе точек, которые в МКЭ принято называть узлами.
Система разбивается на простые конечные элементы (КЭ) напряженно-деформированное состояние которых исследуется заранее.
Так стержневые системы могут быть разбиты на элементы в виде прямолинейных или криволинейных стержней (например, для расчета арок) с различными условиями соединения элементов в узлах. В этом случае дискретная модель является точной копией исходной конструкции (с учетом принятых технических гипотез).
В расчетах пластин наибольшее распространение получили прямоугольные и треугольные конечные элементы. Здесь дискретная модель лишь приближенно отражает поведение исходной конструкции.
Заметим, что даже при одном и том же числе узловых точек различные схемы дискретизации исходной конструкции порождают разницу в окончательных результатах расчета. К сожалению, заранее сказать, какая из возможных схем дискретизации приведет к наименьшей погрешности расчета, невозможно.
Число степеней свободы КЭ, а в конечном итоге число неизвестных МКЭ, определяется количеством наложенных в узлах дополнительных связей.
Условия
равновесия и совместности деформаций
выполняются только в узловых точках -
точках соединения КЭ. Однако это не
значит, что общая жесткость пластины
при этом резко уменьшается, поскольку
зависимость между узловыми усилиями и
деформациями каждого элемента
рассматривается с учетом некоторых
внутренних связей.
Каждый элемент является частью заменяемой среды, т.е. сплошное тело лишь условно делится на отдельные элементы конечных размеров. Выделенный элемент имеет те же физические свойства и геометрические характеристики, что и рассматриваемая конструкция в месте расположения элемента.
Все внешние силы считаются приложенными в узлах, по направлению их возможных перемещений. Вне узловые нагрузки предварительно приводятся к узловым.
При реализации МКЭ наибольшее распространение получили идеи метода перемещений, хотя имеются работы, где рассматривается метод сил и смешанный метод. Предпочтение методу перемещений отдано в основном из-за простоты выбора основной системы, составления матрицы жесткости и формирования вектора внешних нагрузок.
Разрешающее уравнение МКЭ, которое представляет собой матричную форму канонических уравнений метода перемещений, имеет вид:
[r]{Z}={P},
где: [r] - матрица жесткости сооружения в целом,
{Z}- вектор перемещений узловых точек сооружения,
{P}- вектор внешних нагрузок.
Подход к решению задачи МКЭ является единым, как для стержневых систем, так и для пластин, оболочек и объемных тел.
Дальнейшее рассмотрение МКЭ будем проводить на примере плоских стержневых систем.
17.3. Метод конечных элементов в расчетах плоских стержневых систем
Рассмотрение упругих систем вообще и плоских стержневых систем в частности с позиций МКЭ есть представление упругих систем в виде набора элементов с конечным числом степеней свободы, которые соединяются между собой в узловых точках (узлах). Такое представление заданной системы в виде дискретной модели приводит к полной формализации всех этапов расчета. Подход к решению задачи является единым, как для стержневых систем, так и для пластин, оболочек, объемных тел и т.п.
Рассматривать будем МКЭ разработанный на базе метода перемещений, применительно к расчету плоских стержневых систем.
При расчете плоских стержневых систем в МКЭ приняты те же гипотезы, что и в обычном методе перемещений. Несколько уточняется только одна гипотеза: в МКЭ будем учитывать влияние не только изгибных, но и влияние продольных деформаций на перемещения узловых точек сооружения. Т.е. длина стержня в результате деформаций растяжения-сжатия может изменяться. Это положение позволяет в большей степени формализовать выбор основной системы МКЭ и получить результаты расчета более точные, чем в обычном методе перемещений.
Расчет стержневых систем, как и любых других, в МКЭ начинают с разбиения заданной системы на отдельные элементы.
В качестве конечных элементов (КЭ) мы будем рассматривать прямолинейные стержни, имеющие постоянную жесткость по длине. Между собой КЭ могут соединяться жестко или с помощью шарнира. Точки соединения элементов в МКЭ называют узловыми или узлами.
Т.е. основную систему (дискретную модель) МКЭ получают, разбивая заданную систему на отдельные прямолинейные элементы, имеющие постоянную жесткость по длине. При наличии в системе криволинейных стержней или стержней с переменной жесткостью, их, с достаточной степенью точности, разбивают на участки, в пределах которых стержни считают прямолинейными, с усредненной постоянной жесткостью.
Кроме того, алгоритм МКЭ требует, чтобы все внешние нагрузки, действующие на сооружение, были приложены к узловым точкам ее дискретной модели. Поэтому, точки приложения сосредоточенных сил делают узловыми, а нагрузки распределенные по длине стержня, преобразуют к узловым.
Д
ля
преобразования вне узловой нагрузки к
узловой используют таблицы метода
перемещений (см. рис. б).
В узлах, где отдельные элементы соединяются между собой жестко, имеется три неизвестных перемещения, в шарнирных узлах – два. Следовательно, количество неизвестных МКЭ можно определить:
n = 3nж.уз. + 2nш.уз .
П
оложительные
направления и порядок нумерации
неизвестных принимаем следующий
Д
искретная
модель сооружения в целом, которая
связывается с общей системой осей
координат, характеризуется n
параметрами перемещений Zi
и узловых силовых воздействий Pi
, составляющих векторы
Каждый конечный элемент связан с местной
системой осей координат и характеризуется
своими параметрами узловых перемещений
{V} и соответствующими
узловыми усилиями {S}’.
Кроме того, для элементов, в пределах
которых вне узловая нагрузка преобразуется
к узловой, формируют векторы грузовых
усилий {So}.
В разрешающем уравнении МКЭ
[r]{Z}={P},
матрица [r], которая называется матрицей жесткости сооружения в целом, формируется из матриц жесткости отдельных элементов.
17.4 Матрица жесткости КЭ в местной системе осей координат
Матрица жесткости КЭ в местной системе осей координат [r]j’ устанавливает связь между единичными перемещениями и усилиями в узловых точках, для горизонтально расположенного элемента.
Как уже отмечалось, основная система МКЭ представляет собой совокупность отдельных стержней, которые соединяются между собой жестко или с помощью шарнира.
Следовательно, все многообразие КЭ можно свести к трем типам элементов:
а
)
элемент с двух сторон ограниченный
жесткими узлами, такой элемент имеет
шесть степеней свободы, т.е. неизвестными
являются шесть возможных перемещений
узловых точек элемента
б
)
элемент с одной стороны ограниченный
жестким узлом, а с другой шарнирным,
такой элемент имеет пять степеней
свободы
г) элемент с двух сторон ограниченный шарнирными узлами, такой элемент имеет четыре степени свободы
Рассмотрим формирование матриц жесткости для показанных выше трех типов КЭ, в местной системе осей координат.
Для элемента с шестью степенями свободы (тип а) матрица жесткости будет иметь вид
[r]j’=
где r12’ - реактивное усилие в связи 1 от единичного перемещения связи 2 и т.д.
В блочном виде матрица [r]jэ будет иметь вид
[r]j’=
где rнк’ - блок реакций в связях, наложенных на начальный узел, от единичных перемещений связей, наложенных на конечный узел КЭ и т.д.
Численные значения реакций будем определять используя таблицы метода перемещений. Знак реакций будет положительным, если направление реакций совпадает с положительным направлением перемещений Vj.
1
.Формирование
матрицы жесткости для элемента с шестью
степенями свободы (тип а)
З
ададим
последовательно единичные перемещения
связям V1=1,
V2=1...V6=1.
1) V1=1
В соответствие с законом Гука
l
=
l
= V1 = 1,
N =
2
)
V2=1
3
)
V3=1
и т.д. задаем V4= 1, V5 = 1, V6 = 1
Расставляя реакции в соответствующие ячейки, получим матрицу жесткости для элемента с шестью степенями свободы
[r]j’
=
Запишем матрицы жесткости для двух других типов элементов, формируя их аналогично предыдущей.
-
Матрица жесткости для элемента с пятью степенями свободы (тип б) будет иметь вид
[r]j’
=
3. Матрица жесткости для элемента с четырьмя степенями свободы (тип в) запишется
[r]j’ =
17.5 Матрица преобразований (направляющих косинусов)
Рассмотрим теперь КЭ в составе рамы,
расположенный под углом
к оси Х в общей системе осей координат.
Необходимо перейти от матрицы реакций [r]’j в местной системе осей координат к матрице [r] в общей системе координат.
З
адачу
решаем следующим образом. В начале
построим матрицу [c]j,
которая преобразует перемещения КЭ
{z}j в
общей системе осей координат в перемещения
{v}j,
по выражению
{v}j= [c]j {z}j
V1= Z1cos + Z2sin
V2 = -Z1 sin + Z2 cos
V3 = Z3
V6 = Z6
V4 = Z4 cos + Z5 sin
V5= -Z4sin + Z5cos
В матричной форме приведенная выше запись будет иметь вид
или в блочной форме
[c]j
=
где для жесткого узла
[c]н,к =
для шарнирного узла
[c]н,к=
Так как мы рассматриваем плоские упругие системы, то векторы узловых усилий и узловых перемещений, как для отдельного элемента, так и для сооружения в целом, связаны между собой линейно
{S}j’ = [r]j’{V}j - в местной системе осей координат.
{S}j = [r]j{Z}j - в общей системе осей координат
Кроме того
{V}j = [c]j{Z}j,
Аналогично
{S}j’ = [c]j{S}j,
где {S}’,{S}-узловые усилия КЭ соответственно, в местной и общей системах осей координат.
Тогда
{S}j = [c]j-1 {S}j’ = [c]j-1[r]j’{V}j = [c]j-1[r]j’[c]j{Z}j.
Для матрицы направляющих косинусов выполняется равенство
[c]j-1 = [c]jT,
Тогда {S}j = [c]jT [r]j’ [c]j{Z}.
Обозначим
[r]j = [c]jT [r]j’ [c]j - это выражение и является формулой для вычисления матрицы жесткости КЭ в общей системе осей координат.
При формировании матриц жесткости отдельных элементов [r]j’ должны быть зафиксированы начало и конец каждого стержня, так как от этого зависит знак угла , определяющего ориентацию стержня в общей системе осей координат ХОY.
17.6 Матрица жесткости для сооружения в целом
Рассмотрим теперь как формируется матрица жесткости для сооружения в целом, когда имеются матрицы жесткости для всех отдельных элементов.
П
усть
задана какая-либо стержневая система
В
се
узлы будем считать жесткими, т.е. с каждым
из них связано по 3 возможных перемещения.
Матрицу жесткости для всего сооружения покажем в блочном виде, с размерами блоков 3x3, т.к. с каждым узлом связано по 3 возможных перемещения (горизонтальное, вертикальное и поворот узла).