Новости

Моделирование и автоматическое проектирование устройств связи. Моделирование теплового поля однородного стержня

Работа добавлена:






Моделирование и автоматическое проектирование устройств связи. Моделирование теплового поля однородного стержня на http://mirrorref.ru

Министерство образования и науки Российской Федерации

Новосибирский государственный  технический  университет

Факультет радиотехники, электроники и физики

Кафедра КТРС

Расчетно-графическая работа по дисциплине

«Моделирование и автоматическое проектирование устройств связи»

«Моделирование теплового поля однородного стержня»

Вариант 11

Студент: Цвик П.Ю.                                                Проверил: Девятков Г.Н.

Группа: РКС10-31                                                    Оценка:

      Факультет: РЭФ

Дата сдачи: 06.10.16

Новосибирск

2016

  1. Исходные данные

Последние цифры студенческого шифра

11

Метод решения задачи

Метод явных разностных схем

5

Длина стержня,м

9

,°С

80

,°С

0

Коэффициент теплоёмкости С, МДж/м³К

2.4

Коэффициент теплопроводимости К, Дж/мКс

200

  1. Постановка задачи

Имеется однородный стержень длинойL, один конец которого соединен с идеальным теплопроводом. В момент времени  к свободному концу стержня прикладывается «тепловая ступенька». Мощность источника тепла достаточна для поддержания постоянной температуры  на свободном конце стержня, а отвод тепла происходит только за счет теплопроводности стержня в продольном направлении.

Требуется формализовать задачу, получив уравнения теплопроводности непосредственно в конечных разностях, и решить ее методом явных разностных схем, определив температурный профиль по длине   стержня в заданный момент времени .

Рис.1. Однородный стержень

  1. Формализация задачи

Данная задача - одномерная, так как по условию отвод тепла происходит только за счет теплопроводности стержня в продольном направлении, то есть тепловые потоки в поперечных направлениях стержня равны нулю. Выделив из стержня элементарный объем, можно записать уравнение баланса количества теплоты, предварительно разделив обе части уравнения на объем элемента и рассматриваемый период времениτ:

                                            ,                                      (1)

гдеJx+,Jx-- удельные плотности входящего и выходящего тепловых потоков;

hx - длина элементарного объема;

θis+1- θis - приращение температуры элементарного объемаi за время  τ;

s-номер шага по времени, т.е. номер временного слоя.

Рис.2. Элементарный объем

Считая, что свойства среды линейны, на основании закона Фурье рассматриваемые тепловые потоки можно выразить через разности температур соседних объёмов и получить уравнение в явной форме:

                                           ,                                (2)

Уравнение (2) содержит одну неизвестную - .

  1. Алгоритм численного решения и его описание

Данное уравнение необходимо преобразовать относительно :

                                            ,                               (3)

Полученному уравнению соответствует форма расчётной ячейки, представленная нарис.3. Расчётная ячейка позволяет наглядно представить, значения температур каких узлов следует подставлять в уравнение (3) при вычислении неизвестной . Для минимизации расчетов принимаем .

Известны начальные и граничные условия: температура свободного конца всегда равна θ1,а закреплённого на теплоотводе -θср.В начальный момент времени  температура стержня во всех узлах равнаθср.

Таким образом, известна температура всех узлов стержня на нулевом временном слое, а также температуры крайних точек стержня. Это позволяет, начав движение расчётной ячейки в нулевом слое, найти значения температур точек стержня на последующих временных слояхs+1.

Рис.3. Схематичная форма расчетной ячейки в случае явной разностной схемы и ее начальное положение в момент времени .

Параметр  для явной разностной схемы выбирается из условия устойчивости вычислительного процесса:

                                                                                                                   (5)

Параметрhx выбирается из условия заданной точности решения. Приемлемая точность может быть получена при разделении стержня по длине на восемь частей.

  1. Вычисление

Выбираем шагиτ иhxв соответствии с формулой (5) и рекомендациям к решению:

Из заданного критерия точности выбираем значениеτ удобным для дальнейших математических операций, а также обеспечиваем целочисленное значение количества уровнейs:

При вычислении температуры элементов стержня в различных температурных слоях используется постоянный коэффициент А:

Так как выше было обговорено, что для обеспечения заданной точности решения по длине стрежня взято 8 узлов, поэтому:

Таблица 1. Решение задачи с помощью явной разностной схемы

Счетчики циклов

Вычисления

S

(0,…,Smax-1)

i

(L/hx-1,…,1)

(oC)

s=0

i= L/hx-1

= 0

i= L/hx-2

=0

i= L/hx-3

=0

i= L/hx-4

=0

i= L/hx-5

=0

i= L/hx-6

=0

i= L/hx-7

=0

s=1

i= L/hx-1

= 26.337

i= L/hx-2

= 0

i= L/hx-3

=0

i= L/hx-4

=0

i= L/hx-5

=0

i= L/hx-6

=0

i= L/hx-7

=0

s=2

i= L/hx-1

= 35.333

i= L/hx-2

=8.671

i= L/hx-3

= 0

i= L/hx-4

=0

i= L/hx-5

=0

i= L/hx-6

=0

i= L/hx-7

=0

s=3

i= L/hx-1

= 41.261

i= L/hx-2

= 14.594

i= L/hx-3

= 2.855

i= L/hx-4

= 0

i= L/hx-5

=0

i= L/hx-6

=0

i= L/hx-7

=0

s=4

i= L/hx-1

= 45.235

i= L/hx-2

= 19.508

i= L/hx-3

= 5.78

i= L/hx-4

= 0.94

i= L/hx-5

= 0

i= L/hx-6

