10380

- 90 -

Решение.

1.Отметим, что заданная статически неопределимая система имеет две лишние связи, и при ее расчете методом сил число неизвестных равнялось двум. При решении той же задачи методом перемещений число неизвестных, равное в данном случае числу незакрепленных жестких узлов, будет равно только единице, поэтому в этом примере МП будет эффективнее метода сил.

2.Основную систему МП получаем, вводя моментную связь в этом свободном узле (рис. 6.5, б).

3.Каноническое уравнение метода перемещений имеет вид:

4. С помощью стандартных готовых решений (рис. 6.4) строим эпюры изгибающих моментов от единичного значения Z1 и от заданной нагрузки (рис.

6.5, в, г):

2 3

1

Рис. 6.5

5. Вычисляем r11 и R1p0 , рассматривая равновесие вырезанного второго узла рамы:

r11 = 7i, R1p0 = ql2/12.

Полагая для удобства EJ = 2, получим i = EJ/l = 1, откуда r11 = 7, R1p0 = 1/3. 6. Решая уравнение (а) найдем:

Z1 = R1p0/ r11 = – 1/ 21.

-91 -

7.Искомую эпюру изгибающих моментов (рис. 6.5, е) можно построить по формуле (6.2):

Mp = Mp0 + Mi0Zi.

Нетрудно заметить, что она совпадает с эпюрой, полученной ранее в примере (4.3) с помощью метода сил (рис. 4.5, и).

Примечания

1. Метод перемещений в отличие от метода сил не требует проведения кинематической проверки – достаточно убедиться в равновесии узлов построенной эпюры Mp.

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

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

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

Пусть рама загружена произвольной нагрузкой (рис. 6.6, а), а соответствующая ей основная система МП образована введением двух связей: моментной

– i и линейной – j (рис. 6.6, б).

Рассмотрим два состояния этой системы, соответствующие единичным смещениям введенных связей, и обозначим через Mi0 и Mj0 соответствующие им эпюры изгибающих моментов (рис. 6.6, в, г).

Вычислим работу внешних сил первого состояния системы на перемещениях второго состояния:

A12 = rii· θij + rji·δjj = rii· 0 + rji·1 = rji .

Учитывая, что в силу (3.15)

A12 = – W12 = ( Mi0 · Mj0/EJ ) ds,

получим отсюда искомое выражение для определения удельных реакций:

Последнее выражение напоминает формулу (4.4) для вычисления коэффициентов канонических уравнений в методе сил:

- 92 -

ij = ( Mi0 Mj0) = ( Mi0 Mj0 /EJ)ds,

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

ip0= ( Mi0 Mp0) = ( Mi0 Mp0 /EJ)ds.

В действительности это не так: Rip0 ≠ ( Mi0 Mp0), а ( Mi0 Mp0) = 0.

Рассмотрим снова два состояния основной системы МП, и пусть первое по-прежнему соответствует единичному смещению i-й связи (рис. 6.6, в), а второе – единичной силе, приложенной в k-й точке этой системы (рис. 6.6, д).

Рис. 6.6

Работа внешних сил первого состояния системы на перемещениях второго состояния равна нулю:

A12 = rii · θip+ rji · δjp = rii · 0 + rji · 0 = 0,

а поскольку A12 = A21 , то

A21 = P · pi + rip· θii = 0,

откуда

rip = – pi .

- 93 -

То есть реакция в i-й связи основной системы от единичной силы, приложенной в точке k , равна взятому со знаком минус перемещению точки приложения силы от единичного смещения этой связи.

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

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

Таким образом, реакция в i-ой связи ОС МП от заданной нагрузки равна взятой со знаком минус работе всех сил, приложенных к системе, на перемещениях, вызванных единичным смещением этой связи.

Примечание

Поскольку A12 = 0, а A12 = – W12 , то действительно: ( M 0 M 0) = – W12 = 0. При этом

i p

нетрудно доказать, что искомую реакцию можно вычислить по формуле:

Rip0 = ( M 0i M p0) = – Σ( M 0i M p0/ EJ) ds,

где Mi0 – это по-прежнему эпюра моментов в ОС МП от единичного смещения i-й связи, а Mp 0 – эпюра моментов в любой ОС МС от заданной нагрузки. Такой, например, является

