Тогда уравнение МВН в окончательном виде может быть записано
Используя метод Галеркина Wn = Nn и выбирая базисные функции в виде Nm =xmможно получить:
При М=2 из решения системы получаются следующие значения параметров am : a1=11.758 a2=3.458.
Cоответствующие значения
Таблица 2.3
Результаты решения задачи на основе точного решения (ТР) и метода Галеркина (МГ) при М=2
Функция | МГ | ТР |
| 6.743 | 6.754 |
| 15.216 | 15.232 |
| 18.67 | 20.0 |
2.6. Общая формулировка естественных краевых условий для задач теплопроводности
Рассмотрим уравнение теплопроводности
T=T(xi), Q=Q(xi) – заданная функция координат xi, при граничных условиях
Представим аппроксимацию решения этого уравнения в виде разложения по базисным функциям
где
Для определения констант amиспользуем общую формулировку уравнения МВН
для
Первое слагаемое преобразуем с помощью формулы Остроградского-Гаусса
В результате уравнение МВН примет вид
Подставляя в (2.29) представление
[K]{
2.7. Применение МВН в задачах теории упругости
Рассмотрим применение МВН для решения задач теории упругости. Полный комплект уравнений, необходимых для решения задач теории упругости в перемещениях рассмотрен в первой главе и представлен соотношениями (1.1)−(1.21).
В отличие от рассмотренных выше задач теплопроводности для скалярных функций
Рассмотренные в первой главе уравнения можно непосредственно использовать для решения задачи в рамках МВН, однако оказывается более удобным преобразовать их к виду, приводящему к слабой формулировке МВН.
Запишем уравнения равновесия (1.12) и граниченые условия (1.1 и 1.13) для задач теории упругости в форме (2.24)МВН:
для
Первое слагаемое в этом уравнении можно привести к виду:
Положим
Тогда
т.е. условия на
В результате уравнение МВН в слабой формулировке может быть записано
В матричной форме это уравнение можно представить в виде
В соответствии с обозначениями, принятыми в § 1.1 первой главы
уравнение (2.33) примет вид
Конкретный вид матриц
По аналогии с (2.27) положим
Тогда уравнение МВН в окончательном виде может быть записано
Таким образом, уравнение МВН для задач теории упругости может быть сведено к системе алгебраических уравнений относительно параметров аппроксимации искомых функций