Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

2545

.pdf
Скачиваний:
4
Добавлен:
22.06.2024
Размер:
13.44 Mб
Скачать

После определения перемещений Z внутренние усилия S могут быть найдены по формуле (9.41), в которую нужно подставить значение E из выражения (9.40):

S [K]

 

[K]AT Z .

(9.46)

E

Заметим, что в стержнях, загруженных внеузловой нагрузкой, к усилиям (9.41) нужно добавить те усилия S0, которые возникают в стержне от уравновешенной системы сил (см. рис. 9.7, б). Таким образом, для загруженных стержней будем иметь:

S [K]AT Z S0 .

(9.47)

В методе перемещений достаточной проверкой правильности полученного решения является статическая проверка, заключающаяся в проверке выполнения условий равновесия узлов системы. Поэтому для проверки результатов расчета нужно подставить полученный вектор S усилий в исходное уравнение равновесия (9.39), т.е. умножить его слева на статическую матрицу A и результат сравнить с вектором узловых сил P, на действие которых производится расчет:

AS P .

(9.48)

Втом случае, если конструкция рассчитывается на s вариантов загружения, то P и S будут уже не векторами, а матрицами, содержащими

по s столбцов.

Заметим, что получение обратной матрицы K-1 требует большего числа вычислений, чем прямое решение системы уравнений (9.44). Доказано, что имеет смысл обращать матрицу A лишь при большом числе вариантов загружения, т.е. когда отношение числа загружений s к порядку п системы уравнений (9.44) превышает 0,6: s / n > 0,6.

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

Заметим также еще раз, что матричное уравнение (9.44) полностью тождественно системе канонических уравнений классического варианта

метода перемещений: каждый элемент kij матрицы жесткости системы K равен коэффициенту rij, представляющему собой реактивное усилие в связи i по направлению перемещения zi от смещения связи j на единицу.

9.7. Примеры расчета

Пример №1. В качестве первого примера рассмотрим расчет рамы (рис. 9.13) на действие заданной узловой нагрузки.

241

Рис.9.13

Рассматриваемая рама в общем случае (с учетом продольных перемещений) пять раз кинематически неопределима, т.е. степень свободы ее п = 5 (см. рис.9.15). Стержни рамы содержат семь искомых внутренних усилий (рис. 9.14), которыми являются изгибающие моменты в жестких узлах и продольные силы, в связи с чем т = 7. m n 2 , следовательно, данная рама 2 раза статически неопределима.

Рис. 9.14

242

Рис. 9.15

На рис.9.15,а и рис. 9.14 приведены диаграммы P–z и S–e, на которых показаны номера и принятые за положительные направления перемещения узлов рамы, а также сечения, в которых вычисляются моменты Sj, и их положительные направления.

Составим пять уравнений равновесия, соответствующие пяти возможным перемещениям узлов рамы (рис. 9.15,б):

Узел I. FKX 0, Q1

N2

P 0;

M H 1

N2 P 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4a

 

 

 

 

 

 

 

FKY 0, N1 Q2

0; N1

 

M K 2

 

0.

 

 

 

 

3a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Узел II. FKX

0,

N2 Q3 cos N3

sin 0,

 

 

 

N2

M K 3 M H 3 0.8 N3 0.6 0.

 

 

 

 

 

 

 

 

5a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

FKY

0, Q2

N3 cos Q3 sin 0,

 

 

 

 

M K 2 N3 0.8 M K 3 M H 3

0.6 0.

 

 

 

 

 

 

 

3a

 

 

 

 

5a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M 0, M K 2 M H 3 M 0.

 

 

 

 

 

 

 

 

 

 

Запишем уравнение равновесия в матричной форме:

 

 

 

 

 

 

 

 

AS P 0,

 

 

 

 

 

 

 

 

 

(*)

 

 

 

 

 

S s1T s2T s3T ;

 

 

 

 

 

 

 

s N M

T ; s

2

N

M

K 2

T ; s

N

M

H 3

M

K 3

T ;

1

1

H1

 

2

 

 

 

3

2

 

 

 

P P 0 0 0 M T .

243

Матрица коэффициентов уравнений равновесия, или статическая матрица:

 

0

0.25

/ a

1

0

0

0

 

0

 

 

1

0

 

0

0.333 / a

0

0

 

0

 

 

 

 

 

0

0

 

1

0

0.6

0.16

/ a

0.16

A

 

/ a .

 

0

0

 

0

0.333 / a

0.8

0.12

/ a

0.12 / a

 

0

0

 

0

1

0

1

 

0

 

 

 

 

 

Для формирования матрицы внутренней жесткости запишем матрицы жесткости каждого стержня рамы:

EF k1 4a0

0

3EF

4a

 

 

 

EF

 

,

k2

 

3a

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

EF

0

0

 

 

 

 

 

 

5a

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

k3

 

0

