2.2. Метод Бубнова-Галеркина Дифференциальное уравнение упругого равновесия растяжения-сжатия стержня имеет вид:
(2.2.1)
Запишем аппроксимирующую функцию
,
(2.2.2)
где:
ui
– неизвестные параметры, подлежащие
определению;
-координатные функции, удовлетворяющие
граничным условиям.
Внесем аппроксимирующую функцию в дифференциальное уравнение
,
(2.2.3)
при этом в правой части появляется функция невязки e(x), представляющая по физическому смыслу неуравновешенную нагрузку.
Неизвестные
параметры будем определять из условия
равенства нулю работы невязки на
возможных перемещениях, для чего умножим
уравнение почленно на возможное
перемещение
и проинтегрируем по длине стержня.
.
(2.2.4)
В силу ортогональности тригонометрических функций
(2.2.5)
Запишем остальные интегралы (последние два интеграла брались по частям):
(2.2.6)
Внося значения интегралов в выражение для работы невязки на возможных перемещениях, приходим к равенству:
,(2.2.7)
откуда
.
(2.2.8)
Выражение для перемещения с учетом размерности в точном решении принимает вид:
(2.2.9)
Дифференцируя по x и умножая на жесткость растяжения-сжатия EA, приходим к выражению для продольного усилия:
.(2.2.10)
Результаты счета приведены в таблице 1 и на рис. 2. В таблице приведены результаты расчета при удержании одного, десяти и ста членов рядов (2.2.9) и (2.2.10), которые сходятся достаточно быстро для перемещений и медленнее для усилий.
Таблица 3
|
0 |
0,25 |
0,5 |
0,75 |
1,0 |
|
|
1 10 100 точн. р. |
0 0 0 (0) |
0,476 0,477 0,478 (0,478) |
0,674 0,652 0,651 (0,651) |
0,476 0,504 0,504 (0,504) |
0 0 0 (0) |
|
1 10 100 точн. р. |
4,234 4,836 5,018 (5.038) |
2,994 2,640 2,613 (2,613) |
0 0,144 0,141 (-0,141) |
-2,994 -2,581 -2,548 (-2,548) |
-4,234 -5,317 -5,548 (-5,574) |
Рис. 6. Изменение продольного перемещения и нормального усилия по длине стержня (метод Бубнова-Галеркина, совпадение уже при 10-ти членах ряда).
Рис.
7. Изменение нормального усилия по длине
стержня (метод Бубнова-Галеркина,
совпадение уже при 10-ти членах ряда).
2.3. Метод Ритца-Тимошенко
Функционал полной потенциальной энергии системы имеет вид:
.
(2.3.1)
Внесем аппроксимирующую функцию (2.2.2) в выражение (2.3.1):
(2.3.2)
Запишем условие минимума функционала
(2.3.3)
для (2.3.2)
,
(2.3.4)
откуда
.
(2.3.5)
Внося значения интегралов (2.1.5) в (2.3.5), приходим к выражению для i–го неизвестного параметра
, (2.3.6)
которое полностью совпадает с (2.2.8). Совпадение решений, полученных методом Бубнова-Галеркина и методом Ритца-Тимошенко, объясняется тем, что аппроксимирующая функция бралась одна и та же.