- •Методические указания для выполнения расчетно-графических работ по курсу «метод конечных элементов»
- •Общие указания
- •Расчет стержня на растяжение-сжатие
- •Исходные данные
- •Определение неизвестных перемещений узлов
- •Определение реакций опор
- •Определение напряжений в элементах
- •Расчет фермы
- •Варианты схем
- •Пример расчета
- •Составление матриц жесткости конечных элементов
- •Составление общей матрицы жесткости
- •Составление уравнения равновесия системы
- •Исходные данные
- •Пример расчета
- •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 |
схема 7 |
схема 8 |
схема 9 |
схема 10 |
схема 11 |
схема 12 |
схема 13 |
схема 14 |
схема 15 |
схема 16 |
схема 17 |
схема 18 |
схема 19 |
схема 20 |
Пример расчета
Д
ля
заданной балки определить перемещения
и углы поворота в узлах, а также опорные
реакции.
Дано: L1=1 м; L2=2 м; L3=3 м; q=5кН/м; М=20кНм; Р=10кН; Е=2,1·1011Н/м2; I=2000см4.
Решение
1. Составление матрицы жесткости отдельных элементов
Элемент 1 (1-2)
Элемент 2 (2-3)
Элемент 3 (3-4)
2. Составление глобального уравнения равновесия
Запишем глобальную матрицу жесткости:
Для того, что бы записать вектор нагрузок в уравнении равновесия, приведем заданную распределенную нагрузку к эквивалентным силам, сосредоточенным в узлах схемы, согласно изображенному правилу.
Заданная схема примет вид:
В общем виде глобальное уравнение равновесия имеет вид:
С учетом внешней нагрузки, и закреплений опорных узлов, это уравнение примет вид:
(1)
3. Определение перемещений в узлах
Из этого уравнения следует вычеркнуть строки и столбцы, соответствующие нулевым перемещениям узлов, т.е. номер 1, 2 и 3. Полученная матрица должна быть симметричной относительно главной диагонали:
Решая эту систему уравнений, получаем вектор неизвестных перемещений.
4. Определение реакций в опорах
Умножим найденный вектор перемещений на 1-ю, 2-ю и 3-ю строку глобального уравнения равновесия (1). Запишем это произведение в матричном виде:
,
где
.
Откуда
(2)
Описанные выше расчеты рекомендуется выполнить в программном комплексе MathCad.
В соответствии с (2) окончательно найдем реакции:
5. Проверка полученных результатов
Для проверки составим уравнение статики.
.
Так как условие равновесия сил соблюдается задача решена верно.
Приложения приложение а
Постановка задачи
Для заданного стержня необходимо определить перемещения узлов, реакции опор и осевые усилия. Сравнить полученные результаты со значениями найденными вручную.
Дано A1=15см2; A2 =10см2; A3 =5см2; l1 = 1м; l2=2м; l3=3м; P=15 кН; q=10 кН/м; E=2,1·1011Н/м2.
В отличие от предыдущей работы мы не будем отдельно создавать геометрическую модель, а затем разбивать ее на конечные элементы. В данном случае гораздо проще сразу создать узлы и объединить их в элементы.
Шаг 1: Задание типа и параметров конечных элементов. В этой задаче мы будем использовать стержневой конечный элемент Link1. Поскольку условие предполагает три ступени, для их описания воспользуемся тремя конечными элементами одного типа и разных диаметров.
|
|
|
|
|
|
Шаг 2: Выбор действительных параметров элементов.
На этом этапе нам необходимо задать три площади поперечного сечения балки.
Main Menu > Preprocessor > Real Constants (действительные параметры) > Add/Edit/Delete
Нажмите Add (добавить к существующему списку наборов параметров).
Cross-sectional area AREA 15e-4 (Площадь поперечного сечения равна А=15·10-4м2).
|
|
Cross-sectional area AREA 10e-4 (Площадь поперечного сечения равна А=10·10-4м2).
|
|
Cross-sectional area AREA 5e-4 (Площадь поперечного сечения равна А=5·10-4м2).
|
|
Шаг 3: Установка свойств материала.
В данном примере используется изотропный материал с постоянными свойствами.
M
ain
Menu
> Preprocessor
(Препроцессор)
> Material
Props
(Свойства материала) >
Material
Models
(модель материала)Дважды щелкните на Structural (Структурный), Linear (Линейный), Elastic (Упругий), Isotropic (Изотропный).
Введите 2.1e11 для EX (Модуль упругости).
Введите 0.3 для PRXY (Коэффициент Пуассона).
Нажмите OK.
Save_DB.
Шаг 4: Создание узлов.
Main Menu (Главное Меню) > Preprocessor (Препроцессор)> -Modeling (Моделирование)- Create (Создать) > - Nodes (Узлы)- In Active CS (В текущей координатной системе).
Заполните поля появившегося диалогового окна
Node number (Номер узла) 1
X, Y, Z Location in active CS 0 0 0
(Координаты 1-го узла Х=0; Y=0; Z=0)
Apply (Применить).
Node number (Номер узла) 2
X, Y, Z Location in active CS 1 0 0
(Координаты 2-го узла Х=1м; Y=0; Z=0)
Apply (Применить).
Node number (Номер узла) 3
X, Y, Z Location in active CS 3 0 0
(Координаты 3-го узла Х=3м; Y=0; Z=0)
Apply (Применить).
Node number (Номер узла) 4
X, Y, Z Location in active CS 6 0 0
(Координаты 4-го узла Х=6м; Y=0; Z=0)
OK.
Save_DB.
Шаг 5: Создание конечных элементов.
При создании конечных элементов будем учитывать то, что каждый элемент имеет свою площадь поперечного сечения. Эта площадь задается при помощи параметра Real Constants.
|
|
|
|
|
|
Шаг 6: Задание граничных перемещений.
Н
а
этом этапе следует задать нулевое
вертикальное и горизонтальное перемещение
1-го и 4-го узла.
Main Menu (Главное меню) > Solution (Решение) > -Loads (Нагрузки) - Apply (Приложенные) > -Structural (Структурный) - Displacement (Перемещения) > On Nodes (В узлах)
Выберите узел 1 и 4.
Нажмите OK в Picking Menu (Меню выбора).
Выберите All DOF (перемещения по всем осям).
Введите 0 в Value (нулевые перемещения).
Нажмите ОК.
SAVE_DB.
Шаг 7: Приложение внешней нагрузки.
Распределенную нагрузку приложим как две сосредоточенные, действующие в узлах элемента и равные половине общей распределенной нагрузки. При этом во втором узле общая нагрузка будет складываться из половины распределенной и сосредоточенной.
Main Menu (Главное Меню) > Solution (Решение) > -Loads (Нагрузки) - Apply (Приложенные) > -Structural (Структурный) – Force/Moment (Сила/Момент) > On Nodes (К узлам).
Выберете 1-й узел.
Нажмите ОК.
Из списка Direction of force/mom (Направление силы/момента) выберите FX.
В поле Force/moment value (Значение силы/момента) введите заданную силу -10*1000*1/2 (сила направленная вправо считается положительной).
Нажмите Apply.
Выберете 2-й узел.
Нажмите ОК.
Из списка Direction of force/mom (Направление силы/момента) выберите FХ.
В поле Force/moment value (Значение силы/момента) введите заданную силу 10*1000.
Нажмите ОК.
SAVE_DB.
Шаг 8: Расчет.
Main Menu (Главное меню) > Solution (Решение) > -Solve (Счет) - Current LS (Текущая задача).
Изучите содержание появившегося окна и закройте его File > Close.
Нажмите OK для запуска счета.
Нажмите кнопку Close в появившемся желтом информационном окне Solution is done (Расчет выполнен).
Шаг 9: Просмотр перемещений в узлах.
M
ain
Menu > General Postproc > List Results
(Список
результатов)
> Nodal
Solution . . . (Расчет
в
узлах)Выберите DOF solution (Расчет степеней свободы) в левом scroll-меню.
Выберите ALL DOFs (Все степени свободы) в правом scroll-меню.
Нажмите OK.
Сравните полученные перемещения со значениями, найденными вручную.
Нажмите File > Close для закрытия окна.
Шаг 10: Просмотр реакций в закреплениях.
Main Menu > General Postproc > List Results (Список результатов) > Reaction Solu (Расчет реакций).
Нажмите OK.
Сравните полученные реакции со значениями, найденными вручную.
Нажмите File > Close для закрытия окна.
Шаг 11: Просмотр напряжений в элементах.
Main Menu > General Postproc > List Results (Список результатов) > Element Solution.
В открывшемся окне выберите LineElem Results.
Нажмите OK.
Величина SAXL показывает напряжения, возникающие в элементе.
Сравните полученные реакции со значениями, найденными вручную.
Нажмите File > Close для закрытия окна.
