Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
31
Добавлен:
18.07.2022
Размер:
1.3 Mб
Скачать

Выполнение лабораторной работы №4

по дисциплине «Моделирование систем управления»

на тему «Моделирование системы с распределенными параметрами методом конечных элементов»

Параметры по варианту:

Ширина

Высота

Высота

Материал

Материал

Материал

Мощность

дома, м

стен,

м

крыши, м

стен

крыши

печи

печи, кВт

 

 

 

 

 

 

 

 

5

2.4

 

3.6

Пенобетон

Сталь

Сталь

30

1.Провести моделирование объекта теплопроводности в MATLAB:

1.1.Создать двумерную геометрическую модель тела с требуемыми размерами и

формой.

1

1.2. Задать параметры материалов, источников тепла, начальные и граничные условия для решения уравнения теплопроводности.

Задание граничных условий для всех границ, в том числе внешних границ грунта, кроме левой, правой нижней границ грунта:

Для левой, правой, нижней границы грунта:

Задание параметров различных частей и источников тепла:

Снег:

Крыша:

2

Воздух:

Перекрытие:

Стены:

Печь:

3

Фундамент:

Грунт:

1.3. Провести моделирование и визуализировать распределение температуры, а также сетку конечных элементов.

Задание параметров отображения результата моделирования:

4

1.3.1. Моделирование для 200 часов:

Параметры для решения уравнения теплопроводности:

Сетка конечных элементов:

Распределение температуры:

2.Увеличив время моделирования, добиться достижения целевых температур в заданных частях тела.

5

Моделирование для случая, когда снег на крыше начнет таять:

3.Исследовать точность моделирования, сравнив поле температур при различных сетках конечных элементов.

Примечание для пункта 3: нужно взять значения x, y, которые не меньше участка окна в pdetool, где отображается дом. Шаг как можно мельче (если будет слишком мелкий программа не посчитает, выведет ошибку, тогда взять другой). Когда строим contour, нужно, чтобы весь дом отображался на картинке, если не отображается, то раздвинуть границы x, y при разбиении и заново построить рисунки для нового uxy.

Для решения, полученного в пункте 2:

uxy1=tri2grid(p,t,u(:,61),[0:0.05:9],[-1.2:0.05:8])

uxy=uxy1

uxy(isnan(uxy))=0

surf(uxy)

6

contour(uxy)

7

При мелкой сетке:

uxy2=tri2grid(p,t,u(:,61),[0:0.05:9],[-1.2:0.05:8]) uxy=uxy2

uxy(isnan(uxy))=0

surf(uxy)

8

contour(uxy)

Построение графиков разности uxy1-uxy2: delta_12=uxy1-uxy2

surf(delta_12)

9

contour(delta_12)

Вычисление среднеквадратичной и максимальной ошибок: delta_12(isnan(delta_12))=0 sr_err12=norm(delta_12)/numel(delta_12)

Результат: 0.0700 max_err12=max(max(abs(delta_12)))

Результат: 261.3207

10

Соседние файлы в папке Лабораторные работы