Расчёт ферм методом конечных элементов. Учебное пособие по дисциплине Строительная механика
.pdf
МИНИСТЕРСТВО ТРАНСПОРТА РФ
МОСКОВСКАЯ ГОСУДАРСТВЕННАЯ АКАДЕМИЯ ВОДНОГО ТРАНСПОРТА



МГАВТ



М.А. Галабурда
Расчет ферм методом конечных элементов
«Альтаир» - МГАВТ Москва 2006
МИНИСТЕРСТВО ТРАНСПОРТА РФ
МОСКОВСКАЯ ГОСУДАРСТВЕННАЯ АКАДЕМИЯ ВОДНОГО ТРАНСПОРТА
М.А. Галабурда
Расчет ферм методом конечных элементов
Учебное пособие по дисциплине
Строительная механика
«Альтаир» - МГАВТ Москва 2006
1
УДК 621.86
М.А. Галабурда Расчѐт ферм методом конечных элементов. Учебное пособие по дисциплине
Строительная механика -М.: «Альтаир» - МГАВТ, 2006 - 25 с., ил.
В Учебном пособии «Расчѐт ферм методом конечных элементов» даѐтся изложение идеи метода конечных элементов и возможность применения его конкретно к шарнирно-стержневым системам.
Приведѐн пример расчѐта фермы с помощью МКЭ.
Предназначено для студентов, обучающихся по специальности 190602.65 (150900) «Эксплуатация перегрузочного оборудования портов и транспортных терминалов». Кроме того, оно может быть полезным для аспирантов и инженерно-технических работников, занимающихся вопросами прочности стержневых металлоконструкций.
Рецензент – к.т.н. доцент кафедры ППТМиР А.И. Вашунин
Рассмотрены и рекомендованы к изданию кафедрой Портовых подъѐмнотранспортных машин и робототехники МГАВТ (протокол №6 от 15.05 2006 г)
Рекомендованы к изданию учебно-методическим советом МГАВТ (протокол № 7 от 18.05.2006)
Ответственность за оформление и содержание представленных к изданию материалов несут автор и кафедра.
© МГАВТ, 2006г.
© М.А. Галабурда, 2006г.
2
Введение.
В последние годы исключительно широкое использование в прочностных расчетах строительных, авиационных и судовых конструкций получил метод конечных элементов (МКЭ).
Этот метод нашел применение не только в задачах прочности, но и при решении задач прочности, теплопроводности, термоупругости и термопластичности, волновых процессов, гидромеханики, гидроупругости и многих других. В настоящее время МКЭ является одним из наиболее эффективных методов расчета, применяемых в строительной механике. Если сравнительно недавно в строительной механике стержневых статически неопределимых систем главное внимание уделялось поиску таких основных систем, которые бы позволили максимально упростить матрицу единичных коэффициентов, а в идеале свести ее к диагональной, то в настоящее время не число лишних неизвестных определяет трудоемкость расчетов, а трудности, связанные с общением человека с машиной, с созданием соответствующего математического обеспечения.
§1. Идея метода.
Название этого метода в какой-то степени определяет его сущность, так как рассчитываемая конструкция представляется в виде совокупности отдельных элементов, соединенных между собой в конечном числе узловых точек. Эти элементы могут представлять собой как отдельные стержни шарнирно-стержневой (Рис 1) или рамной конструкции (Рис 2), так и плоские элементы треугольной или прямоугольной формы, если речь идет о расчете плоской (Рис 3,4) или оболочечной конструкции (Рис 5). При расчете трехмерных тел их дискретизацию осуществляют с помощью объемных элементов различной формы (Рис 6).
.
3
Рис. 1. Плоская ферма
Рис. 2. Идеализация рамной конструкции совокупностью стержневых элементов
Рис. 3. Идеализация поперечной рамы супертанкера совокупностью плоских конечных элементов
4
Отсюда следует, что членение системы на элементы заданной конфигурации отнюдь не связано с получением отдельных частей, на которые разрезается конструкция. Фактически конечные элементы представляют собой упругие элементы особого типа, на деформации которых наложены связи, заставляющие их изменяться по определенной форме так, чтобы по возможности сохранилась непрерывность деформации расчетной модели.
Повысить точность расчета можно за счет увеличения в расчетной модели числа конечных элементов, что приводит к точному удовлетворению условий неразрывности перемещений в большем числе точек.
Каждому узловому перемещению (линейному или угловому) соответствует узловое усилие. Совокупность этих усилий определяет влияние смежных элементов конструкции на рассматриваемый элемент.
Упругие свойства отдельных элементов (или стержней), на которые разбивается конструкция, определяются их матрицей податливости или матрицей жесткости. Эти матрицы определяют связь между узловыми усилиями и узловыми перемещениями рассматриваемого конечного элемента.
При использовании МКЭ представляется возможным, как это сделано в теории упругости, сформулировать три группы уравнений:
1.Статические, описывающие равновесное состояние системы.
2.Геометрические, устанавливающие связь между деформациями
иперемещениями.
3.Физические, которые связывают между собой силы и деформации.
В соответствии с тем, какие величины принимаются за неизвестные, различают три классических метода решения: метод сил, метод перемещений и смешанный метод.
При рассмотрении стержневых систем, у которых в рамках допущений, принимаемых в сопротивлении материалов, всегда можно точно установить напряженно-деформированное состояние, МКЭ является точным методом.
Иное дело, когда рассматриваются пластинки, плиты, оболочки, или массивные тела трех измерений (см. рис. 4, 5, 6). Здесь, при рассмотрении самих конечных элементов, а также при составлении соответствующих условий неразрывности, не представляется возможным получить точные соотношения, в связи с чем МКЭ
5
приобретает вариационный смысл.
2б среднее
Рис.4. Растяжение полосы с отверстием (плоское напряженное состояние)
Рис. 5 Идеализация оболочки совокупностью треугольных конечных элементов.
6
Рис.6.Трехмерная задача о сосуде высокого давления (96 элементов, 707 узлов, 2121 степень свободы)
Сочленение элементов осуществляется только в узлах. В них полностью удовлетворяются условия равновесия и неразрывности перемещений. Разрезание рассматриваемой области проводится к кажущемуся нарушению условий неразрывности перемещений на участках между узлами, в значительной степени компенсируемому предположением о линейном законе изменения напряжений в любом сечении конечного элемента. Это обуславливает наложение на деформации элемента сильно ограничивающих их связей, которые, с одной стороны, улучшают условия соблюдения неразрывности деформаций, а с другой, - не вызывают концентрации напряжений в узловых точках.
В этих случаях к проблеме подходят следующим образом. 1.Рассматриваемый объект разбивают линиями или
поверхностями на некоторое количество конечных элементов.
2. Предполагают, что элементы связаны между собой только в узловых точках. Основными неизвестными являются перемещения этих узлов.
7
3.Выбирают функцию (или функции), однозначно определяющие перемещения внутри каждого конечного элемента через перемещения узлов этих элементов.
4.Функции смещения, единственным образом определяющие деформированное состояние элемента через узловые перемещения, совместно с граничными условиями и упругими свойствами материала, определяют напряженное состояние как внутри элемента так и на его границах.
5.Определяют систему сил, сосредоточенных в узлах и уравновешивающих напряжения на границе, в результате чего получают жесткостную связь с перемещениями.
После выполнения всех этих пунктов дальнейшее решение задачи может быть выполнено обычным методом матричного решения задач строительной механики. Очевидно, что такой подход является приближенным. Во-первых, не всегда удается добиться, чтобы выбранные нами функции перемещений удовлетворяли условиям непрерывности перемещений между соседними элементами. Поэтому условия совместности перемещений на границах элементов может быть нарушено (хотя внутри каждого элемента это условие выполняется в связи с однозначностью смещения, следующей из функционального представления).
Во-вторых, не исключается местное нарушение условий равновесия внутри каждого элемента и на его границах.
Выбор формы элемента и функций перемещений для конкретных задач зависит от интуиции и квалификации расчетчика. Точность приближенного решения в значительной степени зависит от этого выбора. Подход к решению задачи, описанный здесь, есть не что иное, как формулирование задачи в перемещениях и известен в строительной механике, как метод перемещений.
С тем чтобы придать изложению МКЭ большую ясность, а восприятию большую простоту, ограничимся рассмотрением решения с помощью МКЭ наиболее простых задач, таких как шарнирностержневые системы.
8
§2. Расчет по МКЭ шарнирно-стержневых систем.
Рассмотрим произвольную плоскую шарнирно-стержневую систему, нагруженную только в узлах (рис. 7)
Y |
F2 |
2 |
F3 |
|
|
2 |
3 |
|
|
|
|
|
||
|
1 |
|
|
3 |
|
6 |
|
7 |
|
1 |
|
|
|
4 |
|
5 |
5 |
4 |
X |
|
|
Рис. 7.Плоская шарнирно-стержневая система
Из теории ферм известно, что напряженно-деформированное состояние стержней является однородным. Поэтому, в данном случае, в качестве конечных элементов, на которые разбивается конструкция, являются стержни фермы. Число узлов в ферме будем обозначать через m, а число стержней через s.
Прежде, чем переходить к рассмотрению равновесия |
шарнирно- |
стержневой системы, изложим прием ее аналитического |
описания. |
Под последним следует понимать информацию о геометрии системы, об условиях сопряжения элементов (стержней) в узлах и, соответственно, о граничных условиях.
В первую очередь, дадим описание структуры системы. Для этого введем понятие структурной матрицы SC , которая представляет собой прямоугольную матрицу, число строк которой равно числу узлов m, а число столбцов – числу стержней (элементов) s.
Осуществив произвольную нумерацию узлов и стержней, следует построить структурную матрицу так, чтобы в каждом ее столбце находились лишь два числа: 1 и -1. Номер столбца совпадает с номером стержня. При этом 1 располагается в строке, номер которой совпадает с номером узла, соответствующего началу стержня, а -1 в строке, номер которой совпадает с номером узла, к которому примыкает конец стержня (элемента).
За начало стержня всегда будем принимать то его сечение, которое
9
