Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Методы расчета элементов конструкций подвижного состава. Учебно-методическое пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
970 Кб
Скачать

Рисунок 4.7 – Решение МКЭ для стержня с различными по длине характеристиками

4.3 Задачи для самостоятельного решения

1.Для 6, 8 и 14 участков произвести расчет МКЭ и построить графики изменения перемещений для трех вариантов граничных условий.

2.Для схемы, приведенной на рисунке 2.12, произвести расчет

МКЭ и построить графики изменения перемещений и по длине

стержня. Привести решения для 6, 8 и 14 участков разбиения.

41

5 Метод Бубнова-Галеркина

В отличие от метода конечных разностей, где использовалась линейная аппроксимация для участков разбиения, в методе Бубнова -

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

u x ui fi x , i 1

(5.1)

где ui

- неизвестные параметры, подлежащие определению;

 

f

x - координатная функция, выбираемая таким образом,

 

i

 

 

 

 

 

чтобы она удовлетворяла граничным условиям.

 

 

При

внесении

аппроксимирующей

функции

в

дифференциальное уравнение, т.к. она не является точным решением,

появляется функция невязки x , которая по физической сути представляет собой неуравновешенную нагрузку.

d

2

u x

 

n

 

 

 

 

x

n

x

 

 

u

f

 

 

 

 

2

 

i

 

 

 

dx

 

EF

i

 

 

 

EF

 

 

 

 

i 1

 

 

 

 

 

(5.2)

f

j

x

 

 

Умножая почленно уравнение (5.2) на возможное перемещение и интегрируя по области, приходим к соотношению:

l

 

 

 

 

l

n

 

l

ui

fi

x f j

x dx

f j

x dx x f j x dx

 

EF

0

i 1

 

 

 

0

 

0

 

 

 

 

 

(5.3)

Первое слагаемое представляет собой работу внутренних сил на возможных перемещениях, второе – работу внешних сил, правая часть

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

42

l

 

x f j

x dx 0.

0

 

(5.4)

Определенные из (5.3) с учетом (5.4) после взятия интегралов параметры ui вносятся в выражение для аппроксимирующей функции

(5.1).

Покажем ход решения методом Бубнова-Галеркина разных вариантов граничных условий.

5.1 Геометрические граничные условия

В этом случае (рисунок 2.1) аппроксимирующую

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

 

i x

u x ui

sin .

i 1

l

функцию

(5.5)

Она действительно удовлетворяет граничным условиям, так как

 

sin0 sin i 0 .

 

 

Внесем соотношение (5.5) в уравнение упругого равновесия

(1.5)

 

 

 

 

 

 

 

 

 

n

i

2

i x

 

n

x

 

 

 

 

ui

 

sin

 

 

 

(5.6)

 

l

EF

i 1

l

 

 

 

 

 

Умножим почленно на sin jl x и проинтегрируем от 0 до l:

 

l

n

 

i

2

 

i x

 

j x

 

 

l

n

 

j x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

i

 

 

sin

 

sin

 

 

dx

 

 

sin

 

dx 0.

(5.7)

 

0

i 1

 

l

 

 

 

l

 

 

l

 

 

0

EF

 

l

 

 

Здесь мы положили, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

j x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x sin

l

 

dx 0.

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

Первый интеграл, входящий в (6.7) в силу ортогональности:

43

l

i x

 

j x

0,

если i

j;

 

 

 

 

sin

l

sin

l

dx l

, если i

j.

0

 

 

 

 

 

 

2

 

 

Второй интеграл с учетом, что j=i.

l

n

 

i x

 

n

l

cos i 1 .

 

sin

dx

EF

l

EF i

0

 

 

 

Выражение (5.7) с учетом (5.8), (5.9) приобретает вид:

i

 

2

l

 

 

 

n

 

l

 

cos i 1 0,

 

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

i

l

 

 

2

 

 

 

EF

i

 

 

 

 

 

 

 

 

 

 

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2nl

2

3 cos i 1 .

 

 

ui

 

 

 

 

 

 

