Основные определения метода конечных элементов

Построение объемлющей эпюры изгибающих моментов

Пример расчета стержневой системы методом конечных элементов

Формирование и решение системы уравнений МКЭ. Определение внутренних усилий в элементах

Примеры построения матрицы жесткости конечного элемента

Преобразование матрицы жесткости при переходе от одной системы координат к другой

Конечный элемент. Матрица жесткости конечного элемента

Основные определения метода конечных элементов

Лекция 17

III желудочек

III желудочек ( ventriculus tertius) –представляет собой узкую вертикальную щель, расположенную между двумя зрительными буграми, имеет следующие стенки:

1) боковые - образованы медиальными поверхностями двух зрительных бугров;

2) верхняя - образована мягкой мозговой оболочкой, в состав которой входит сосудистое сплетение.

3) нижняя - ее формируют все образования гипоталамуса. На нижней стенке имеются два углубления - карман воронки и зрительный карман.

4) передняя - образована столбами свода, передней мозговой спайкой и конечной пластинкой. На передней стенке имеются два отверстия, ведущие в боковые желудочки - межжелудочковые отверстия.

5) задняя - здесь располагаются все образования эпиталамуса с шишковидным углублением, задняя мозговая спайка, а ниже открывается водопровод мозга, ведущий в IV желудочек.

Расчет статически неопределимых стержневых систем методом конечных элементов

Учебные вопросы:

-----------------------------------------------------------------------------------------------------------------

Метод конечных элементов (МКЭ) - основной метод современной строительной механики, лежащий в основе подавляющего большинства современных программных комплексов, предназначенных для выполнения расчетов строительных конструкций на ЭВМ. МКЭ также используется для решения других разнообразных задач как в области прочностных расчетов, так и во многих других сферах, например задачах гидродинамики, электромагнетизма, теплопроводности и многих других. Математические основы метода были впервые сформулированы Р.Курантом в 1943г., а термин “конечный элемент” впервые был введен Р.Клафом в 1960 г.

Метод конечных элементов позволяет практически полностью автоматизировать расчет стержневых систем, хотя, как правило, требует выполнения значительно большего числа вычислительных операций по сравнению с классическими методами строительной механики. Однако, в современных условиях большой объем вычислений не является серьезной проблемой, и, в связи с этим, при внедрении ЭВМ в инженерную практику МКЭ получил широчайшее распространение. Поэтому, знание основ метода конечных элементов и современных программных средств, позволяющих на его основе решать разнообразные задачи, в наше время для инженера является абсолютно необходимым.

Рис.17.1

В МКЭ стержневая система мысленно разбивается на отдельные части - конечные элементы, соединяющиеся между собой в узлах (рис.17.1). Узлы могут быть жесткими и шарнирными. Совокупность соединенных между собой и прикрепленных к основанию конечных элементов образует расчетную схему метода, называемую конечно-элементной схемой или конечно-элементной моделью или просто системой элементов. Элементы и узлы конечно-элементной схемы нумеруются.

Внешняя нагрузка считается приложенной только в узлах конечно-элементной схемы. В общем случае переход от заданной нагрузки к узловой осуществляется следующим образом. На основании принципа суперпозиций рассматриваемое состояние стержневой системы может быть представлено как сумма двух состояний (рис.17.2). В первом состоянии (задача 1) вводятся связи, препятствующие всем возможным смещениям узлов системы, аналогично тому, как образуется основная система в методе перемещений. При этом, однако, продольными деформациями стержней не пренебрегают. От действия заданных нагрузок во введенных связях возникают реакции. Во втором состоянии (задача 2) узлы конечно-элементной схемы не закреплены от смещений, но к ним прикладываются усилия равные по модулю реакциям в связях, определенным в первом состоянии, но противоположные им по направлению (рис.9.2). Расчет системы в первом состоянии не представляет труда. В частности, если конечно-элементная схема создается таким образом, чтобы элементы представляли собой отдельные стержни (элементы 1, 2 и 3 на рис.17.2), то для каждого из таких элементов имеется табличное решение, позволяющее определить реакции в связях и построить эпюры внутренних усилий по их длине. Для расчета же системы во втором состоянии, т.е. для решения задачи 2, и применяется метод конечных элементов. Окончательное решение задачи будет представлять собой сумму решений этих двух задач.

Рис.17.2

Рис.17.3

