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

1846

.pdf
Скачиваний:
2
Добавлен:
16.06.2024
Размер:
3.28 Mб
Скачать

Решив матричное уравнение (1), найдём вектор неизвестных:

 

 

 

 

 

X

A 1

 

p .

 

 

 

 

(4)

Учитывая соотношения (2) и (3), уравнение (4) преобразуем к виду:

 

 

 

0

T

 

 

0

1

0

T

0

 

 

 

 

X Lm

B Lm

Lm

B Lmp P .

(5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вектор изгибающих моментов в пронумерованных сечениях заданной статически неопределимой стержневой системы найдём, используя соотношение:

 

L0mp

 

L0m

 

.

(6)

M

P

X

Таким образом, для численной реализации на ЭВМ матричного метода сил необходимо составить четыре исходные матрицы: L0m , B, L0mp , P .

Окончательная эпюра изгибающих моментов должна удовлетворять как уравнениям равновесия, так и всем кинематическим условиям. В матричной форме это равносильно отысканию полного вектора перемещений Z , значение которого должно быть равно нулю:

 

L0m T B

 

0.

(7)

Z

M

Формула (7) используется для проверки правильности построенного решения.

При действии на статически неопределимую стержневую систему раздельно "t" вариантов внешних воздействий система канонических уравнений имеет вид:

 

 

A X

0 .

 

 

 

(8)

Здесь

X

матрица неизвестных, имеющая "t"

столбцов;

 

A

матрица перемещений

в основной системе от

L0m T B A0mp P

единичных сил X1 X 2

... X n 1;

 

матрица перемещений

от всех "t" вариантов

 

 

нагрузок в основной системе;

 

 

 

 

A0mp

матрица влияния изгибающих моментов в основ-

 

 

ной системе от всех "t" вариантов внешних

 

 

единичных

нагрузок

F 1 , F 1 ,

..., F 2 , F 2 ,

 

 

 

 

1

2

1

2

 

 

..., F t 1, приложенных раздельно;

 

 

 

 

k

 

 

 

 

 

 

P

матрица всех "t" вариантов

внешних

воздей-

 

 

ствий.

 

 

 

 

 

Решение системы (8) имеет вид:

 

 

 

 

 

 

 

X A 1 .

 

 

 

(9)

91

Окончательные значения изгибающих моментов от всех "t" вариантов внешних воздействий получим в виде матрицы:

M 1

M 2

...

M t

 

 

 

1

1

 

1

 

 

 

 

M 1

M 2

...

M t

0

0

(10)

M 2

2

...

2

 

Amp P Lm X.

... ...

...

 

 

 

 

1

2

...

t

 

 

 

M m

M m

M m

 

 

 

Таким образом, при расчёте системы методом сил в матричной форме необходимо:

a)выбрать основную систему метода сил;

b)построить единичные эпюры моментов от единичных воздействий

X1 X 2 ... X n 1;

c)построить единичные грузовые эпюры моментов от грузовой

нагрузки, принимаемой равной единице: F1 1 , F21 ,..., F1 2 , F2 2 ,..., Fk t 1; d) составить матрицы L0m , B, L0mp , P и выполнить их перемножение.

П р и м е р .

Построить эпюры внутренних усилий в статически неопределимой раме (рис. 24) методом сил в матричной форме [28].

Исходные данные: h 3 м, l 6 м, F 3 kH, q 2 kHм.

Р е ш е н и е .

Рама является три раза статически неопределимой. Для получения основной системы отбросим правую жёсткую опору и заменим её действие тремя неизвестными силами X1, X 2 , X3 (рис. 25).

Рис. 24

Рис. 25

Разобьем раму на три участка I, II, III. На левой стойке наметим два сечения: 1, 2; на ригеле – три сечения: 2, 3, 4; на правой стойке – два сечения: 4, 5 (рис. 25).

92

Строим единичные эпюры моментов в основной системе (рис. 26):

Рис. 26

Строим эпюры моментов от единичной грузовой нагрузки (рис. 27):

18

+

 

4,5

q 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

Mq

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 27

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Формируем расчётные матрицы:

 

 

 

 

 

6

0

1

 

 

 

18

3

 

 

6

3

 

 

 

 

 

18

0

 

 

 

1

 

 

 

 

 

 

 

6

3

1

 

 

 

 

18

0

 

0

 

3

3

 

0

 

 

 

4,5

0

 

Lm

1 ;

Lmp

 

;

 

 

0

3

 

 

 

 

 

0

0

 

 

 

1

 

 

 

 

 

 

 

0

3

 

 

 

 

 

0

0

 

 

 

1

 

 

 

 

 

 

 

0

0

1

 

 

 

 

0

0

 

Матрица податливости:

 

 

 

 

 

 

 

 

 

 

 

 

b

1

 

0

0

 

 

 

 

 

 

 

 

 

b2

0

 

 

 

 

 

 

 

B

0

 

 

 

 

 

 

 

 

0

 

0

 

 

 

 

 

 

 

 

 

 

b3 .

F 1

M F

 

q

 

2

 

 

P

 

 

 

.

 

 

 

 

3

 

 

