- •1 Матриці і дії над ними
- •1.1 Основні поняття
- •1.2 Дії над матрицями
- •1.2.1 Складання і віднімання матриць
- •1.2.2 Множення матриці на число
- •1.2.3 Множення вектора на матрицю
- •1.2.4 Перемноження векторів
- •1.2.5 Множення матриці на матрицю
- •1.3 Системи лінійних рівнянь алгебри
- •1.3.1 Загальні відомості про системи лінійних рівнянь
- •1.3.2 Матричний метод розв’язання систем лінійних рівнянь
- •1.3.3 Метод Гауса
- •1.4 Лабораторна робота 1
- •2 Визначення внутрішніх зусиль та напружень у конструкціях, що перебувають в одноосному напруженому стані
- •2.1 Загальне уявлення про метод кінцевих елементів
- •2.2 Розтягання (стикання) призматичних стрижнів. Визначення напружень. Розрахунок на міцність
- •2.3 Лабораторна робота 2
- •2.4 Лабораторна робота 3
- •2.5 Поняття про кручення. Побудова епюри крутних моментів. Напруження і деформації при крученні круглого вала
- •2.6 Лабораторна робота 4
- •3 Згинання. Побудова епюр
- •3.1 Поняття про згинання балки. Види опор й опорні реакції. Внутрішні зусилля в балці, їх визначення і правила знаків
- •3.2 Лабораторна робота 5
- •3.3 Лабораторна робота 6
- •4 Складний опір
- •4.1 Поняття про складний опір
- •4.1.1 Складне і косе згинання
- •4.1.2 Згинання з крученням круглих валів
- •4.2 Лабораторна робота 7
- •4.3 Лабораторна робота 8
- •5.1.3 Плоский деформований стан
- •5.1.4 Зв'язок між деформаціями і переміщеннями
- •5.1.5 Зв'язок між напруженнями і переміщеннями
- •5.2 Особливості осесимметричної задачі теорії пружності
- •5.3 Лабораторна робота 9
- •5.4 Лабораторна робота 10
- •6 Розрахунки в спеціальних cae системах на прикладі пакету cosmos/m
- •6.1 Загальні відомості
- •6.1.1 Основний екран і головне меню
- •6.1.2 Алгоритм ке-розрахунку в cosmos/m
- •6.1.3 Геометричні примітиви в geostar
- •6.1.4 Властивості елементів
- •6.1.5 Параметрична генерація кe-сітки
- •6.1.6 Автоматична генерація одно- і двомірних ке-сіток
- •6.2 Команди cosmos/m
- •6.2.1 Меню geometry
- •6.2.2 Меню Meshing
- •6.2.3 Меню Propsets
- •6.2.4 Меню loadsbc
- •6.2.5 Меню Analysis
- •6.3 Лабораторна робота 11
- •6.4 Лабораторна робота 12
- •6.5 Лабораторна робота 13
- •6.6 Лабораторна робота 14
- •6.7 Лабораторна робота 15
- •Глава 7 обчислювальна гідродинаміка
- •7.1. Поняття про рівняння руху в’язких середовищ
- •7.1.1. Рівняння імпульсу в консервативній формі
- •7.1.2 Формулювання початкових і граничних умов
- •7.1.3 Рівняння перенесення вихору
- •7.1.4 Граничні умови для вихорового рішення
- •7.1.5 Вихоровий розв’язок для закругленого штампу
- •7.1.6 Метод маркерів
- •7.6 Лабораторна робота 16
- •7.4 Лабораторна робота 17
- •139/2010. Підп. До друку Формат 60х84/16.
- •84313, М. Краматорськ, вул. Шкадінова, 72
3.2 Лабораторна робота 5
Для двоопорної балки побудувати епюри поперечних сил і згинальних моментів, обрати прямокутний поперечний переріз із співвідношенням сторін h=1,5b, якщо []=160 МПа.
Вихідні дані вибрати з таблиць 3.1 та 3.2 і рисунка 3.10.
|
|
|
|
|
|
|
|
|
|
Рисунок 3.10
|
|
|
|
|
|
|
|
|
|
Рисунок 3.10, аркуш 2
|
|
|
|
|
|
|
|
|
|
Рисунок 3.10, аркуш 3
Таблиця 3.1
Величина |
Значення величини відповідно до першої цифри номера варіанта |
|||||||||
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
|
M, кН∙м |
12 |
14 |
16 |
18 |
20 |
18 |
16 |
14 |
12 |
16 |
P, кН |
5 |
6 |
7 |
8 |
5 |
6 |
7 |
8 |
7 |
6 |
q, кН/м |
3 |
2 |
4 |
5 |
4 |
3 |
2 |
5 |
3 |
4 |
Таблиця 3.2
Величина, м |
Значення величини відповідно до другої цифри номера варіанта |
|||||||||
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
|
a |
2,0 |
3,0 |
2,4 |
2,8 |
1,6 |
1,8 |
2,4 |
1,8 |
1,6 |
1,4 |
b |
1,6 |
1,8 |
1,2 |
3,2 |
3,4 |
3,6 |
3,8 |
4,0 |
2,6 |
2,4 |
c |
1,2 |
1,4 |
1,6 |
1,8 |
1,5 |
2,2 |
2,4 |
2,6 |
2,5 |
2,7 |
Указівки до виконання лабораторної роботи 5
Поділяємо стрижень на кінцеві елементи. Нумеруємо отримані елементи і вузли (зручно розділяти таким чином, щоб кінцеві елементи співпадали з ділянками стрижня).
Обчислюємо матрицю жорсткості кожного з елементів, де Jхi – момент інерції перерізу i-го елемента, i – довжина i-го елемента, Е – модуль пружності.
Формуємо матрицю жорсткості конструкції як ансамбль матриць жорсткості кінцевих елементів, з яких складається конструкція. Її порядок дорівнює кількості вузлів, яка помножена на 2.
Записуємо систему рівнянь виду:
де К – матриця жорсткості конструкції в цілому;
– вектор переміщень (лінійних та кутових) усіх вузлів;
– вектор вузлових навантажень.
На ділянках, де діє розподілене навантаження, вважається, що його дія рівносильна впливу вузлових навантажень: у лівому вузлі Pi = ql/2, Mi = ql2/12, а у правому PJ = ql/2, MJ = – ql2/12 (рис. 3.11, б).
Рисунок 3.11
5 З розв’язання системи рівнянь знаходяться переміщення і кути повороту вузлів.
6 Визначити величини поперечних сил Рi і згинальних моментів Mi:
Побудувавши графік зміни поперечних сил і згинальних моментів уздовж стрижня, отримаємо епюри поперечних сил і згинальних моментів.
З умови міцності визначити необхідний момент опору перерізу:
,
де
– осьовий момент опору перерізу елемента,
[] – допустимі нормальні напруження.
Визначити необхідні розміри поперечного перерізу, ураховуючи, що
.
Приклад виконання лабораторної роботи 5
Для балки, показаної на рисунку 3.12, а, побудувати епюри поперечних сил і згинальних моментів. Підібрати прямокутний поперечний переріз із співвідношенням сторін h=2 b, якщо []=160 МПа.
Розв’язання. Нумеруємо вузли і елементи (номера елементів обведені колом). Обчислюємо матриці жорсткості елементів. До елемента 1 входять вузли 1 і 2, його довжина – 3 м:
.
Елемент 2: входять вузли 2 і 3, його довжина – 2 м:
.
Тепер формуємо матрицю жорсткості конструкції. Оскільки кожен вузол має два степені вільності, а вузлів 3, то розмірність матриці жорсткості конструкції 6х6:
При записі вектора навантажень враховуємо, що на перший елемент діє розподілене навантаження. Його дію замінюємо вузловими навантаженнями, як було сказано вище. У даному випадку знак розподіленого навантаження беремо з мінусом, тому що воно направлене вниз.
Ураховуємо граничні умови:
V1=0, 1=0, V3=0.
Після їх урахування отримаємо систему рівнянь:
Після розв’язання одержимо:
Визначаємо зусилля на кінцях елементів.
Розглядаємо 1-й елемент:
Прикладаємо отримані вузлові сили до вузлів 1-го елемента (рис. 3.13,а).
Для контролю необхідно виконати перевірки:
∑Pу=0
20,84+9,16-10·3=0,
∑Ма=0
14,5-20,84·3+10·3·1,5+3,2=0,
14,5-62,52+45+3,2=0.
А,б,в вінести перед названием!
Розглядаємо 2-й елемент:
Прикладаємо отримані вузлові сили до вузлів 2-го елемента (рис. 3.14,а).
А,б,в вінести перед названием!
Зістикувавши епюри по ділянках, отримуємо остаточні епюри поперечних сил (рис. 3.12, б) і згинальних моментів (рис. 3.12, в):
,
враховуючи, що
.
Звідки
.
