Материал: 442

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

2. Форма занятия (семинар, практическое занятие и др.) Лекция.

3.Контингент (факультет, курс, группа). Факультет АДМ, 3-й курс, группа 3МТб.

4.Тема занятия.

Формирование матрицы преобразования координат для перехода от местной системы координат к общей.

5. Цели занятия:

• обучающая:

- студенты должны знать основные инструменты и методы модификаций сложных объектов; - владеть навыками построения и модификации сложных объектов;

• развивающая:

- развитие памяти;

 

 

- развитие коммуникативных способностей учащихся в совмест-

ной деятельности: умение вести диалог, выслушивать оппонен-

та, аргументированно обосновывать своюИточку зрения;

• воспитательная:

 

 

- воспитание дисциплинированности;

 

 

Д

- воспитание умения прислушиваться к мнению своих товари-

щей;

А

 

- воспитание стремлен я для приобретения новых знаний.

6. План занят я:

б

 

 

 

1. Построен е матр цы перехода от локальной системы ко-

ординат к глобальнойив табличном виде.

2. Определение направляющих косинусов.

3. ФормированиеС матрицы перехода.

7. Структурные элементы занятия:

• Организационный этап.

Приветствие учащихся, фиксация отсутствующих.

• Целевая установка.

Постановка целей занятия, сообщение темы.

• Выдача нового материала.

Последовательное изложения учебного материала по намеченному плану. Материал лекции приведен ниже.

• Контроль знаний.

Метод устного контроля. Выявление уровня знаний из смежных дисциплин о типах конструкций стен и способах их усиления.

21

• Подведение итогов.

Оценка успешности достижения целей занятия.

Материал лекции «Общая схема расчета конструкций методом конечных элементов (МКЭ)».

Представим алгоритм МКЭ в форме метода перемещений основными блоками.

1.Задание исходной информации, которая должна включать:

-количество элементов и узлов;

-координаты узловых точек конструкции в общей системе координат;

-взаимное расположение конечных элементов и узлов в местной и общей системе координат;

-внешние узловые, поверхностные и объемные силы;

-условия опирания конструкции (граничныеИусловия);

-геометрические и жесткостные параметры каждого элемента конструкции. Д

2.Построение матрицы жесткости конечного элемента в местной системе координат. А

3.Формирование матрицы преобразования координат для перехода от местной системы координатбк общей.

4.Определение матрицы жесткости конструкции в общей системе координат. и

5.Приведение поверхностных и объемных внешних сил, дейст-

вующих на элемент, к стат чески эквивалентным узловым силам (формированиеСлокального вектора узловых сил).

6.Определение суммарных узловых сил в каждом узле всей конструкции (формирование глобального вектора узловых сил).

7.Учёт граничных условий.

8.Определение узловых перемещений конструкции в общей системе координат.

9.Определение внутренних узловых усилий и напряжений в отдельных элементах конструкции.

Решение приведенной ниже задачи выполняется по описному выше алгоритму.

Рассмотрим простую шарнирно опертую балку, загруженную равномерно распределенной нагрузкой: q = 1 кН/м, l = 10 м, ЕJ =

=const (рис. П.5.2).

22

Рис. П.5.2. Расчетная схема балки

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

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

 

 

 

 

 

 

 

 

 

И

 

 

 

 

 

 

Д

 

Рис. П.5.3. Расчетная схема балочного конечного элемента

2. Локальная матрица жесткости имеет размерность 4×4 (коли-

 

 

 

 

А

 

 

 

чество узлов элемента умножается на количество степеней свободы в

узле):

и

 

 

 

 

 

 

 

 

 

 

k11

k12

k13

k14

 

[k]= k21

k22

k23

k24

 

С

k

 

k

 

k

33

k

34

 

 

 

б31 32

 

 

 

 

 

k

41

k

42

k

43

k

44

 

 

 

 

 

 

 

Коэффициенты матрицы жесткости определяются по формуле

kij = 0l ЕJ ×Ni" (x)×N"j (х)dx ,

где Ni" (x), N"j (х) ‒ вторые производные функций формы в i-м и j

единичных состояниях (i, j =1,…,4). Функция прогибов задается выражением

w(x) = N1(x) 1 + N2(x) 2 + N3(x) 3 + N4(x) 4 ,

23

где Ni(x) –координатные функции; i – неизвестные узловые пере-

мещения (i = 1,…,4).

Рассмотрим вычисление одного из коэффициентов локальной матрицы жесткости:

l

l

6

 

2x

 

2

 

12 EJ

k11 = 0

N1

(x) ×EJ ×N1

(x)dx = EJ

 

 

(

 

1)

 

dx =

 

3 .

 

2

l

 

l

 

 

 

0

l

 

 

 

 

 

 

Полностью локальная матрица жесткости балочного конечного элемента имеет вид

 

 

 

6

3l

6

3l

 

 

[k ]=

2EJ

 

3l

2l

2

3l

l

2

 

 

 

 

 

 

 

 

l 3

6 3l

6

3l

 

 

 

 

3l

 

Д

 

 

 

 

 

l2

3l

2l2

 

 

 

 

А

 

 

 

 

 

Подставляя численное значение

l =10 мИв выражение , записан-

ное для матрицы жесткости балочного конечного элемента, получим

 

 

0.012

0.06

0.012

0.06

 

[k ]= EJ

 

0.06

0.04

0.06

0.2

 

 

 

0.06

0.012

 

 

0.012

0.06 .

 

 

б

0.06

0.04

 

 

 

0.06

0.2

 

и

 

 

 

 

3. Данный пункт пропускаем, так как расчетная схема конструк-

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

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

[K ]= [k ].

5.Вектор узловых сил имеет размерность 1×4 (совпадает с числом степеней свободы балочного элемента):

24

F1 {F}= FF2 ,3F4

где Fi – узловые силы, определяемые по формуле

 

Fi

= lq(x)Ni (x)dx,

i =1,...,4 .

 

 

 

0

 

 

 

 

 

 

 

 

Рассмотрим вычисление одного компонента вектора{F}:

 

 

l

 

 

 

 

 

И

 

F1 = ql (1+

2x33

3x22 )dx =

ql .

 

 

0

l

 

l

 

 

2

 

 

 

2

 

Д

2

 

Аналогично определяем остальные компоненты локального

вектора:

 

А

 

 

 

 

 

б

 

 

 

 

 

 

 

 

F2 = ql

; F3 = ql

; F4 = − ql .

 

и

 

 

2

 

 

12

 

 

12

 

 

 

 

Полностью локальный вектор узловых сил можно представить

так:

С

 

 

 

 

 

 

 

 

 

 

 

 

 

ql

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

ql

 

 

 

 

 

{F}

=

 

12

 

 

 

 

 

 

 

 

 

 

 

 

 

ql

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ql 2

 

 

 

 

 

 

 

12

 

 

 

 

 

 

 

 

 

 

Подставив численные значения q=1 кН/м, l=10 м в выражение для вектора {F}, получим

25

Источник: https://studfile.net/preview/16406499/