F

 

 

93

Здесь

b

 

 

 

h

2

1

 

3 4

2

1

 

 

1

 

4

 

2

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6EI1

 

 

 

 

6EI

 

1

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

1

2

 

 

2

 

EI

 

4

 

 

 

 

 

 

 

 

 

 

 

 

1

0

0

 

 

 

1

0

 

0

 

 

1

0

0

b

 

 

l

 

 

0

4

0

 

 

 

6

 

0

4

 

0

 

 

1

 

0

4

0

;

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6EI2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

EI

 

 

 

 

 

 

 

 

 

 

0

0

1

 

 

6EI

0

0

 

1

 

 

 

0

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

h

2

1

 

3 2

2

1

 

1

2

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6EI3

 

1

2

 

 

6EI

 

1

2

 

 

 

 

1

 

 

 

 

 

 

 

 

 

EI

То есть матрица податливости получает вид:

 

 

 

4

2

0

0

0

0

 

 

 

2

4

0

0

0

0

 

 

 

 

1

 

0

0

1

0

0

0

B

 

0

0

0

4

0

0

 

 

 

 

EI

0

0

0

0

1

0

 

 

 

 

 

 

0

0

0

0

0

2

 

 

 

 

 

 

0

0

0

0

0

1

Перемножаем матрицы:

A L0m T B L0m

1.

0

0

0

0

0

1

2 .

 

 

 

 

 

 

 

 

 

4

2

0

0

0

0

0

6

0

1

 

 

 

 

 

 

 

 

 

 

 

2

4

0

0

0

0

0

 

 

6

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

6

6

6

3

0 0 0

 

0

0

1

0

0

0

0

 

 

6

3

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

3

3

3 3 3 0

 

 

0

0

0

4

0

0

0

 

 

3

3

1

 

 

 

EI

1

1

1

1 1 1 1

 

 

0

0

0

0

1

0

0

 

 

0

3

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

0

0

0

0

0

2

1

 

 

0

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

0

0

0

0

0

1

2

 

 

0

0

1

 

 

 

 

504

162

90

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

162

108

45

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

EI

90

 

45

24

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

94

 

 

 

 

 

 

 

 

1 p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L0m T B L0mp

 

 

 

 

 

 

 

 

 

 

 

p

 

2 p

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

2

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

4

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

6

 

 

6

6

 

3 0 0

0

 

0 0 1 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

3

 

 

3

 

 

3 3 3 0

 

 

0 0 0 4

 

EI

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

1

 

 

1 1 1 1

 

 

0 0 0 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

0

0

 

 

 

18

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

18

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2592

 

 

 

 

 

 

 

 

 

 

18

0

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

810 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4,5

0

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

EI

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

 

 

 

 

 

 

 

 

 

450

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Далее находим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

EI

 

567

162

2430

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

162

3996

8100

.

 

 

 

 

 

 

 

 

 

 

 

 

 

93312

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2430

8100

28188

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теперь находим вектор неизвестных:

0

0

0

 

 

0

0

0

 

 

 

 

0

0

0

 

 

0

0

0

 

 

 

1

0

0

 

 

 

 

0

2

1

 

 

 

 

0

1

2

 

 

 

X1

 

 

 

 

 

 

 

 

 

1

 

 

p

 

 

 

X X2

 

A

 

 

 

X3

 

 

 

 

 

 

 

 

 

 

 

 

 

EI

 

567

162

2430

 

 

1

 

2592

 

 

5,44

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

162

3996

8100

 

 

 

 

810

 

 

0,13

.

93312

 

 

 

2430

8100

28188

 

 

EI

450

 

 

1,88

 

 

 

 

 

 

 

 

 

 

 

95

Подсчитываем значения изгибающих моментов в пронумерованных сечениях заданной стержневой системы:

 

M

1

 

 

 

 

 

 

 

 

 

18

3

 

 

 

 

 

6 0

 

 

 

 

 

 

 

 

 

 

 

 

18

0

 

 

 

 

 

 

 

6 3

 

M 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

2

 

 

 

 

 

 

 

 

 

 

18

0

 

 

 

 

 

 

 

6 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

L0

P L0

 

 

4,5

0

 

3 3

M M

3

 

X

 

 

 

 

 

 

 

 

mp

 

 

m

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

 

 

 

 

0 3

 

M 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

4

 

 

 

 

 

 

 

 

 

 

0

0

 

 

 

 

 

 

 

0 3

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

 

 

 

 

 

 

 

0 0

 

M5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выполняем проверку:

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

Lm

B M

 

 

 

0

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

EI

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Строим окончательную эпюру моментов (рис. 28):

1

 

 

 

 

3,75

 

 

 

 

 

4,88

 

1

 

 

 

 

 

1

 

5,44

 

 

4,88

 

 

 

 

 

 

 

 

1

 

0,13

 

5,81

 

 

1,88

 

 

1,5

 

1

 

 

 

 

 

 

 

 

 

1,5

 

1

 

 

 

 

 

1

 

 

 

 

1,88

 

Рис. 28

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

§13. Расчёт стержневых систем методом перемещений

в матричной форме

