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

4. Расчет температуры в холодильной камере

Расчет будем производить на основании [7]. В инженерной практики принято промышленные холодильные камеры описывать линейным дифференциальным уравнением 1‑го порядка с постоянными коэффициентами. Камеры являются весьма инерционными объектами. Так, например, постоянная времени Т рассматриваемой холодильной камеры равна 100 ч.

Однако промышленные холодильные камеры фактически являются многоемкостными объектами и более точно их следует описывать дифуравнениями выше первого порядка с тем, чтобы проверить насколько целесообразна их апроксимация дифуравнениями первого порядка.

В [7] предлагается описывать холодильную камеру линейным уравнением второго порядка с постоянными коэффициентами:

d²Dt dDt

Т₁ * Т₂ ¾¾ + (Т₁ + Т₂) ¾¾ + Dt = Dt_y

dt² dt

Используя уравнение (1) и, пренебрегая запаздыванием объекта, проводили исследования двухпозиционной системы регулирования в холодильной камере. Расчет проводили методом Рунге – Кутта (исходный текст программы приведен в приложении).

Метод Рунге – Кутта предназначен для дифференциального уравнения второго порядка вида (c учетом того, что Dt_y изменяет свое значение в зависимости от того работают компрессоры или происходит нагрев воздуха в камере за счет естественного притока тепла, или же при отрицательных температурах окружающей среды когда воздух в камере исскуственно подогревается за счет нагревательных элементов или температура в камере падает за счет естественного оттока тепла через стены камеры)

d²Dt

Dt»= ¾¾ = F (t, Dt, Dt’, Dt_y),

dt²

имеющий погрешность R~(h⁵), реализовался с помощью следующих формул [8]:

К₁=h * F(t_i; Dt _i; Dt’_i; Dt_yi);

К₂=h * F(t_i +(h/2); Dt _i +(h/2)* Dt’_i+(h/8)* К₁;Dt’_i + (К₁/2); Dt_yi);

К₃=h * F(t_i +(h/2); Dt _i +(h/2)* Dt’_i+(h/8)* К₁;Dt’_i + (К₂/2); Dt_yi);

К₄=h * F(t_i +h; Dt _i +h* Dt’_i+(h/2)* К₃;Dt’_i + К₃; Dt_yi);

Dt _i+1=Dt _i+ h*[Dt’_i +(К₁+ К₂ + К₃)/6];

Dt’_i+1 =Dt’_i + (К₁+ 2*К₂ + 2*К₃ + К₄)/6

Расчет проводился на участке от 0 ч до 200 ч при следующих начальных условиях:

Dt ₀= Dt0;

Dt’₀ = 0.

Вариант 1. Т₁ =100 ч, Т₂ =10 ч, заданный диапазон 0,5 – 1 °С, установившаяся температура при ее росте 10 °С и установившаяся температура при ее снижении минус 3 °С. При этом были получены следующие результаты: фактический диапазон поддержания температуры составил 0,45 – 1,25 °С, а период колебаний 54,2 часа. График переходного процесса и протокол работы приведен в приложении.

При описании холодильной камеры линейным дифуравнением первого порядка следующего вида:

dDt

Т ¾¾ + Dt = Dt_y

dt

провели аналогичные исследования системы двухпозиционного регулирования, т.е. полагали Т= Т₁ +2* Т₂=120 ч, а остальные данные были такими же, как и в варианте 1. При этом температура поддерживалась в заданном диапазоне (запаздыванием пренебрегали), а период колебаний составил 22,5 ч.

Из приведенных данных следует, что фактический диапазон поддержания температуры при более точном математическом описании холодильной камеры увеличивается в 1,6 раза а период колебаний возрастает в 2,5 раза. Следовательно для приведенных исходных данных рассматривать камеру в упрощенном варианте не следует.

Вариант 2. Т₂ = 0,5 ч, а остальные данные аналогичны варианту 1. По данному варианту получили, что температура поддерживается в заданном диапазоне, а период колебаний составил 21,3 ч. Исследования в упрощенном объекте (Т = 101 ч) показало, что период колебаний получился равным 19 ч. Как видим, для варианта 2 апроксимация холодильной камеры апериодическим звеном первого порядка вполне допустима.

Вариант 3. Поддержание рабочей температуры в камере происходит за счет работы электронагревателей при Т₁ =100 ч, Т₂ =15 ч, заданный диапазон 0,5 – 1 °С, установившаяся температура при ее росте 4 °С и установившаяся температура при ее снижении минус 5 °С. При этом были получены следующие результаты: фактический диапазон поддержания температуры составил 0,307 – 1,082 °С, а период колебаний 73 часа. График переходного процесса и протокол работы приведен в приложении.

Вариант 4. Т₂ =1.5 ч, а остальные данные аналогичны варианту 3. По данному варианту получили, что температура поддерживается в заданном диапазоне, а период колебаний составил 30,3 ч. Исследования в упрощенном объекте (Т = 103 ч) показало, что период колебаний получился равным 29 ч. Как видим, для варианта апроксимация холодильной камеры апериодическим звеном первого порядка вполне допустима.

Как мы можем видеть из рассмотренного выше целесообразно производить апроксимацию холодильной камеры апереодическим звеном первого порядка только в тех случаях когда постоянная времени Т₂ составляет не более чем 0,01…0,025 Т₁ то есть ее влияние на качество переходного процесса – несущественно.

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

На основании вышеизложенного можно сделать следующий вывод: в современной инженерной практике при использовании средств вычислительной техники необходимо для повышения точности расчетов рассматривать промышленную холодильную камеру, как апериодическое звено второго порядка (при Т₂ > 0,01…0,025 Т₁₎

Похожие работы

Проектирование холодильной установки

Скачать

35552

0

0

... батареи и воздухоохладители)? испарители для охлаждения промежуточного хладоносителя; отделители жидкости или защитные ресиверы (при надобности). Автоматизация испарительной системы холодильной установки с промежуточным хладоносителем (рассолом) предусматривает: а) автоматическое регулирование температуры воздуха в охлаждаемых помещениях; б) автоматическое регулирование температуры ...

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

Скачать

79391

6

0

... камере понижается ниже допустимой , то в данном случае включаются электронагреватели , которые встроены в ВО. Включением и выключением поддерживают заданную температуру в камере. 3.2 Работа узлов функциональной схемы автоматизации холодильного модуля. Основной регулируемой величиной в данной схеме есть температура воздуха в холодильной камере. Ее регулируют включением и выключением КМ , а зимой ...

Холодильное оборудование

Скачать

23500

0

0

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

Холодильники и холодильное оборудование: устройство и эксплуатация

Скачать

33925

1

3

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