ПОИСК Статьи Чертежи Таблицы Соотношения МКЭ для треугольного конечного элемента в плоской задаче теплопроводности из "Метод Конечных элементов в расчетах деталей тепловых двигателей " Рассмотрим блок матрицы теплопроводности, который имеет размерность (1 X 1), т. е. выражается в скаляр, поскольку в каждом узле имеется одна степень свободы. [c.61] Выражение (3.29) легко программируется, но прежде отметим следуюш,ее. Матрица геплопроводности Р симметрична, поэтому будем придерживаться принятой ранее стратегии программирования, т. е. вычислять симметричную часть матрицы и хранить ее в виде одномерного массива (рис. 3.1). Откажемся здесь от использования подпрограммы сортировки элементов матрицы в одномерный массив, поскольку речь идет о сортировке всего лишь одного числа, а не блока матрицы, и это прош,е сделать сразу вслед за вычислением элемента матрицы соответствуюш,им оператором Фортрана. Введем в рассмотрение одномерный массив С, в котором хранятся коэффициенты теплопроводности материала. [c.61] Параметрам NPE, NDM, DET присваиваются значения соответственно 3, 2 и NPE — число узлов конечного элемента, двухмерный массив ВЕ был определен ранее в (2.70). [c.62] Очевидно, что эта часть программы должна быть окаймлена внешним циклом по конечным элементам. Здесь также используется суммирование с накоплением результата в массиве SE, необходимость которого будет видна позже, при рассмотрении конечных элементов высокого порядка. [c.62] Выражения (3.36) и (3.37) мо-I ут быть обобщены в предположении линейной зависимости для q, h, Too на стороне элемента. [c.63] Вернуться к основной статье