=0

i= L/hx-7

=0

s=5

i= L/hx-1

= 48.211

i= L/hx-2

=23.458

i= L/hx-3

=8.706

i= L/hx-4

=2.224

i= L/hx-5

=0.309

i= L/hx-6

=0

i= L/hx-7

=0

s=6

i= L/hx-1

= 50.527

i= L/hx-2

= 26.751

i= L/hx-3

=11.429

i= L/hx-4

=3.728

i= L/hx-5

=0.838

i= L/hx-6

=0.102

i= L/hx-7

=0

s=7

i= L/hx-1

= 52.403

i= L/hx-2

=29.534

i= L/hx-3

= 13.938

i= L/hx-4

=5.312

i= L/hx-5

=1.547

i= L/hx-6

=0.311

i= L/hx-7

=0.034

s=8

i= L/hx-1

= 53.959

i= L/hx-2

= 31.928

i= L/hx-3

= 16.232

i= L/hx-4

= 6.912

i= L/hx-5

=2.379

i= L/hx-6

=0.626

i= L/hx-7

=0.114

s=9

i= L/hx-1

= 55.279

i= L/hx-2

= 34.014

i= L/hx-3

= 18.331

i= L/hx-4

= 8.488

i= L/hx-5

= 3.294

i= L/hx-6

=1.035

i= L/hx-7

=0.245

s=10

i= L/hx-1

= 56.417

i= L/hx-2

= 35.852

i= L/hx-3

= 20.254

i= L/hx-4

= 10.019

i= L/hx-5

= 4.26

i= L/hx-6

=1.519

i= L/hx-7

=0.424

Таким образом, получены значения температуры рассматриваемого стержня в восьми его точках (узлах) в десяти временных слоях. Согласно заданию, требуется получить температурный профиль стержня в момент времени .

Данному профилю соответствует временной слой  под номером десять,так как

Рис.4. Распределение температуры по стержню

Теперь проведем расчет с учетом нарушенного условия спектральной устойчивости, т.е. формулы (5).

Примем и выполним перерасчет необходимых параметров:

Подставим данные в уравнение (3) и выполним математические вычисления, построим необходимые графики.

Полученные значения представим в виде таблицы.

Таблица 2. Проверка условия спектральной устойчивости

Счетчики циклов

Вычисления

S

(0,…,Smax-1)

i

(L/hx-1,…,1)

(oC)

s=0

i= L/hx-1

= 0

i= L/hx-2

=0

i= L/hx-3

=0

i= L/hx-4

=0

i= L/hx-5

=0

i= L/hx-6

=0

i= L/hx-7

=0

s=1

i= L/hx-1

=131.687

i= L/hx-2

=0

i= L/hx-3

=0

i= L/hx-4

=0

i= L/hx-5

=0

i= L/hx-6

=0

i= L/hx-7

=0

s=2

i= L/hx-1

=71.317

i= L/hx-2

=43.354

i= L/hx-3

=0

i= L/hx-4

=0

i= L/hx-5

=0

i= L/hx-6

=0

i= L/hx-7

=0

Рис.5. График температурного профиля по длине стержня в момент времени

Заключение

Метод явных разностных схем благодаря простоте и наглядности позволяет получать решение путем определения значений только одного неизвестного параметра. Требуемые на вычисления относительно большие затраты времени незначительны с применением автоматизированных средств расчета. С помощью данного метода был произведён расчет теплового поля однородного стержня, с одного конца закрепленного на теплоотводе и нагреваемого с другого.

В процессе решения получается полная картина распределения температуры по всей области стержня. Имеется возможность моделирования решения на персональном компьютере с последующим дополнением и усложнением модели. С помощью данного метода был произведён расчет теплового поля однородного стержня, с одного конца закрепленного на теплоотводе и нагреваемого с другого.

По длине стержня температура убывает практически экспоненциально, площадь под кривой распределения температуры возрастает с ростом номера временного слоя.

Нарушение условия устойчивости приводит к получению некорректных результатов.

Список литературы

  1. Маквецов Е.Н. Модели из кубиков  – М.: Сов. радио, 1978, с.52-62,181-184.

  1. Малика А.С. Автоматизация конструирования РЭА. – М.:

 Высшая школа, 1980, с.180-188, 190, 205-206.

Моделирование и автоматическое проектирование устройств связи. Моделирование теплового поля однородного стержня на http://mirrorref.ru


Похожие рефераты, которые будут Вам интерестны.

1. Реферат МОДЕЛИРОВАНИЕ ЭЛЕКТРОСТАТИЧЕСКОГО ПОЛЯ

2. Реферат Моделирование на основе парной регрессии и корреляции Моделирование одномерных временных рядов

3. Реферат ПРАКТИЧЕСКОЕ РУКОВОДСТВО ДЛЯ ВЫПОЛНЕНИЯ ЛАБОРАТОРНЫХ РАБОТ ПО КУРСУ МОДЕЛИРОВАНИЕ ЭЛЕКТРОННЫХ УСТРОЙСТВ

4. Реферат Классификация моделей по способу построения. Интуитивное и научное моделирование. Натуральное и аналоговое моделирование

5. Реферат Статистическое моделирование связи методом регрессионно-корреляционного анализа

6. Реферат Моделирование технологических процессов как инструмент проектирование терминалов

7. Реферат Модели и моделирование: основные понятия (объект, модель, моделирование), основные цели моделирования

8. Реферат Устройства связи с объектом управления. Классификация. Примеры устройств связи с объектом управления

9. Реферат Исследование влияния обратной связи на показатели усилительных устройств

10. Реферат Моделирование организаций