EF i

 

 

 

 

 

 

 

 

 

 

 

 

Внесем (5.11) в (5.5):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u x

2nl

2

 

n

1

cos i 1 sin

i x

 

 

 

 

.

 

 

 

 

3

 

3

 

 

 

EF

 

 

i

 

 

l

 

 

 

 

i 1

 

 

 

 

(5.8)

(5.9)

(5.10)

(5.11)

Иногда решение методом Галеркина называют точным,

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

u x

4nl

2

n

1

 

i x

 

 

 

 

sin

.

EF

3

i

3

l

 

 

 

 

 

 

i 1

 

 

 

Решение в Excel представлено на рисунке 5.1.

(5.12)

5.2 Смешанные граничные условия

Здесь (рисунок 2.6) аппроксимирующая функция с учетом

граничных условий имеет вид:

n

 

i x

 

 

u x ui

sin

Cx.

(5.13)

2l

i 1

 

 

 

44

Рисунок 5.1 - Решение методом Бубнова-Галеркина для геометрических граничных условий

Здесь второе слагаемое поддерживает силовое граничное условие

EF

du

P

 

 

dx x l

45

Найдем исходя из него параметр С

 

 

 

 

EF

du

 

 

n

u

i

cos

i x

 

 

P

 

 

EF

 

 

 

 

C

 

dx

 

 

i

2l

 

2l

 

 

 

 

x l

i 1

 

 

x l

 

 

Принимая i только нечетными (при этом cos обнуляется) будем

иметь

C

Р

EF

 

Тогда (6.13) можно записать в виде

n

 

i x

 

P

 

u x ui

sin

2l

 

EF

x.

i 1

 

 

 

(5.14)

Найдем функцию невязки, внеся выражение (5.13) в

дифференциальное уравнение (1.5) и взяв вторую производную

 

 

 

 

 

n

i

 

2

 

 

i x

 

n

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.15)

 

 

 

 

 

ui

 

 

 

 

sin

 

 

 

 

 

 

 

 

 

 

 

 

 

2l

 

 

2l

EF

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

Умножим почленно на

sin

j x

 

и проинтегрируем от 0 до l

 

2l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l n

 

i

 

2

i x

 

 

 

j x

 

 

 

l

 

n

 

j x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

i

 

sin

 

 

sin

 

 

 

dx

 

 

 

sin

 

dx 0.

(5.16)

 

0 i 1

 

2l

 

 

2l

 

 

 

 

2l

 

 

 

 

0

EF

 

2l

 

 

Здесь мы положили, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

x sin

j x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

2l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Первый интеграл, входящий в (6.16) в силу ортогональности:

l

i x

 

j x

0,

если i

j;

 

 

 

 

sin

2l

sin

2l

dx l

, если i

j .

0

 

 

 

 

 

 

2

 

 

(5.17)

Второй интеграл возьмем с учетом, что j=i. При удержании нечетных гармоник ряда

46

l

i x

 

2l

 

i

 

 

2l

 

sin

dx

 

.

 

i

cos

 

1

i

0

2l

 

 

2

 

 

 

Выражение (5.16) с учетом (5.17), (5.18) приобретает вид:

i

 

2

l

 

n

2l

 

 

 

0,

u

 

 

 

 

 

i

 

 

2

 

EF i

 

2l

 

 

 

откуда

 

 

 

 

 

 

 

 

 

 

16nl

2

 

 

ui

 

 

3 .

 

 

 

 

 

 

 

 

 

EF i

 

 

(5.18)

(5.19)

(5.20)

Внесем (5.20) в (5.13):

u x

16nl

2

n

1

 

i x

 

P

 

 

 

sin

 

x.

EF

3

i

3

2l

EF

 

 

i 1

 

 

 

 

Решение в Excel представлено на рисунке 5.2.

5.3 Задачи для самостоятельного решения.

(5.21)

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

47

Рисунок 5.2 - Решение методом Бубнова-Галеркина для смешанных

граничных условий

48

6 Метод Ритца-Тимошенко

В методе Ритца-Тимошенко потенциальной энергии, который в имеет вид:

 

