Курсовая работа: Расчет стационарного теплового поля в двумерной пластине

, (10)

где глобальная матрица теплопроводности K и глобальный вектор нагрузки F определяются по формулам

, . (11)

Для решения задачи (10) применялся следующий алгоритм:

· Вычисление разложения матрицы ().

· Оценка числа обусловленности. Если число обусловленности больше ( определяется точностью вычислительной машины), то выдается предупреждение, так как малые отклонения в коэффициентах матрицы могут привести к большим отклонениям в решении.

· . .

Реализация описанного выше метода проводилась на языке программирования С++ и FORTRAN в среде интегрированной среде разработки Microsoft Visual C++ 6.0. Конечные результаты данной работы приведены на рис.4 - 7.

Рис. 4

Рис.5

Рис.6

Рис.7

Список литературы:

1. Амосов А.А, Дубинский Ю.А, Копченова Н.В. Вычислительные методы для инженеров: Учеб. пособие. – М.: Высш. шк., 1994. – 544 с.

2. Сегерлинд Л. Применение метода конечных элементов. – М.: Мир, 1979. – 392 с.

3. Станкевич И. В. Сеточные методы (лекции и семинары 2002 года).