2286
.pdf
Задачи для самостоятельного решения
Рассчитать раму (рис. 7.12) методом перемещений.
Рис. 7.12. Задачи для самостоятельного решения
81
8. МАТРИЧНАЯ ФОРМА РАСЧЕТА СТЕРЖНЕВЫХ СИСТЕМ МЕТОДОМ ПЕРЕМЕЩЕНИЙ
8.1. Расчет рамы с узловой нагрузкой матричным методом перемещений
Пример№1.
Построить эпюру изгибающих моментов для рамы с узловой нагрузкой (рис. 8.1а), используя матричный метод перемещений. Характеристики материала и поперечного сечения стержней: модуль упругости Е=2·1011 Н/м2, площадь поперечного сечения F=20·10-4 м2, момент инерции I=11·10-6 м4.
Рис. 8.1. Расчетная схема рамы и схема узловых перемещений
Рассматриваемая рама в общем случае (с учетом продольных перемещений) пять раз кинематически неопределима, т.е. степень ее свободы п = 5 (рис. 8.1, б). Схемы внутренних усилий в стержнях и узлах рамы показаны на рис. 8.2, 8.3.
Рис. 8.2. Схема внутренних усилий в стержнях фермы
82
Рис. 8.3. Вырезанные узлы
Составим уравнения равновесия вырезанных узлов (рис. 8.3). Узел А:
Fkx 0 : Q1 N2 10000 0;
Fky 0 : N1 Q2 0 .
Узел B:
Fkx 0 : N2 Q3 0 ;
Fky 0 : Q2 N3 0;
M узл 0 : M2 M3н 5000 0 .
Выразим поперечные силы через изгибающие моменты, составляя для каждого стержня уравнение моментов относительно точки начала (рис. 8.2):
M |
1 |
Q 4 0 |
|
Q |
M1 |
; |
|
||||||||
|
|
|
|
||||||||||||
|
|
|
1 |
|
|
|
|
1 |
4 |
|
|
||||
M |
2 |
Q 3 0 |
Q |
M2 |
; |
|
|
||||||||
|
|
|
|||||||||||||
|
|
|
2 |
|
|
|
2 |
|
3 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
M3к M3н . |
|||
M |
3н |
M |
3к |
Q |
4 0 |
|
Q |
|
|||||||
|
|
|
|
3 |
|
|
3 |
|
|||||||
4
Перепишем уравнения равновесия вырезанных узлов с учетом полученных выражений:
M41 N2 10000 0 ;
N1 M32 0 .
N2 M3к M3н 0 ; 4
M32 N3 0 ;
M2 M3н 5000 0 .
Уравнения равновесия в матричной форме имеют вид:
AS P 0.
Тогда вектор внутренних усилий
S N1 |
M1 N2 M2 N3 M3н M3к T ; |
|
83 |
вектор внешней нагрузки
P 10000 0 0 0 5000 T ;
статическая матрица (матрица коэффициентов при неизвестных)
N1 |
M1 |
N2 |
M 2 |
N3 |
|
M3н M3к |
||||
|
0 |
1 |
1 |
0 |
0 |
0 |
0 |
|
1 |
|
|
4 |
|
||||||||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
1 |
0 |
0 |
0 |
0 |
0 |
|
2 |
||
|
|
|
|
3 |
|
|
|
|
|
|
A |
0 |
0 |
1 |
0 |
0 |
1 |
|
1 |
|
|
|
4 |
4 |
|
3 |
||||||
|
|
|
|
|
|
|
|
|||
|
0 |
0 |
0 |
1 |
1 |
0 |
0 |
|
|
|
|
|
|
||||||||
|
0 |
0 |
0 |
3 |
0 |
1 |
0 |
|
4 |
|
|
1 |
|
|
|||||||
|
|
|
|
|
|
|
|
|
5 |
|
Для формирования матрицы внутренней жесткости запишем матрицы жесткости каждого стержня рамы:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
EF |
|
0 |
|
|
EF |
|
|
|
|
|
EF |
|
|
|
|
|
|
4 |
|
|||
|
0 |
|
|
|
0 |
|
|
|
|
|
||||||||
k |
|
4 |
|
, |
k |
|
|
3 |
|
, |
k |
|
|
0 |
4EI |
|||
|
|
|
|
|
|
|
|
|
||||||||||
1 |
|
0 |
3EI |
|
|
2 |
|
0 |
3EI |
|
|
3 |
|
|
|
4 |
||
|
|
4 |
|
|
|
|
|
3 |
|
|
|
|
|
|
|
2EI |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Тогда матрица внутренней жесткости всей системы:
2EI4 .
4EI
|
EF |
0 |
0 |
0 |
0 |
|
0 |
|
0 |
|
|
4 |
|
|
|
||||||
|
3EI |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
0 |
0 |
0 |
0 |
|
0 |
|
0 |
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
EF |
|
|
|
|
|
|
|
|
0 |
0 |
0 |
0 |
|
0 |
|
0 |
|
|
3 |
|
|
||||||||
|
|
|
|
|||||||
|
0 |
0 |
0 |
3EI |
0 |
|
0 |
|
0 |
|
k |
3 |
|
|
|
||||||
|
|
|
|
EF |
|
|
|
|
|
|
|
0 |
0 |
0 |
0 |
|
0 |
|
0 |
|
|
|
4 |
|
|
|
||||||
|
|
|
|
|
4EI |
|
|
|
||
|
|
|
|
|
|
|
||||
|
0 |
0 |
0 |
0 |
0 |
|
2EI |
|||
|
|
4 |
4 |
|
||||||
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
2EI |
4EI |
|
|
|
0 |
0 |
0 |
0 |
0 |
|
|
|||
|
4 |
|
4 |
|
||||||
|
|
|
|
|
|
|
|
|
||
84
С помощью системы компьютерной алгебры MathCad определим перемещения и внутренние усилия:
Z A k AT 1 P ; S k AT Z .
Получим следующие значения:
|
|
|
|
|
|
2613,216 |
|
|
|
0,0252 |
|
|
|
10381,719 |
|
|
|
2,6132 10 5 |
|
|
|
|
|
|
|
|
|
|
7404,57 |
|
|
Z |
|
0,0251 |
|
; |
|
7839,647 |
|
|
|
S |
. |
||||
|
2,6132 10 5 |
|
|
|
2613,216 |
|
|
|
|
3,5809 10 3 |
|
|
|
12839,647 |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
16778,634 |
|
|
|
|
|
|
|
|
|
Пополученнымзначениямстроимэпюруизгибающихмоментов(рис. 8.4).
Рис. 8.4. Эпюра изгибающих моментов
8.2. Расчет рамы с пролетной нагрузкой матричным методом перемещений
Пример №1.
Построить эпюру изгибающих моментов для рамы с пролетной нагрузкой (рис. 8.5), используя матричный метод перемещений. Характеристики материала и поперечного сечения стержней: модуль упругости
Е=2·1011 Н/м2, площадь поперечного сечения F=20·10-4 м2, момент инерции
I=11·10-6 м4.
С помощью таблицы метода перемещений строим эпюру S0 от внешней нагрузки (т.е. грузовую эпюру) и приводим пролетную нагрузку к узловой
(рис. 8.6).
На рис. 8.7 представлена схема узловых перемещений системы, схемы внутренних усилий в стержнях и узлах рамы показаны на рис. 8.8, 8.9.
85
Рис. 8.5. Расчетная схема рамы
Рис. 8.6. Эпюра моментов от пролетной нагрузки
Рис. 8.7. Схема узловых перемещений
86
Рис. 8.8. Схема внутренних усилий в стержнях фермы
Рис. 8.9. Вырезанные узлы
Составим шесть уравнений равновесия (рис. 8.9), соответствующих шести возможным перемещениям рамы (рис. 8.7):
Узел А:
Fkx 0 : N1 0 .
Узел В:
Fkx 0 : N1 Q2 N3 0 ;
Fky 0 : Q1 N2 Q3 15000 0 ;
M узл 0 : M1 M2к M3 18000 0 .
Узел С:
Fkx 0 : N3 Q4 5000 0 ;
Fky 0 : Q3 N4 9000 0 .
Выразим поперечные силы через изгибающие моменты, составляя для каждогостержняуравнениемоментовотносительноточкиначала(рис. 8.8):
M |
Q 4 0 |
Q |
M1 |
; |
|
||||
1 |
1 |
1 |
4 |
|
|
|
|
|
87
M |
2н |
M |
2к |
Q |
6 0 |
|
Q |
|
M2к M2н ; |
||||
|
|
2 |
|
|
2 |
|
6 |
||||||
|
|
|
|
|
|
|
|
|
M3 |
|
|
||
M |
3 |
Q 6 0 |
Q |
; |
|
||||||||
|
|
|
|||||||||||
|
|
3 |
|
|
|
3 |
6 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|||
M |
4 |
Q 6 0 |
Q |
M4 |
. |
|
|||||||
|
|
||||||||||||
|
|
4 |
|
|
|
4 |
6 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|||
Перепишем уравнения равновесия вырезанных узлов с учетом полученных выражений:
N1 0 ;
N1 M2к M2н N3 0 ; 6
M41 N2 M63 15000 0;
M1 M2к M3 18000 0
N3 M64 5000 0 ;
M63 N4 9000 0 .
Уравнения равновесия в матричной форме имеют вид:
AS P 0.
Тогда вектор внутренних усилий
S N1 M1 N2 M2н M2к |
N3 M3 N4 M4 T ; |
вектор внешней нагрузки
P 0 0 15000 18000 5000 9000 T ;
статическая матрица (матрица коэффициентов при неизвестных)
|
N1 |
M1 |
N2 M 2н M 2к |
N3 M3 |
|
N4 |
|
M 4 |
||||
1 0 0 0 0 |
0 |
0 |
0 |
0 1 |
||||||||
|
|
|
|
|
1 |
1 |
|
|
|
|
|
|
1 |
|
0 |
0 |
1 |
0 |
0 |
0 |
2 |
||||
|
|
|
|
|
6 |
6 |
|
|
|
|
|
|
|
0 |
|
1 |
1 0 0 |
0 |
1 |
0 |
0 |
|
3 |
||
|
|
|
||||||||||
A |
0 |
|
4 |
0 |
0 |
1 0 |
6 |
0 |
0 |
|
|
|
|
|
1 |
1 |
|
4 |
|||||||
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
0 |
|
0 |
0 |
0 |
0 |
1 |
0 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
6 |
|
5 |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
0 0 |
0 |
0 |
0 |
6 |
1 |
0 |
|
||||
|
6 |
|||||||||||
88
Для формирования матрицы внутренней жесткости запишем матрицы жесткости каждого стержня рамы:
|
|
|
|
|
|
|
|
|
|
EF |
|
0 |
|
0 |
|
|
|||
|
|
3EF |
|
|
|
|
|
|
6 |
|
|
|
|
|
|||||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|||||||
k |
|
|
4 |
|
|
, k |
|
|
0 |
|
4EI |
|
|
, |
|||||
|
|
|
|
|
|
|
|
2EI |
|||||||||||
1 |
|
|
0 |
|
3 3EI |
|
2 |
|
|
|
|
|
|
6 |
|
6 |
|
|
|
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
2EI |
4EI |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
0 |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
6 |
|
6 |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
3EF |
|
0 |
|
|
|
|
|
EF |
0 |
|
|
|
||
|
|
k |
|
|
6 |
|
|
|
|
k |
|
|
|
6 |
|
|
|
||
|
|
3 |
|
3 |
|
|
, |
4 |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
0 |
3EI |
|
|
|
|
0 |
3EI |
|
|
||||||
|
|
|
|
|
|
6 |
|
|
|
|
|
|
|
6 |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Тогда матрица внутренней жесткости всей системы:
|
3EF |
0 |
0 |
|
0 |
|
0 |
0 |
0 |
0 |
|
4 |
|
|
|||||||
|
|
|
|
|
|
|
|
|
||
|
|
9EI |
|
|
|
|
|
|
|
|
|
0 |
0 |
|
0 |
|
0 |
0 |
0 |
0 |
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
EF |
|
|
|
|
|
|
|
|
0 |
0 |
|
0 |
|
0 |
0 |
0 |
0 |
|
6 |
|
|
||||||||
|
|
|
||||||||
|
0 |
0 |
0 |
4EI |
|
2EI |
0 |
0 |
0 |
|
|
||||||||||
|
|
|
|
|
6 |
|
6 |
|
|
|
|
0 |
0 |
0 |
|
2EI |
4EI |
0 |
0 |
0 |
|
k |
6 |
|
6 |
|||||||
|
|
|
|
|
|
3EF |
|
|
||
|
|
|
|
|
|
|
|
|
||
|
0 |
0 |
0 |
|
0 |
|
0 |
0 |
0 |
|
|
|
|
|
|
|
|
|
6 |
|
|
|
|
|
|
|
|
|
|
|
9EI |
|
|
0 |
0 |
0 |
|
0 |
|
0 |
0 |
0 |
|
|
|
6 |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
EF |
|
0 |
0 |
0 |
|
0 |
|
0 |
0 |
0 |
|
|
|
6 |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
0 |
0 |
0 |
|
0 |
|
0 |
0 |
0 |
0 |
|
|
|
||||||||
С помощью системы компьютерной алгебры MathCad перемещения и внутренние усилия:
Z A k AT 1 P ; S k AT Z .
0
0
0
0
0
0
0
0
3EI
6
определим
89
Получим следующие значения:
|
|
|
|
|
|
0 |
|
|
0,031 |
|
|
|
|
3208,281 |
|
|
|
|
|
||||
|
0,031 |
|
|
|
|
15405,558 |
|
|
|
|
|
|
11894,981 |
|
|
|
2,3108 10 4 |
; |
|
|
|||
Z |
7,0591 10 4 |
|
S 12412,647 |
. |
|||
|
|
|
|
4051,271 |
|
||
|
0,031 |
|
|
|
|
|
|
|
|
|
|
|
2379,072 |
|
|
|
1,4095 10 |
4 |
|
|
|
9396,512 |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
5692,371 |
|
|
|
|
|
|
|
|
|
По полученным значениям строим эпюру изгибающих моментов S
(рис. 8.10).
Рис. 8.10. Эпюра изгибающих моментов от узловой нагрузки
Окончательная эпюра моментов находится путем сложения двух эпюр: Mок S0 S . Эта эпюра показана на рис. 8.11.
Рис. 8.11. Окончательная эпюра моментов
90
