|
|||||||
АвтоАвтоматизацияАрхитектураАстрономияАудитБиологияБухгалтерияВоенное делоГенетикаГеографияГеологияГосударствоДомДругоеЖурналистика и СМИИзобретательствоИностранные языкиИнформатикаИскусствоИсторияКомпьютерыКулинарияКультураЛексикологияЛитератураЛогикаМаркетингМатематикаМашиностроениеМедицинаМенеджментМеталлы и СваркаМеханикаМузыкаНаселениеОбразованиеОхрана безопасности жизниОхрана ТрудаПедагогикаПолитикаПравоПриборостроениеПрограммированиеПроизводствоПромышленностьПсихологияРадиоРегилияСвязьСоциологияСпортСтандартизацияСтроительствоТехнологииТорговляТуризмФизикаФизиологияФилософияФинансыХимияХозяйствоЦеннообразованиеЧерчениеЭкологияЭконометрикаЭкономикаЭлектроникаЮриспунденкция |
Алгоритм. 1. Считываем сетку, элементы, материал из файла - формируем массивы mesh, material, elements1. Считываем сетку, элементы, материал из файла - формируем массивы mesh, material, elements. 2. В цикле по элементам: для i = 1 до m 2.1 Определяем номера узлов текущего элемента: Для номера узлов текущего элемента. 2.2 Формируем элементную матрицу (формирование матрицы для задачи теплопроводности см. ниже). 2.3 Проводим процесс сборки: формируем глобальную матрицу жесткости. Для Для 2.4 Формируем элементный вектор правых частей . 2.5 Проводим процесс сборки – формируем глобальный вектор правых частей. Для 3. Вносим граничные условия Дирихле в матрицу и вектор правых частей: 3.1 3.2 Для 4. Решаем систему 5. Выводим результат 6. Определение расчетных величин в элементах.
Процедура формирования локальной матрицы жесткости для задачи теплопроводности (п.2.2 алгоритма, формулы (10), (11), (12), (14)). 1. Определяем координаты узлов текущего элемента Для
2. Формируем
;
3. Вычисляем площадь треугольника 4. Формируем матрицу по формуле (12). 5. Если на одной из сторон задан конвективный теплообмен, например на стороне 12, то вычисляем длину стороны и матрицу 6. 7. Формируем матрицу 8. Если на одной из границ задан тепловой поток, его надо учесть в правой части. Пусть, например, поток задан на стороне 23. Вычисляем длину стороны К вектору правых частей добавляем вектор: 9. Если на одной из границ задан конвективный теплообмен, его надо учесть в правой части. Пусть, например, теплообмен задан на стороне 13. Вычисляем длину стороны . Тогда к вектору правых частей добавляем вектор Определение расчетных величин в элементах (пункт 6). Если есть необходимость вычислить градиенты температур в узлах элементов
В цикле по элементам 1. Определяем координаты узлов текущего элемента и значения решений Для
2. Формируем
; 3. Вычисляем площадь треугольника 4. Вычисляем градиенты температур на всем элементе (они будут постоянны по элементу):
Поскольку градиенты получаются кусочно-постоянными на элементах, можно провести сглаживание значений. Например, если один узел принадлежит четырем элементам, то значению в узле приписывают среднее по этим четырем элементам. Поиск по сайту: |
Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. Студалл.Орг (0.005 сек.) |