Физика / Реферат: Методы оценки температурного состояния

Реферат: Методы оценки температурного состояния

Рис.1.2 Пример нерегулярной согласованной разностной сетки.

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

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

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

Построение разностных схем может осуществляться на основе метода конечных элементов. Для построения конечномерного подпространства исходная расчетная область разбивается на некоторые элементарные ячейки. В двумерном случае в качестве таковых наиболее подходящими являются треугольники, причем внутри таких ячеек приближенное решение является линейной функцией. Такого вида сетки выбраны в связи с возможностью решения задач в областях достаточно произвольной формы. На рис.1.3 показана равномерная сетка с шагом .

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

В рассматриваемой области фиксируется конечное число точек, которые называют узлами (узловыми точками). Непрерывная величина аппроксимируется моделью, состоящей из отдельных элементов. На каждом из этих элементов исследуемая непрерывная величина аппроксимируется кусочно-непрерывной функцией. Выбираются аппроксимирующие базисные функции в виде кусочных полиномов малой степени или полиномов более высоких степеней (квадратичных, кубических и др.) Полином, связанный с данным элементом называют функцией элемента. Далее строится разностная схема.

Разностные схемы для нестационарных задач.

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

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

У данной задачи есть и начальные, и граничные условия, поэтому задача является нестационарной (смешанной) краевой. Задача имеет нелинейный характер, т.е. теплофизические свойства среды зависят от температуры и граничные условия нелинейны.

Метод конечных разностей.

В качестве метода решения системы дифференциальных уравнений выбирается численный математический метод конечных разностей - широко известный и простейший метод интерполяции. Метод конечных разностей означает по сути обратный переход от дифференциальной модели к интегральной. При осуществлении данной методики строится конечно-разностная сетка и записываются конечно-разностные аналоги дифференциальных уравнений теплопроводности (разностная схема). Для аппроксимации дифференциальных уравнений теплопроводности применяется неявная консервативная итерационная разностная схема, реализуемая градиентным методом покоординатного спуска (Гаусса-Зейделя), являющимся классическим итерационным методом решения системы линейных уравнений. Неявной она является потому, что содержит несколько неизвестных значений функции на новом слое. Это означает, что значение функции нельзя явно выразить через значения функции на данном слое. Такая схема является безусловно устойчивой. Неявность разностной схемы достигается применением итерационной процедуры на каждом временном слое. Решение в узлах сетки получается приближенным (разностным). Поскольку одна из переменных имеет физический смысл времени t , то сетка строится так, чтобы среди ее линий были линии t= t_m _, где m - номер индекса дискретного момента времени. То есть переменная t не непрерывна, а увеличивается на дискретное значение. Решение численной задачи получается в виде таблицы.

Экономичные разностные схемы нестационарной теплопроводности.

Поскольку при использовании неявных схем вычислительные затраты высоки, применяют методы реализации разностных схем, которые по вычислительной реализации были бы аналогичны явным схемам. К таким методам относятся явный итерационный метод, метод переменных направлений, попеременно-треугольный метод, итерационный метод с эллиптическим оператором B . Для явных схем число арифметических операций, приходящихся на один узел сетки не зависит от общего числа узлов. Разностные схемы метода переменных направлений основываются на представлении оператора по пространственным переменным в виде суммы операторов, каждый из которых является одномерным [1], [2].

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

2. Постановка нестационарной краевой задачи теплопроводности в

системе, включающей прошивную оправку

2.1 Условия однозначности или краевые условия задачи

Геометрические условия.

Оправка - это сплошное тело сложной формы (при решении задачи термоупругости не рассматривается возможное наличие в оправке специальных каналов для подачи охлаждающей жидкости, хотя они часто применяются на практике). Диаметр оправки зависит от внутреннего диаметра гильзы. Оправка подразделяется на участки различной геометрической формы: сферическую часть, коническую часть до пережима, коническую часть после пережима и часть штока, примыкающую к оправке. Длины этих участков рассчитываются по известным формулам.

Постановка краевой задачи зависит от выбора системы координат. Простейший подход к решению задач в нерегулярных областях состоит в использовании криволинейных координат, в которых расчетная область становится регулярной (понятия регулярной и нерегулярной областей были рассмотрены в разделе 1). Для сферического участка I принята сферическая система координат. Для участков II, III, IV принята цилиндрическая система координат.

Диаметр оправки на третьем участке равен:

,(2.1)