Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
1360.pdf
Скачиваний:
35
Добавлен:
07.01.2021
Размер:
1.08 Mб
Скачать

4.2.2.Расчет балки, разбитой на два конечных элемента

1.Рассмотрим задачу, приведенную на рис. 4.1, разбив балку на два конечных элемента с двумя степенями свободы в узле и с заданным направлением глобальных осей (рис. 4.3).

Рис. 4.3. Конечные элементы балки

2. Локальные матрицы жесткости [k1] и [k2 ] определим, подста-

вив численные значения l1 =

 

l2 = 5 м в выражение (2.45):

 

 

 

 

 

0,096

0,24

0,096

0,24

 

 

 

 

 

 

 

 

 

 

 

И

 

 

 

 

[k1

 

 

 

0,24

0,8

 

0,24

0,4

 

] = [k 2 ] = EJ

 

 

 

.

 

 

 

 

0,096

0,24

0,096

0,24

 

 

 

 

 

 

 

 

 

Д

 

 

 

 

 

 

 

 

0,24

0,4

 

0,24

0,8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

 

 

 

 

 

 

3. Взаимосвязь между локальной и глобальной нумерацией ко-

нечных элементов, узлов бстепеней свободы представим в табличной

форме.

 

и

 

 

 

 

 

 

Таблица 4.1

 

 

 

 

 

 

 

 

 

 

 

Взаимосвязь локальной и глобальной нумераций

 

 

 

 

С

 

 

 

 

 

 

Номер

 

Номера узлов

 

 

Номера степеней свободы

 

элемента

 

локальные

 

 

 

глобальные

 

локальные

 

глобальные

 

 

I

 

1

 

 

 

1

 

 

1,2

 

 

1,2

 

 

 

2

 

 

 

2

 

 

3,4

 

 

3,4

 

 

 

 

 

 

 

 

 

 

 

 

 

II

 

1

 

 

 

2

 

 

1,2

 

 

3,4

 

 

 

2

 

 

 

3

 

 

3,4

 

 

5,6

 

 

 

 

 

 

 

 

 

 

 

 

4. Глобальная матрица жесткости формируется путем объединения общих степеней свободы соседних элементов. Глобальная матрица жесткости имеет размерность 6×6 и имеет вид

53

 

 

0,096

0,24

0,096

0,24

0

0

 

 

 

0,24

0,8

0,24

0,40

0

0

 

[K]= EJ

 

 

0,096

0,24

0,096

+0,096

0,024

0,024

0,24

0,24

.

 

 

0,24

0,4

0,24

0,24

0,8

+0,8

0,24

0,24

 

 

 

 

 

 

0

0

0,096

0,24

0,096

0,024

 

 

0

0

0,24

0,4

0,024

0,8

 

 

 

 

5. Для построения локальных векторов узловых сил используем выражение (4.2), в котором произведем замену l на l1 = l2 = 5 м:

 

2,5

 

 

2,083

 

{FI } = {FII } =

.

 

2,5

 

Д

 

2,083

 

 

 

А

 

 

6. Формирование глобального вектораИсвободных членов осу-

ществляется по принципу объединения компонентов локальных векторов соседних конечных элементов, соответствующих общим степеням свободы. Глобальный вектор свободных членов {P } имеет раз-

мерность 1×6 и состоит из локальных векторов первого и второго

элементов {FI }и {FII }:

бq l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С

 

 

 

 

q l12

 

 

 

 

 

2,5

 

 

 

q l

 

 

 

12q l

 

 

 

 

 

2,083

 

 

 

1

2

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5,0

 

{P } =

 

 

2

 

 

 

 

 

2

 

 

=

 

 

 

 

 

 

 

 

 

 

0,0

.

 

 

q l12

 

q l 22

 

 

 

 

 

 

12

+

12

 

 

 

2,5

 

 

 

 

 

 

 

 

 

 

 

q l2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,083

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

q l

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

54

 

 

 

 

 

 

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

 

 

 

0,096

 

0,24

0,096

 

0,24

 

0

 

0

 

 

 

 

0,24

 

0,8

 

0,24

 

0,40

 

0

 

0

 

[K ]= EJ

 

 

 

 

 

 

 

0,096

0,24

0,192

 

 

0

 

 

0,24

 

0,24

.

 

 

 

0,24

 

0,4

 

 

0

 

 

 

1,6

 

0,24

 

0,4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

0

0,096

0,24

 

0,096

 

0,024

 

 

 

