19965

Решение задачи о поле температуры в облучательном устройстве при отсутствии утечек тепла в торцы

Лекция

Физика

Поставить и решить вспомогательную задачу Б и закончить рассмотрение задачи о радиальном распределении температуры в облучательном устройстве при отсутствии утечек тепла в торцы. Обосновать необходимость использования метода конечных элементов (МКЭ) для расчета полей температуры в облучаемых образцах. Приступить к постановке задачи расчета поля температуры МКЭ для цилиндрического образца.

Русский

2013-08-13

39.33 KB

0 чел.

Конспект занятия 13.

Цель.

    Поставить и решить вспомогательную задачу Б и закончить рассмотрение задачи о  радиальном распределении температуры в облучательном устройстве при отсутствии утечек тепла в торцы. Обосновать необходимость использования метода конечных элементов  (МКЭ) для расчета полей температуры в облучаемых образцах. Приступить к постановке задачи расчета поля температуры МКЭ для цилиндрического образца.

План.

1. Постановка и решение вспомогательной задачи Б.

2.Решение задачи о поле  температуры в облучательном устройстве при отсутствии утечек тепла в торцы.      

3. Постановка задачи расчета поля температуры МКЭ для цилиндрического образца.

   Поле температуры в образце (задача Б)

    На поверхности цилиндра с коэффициентом теплопроводности λ0,1   задана температура Т1   , внутри цилиндра действуют внутренние источники тепла qv01 , в центре цилиндра температура имеет экстремум.

Граничные условия:

 dT/dr | r= 0                                                                                           (19)

T | r= R1= Т1                                                                                          (20)

    Поле температуры описывается уравнением (13) и (14).

Из (19)  C1 = 0, тогда  из (20) определяем:

С2 = Т1+ qv,0,1 R21/4 λ0,1

Поле температуры в цилиндре (образце) имеет вид:

Рис.3.7 Модель расчета поля температуры при отсутствии утечек тепла в торцы.

λk.k+1

Rk

Rn

3

1

2

K

K+1

K+2

K+3

n-1

n

qv01v01

qv23

λ23

qvk.k+1

Q

α.Tcp

T=T(r)

Qk

λk,k+1

qvk,k+1

Rk

Tk

Tk+1

А

dT/dr|r=0=0

T=T(r)

λ01

qv01

R1

T1

Б

Т=Т1+ qv,0,1 (R21-r2)/4 λк, к+1

Поток тепла с поверхности цилиндра:

Qk = - 2π λ0,1  R1 dT/dr | r= R1 = πqv,0,1R21 = πqv,0,1(R21 R20),

где R0 = 0

Определяем потоки тепла, пользуясь результатами задач, рассмотренных выше:

        n

Qn = Σ πqv,k,k+1(R2k+1 R2k) при R0 = 0.

       k=0

     Используя краевое условие (6), имеем:

Tn -Tc = Qn/2 παRn

     Определяем перепады температуры:

- на оболочке

Тn-1 - Тn = Av,n-1,n+ An-1,n ,

- в газовой прослойке:

Tn -Tc = Qn-2,n-1 /hn-2,n-1

- на к-ом экране:

Тk - Тk+1 = Av,k,k+1+ Ak,k+1 ,

- в к-1 прослойке:

Тk-1 - Тk = Qk-1,k /hk-1,k

-в экране с радиусами R2 и  R3 :

T2T3 = Av,2,3+ A2,3 ,

- в прослойке с радиусами R1 и  R2:

Т1 - Т2 = Q1,2 /h1,2 ,

- в образце:

Т0 1= qv,0,1 R21/4 λ0,1

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

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

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

    Для расчета температурных полей в образце реакторной установки целесообразно воспользоваться методом конечных элементов.

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

1.Геометрические условия задают цилиндрический осе симметричный образец.

2.Физические   условия   задают    распределение     источников

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

3.Временные условия рассматривают стационарную задачу:

dT/ =0                                                                                          (21)

4.Граничные условия.

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

- условия первого рода:

T|z=0, 0≤  rR = T (0, r)                                                                         (22)

T|z=H, 0≤  rR = T (H, r)                                                                       (23)

- условия третьего рода:

- λ dT/dr |z=0, 0≤  rR = α (0,r) [T (0, r) – Tc0]                                    (24)

- λ dT/dr |z=H, 0≤  rR = α (H,r) [T (H, r) – TcH]                                (25)

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

Q= 2πRα(z,r) [T (z, r) – Tcr]                                                             (26)

    Решение задачи методом конечных элементов.

    Дискретизация геометрической области проводится по схеме представленной в верхней части рис. 3.8.

    Определение стационарных двумерных полей температуры основано на простейшем варианте метода конечных элементов. Ищется решение стационарного уравнения теплопроводности:

                            

div [ λ(T) grad T( r )] +  qv(r) =0 ,                                                    (27)         

где

Т(r) - температура образца;

λ(Т) - коэффициент теплопроводности в общем случае, зависящий от температуры;

qv(r)- плотность внутренних источников тепла может быть функцией координат.

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

элементов и с учетом симметрии задачи цилиндрический образец разбивается на N кольцевых элементов и  М элементов по высоте.  

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

    В рассматриваемом случае уравнения теплового баланса    элементов принимают вид:

N(i)

Σ γ(i,j)[T(i)-T(j)] + qv(i)S(i)+QL(i) = 0                                           (28)

j=1

где

S(i)- площадь получаемого  при таком разбиении элемента;

Т(i)- температура элемента;

qv(i) плотность внутренних источников тепла;

QL(i)- поток тепла в элемент из внешней среды;

