- •Методические указания для выполнения расчетно-графических работ по курсу «метод конечных элементов»
- •Общие указания
- •Расчет стержня на растяжение-сжатие
- •Исходные данные
- •Варианты схем
- •Пример расчета
- •Составление матриц жесткости конечных элементов
- •Составление общей матрицы жесткости
- •Составление уравнения равновесия системы
- •Определение неизвестных перемещений узлов
- •Определение реакций опор
- •Определение напряжений в элементах
- •Расчет фермы
- •Исходные данные
- •Пример расчета
- •1. Составление матриц жесткости конечных элементов
- •2. Составление общего уравнения равновесия
- •3. Решение системы линейных алгебраических уравнений
- •4. Определение реакций в опорах
- •5. Проверка
- •6. Определение напряжений в элементах
- •Расчет балки
- •Варианты схем
- •Пример расчета
- •1. Составление матрицы жесткости отдельных элементов
- •2. Составление глобального уравнения равновесия
- •3. Определение перемещений в узлах
- •4. Определение реакций в опорах
- •5. Проверка полученных результатов
- •Приложения приложение а
- •Приложение б решение систем линейных алгебраических уравнений на базе програмного комплекса mathcad
- •1.Знакомство с программным комплексом MathCad
- •1.1 Основные элементы интерфейса
- •1.2 Основные принципы работы с программой
- •2. Решение систем линейных алгебраических уравнений
- •Приложение в расчет балки методом конечных элементов в программном комплексе ansys
- •Рекомендуемая литература
Варианты схем
1 |
|
2 |
|
3 |
|
4 |
|
5 |
|
6 |
|
Пример расчета
Дано
A1=15см2; A2 =10см2; A3 =5см2; l1 = 1м; l2=2м; l3=3м; P=15 кН; q=10 кН/м; E=2,1·1011Н/м2
Решение
Составление матриц жесткости конечных элементов
В общем виде матрица жесткости для стержневого конечного элемента имеет вид:
где Е – модуль упругости; А – площадь поперечного сечения элемента; L – длина элемента.
Выполним некоторые подготовительные операции.
Пусть А=5см2, L=1м;
тогда A1=3А; A2 =2А; A3 =А;
L1 = L; L2=2L; L3=3L.
Элемент 1 (1-2)
Элемент 2 (2-3)
Элемент 3 (3-4)
Составление общей матрицы жесткости
Объединим составленные матрицы жесткости отдельных элементов в общую матрицу жесткости.
Правильно составленная матрица жесткости должна быть симметричной относительно главной диагонали.
Составление уравнения равновесия системы
Уравнение равновесия имеет вид:
K·u=F,
где u – вектор перемещений узлов; F – вектор внешних нагрузок, приложенных к каждому узлу.
,
т.к. 1-й и 4-й узел находятся в опорах, их
перемещения равны нулю, тогда вектор
перемещений примет вид
.
Вектор
внешних нагрузок в общем виде равен:
.
Для корректного использования, заданную нагрузку необходимо преобразовать следующим образом:
в опорах действуют реакции R1 и R4, которые тоже являются внешней нагрузкой;
распределенная нагрузка действует на протяжении всего 1-го элемента, однако МКЭ учитывает только нагрузку, приложенную в узлах, поэтому всю нагрузку q следует равномерно распределить по узлам этого элемента.
С учетом замечаний схема нагружения имеет вид:
Вектор
внешних нагрузок запишется следующим
образом:
Запишем окончательно уравнение равновесия
.
Определение неизвестных перемещений узлов
Поскольку перемещение первого и четвертого узлов равны нулю, можно вычеркнуть соответствующие строчки и столбцы матрицы жесткости:
.
Или в виде системы:
Подставим численные значения параметров стержня и приложенной нагрузки (все значения подставляются в СИ).
,
.
Разделим левую и правую часть каждого уравнения системы на жесткость EA/L.
Умножим второе уравнение на 4.
И сложим эти уравнения
4,32u3=0,95·10-4.
Отсюда найдем u3
.
Перемещение u2 найдем из второго уравнения
.
Определение реакций опор
Подставив найденные перемещения в первое уравнение равновесия, найдем неизвестную реакцию R1.
,
откуда
.
Аналогично найдем реакцию в 4-м узле R4.
.
В качестве проверки воспользуемся уравнением статики и найдем сумму проекций сил на ось Х.
ΣХ=0; R1 – ql1 + P + R4= -4,23 – 10·1 + 15 – 0,77=0.
Определение напряжений в элементах
;
;
.
Расчет фермы
Для заданной фермы необходимо определить перемещения узлов, реакции опор, и напряжения в стержнях.
Исходные данные принимаются по двум последним цифрам номера зачетной книжки. Для всех вариантов принять Е=2,1·1011Н/м2; А=20см2.
1 |
|
|
2 |
|
3 |
|
|
4 |
|
5 |
|
|
6 |
|
7 |
|
|
8 |
|
9 |
|
|
0 |
|
