Курсовая работа: Расчет стационарного теплового поля в двумерной пластине
, (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 года).