|
|||||||
АвтоАвтоматизацияАрхитектураАстрономияАудитБиологияБухгалтерияВоенное делоГенетикаГеографияГеологияГосударствоДомДругоеЖурналистика и СМИИзобретательствоИностранные языкиИнформатикаИскусствоИсторияКомпьютерыКулинарияКультураЛексикологияЛитератураЛогикаМаркетингМатематикаМашиностроениеМедицинаМенеджментМеталлы и СваркаМеханикаМузыкаНаселениеОбразованиеОхрана безопасности жизниОхрана ТрудаПедагогикаПолитикаПравоПриборостроениеПрограммированиеПроизводствоПромышленностьПсихологияРадиоРегилияСвязьСоциологияСпортСтандартизацияСтроительствоТехнологииТорговляТуризмФизикаФизиологияФилософияФинансыХимияХозяйствоЦеннообразованиеЧерчениеЭкологияЭконометрикаЭкономикаЭлектроникаЮриспунденкция |
Алгоритм. 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. Формируем матрицу 5. Если на одной из сторон задан конвективный теплообмен, например на стороне 12, то вычисляем длину стороны 6. 7. Формируем матрицу 8. Если на одной из границ задан тепловой поток, его надо учесть в правой части. Пусть, например, поток задан на стороне 23. Вычисляем длину стороны 9. Если на одной из границ задан конвективный теплообмен, его надо учесть в правой части. Пусть, например, теплообмен задан на стороне 13. Вычисляем длину стороны Определение расчетных величин в элементах (пункт 6). Если есть необходимость вычислить градиенты температур в узлах элементов
В цикле по элементам 1. Определяем координаты узлов текущего элемента и значения решений Для
2. Формируем
3. Вычисляем площадь треугольника 4. Вычисляем градиенты температур на всем элементе (они будут постоянны по элементу):
Поскольку градиенты получаются кусочно-постоянными на элементах, можно провести сглаживание значений. Например, если один узел принадлежит четырем элементам, то значению в узле приписывают среднее по этим четырем элементам. Поиск по сайту: |
Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. Студалл.Орг (0.004 сек.) |