АвтоАвтоматизацияАрхитектураАстрономияАудитБиологияБухгалтерияВоенное делоГенетикаГеографияГеологияГосударствоДомДругоеЖурналистика и СМИИзобретательствоИностранные языкиИнформатикаИскусствоИсторияКомпьютерыКулинарияКультураЛексикологияЛитератураЛогикаМаркетингМатематикаМашиностроениеМедицинаМенеджментМеталлы и СваркаМеханикаМузыкаНаселениеОбразованиеОхрана безопасности жизниОхрана ТрудаПедагогикаПолитикаПравоПриборостроениеПрограммированиеПроизводствоПромышленностьПсихологияРадиоРегилияСвязьСоциологияСпортСтандартизацияСтроительствоТехнологииТорговляТуризмФизикаФизиологияФилософияФинансыХимияХозяйствоЦеннообразованиеЧерчениеЭкологияЭконометрикаЭкономикаЭлектроникаЮриспунденкция

Алгоритм. 1. Считываем сетку, элементы, материал из файла - формируем массивы mesh, material, elements

Читайте также:
  1. IV. Алгоритм действий командира (начальника) при увольнении военнослужащего в связи с невыполнением им условий контракта
  2. LZW-модификация алгоритма Лемпеля-Зива
  3. Zip–модификация алгоритма Лемпеля-Зива
  4. А.3.3. Алгоритм медикаментозного лікування
  5. Алгоритм
  6. Алгоритм
  7. АЛГОРИТМ
  8. Алгоритм 1.11. Пошук невідкладних дій (перша медична допомога) симптоматичної допомоги при гострих струєннях.
  9. Алгоритм 1.4. Діагностичний і лікувальний (перша медична допомога) пошук при гіпертонічній кризі
  10. Алгоритм 2.4. Транспортна іммобілізація
  11. Алгоритм First Come First Served (FCFS)

1. Считываем сетку, элементы, материал из файла - формируем массивы 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. Вычисляем градиенты температур на всем элементе (они будут постоянны по элементу):

 

Поскольку градиенты получаются кусочно-постоянными на элементах, можно провести сглаживание значений. Например, если один узел принадлежит четырем элементам, то значению в узле приписывают среднее по этим четырем элементам.


1 | 2 | 3 | 4 |

Поиск по сайту:

Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. Студалл.Орг (0.004 сек.)