1

l

du x

Э

 

2

EF

dx

 

0

 

 

 

 

 

записывается функционал полной случае растяжения-сжатия стержня

 

2

l

 

 

 

 

 

 

dx n x u x dx.

(6.1)

 

 

0

 

 

 

 

Для постоянной по длине нагрузки после внесения аппроксимирующей функции в функционал (6.1) он приобретает вид:

 

 

 

 

 

 

 

1

l

 

n

df x

2

l

 

n

fi x dx.

 

 

 

 

 

 

Э

 

 

 

 

 

 

 

 

 

EF

ui

i

 

dx

n ui

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 0

i 1

dx

 

 

0

i 1

 

 

 

 

 

1

l

n

2 dfi x

2

n

dfi x

 

df

 

x

l

n

 

 

 

 

 

 

 

 

j

 

 

 

 

 

2

EA

 

 

ui

 

dx

 

2 ui

dx

u j

dx

 

dx n ui

 

 

 

 

 

 

 

 

j i

 

 

 

i 1

 

 

0

i 1

 

 

 

 

 

 

 

 

 

 

 

 

0

Неизвестные параметры ui находятся из условия функционала полной потенциальной энергии

Э

0.

u

 

i

 

(6.2)

f

x dx.

i

 

минимума

(6.3)

Внося полученные соотношения для ui в соответствующую аппроксимирующую функцию, получаем выражение для перемещения.

В случае одинаковых аппроксимирующих функций решение методом Ритца-Тимошенко полностью совпадает с решением методом Бубнова-

Галеркина.

6.1 Геометрические граничные условия

Для геометрических граничных условий (рисунок 2.1) при внесении аппроксимирующей функции

u ui sin i x

n

 

i 1

l

функционал приобретает вид:

49

 

 

 

 

EF

l

 

 

n

i

 

 

i x

2

 

 

 

l

n

 

 

 

i x

 

 

 

 

 

Э

 

2

 

 

 

ui

l

cos

 

l

 

dx n ui

sin

 

l

dx

 

 

 

 

 

 

 

0

 

i 1

 

 

 

 

 

 

 

 

0

i 1

 

 

 

 

 

 

 

 

EF

l

n

 

 

2 i

 

2

2 i x

 

 

 

 

n

 

 

 

i

j

 

 

i x

 

 

j x

 

 

 

 

i

 

 

 

 

 

 

 

 

 

i

 

j

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

u

 

 

l

 

cos

 

l

2

 

 

u

u

 

 

l

l

cos

l

 

cos

l

dx

 

0

i 1

 

 

 

 

 

 

 

 

j i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

n

 

 

 

 

 

i x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n ui

sin

l

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Первый интеграл равен

(6.4)

l

 

i x

 

1

l

 

 

2i x

1

l

 

2i x

cos

2

dx

 

cos

sin

 

 

 

1

 

dx

x

2i

 

 

0

 

l

 

2

0

 

 

l

 

2

 

l

 

В силу ортогональности второй интеграл

l

0

l . 2

(6.5)

l

i x

 

j x

 

 

 

 

 

 

 

cos

l

cos

l

dx 0

при i j,

0

 

 

 

 

 

 

 

 

А третий интеграл уже рассматривался в методе Галеркина. Для нечетных гармоник будем иметь

 

l

l

 

i x

 

l

cos i 1 2nui

l

 

nui

 

sin

dx nui

.

i

l

i

i

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда выражение (6.4) приобретает вид

 

EF

2

i

2

l

n

l

 

 

n

 

 

 

 

 

 

Э

2

ui

 

l

 

 

2

2n ui

i

 

i 1

 

 

 

i 1

Из условия минимума функционала

(6.6)

Бубнова-

(6.7)

Э

i

2

l

 

 

 

2nl

 

 

u

 

 

0

u

EF

 

 

 

 

i

i

 

l

 

 

2

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

следует выражение для параметра

ui

4nl 2

 

,

EF i 3

Тогда функция перемещений полностью совпадает с (5.12)

(6.8)

(6.9)

50

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