В задаче 2 усилия, действующие на любой элемент приложены исключительно в узлах. В этом случае перемещения узлов любого элемента, взятого в отдельности (рис.17.3), однозначно определяют усилия и перемещения в любой точке этого элемента. Как известно, для стержневых систем решение такой задачи может быть найдено точно.

Каждый, взятый отдельно от системы, конечный элемент должен быть достаточно простым, чтобы имелась возможность легко определить перемещения и усилия в любом сечении стержней элемента по заданным перемещениям его узлов. Связь между перемещениями узлов элемента и усилиями в них задается при помощи матрицы жесткости элемента. Количество перемещений узлов элемента, которые однозначно определяют состояние данного элемента называют числом степеней свободы элемента. Оно определяется по формуле:

(17.1)

где -число шарнирных узлов в элементе, а - число жестких узлов в элементе. Действительно, если узел представляет собой шарнир, то его положение на плоскости можно охарактеризовать двумя линейными перемещениями, например в вертикальном и горизонтальном направлениях. В случае жесткого узла необходимо еще дополнительно к линейным смещениям задать его поворот.

Рис.17.4

На рис.17.4 первый элемент характеризуется четырьмя степенями свободы, т.к. он содержит два шарнирных узла. При отсутствии нагрузки, кроме приложенной в самих узлах, положение на плоскости любой точки этого элемента определяется четырьмя параметрами - двумя вертикальными и двумя горизонтальными перемещениями узлов элемента. У второго элемента на рис.17.4 - пять степеней свободы - к четырем линейным смещениям добавляется поворот в одном из узлов. У третьего элемента - шесть степеней свободы, которым соответствуют четыре линейных и два угловых перемещения.

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

Рис.17.5

Например, конечно-элементная схема висячей системы, изображенной на рис.17.5. содержит в себе 28 шарнирных узлов, следовательно характеризуется 56 степенями свободы.

Рис.17.6

В конечно-элементной схеме балки (рис.17.6) используется один жесткий и три шарнирных узла. Следовательно, в соответствии с (17.1) эта схема характеризуется 3×2+1×3=9 степенями свободы.

Для всех элементов, из которых состоит конечно-элементная схема, должны быть построены матрицы жесткости элементов. В программных комплексах, реализующих алгоритм метода конечных элементов, хранятся готовые матрицы жесткости для элементов различных типов.

На практике, при расчете плоских стержневых систем используют готовые матрицы жесткости для элементов только трех типов: простых стержней с двумя жесткими узлами, двумя шарнирными узлами, одним жестким и одним шарнирным узлом (рис.9.4). В этом случае при разбивке стержневой системы на элементы узлы вводятся в местах соединения и изломов стержней, в опорах, шарнирах и на свободных концах консольных стержней. В принципе узел может быть введен и в любых других точках, например, в точках приложения сосредоточенных сил.

В учебных целях могут использоваться и элементы других типов (рис.17.7), в том числе и включающие в себя опорные закрепления.

Рис.17.7

Из построенных матриц жесткости элементов формируется матрица жесткости системы. Для этого все матрицы жесткости элементов и матрица жесткости системы должны быть сформированы в единой системе осей координат, называемой глобальной системой осей координат. При расчете плоских стержневых систем традиционно используется следующая глобальная система осей координат (рис.17.8): ось 1 направлена вправо, ось 2 - вверх, ось 3 - против часовой стрелки.

Рис.17.8

Матрицы жесткости элементов могут формироваться и храниться в памяти ЭВМ в своих, локальных системах осей координат, в общем случае отличных от глобальной системы осей координат. В данной ситуации при помощи специальной процедуры эти матрицы должны быть перестроены для глобальной системы осей координат.

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

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

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

Итак, расчет стержневой системы методом конечных элементов в форме метода перемещений состоит из следующих этапов:

Создание конечно-элементной схемы (разбивка системы на элементы и их нумерация).

Сведение заданной внешней нагрузки к узловой.

Формирование матриц жесткости всех элементов системы в локальных системах координат и их преобразование в глобальную систему координат.

Формирование глобальной матрицы жесткости, системы уравнений метода конечных элементов и ее решение.

Определение усилий в элементах от действия узловой нагрузки.

Определение окончательных значений усилий в элементах путем сложения решений задач 1 и 2.

Далее подробно рассмотрим все эти этапы.

Дата добавления: 2014-01-04