эпюра, приведенная на рис. 6.6, е.

Легко убедиться, что в этом случае значение R 0 = – 3Pl/16 совпадает с табличным

ip

значением, указанным на рис. 6.4.

- 94 -

ГЛАВА 7. ПОНЯТИЕ О РАСЧЕТЕ СНС МЕТОДОМ КОНЕЧНЫХ ЭЛЕМЕНТОВ

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

Метод конечных элементов (МКЭ) является обобщением рассмотренного выше метода перемещений на двух- и трехмерные системы. Он успешно применяется для расчета самых разнообразных строительных конструкций и сооружений, в том числе – тонкостенных систем и массивов.

В отличие от классических форм МС и МП МКЭ является компьютерным методом. Он появился в 60-е годы прошлого века и в настоящее время является самым распространенным методом расчета, реализованным в большинстве систем автоматизированного проектирования.

Особенность метода в том, что он позволяет решать задачи теории упругости методами строительной механики.

Реализация МКЭ включает следующие этапы:

1)разбиение конструкции на конечные элементы с учетом ее геометрии, физических свойств материала и нагрузки;

2)анализ отдельного конечного элемента (КЭ), в ходе которого строится

матрица жесткости КЭ – [Rэ] и определяется вектор приведенной узловой нагрузки – {Pэ};

3) анализ всей конструкции и формирование системы конечно-элементных уравнений:

где [R], {Z} и {P} – соответственно матрица жесткости, вектор неизвестных и вектор приведенной узловой нагрузки всей системы;

4)решение системы уравнений и определение по вектору {Z} напряжений

иусилий в отдельных конечных элементах.

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

Уже к 70-м годам прошлого века была создана хорошо апробированная система конечных элементов, приспособленных для решения самых разнообразных задач.

При решении плоской задачи и задачи изгиба пластин и плит применяют треугольные и прямоугольные КЭ с различным числом степеней свободы. Для решения трехмерных задач применяют КЭ в виде параллелепипеда или тетраэдра. Разработаны конечные элементы для решения осесимметричных задач и целое семейство плоских и криволинейных КЭ для расчета оболочек.

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

- 95 -

телом. Помимо строительной механики и теории упругости МКЭ находит применение в решении широкого круга задач математической физики.

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

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

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

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

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

которые, в отличие от МП, будем обозначать не 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),

- 96 -

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

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

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

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

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

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

(7.4)

z1 =1

z2 =1

z3 =1

z4 =1

Рис. 7.1

- 97 -

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

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

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

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

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

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

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

Z1 = v (0) = 1 + a1x1+ a 2x12+ a 3x13 ; Z2 = v'(0) = a1+2a2x1+ 3a3x12 ;

Z3 = v (l) = 1 + a1x2 + a 2x22 + a3x23 ; Z4 = v'(l) = a1 +2a2x2+ 3a x 2,

3 2

или иначе

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

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

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

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

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

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

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 .

- 98 -

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

Вектор приведенной узловой нагрузки находим по формуле (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 локальной системы координат.

- 99 -

ЛИТЕРАТУРА

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

2.Куликов И.С. Статика твердого тела: учеб. пособие / И.С. Куликов, Г.А. Маковкин. – Н. Новгород: Изд-во ННГАСУ, 2008. – 72 с.

3.Куликов И.С. Статика деформируемого тела: учеб. пособие / И.С. Куликов. – Н. Новгород: Изд-во ННГАСУ, 2008. – 64 с.

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

264 с.

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

дат, 1981. – 191с.

6.Розин Л.А. Расчет статически определимых стержневых систем: учеб.

пособие / Л.А.Розин, И.А.Константинов, В.А.Смелов. Ленинград: Изд-во ЛГУ,

1983. 228 с.

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

8.Сборник задач и упражнений по строительной механике : учеб. пособие. Ч.1 : Статически определимые системы /Лампси Б.Б. и др. Н.Новгород : ННГАСУ, 2015.

9.Сборник задач и упражнений по строительной механике : учеб. пособие. Ч.2 : Статически неопределимые системы /Лампси Б.Б. и др. Н.Новгород : ННГАСУ, 2015.

10.Смирнов В.А., Иванов С.А., Тихонов М.А. Строительная механика:

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

1984. 208 с.

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