0

 

 

0

 

0,24

 

0,4

 

0,024

 

0,8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

И

 

 

 

8. Система разрешающих уравнений (4.3) принимает вид

 

 

 

0,8

0,24

0,4

 

 

0

 

2

 

 

 

 

2,083

 

 

 

 

0,24

0,192

0

 

0,24

 

 

 

 

5

 

 

 

EJ

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

= 0.

 

 

0,4

0

 

1,6

 

0,4

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0,24

 

0,4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,04Д

 

2,083

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

В результате решения этойАсистемы получаем значения узловых

перемещений, составляющ х вектор {}:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

w1

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

и

 

θ

1

 

 

 

 

 

41,665 EJ

 

 

 

 

 

{} =

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

=

w2

 

=

 

 

130,20 EJ

.

 

 

 

 

 

 

 

 

 

θ

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

3

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

41,665 EJ

 

 

 

 

 

 

 

 

 

 

θ3

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

55

 

 

 

 

 

 

 

 

 

 

 

 

 

41,665

 

 

 

 

 

130,2

 

 

 

 

 

 

w

(x ) = N

(x ) 0

+ N

 

(x )

 

 

 

 

 

+ N

 

(x )

 

 

+ N

 

(x ) 0

;

 

 

 

EJ

 

 

EJ

 

1 1

 

1

 

1

 

 

 

2

1

 

 

 

3

 

1

 

 

 

 

 

4

1

 

w

(x

 

 

 

 

(x

 

)

130,2

+ N

 

 

(x

 

) 0 + N

 

(x

 

) 0 N

 

(x

 

)

41,665

 

) = N

 

 

 

EJ

 

 

 

 

 

 

 

 

EJ

,

2

 

2

 

1

 

2

 

 

 

 

2

 

2

 

 

3

 

2

 

 

 

4

 

2

 

 

где w1(х1) – прогиб на первом участке (0 < x1 < l1); w2(х2) – прогиб на втором участке (0 < x2 < l2).

Рис. 4.4. Эпюра прогибов w EI

 

 

Д

 

9. Из курса сопротивления материалов известна зависимость

А

MИ

 

d 2 w

 

(4.5)

d x

2

= −

EJ

.

б

 

 

 

 

 

Для построения эпюры моментов воспользуемся формулой (3.23):

M = − EJ w′′(x) = − EJ [N1′′(x) 1 + N2′′(x) 2 + N3′′(x) 3 + N4′′(x) 4 ],

где Ni(x) – вторые про зводные функций формы, определяемые по

формулам (2.22).

и

Учитывая найденные значения узловых перемещений, запишем

выражения для изгибающих моментов в пределах первого и второго

 

 

 

 

С

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

конечных элементов в следующем виде:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

 

 

 

′′

 

 

+ N

′′

(x )

41,665

+ N

 

′′

(x )

130,2

+ N

′′

(x )0

 

1

(x ) = −ЕJ N (x )0

2

 

 

EJ

 

3

EJ

4

;

 

1

 

1 1

 

 

1

 

 

 

 

 

1

 

 

 

 

1

 

M

 

(x

 

 

′′

 

)

130,2

+ N

′′

(x

 

)0

+ N

′′

(x

 

)0 N

′′

(x

 

 

)

41,665

2

2

) = −ЕJ N (x

2

EJ

 

2

2

3

2

4

2

EJ

.

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как вторые производные координатных функций являются функциями линейными, то и выражения для изгибающих моментов в

56

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

Рис. 4.5. Эпюра изгибающих моментов М, кН м

Изображенная на рис. 4.5 эпюра Мявляется узловой, так как

она является результатом узловых воздействий, то есть перемещений

изгибающих моментов. Эпюра Мр от распределеннойИ нагрузки строится в основной системе метода перемещений (рис. 4.6).

1,…,6. Поскольку балочный элемент испытывает кроме узловых

воздействий еще и воздействие равномерно распределенной нагруз-

ки, то это должно быть учтено при построении окончательной эпюры

Для получения суммарной эпюры моментов ординаты эпюры

Мот узловой нагрузки суммируются с ординатами эпюры Мр от ме-

стной нагрузки:

 

 

 

Д

 

 

M = M + M р.

(4.5)

 

 

 

 

 

А

 

 

 

б

 

 

 

и

 

 

 

С

 

 

 

 

Рис. 4.6. Эпюра Мр от распределенной нагрузки

Окончательная эпюра М представлена на рис. 4.7.

57

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