Применение МКЭ для расчета стержневых систем

Расчет стержневых систем с помощью МКЭ приводит к тем же уравнениям, что и обычный метод перемещений, хотя подходы к их построению несколько отличаются по форме.

При этом КЭ рамы с 6 степенями свободы имеет на концах два узла, в каждом из которых введено по три связи – две линейных и одной моментной.

Матрица жесткости такого КЭ имеет 6 порядок и ее элементами являются реакции rij в шести введенных связях от единичных смещений этих связей.

Компонентами вектора приведеннойузловой нагрузки являются взятые со знаком минус реакции во введенных связях от приложенной к КЭ нагрузки, которые, в отличие от МП будем обозначать не Rip0 , а rip0:

{Pэ}= – [r1p0, r2p0, … , r6p0]Т,

где индексом «т» обозначена операция транспонирования.

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

Анализ КЭ заключается в определении реакций во введенных связях {Sэ}= [S1, S2, S3, S4]Т от кинематических воздействий {Zэ}= [Z1, Z2, Z3, Z4]Т и от действующей местной нагрузки. Первая зависимость имеет вид:

{Sэ} = [Rэ]{Zэ}.

Для построения матрицы жесткости [Rэ] рассмотрим КЭ при единичных кинематических воздействиях (рис. 7.1, б), объединив соответствующие им функции формы в матрицу-строку:

[N] = [N1(x), N2(x), N3(x), N4(x)].

Этим функциям формы соответствуют уже известные эпюры моментов, приведенные на рис. 6.4 и 7.1, в, которые также объединим в вектор:

{M} = [M1(x), M2(x), M3(x), M4(x)]Т.

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

M(x) = – EJ v''(x),

можно записать:

{M} = – EJ [N'']Т. (7.2)

Воспользовавшись соотношением (6.4):

rij = rji = Sò(`Mi0 ·`Mj0/EJ )ds ,

получим с учетом (7.2):

[Rэ] = (1/ EJ) ∫{M}{M}Тdx = EJ ∫ [N'']Т[N''] dx. (7.3)

Для построения вектора приведеннойузловой нагрузки учтем, что на основании принципа суперпозиции уравнение изогнутой оси КЭ можно представить в виде:

v(x) = Σ Ni (x) · Zi = [N]{Zэ}. (7.4)

Применение МКЭ для расчета стержневых систем - student2.ru

Рис. 7.1

Поэтому, дополнив соотношение (6.5) работой распределенной нагрузки, приложенной к КЭ, и сменив обозначения, получим:

rip0 = – [S Pk · Nik + ∫ q(x) · Ni (x) dx].

Таким образом, искомый вектор приведеннойузловой нагрузки равен:

{Pэ}= S Pk ·[Nk]Т + ∫ q(x) [N]Т dx. (7.5)

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

Представим уравнение изогнутой оси КЭ в виде полинома:

v (x) = [H] {a}, (7.6)

где [H] = [1, x, x2, x3], а {a} = [a1, a1, a1, a1]Т.

Приравнивая (7.4) и (7.6) в узловых точках КЭ, то есть x1 = 0 и x2 = l, получим:

Z1 = v (0) = 1 + a1x1+ a2x12+ a3x13 ;

Z2 = v'(0) = a1+2a2x1+ 3a3x12 ;

Z3 = v (l) = 1 + a1x2 + a2x22 + a3x23 ;

Z4 = v'(l) = a1 +2a2x2+ 3a3x22,

или иначе

{Zэ} = [L] {a},

где [L] = [[H(x1)]T, [H'(x1)]T, [H(x2)]T, [H'(x2)]T]Т.

Обратная зависимость

{a} = [L]–1{Zэ}

после подстановки в (7.6) приводит с учетом (7.4) к искомой формуле:

[N] = [H] [L]–1. (7.7)

В скалярной форме последняя зависимость имеет вид:

N1(x) = 1 – 3ξ2 + 2 ξ3;

N2(x) = l (ξ – 2 ξ2 + ξ3);

N3(x) = 1 – 3η 2 + 2η 3;

N4(x) = l (– η + 2η 2 – η 3),

где ξ = x/l , η = (l – x)/l .

Подставляя (7.7) в (7.3), получим искомую матрицу жесткости КЭ балки:

[Rэ] = EJ ∫ [N¢¢]T×[N¢¢]dx = (EJ/l) Применение МКЭ для расчета стержневых систем - student2.ru .

Вектор приведенной узловой нагрузки находим по формуле (7.5). Для постоянной равномерно распределенной нагрузки q(x) = q он имеет вид:

{Pэ} = [P1, P2, P3, P4]Т = (ql/12)[ 6, l, 6, –l ]Т.

Примечание.

Двумерным аналогом балочного КЭ является прямоугольный КЭ с 12 степенями свободы для расчета плит, изогнутая поверхность которого аппроксимируют полиномом

w(x,y) = [H]{a},

где [H]= [1, x, y, x2, xy, y2, x3, x2y, xy2, y3, x3y, xy3].

В каждой из четырех узловых точек такого КЭ, расположенных в его вершинах, вводят по три связи: линейную, препятствующую вертикальным перемещениям в направлении оси Oz, и две моментные в направлениях осей Ox и Oy.

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

ЛИТЕРАТУРА

1. Куликов И.С.Расчет конструкций на деформируемом основании: Учебное пособие / И.С. Куликов. – Горький: Изд-во ГИСИ им. В.П.Чкалова, 1986. – 72 с.

2. Масленников А.М., Егоян А.Г. Основы строительной механики для архитекторов: Учебное пособие / А.М.Масленников, А.Г.Егоян. – Ленинград: Изд-во ЛГУ, 1988. – 264 с.

3. Оксанович Л.В. Невидимый конфликт / Л.В.Оксанович. – М.: Стройиздат, 1981. – 191с.

4. Розин Л.А., Константинов И.А., Смелов В.А.Расчет статически определимых стержневых систем: Учебное пособие / Л.А.Розин, И.А.Константинов, В.А.Смелов. - Ленинград: Изд-во ЛГУ, 1983. - 228 с.

5. Сборник заданий для курсовых работ по теоретической механике: Учебное пособие / Под ред. А.А.Яблонского. – М.: Высш.шк., 1985. – 367 с.

6. Смирнов В.А., Иванов С.А., Тихонов М.А.Строительная механика: Учебник для вузов / В.А.Смирнов, С.А.Иванов, М.А.Тихонов. - М.:Стройиздат, 1984. - 208 с.

7. Rakowski G.Komputerowa mechanika konstrukcji / G. Rakowski. –Warszawa: Wydawnictwa Politechniki Warszawskiej, 1977. – 279 с.

ОГЛАВЛЕНИЕ

Стр.

Предисловие........................................................................................................ 3

Глава 1. Введение.............................................................................................. 4

1.1. Предмет строительной механики и ее задачи.................................... 4

1.2. Кинематический анализ сооружений................................................. 5

1.3. Основные уравнения строительной механики................................. 18

Глава 2. Расчет статически определимых стержневых систем................. 21

2.1. Свойства статически определимых систем...................................... 21

2.2. Внутренние усилия в рамах............................................................. 21

2.3. Расчет плоских ферм........................................................................ 31

2.4. Расчет трехшарнирных арок........................................................... 36

Глава 3. Определение перемещений в СОС................................................. 40

3.1. Работа сил, приложенных к твердому телу..................................... 40

3.2. Работа сил, приложенных к деформируемому телу....................... 42

3.3. Общие теоремы строительной механики......................................... 44

3.4. Работа внутренних сил плоской стержневой системы.................... 48

3.5. Интеграл Мора-Максвелла.............................................................. 50

3.6. Формула Верещагина....................................................................... 51

3.7. Примеры определения перемещений............................................... 54

Глава 4. Расчет статически неопределимых балок и рам методом сил... 58

4.1. Свойства статически неопределимых систем................................... 58

4.2. Суть метода сил. Канонические уравнения МС.............................. 60

4.3. Определение внутренних усилий в МС........................................... 65

4.4. Проверка правильности решения в МС.......................................... 68

4.5. О выборе ОС МС. Признаки ортогональности эпюр..................... 70

4.6. Расчет симметричных систем в МС................................................. 73

4.7. Расчет неразрезных балок................................................................ 76

Глава 5. Расчет статически неопределимых арок и ферм методом сил... 78

5.1. Расчет статически неопределимых ферм......................................... 78

5.2. Расчет статически неопределимых арок.......................................... 79

Глава 6. Расчет СНС методом перемещений............................................... 81

6.1. Суть метода перемещений. Основная система МП......................... 81

6.2. Канонические уравнения метода перемещений............................... 82

6.3. Вычисление коэффициентов канонических уравнений................... 83

6.4. Общий метод вычисление коэффициентов...................................... 86

Глава 7. Понятие о расчете СНС методом конечных элементов.............. 89

7.1. Суть метода конечных элементов.................................................... 89

7.2. Применение МКЭ для расчета стержневых систем......................... 90

Литература........................................................................................................ 93

Наши рекомендации