Для "n" раз кинематически неопределимой системы канонические уравнения метода перемещений в матричной форме имеют вид:

Rr

Z

 

R

p 0 .

(1)

96

Здесь Rr матрица реакций в присоединённых связях в основной системе от единичных перемещений этих связей,

 

 

 

 

r11

r12

r1n

 

 

 

 

 

 

R

 

r

r

r

 

;

 

 

 

 

21

21

2n

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

rn2

 

 

 

 

 

 

 

 

 

rn1

rnn

 

 

Z1, Z2 , , Z4 T ;

 

 

вектор неизвестных перемещений,

 

 

 

 

Z

Z

 

p вектор реактивных

усилий в присоединённых связях от внешней

R

 

нагрузки,

 

p R1 p

R2 p ... Rnp T ;

 

R

rij – реакция в i-й связи от единичного перемещения j-й связи.

Для вычислений компонентов матрицы реакций Rr и вектора реактивных усилий R p воспользуемся способом перемещения эпюр:

 

 

rik MEIi M k dz;

Rip MEIi M P1 dz.

(2)

В силу того что элементы rik , Rip вычисляются так же,

как и ik , ip

метода сил, то и матрицы Rr и RP

можно вычислять аналогично матрицам

A ,

 

p метода сил:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Rr Lom T B L0m ;

 

 

p Lom T B L0mp

 

.

(3)

 

 

R

P

Здесь L0m матрица влияния изгибающих моментов в основной системе

метода перемещений от единичных перемещений присоединённых связей Z1 1, Z2 1, ..., Zn 1. Матрица содержит "n"

столбцов и "m" строк;

n число единичных перемещений;

m число сечений, в которых вычисляются внутренние усилия; B матрица податливости отдельных, несвязанных элементов;

L0mp матрица влияния изгибающих моментов в любой статистически определимой системе, полученной из заданной, от внешних сил,

 

 

равных единице, F1 1, F2 1, ..., Fk

1;

 

 

 

вектор внешних сил,

 

F1 F2

...

Fk .

 

 

P

P

 

Решив матричное уравнение (1), найдём вектор неизвестных:

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

Z

Rr1

R

p L0m T B L0m

 

 

L0m T B L0mp

P

.

(4)

Окончательные значения изгибающих моментов в пронумерованных сечениях заданной стержневой системы определяют по уравнению

 

Lmp

 

L0m

 

.

(5)

M

P

Z

97

Здесь Lmp матрица влияния изгибающих моментов в основной системе

метода перемещений от внешних сил, равных единице.

При действии на статически неопределимую систему "t" вариантов внешних нагрузок векторы Z и R p переходят в прямоугольные матрицы

порядка n t , где n степень кинематической неопределимости системы. Матрица реакций Rr вычисляется по формуле (3).

Таким образом, для расчёта системы методом перемещений в матричной форме необходимо построить матрицы L0m , L0mp , B, Lmp и вектор P . При

этом алгоритм расчёта будет следующим:

а) выбор основной системы метода перемещений; б) построение единичных эпюр моментов от единичных воздействий

Z1 Z2 ... Zn 1;

в) построение единичных грузовых эпюр моментов от грузовой нагрузки, принимаемой равной единице, в статически определимой системе; г) построение единичных грузовых эпюр моментов от грузовой нагрузки, принимаемой равной единице, в основной системе метода пере-

мещений;

д) составление матриц L0m , L0mp , B, Lmp , P и их перемножение.

Проверкой правильности полученного решения является статическая проверка, заключающаяся в проверке выполнения условий равновесия узлов стержневой системы.

П р и м е р .

Построить эпюры внутренних усилий в статически неопределимой раме (рис. 29) методом перемещений в матричной форме [28].

Исходные данные: F 6 кН; q 43 кНм . Длины и жёсткости стержней

показаны на рис. 29. Для упрощения будем полагать, что EI 1. Стержневая системы будет три раза кинематически неопределима n 3 .

Р е ш е н и е .

Разбиваем заданную систему на семь участков, руководствуясь тем, чтобы в пределах каждого участка эпюры изгибающих моментов не имели разрывов и точек перегиба (на основной системе эти участки отмечены арабскими цифрами в кружочках). На каждом участке намечаем расчётные сечения: на участке 1 – сечения 1 и 2; на участке 2 – сечения 3 и 4; на участке 3 – сечения 5 и 6; на участке 4 – сечения 7 и 9; на участке 5 – сечения 8 и 9; на участке 6 – сечения 9, 10 и 11; на участке 7 – сечения 12 и 13.

98

Принимаем рабочее правило для изгибающих моментов: направление положительных моментов отмечено знаком (рис. 30).

Рис. 29

Основная система метода перемещений показана на рис. 30.

3

 

5

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

4

 

Z

 

 

 

 

Z2

 

 

 

 

3

1

4

9

6 10

11

 

Z3

 

2

 

 

 

 

7

 

 

 

12

 

 

 

1

 

 

 

5

 

7

 

 

 

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

13

 

 

Рис. 30

99

Строим единичные эпюры (рис. 31, а,б,в,г,д,е,ж), помня, что EI=1.

100

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