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