8(495)909-90-01
8(964)644-46-00
pro@sio.su
Главная
Системы видеонаблюдения
Охранная сигнализация
Пожарная сигнализация
Система пожаротушения
Система контроля удаленного доступа
Оповещение и эвакуация
Контроль периметра
Система домофонии
Парковочные системы
Проектирование слаботочных сетей
Аварийный
контроль
Раздел: Документация

0 ... 10 11 12 13 14 15 16 ... 87

2.3.1. Матрица жесткости

Для моделирования упругих одномерных элементов конструкций, несущих изгибнуц нагрузку (балки), используют балочный элемент. Характеристиками этого типа конечнц

элементов являются длина элемента L, момец инерции площади поперечного сечения I и % дуль упругости Е. Линейный балочный элемент схематично показанный на рис. 2.12, ограничу двумя узлами i и j, каждый нз которых имеет, щ, минимум, две степени свободы: прогиб V н ущ поворота сечения относительно осн

9 = dv/dx. В узлах действуют перерезываю.]

щая сила F и изгибающий момент М тельно оси Z.

vi>Fi

*---

V

у—

j) X

. L

к.

Mj

Puc2J2

Согласно элементарной теории изгиба:

ОТНОСЦ.

My

(2.»)

(2.22)

Используя результаты теории балок, вычисляем коэффициенты матрицы жесткости) системе координат*, г для элемента с узлами!, j (рис. 2.12):

V,-

12

6L

-12

6L

EI

6L

41?

-6L

21?

Mt

L3

-12

-6L

12

-6L

VJ

FJ

6L

21?

-6L

41?

eJ.

Mj\

(2.23!

Реальные балки воспринимают не только сдвиговые нагрузки и изгиб, но и осевые * грузки. Матрицу жесткости для балочного элемента в этом общем случае получим ко» бинированием матрицы жесткости (2.23) с матрицей жесткости для стержневого элемент (2.17):

23. Балочный элемент


к =

ЕА L

О О

ЕА L

О -О

Обратим внимание на то, что элементы матрицы к имеют различные размерности. Это связано с тем, что н компоненты векторов перемещений, н компоненты сил также имеют различные размерности ([м], [рад] и [И], [Нм] соответственно).

Для вычисления матрицы жесткости в пространственной системе координат она записывается сначала в плоской локальной системе координат, а затем переводится в глобальную пространственную систему координат.

Рассмотрим несколько примеров.

0

0

ЕА L

0

0

12EI

6EI

0

12EI

6EI

L3

L2

L3

L2

6Е1

4EI

0

6EI

2Е1

1?

L

" L2

L

0

0

ЕА L

0

0

12EI

6EI

0

12EI

6EI

I3

" L2

I3

" L2

6EI

2EI

0

6EI

4EI

L2

L

~ I2

L

Y

Р

л

\

L

1-1

2 3 L

1 *

2.3.2. Примеры

Рис.2.13

Шимер 1. Брус (рнс. 2.13) защемлен с двух концов н нагружен посредине сосредоточенной силой Р н изгибающим моментом М. Определить смещение н угол поворота узла 2, а та°ке опорные реакции.

££Щение. Сформируем матрицы жесткости отдельных элементов:

V; в] V2 в2V2 в2 V3 в3

Е1

126L -126L

6L4L2 -6L21?

-12-6L 12-6L

6L

2L2 -6L 4L2

к2~7

6L -12 6L 4L2 -6L 2L2

12 6L

-12 -6L 12 -6L 6L 21? -6L 41?


Запишем глобальное конечно-элементное уравнение

El

0i

V,

о,

0,

" 12

61

-12

6L

0

0

V

6L

41}

-6L

2L1

0

0

0,

м

-12

-6L

24

0

-12

6L

6L

2L2

0

и2

-6L

2L2

М2

0

0

-12

-6L

12

-6L

0

0

6L

2L1

-6L

.

и зададим граничные условия: V/ = v3 = 0, 61=63= 0, F2y = —Р, М 2 = М. i

учетом этого глобальное конечно-элементное уравнение в матричной форме будет выпц. деть следующим образом:

EI

24 0

\V2U

L3

0 8L2

W 1

м

Решением этого уравнения является:

V

-PL2

L

24EI

в2.

ЗМ

Из глобального конечно-элементного уравнения получаем формулы для реакш опор:

F1Y

-12

6L

V

2P + 3M/L

М]

EI

-6L

2L2

1

PL + M

F3Y

-12

-6L

4

2P-3M/L

М3

6L

U2.

-PL + M

Напряжения на концах бруса можно вычислить по формуле:

My

— re — . s

a = aY = -

Задача решена. У

1 111 11111111

Пример 2. Дана консольная балка, нагружен* поперечной равномерно распределенной нагрузи" р (рнс. 2.14). Определить смещение и угол повор та правого конца балки, а также реакции опоры.

х Решение. Сначала рассмотрим вопрос приведен &>1*распределенной нагрузки к узлам. Можно показе

что нагрузка трансформируется к узлам, как гкЛ зано на рнс. 2.15. Проверить правильность э1 Рис.2.14схемы можно, сравнивая величину упругой рабе!

деформации для обеих схем. Применение этого к данной задаче приводит к расчетной схеме, показанной 1

1

рис. 2.16. Здесь: / = pL/2, fit = pL2 jl2.



0 ... 10 11 12 13 14 15 16 ... 87