γ( i,j)- коэффициент, характеризующий перенос тепла между соседними i-ым и j -ым элементами;

N(i)- число элементов, обменивающихся теплом с элементом, равно четырем во внутренней области и трем для элементов, лежащих на границе области.

При составлении системы уравнений (28) предполагалось, что потоки тепла Q между соседними элементами пропорциональны разности температур в этих элементах:

Q = γ( i,j) [T(i)-T(j)]                                                                          (29)

    Выражение, определяющее  γ( i,j)  , может быть получено при рассмотрении  соотношения для потока тепла между   i-ым  и  j -ым  элементами в радиальном направлении:

Q = λ Lij grad T  | ij                                                                            (30)

где λ - коэффициент теплопроводности материала;   L- протяженность границы между элементами; grad T  | ij - градиент температуры на границе между   i-ым и  j -ым элементами.

H 

r

z

r 

Рис. 3.8

0

0

0

0

0

0

0

0

M*N

Рис.3.8. Схема расположения конечных элементов и структура матрицы [B].


 

А также другие работы, которые могут Вас заинтересовать

41706. Программная реализация алгоритмов линейной структуры на языке программирования Turbo Pascal 153.86 KB
  Описание объявление данных содержит упоминание всех объектов используемых в программе и включает в себя: раздел подключаемых библиотек модулей определяется служебным словом USES и содержит имена подключаемых модулей: uses CRT Grph; раздел описания меток: любой оператор в программе может быть помечен меткой. Типы данных. Под типом данных понимается множество допустимых значений этих данных а также совокупность операций над ними. Тип определяет также и формат внутреннего представления данных в памяти компьютера.
41707. Знакомство и приемы работы в интегрированной среде языка Турбо-Паскаль 52.4 KB
  на экране монитора ПК могут формироваться следующие видимые компоненты среды: меню окна диалоговые окна и строка состояния. Меню называется прямоугольный участок экрана содержащий ряд ключевых слов и предназначенный для выбора одного из предлагаемых вариантов работы. Содержащиеся в меню ключевые слова обозначают возможные альтернативы называемые в дальнейшем опциями option выбор. Интегрированная среда Турбо-Паскаля обеспечивает иерархическую систему меню.
41708. ИССЛЕДОВАНИЕ ТИПОВЫХ ЗВЕНЬЕВ ЛИНЕЙНЫХ АСР 165.27 KB
  Переходная характеристика Время сек Импульсная характеристика Время сек Рисунок 1 Переходная и импульсная характеристика звена W1 Переходная характеристика Время сек Импульсная характеристика Время сек Рисунок 2 Переходная и импульсная характеристика звена W2 Переходная характеристика Время сек Импульсная характеристика Время сек Рисунок 3 Переходная и импульсная характеристика звена W3 Переходная характеристика Время сек Импульсная характеристика Время сек Рисунок 4 Переходная и импульсная характеристика звена W4...
41709. Расчет поперечной рамы стального каркаса одноэтажного здания на действие постоянной нагрузки 702.68 KB
  Поперечное сечение в виде 2х уголков с параметрами уголка 160×100×10 Высота низа колонн 8 м Поперечное сечение в виде трубы с параметрами 720×12 Высота верха колонн 6 м Поперечное сечение в виде трубы с параметрами 720×8 q = 24 кН м P = 32 кН ℓ = 32 м Цели лабораторной работы: Произвести расчет поперечной рамы стального каркаса одноэтажного здания с помощью ПК ЛИРА Определить для основных сечений колонн и ригеля максимальные значения нормальных и касательных напряжений Сделать вывод о работе конструкции Ход работы: Создание расчетной...
41710. Команды MSDOS 251.68 KB
  Название команды Синтаксис команды Создание файла с консоли copycon имя файла Удаление файла del имя файла Переименование файла ren имя файла 1 имя файла 2 Редактирование файла edit имя файла Переход на диск имя диска Переход в каталог cd путь Создание каталога md имя каталога Удаление каталога rd имя каталога Очистка экрана Cls Вывод содержимого файла на экран type имя файла Копирование файла copy путь 1 что копируется путь 2 куда копируется Поиск файла filefind имя файла Работа с командной строкой Prompt Информация о команде команда ...
41711. Индекс FTSE 100 109.84 KB
  Индекс начал рассчитываться с 3 января 1984 года с уровня 1000 пунктов. Компании чьи акции учитываются в расчете индекса FTSE 100 должны удовлетворять условиям выставляемым FTSE Group: а должны входить в список Лондонской фондовой биржи; б стоимость акций индекса FTSE 100 выражается в фунтах или евро;в пройти тест на принадлежность к определенному государству;г акции FTSE 100 должны быть в свободном обращении и быть легко ликвидными. Индекс FTSE 100 включает 100 компаний британской экономики хотя на самом деле список состоит из...
41712. Основы работы в ГИС ArcView 1.75 MB
  Интерфейс пользователя rcView Окно программ Окно приветствия Строка меню Строка кнопок Строка инструментов Окно вида Меню работы с файлами Меню редактирования Меню вид Меню тема Меню графика Меню окно Строка кнопок:
41713. Классы. Основные понятия 260.14 KB
  Основные понятия Цель работы: Познакомиться с понятием класса. Продемонстрировать работу с классом создание экземпляра класса изменение значений полей вызов перегруженных методов Выполнить требования к лабораторным работам см. Контрольные вопросы: Что такое метод Что такое возвращаемый тип каким он может быть Что значит void Что такое параметры метода какие параметры бывают Что такое перегрузка методов и для чего она нужна Что такое класс Какие члены класс может содержать Что такое конструктор для чего служит...