4EF

2EF .

3EF

 

 

 

 

5a

5a

 

3a

 

 

 

 

 

2EF

4EF

 

 

 

 

 

0

5a

5a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

244

Матрица внутренней жесткости может быть теперь представлена в виде:

EF

0

0

0

0

0

0

 

 

4a

 

 

 

 

 

 

 

 

 

3EI

 

 

 

 

 

 

0

0

0

0

0

0

 

 

 

4a

 

 

 

 

 

 

 

 

 

EF

 

 

 

 

 

0

0

0

0

0

0

 

3a

 

 

 

 

 

 

 

 

 

 

0

0

3EI

0

0

0

 

k 0

3a

.

 

 

 

 

EF

 

 

 

 

 

0

0

0

0

0

 

0

5a

 

 

 

 

 

 

4EI

 

 

 

 

 

 

 

 

 

0

0

0

0

0

2EI

 

5a

5a

 

 

 

 

 

 

 

 

 

 

 

 

2EI

4EI

 

 

 

0

0

0

0

 

0

5a

5a

 

 

 

 

 

 

 

 

Здесь k [K]..

Уравнения равновесия перепишем в виде:

A[K]AT Z P ,

или

KZ P .

Решим систему уравнений:

 

Z K 1P (A[K ]AT ) 1 P.

 

 

 

 

При а=1 м, Е=21011 Н/м2, Р=10000

Н, М=10000

Н м, F=20 10-4 м2,

I=11 10-6 м4 получим:

 

 

 

 

 

 

 

 

 

1,629 10

5

0,017

0,013 2,032

10

3

T

.

Z 0,017

 

 

 

Внутренние усилия от узловой нагрузки имеют величину:

 

 

7117 8221 4892 6296

 

 

 

S kAT Z 1629

14892 13100 T .

Эпюра изгибающих моментов показана на рис. 9.16.

245

Рис. 9.16

Пример №2. Рассмотрим пример расчета двухпролетной рамы матричным методом перемещений.

Построить эпюру изгибающих моментов, используя матричный метод перемещений. Расчетная схема рамы показана на рис. 9.17.

Рис. 9.17

Пусть Е=2 1011Н/м2, Р=16 кН, q=4000 кН/м, F=20 10-4 м2, I=11 10-6 м4 .

Грузовая эпюра, т.е. эпюра от пролетных нагрузок S0 для отдельных стержней, показана на рис. 9.18.

246

Рис. 9.18

На рис. 9.19 представлена схема узловых перемещений системы (диа-

грамма P z ). При этом вектор P [5000 9000

0 15000 18000 0 0]T .

Рис. 9.19

На рис. 9.20 изображена схема внутренних усилий в стержнях (диаграмма S e ).

Рис. 9.20

247

На рис. 9.19 показаны узлы системы с приложенными к ним силами. Составим уравнения равновесия узлов рамы (рис. 9.21).

Рис. 9.21

Узел А:

Fkx 0, N2 Q1 5 0;

Fky 0, N1 Q2 9 0.

Узел Б:

Fkx 0, N2 N4 Q3 0;

Fky 0, N3 Q2 Q4 15 0;

M 0, M2 M3н M4 18 0.

Узел С:

Fkx 0, N4 Q5 0;

Fky 0, N5 Q4 0.

248

Заменим поперечные силы концевыми моментами (см. рис. 9.20):

Q1 M1 / 6; Q2 M2 / 6;

Q3 (M3к M3н) / 6;

Q4 M4 / 6;

Q5 M5 / 6.

Имеем 11 неизвестных внутренних усилий:

N1, M1, N2 , M2 , N3, M3н, M3k , N4 , M4 , N5, M5.

Расположим усилия S1...S11 в следующем порядке:

1

2

3

4

5

6

7

8

9

10

11

N1

M1

N2

M2

N3

M3н

M3k

N4

M4

N5

N5

Составим таблицу, соответствующую матричному уравнению равно-

весия:

AS P 0.

Или для программного комплекса Mathcad:

 

0

0.16667

1

0

0

0

0

0

0

0

0

 

 

 

1

0

0

0.16667

0

0

0

0

0

0

0

 

 

 

 

 

 

0

0

1

0

0

0.16667

0.16667

1

0

0

0

 

A1

 

0

0

0

0.16667

1

0

0

0

0.25

0

0

 

 

 

0

0

0

1

0

1

0

0

1

0

0

 

 

 

 

 

 

0

0

0

0

0

0

0

1

0

0

0.16667

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

0

0

0

0

0

0

0.25

1

0

 

249

Матрица внутренней жесткости стержневой системы k=[K] имеет вид:

Как известно, в данном примере Е=2 1011Н/м2, F=20 10-4 м2, I=11 10-6 . Решим систему матричных уравнений относительно узловых

перемещений:

Z A k AT 1 P .

250

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]