LS-Sb90926
.pdf
ние почти всегда сопряжено со значительным увеличением времени нагрева, что нецелесообразно. Эта задача не имеет собственного практического значения, если не предусматривать в конструкции нагревателя специальных мер по защите поверхности нагреваемых тел. Вторая задача на максимум точности нагрева также чрезмерно академична и не имеет прикладного значения, хотя ей и посвящено большое количество работ.
Таким образом, наиболее актуальной является первая задача – максимального быстродействия, к тому же при ее решении, как правило, автоматически решается третья задача, т. е. обеспечивается минимальный расход электроэнергии. Рассмотрим данную задачу подробнее.
Время нагрева заготовки tн = tн ( p0 (t ), q0 (t ), f (t )) является функциона-
лом от законов изменения удельной мощности, тепловых потерь и частоты тока. Возможно трехканальное управление:
а) управление частотой тока (практически не используется);
б) управление удельной мощностью тепловых потерь. Значение q0 нарастает к концу нагрева и существенно влияет на температурное поле. Влияние это двоякое. С ростом q0 распределение температуры по сечению
утрачивает монотонность |
и приобретает |
максимум (рис. 2.1): |
Tп < Tmax , |
||
Tmax = = T (βmax ) , βmax ≈ |
|
|
, где ηт – |
|
|
|
ηт |
термический КПД |
индуктора; |
||
Tп - температура на поверхности заготовки. При глубинном нагреве возможен перегрев внутренних слоев. В то же время увеличение тепловых потерь сокращает время нагрева, так как позволяет поднять уровень мощности в заготовке. Очевидно, q0 и, следовательно, толщина тепловой изоляции должны варьироваться с целью минимизации энергозатрат при ограничении по температуре в точке максимума;
|
|
|
T |
Tmax |
|
|
|
|
|
|
|
|
|
Tп |
|
p0 , |
|
|
|
|
|
|
|
p0(t) |
|
|
|
|
|
|
|
T |
|
|
Tц |
|
|
|
|
|
Tп |
T |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Tц |
|
|
|
|
|
|
0 |
|
t |
0 |
|
βн |
1 |
β |
|
Рис. 2.1 |
|
Рис. 2.2 |
|
|
|
|
|
|
|
32 |
|
|
|
|
в) главный канал управления - функция |
p0 (t ) . Из теории теплопровод- |
|||
ности известно, что минимальным время нагрева будет при Tп = max = const . |
||||
Для этого надо, чтобы |
p0 (t ) t=0 ® ¥ (рис. 2.2). Такой режим нереализуем. |
|||
Всегда p0 ≤ pmax из соображений максимально допустимой токовой нагруз- |
||||
ки в обмотке индуктора и из условия |
p0 , |
|
||
ограничения по механическим напряже- |
T |
|
||
ниям во избежание термических трещин |
Tп |
|
||
p0max |
T |
|||
(рис. 2.3). |
|
|
Tц |
|
|
|
p0 |
|
|
Кривые изменения удельной мощно- |
p0(τ) |
|
||
сти, удовлетворяющие первой задаче оп- |
|
|
||
тимизации, |
поддаются |
аналитическому |
|
t |
расчету в линейном приближении в виде |
Рис. 2.3 |
|
||
θ (τ) и Ki (τ) . На практике мощность ре- |
|
|
||
гулируется ступенчато. В такой постановке задача оптимального управления |
||||
сводится к экстремальному уравнению, т. е. к поиску минимума функции |
||||
(n − 1) переменных, где n – число ступеней. |
|
|
||
|
2.1.3. Произвольный закон управления мощностью |
|
||
Ставится задача определения закона управления мощностью при произвольном изменении Ki (τ) . Воспользуемся критериальной формой уравнения теплопроводности
∂θ |
- Ñ2q = Ki (t) w(a, b) . |
(2.1) |
¶t |
|
|
Примем Ki (τ) = Ki0 = const . Тогда, как было показано ранее, решением уравнения (2.1) будет θ (β, τ) = Ki0ϕ (α, β, τ) . Если же Ki (τ) - любая гладкая
функция, то решение уравнения (2.1) можно найти с помощью интеграла наложения:
q(b, t) = ∫τ Ki (t)j¢(t - x)d x . |
(2.2) |
0 |
|
При решении интегрального уравнения (2.2) воспользуемся операторным методом и приведем его основные положения. С помощью операторного метода, или преобразования Лапласа, функция-оригинал f (t ) переводится в
33
функцию-изображение f (s) , или, короче, оригинал f (t ) – в изображение
|
|
(s) . Связь между f |
(t ) и |
|
(s) называется соответствием и обозначается |
||||
|
f |
f |
|||||||
|
следующим образом: f (t ) { |
|
(s )}, или |
|
(s) = |
||||
посредством символа |
|
f |
|||||||
f |
|||||||||
=−1{ f (s)}. Преобразование функций операторным методом осуществляется
спомощью таблицы соответствий (см., например, [1]) и правил (теорем) отображения выполняемых над функциями операций.
Интегральная комбинация функций Ki (τ) и ϕ′ (τ − ξ) , содержащаяся в правой части уравнения (2.2), называется сверткой этих функций и символически обозначается Ki ϕ′ (читается: «функция Ki , свернутая с функцией ϕ′ »). Итак, имеется следующее правило (теорема) о свертке:
τ
∫ Ki(τ)ϕ′(τ − ξ)dξ = Ki ϕ′ K ϕ'.
0
Следовательно, уравнение (2.2) в изображениях примет вид θ = Ki ϕ, откуда
Ki = θ
ϕ′ . Таким образом, поставленная задача решена в изображениях. Чтобы найти решение в оригиналах, необходимо найти изображения θ и
ϕ '. Вид изображения θ определяется заданным режимом нагрева. При определении изображения ϕ ' следует иметь в виду, что его оригиналом является производная ϕ′ . Воспользуемся теоремой дифференцирования для оригинала
ϕ′ sϕ − ϕ(0) и учтем тот факт, что всегда ϕ (0) = 0 . В результате имеем
ϕ′ sϕ и ϕ '= sϕ . Окончательным решением в изображениях поставленной задачи будет
|
= |
|
θ |
|
|
. |
(2.3) |
||
Ki |
|||||||||
|
|
|
|
|
|||||
|
|
sϕ |
|
||||||
2.1.4. Аппроксимация переходных характеристик
Оригиналы переходных характеристик содержат ряды, поэтому получить их изображения практически невозможно. В то же время основная составляющая переходной характеристики – S -функция – по характеру зависимости от τ близка к экспоненте (см. рис. 1.9).
Оказалось возможным представить S -функцию в переходной характе-
34
ристике ϕ (α, β, τ) = τ + S (α, β, τ) в виде комбинации экспонент:
S (α, β, τ) = S p (α, β) 2 |
(1 − e−τ τ01 ) − (1 − e−τ τ02 ) |
, |
(2.4) |
|
|
|
|
где S p (α, β) , как и прежде (см. 1.3.2), характеризует регулярный режим нагрева. Если τ01 = τ02 = τ00 , то
S (α, β, τ) = S p (α, β)(1 − e−kτ ), |
(2.5) |
|
где k = (τ00 )−1 . |
|
|
В горячем режиме нагрева сплошного цилиндра α3 = m = |
|
|
2R . |
||
Допустимость предложенной аппроксимации исследована для S (m,β, τ) при β = 1; 0,9; 0 . Установлено, что переходную характеристику для центра цилиндра (β = 0 ) следует определять по формуле (2.4), когда t01 ¹ t02 , а для поверхности (β = 1; 0,9 ) – по формуле (2.5), когда t01 = t02 . Необходимые данные для расчета этих характеристик приведены в табл. 2.1 (для β = 0 ) и 2.2 (для β = 1 ).
|
|
|
|
|
|
|
|
Таблица 2.1 |
|
|
|
|
|
|
|
|
|
|
|
m |
≤ 1,5 |
3 |
4 |
6 |
8 |
10 |
20 |
40 |
∞ |
Sp ×104 |
−837 |
−886 |
−959 |
−1100 |
−1164 |
−1196 |
−1237 |
−1247 |
−1250 |
|
|
|
|
|
|
|
|
|
|
t01 ×104 |
580 |
584 |
591 |
605 |
616 |
625 |
640 |
642 |
642 |
|
|
|
|
|
|
|
|
|
|
t02 ×104 |
444 |
436 |
427 |
418 |
418 |
422 |
429 |
430 |
430 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 2.2 |
|
|
|
|
|
|
|
|
|
|
|
m |
1,5 |
3 |
4 |
6 |
8 |
10 |
20 |
40 |
∞ |
Sp ×104 4 |
421 |
477 |
562 |
738 |
851 |
923 |
1080 |
1163 |
1250 |
|
|
|
|
|
|
|
|
|
|
t00 ×104 |
597 |
574 |
547 |
508 |
484 |
468 |
430 |
409 |
403 |
|
|
|
|
|
|
|
|
|
|
Таким образом, формулы (2.4), (2.5) и табл. 2.1 и 2.2 могут служить основой для определения закономерностей управления мощностью при различных режимах нагрева.
35
2.2. Ускоренный индукционный нагрев как вариант оптимального управления по быстродействию
2.2.1. Задача на минимум времени нагрева
без ограничений мощности
Необходимо определить режим нагрева, обеспечивающий за минимальное время достижение (при заданной температуре поверхности изделия) за-
Ki(0– ) → ∞ |
|
данного градиента температуры по сече- |
|||
|
нию нагреваемого тела. |
|
|||
Ki(τ) Ki(0+) |
|
Минимальное время нагрева до за- |
|||
|
данного перепада температуры между по- |
||||
|
|
|
|
верхностью и центром тела достигается |
|
|
|
|
|
при мгновенном установлении |
заданной |
|
|
|
температуры на поверхности θп |
и поддер- |
|
|
0 |
|
τ |
||
|
|
Рис. 2.4 |
|
жании ее до конца нагрева, т. е. до момен- |
|
|
|
|
та, когда перепад температуры между по- |
||
|
|
|
|
||
верхностью и центром уменьшится до заданного значения. В этом случае Ki (τ = 0− ) → ∞ , а Ki (τ = 0+ ) – значение ограниченное (рис. 2.4).
Используем полученное ранее решение |
|
= |
|
θ |
|
|
= |
|
|
θ |
|
и определим изоб- |
|||
Ki |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
sϕ |
|
ϕ |
' |
|||||||||
ражения θ и ϕ '. По условию поставленной задачи на поверхности (β = 1 )
оригинал θ (τ) = θп = const , а его изображение будет θп θп . s
Используем переходную характеристику при β = 1 , аппроксимированную экспонентой:
ϕ(τ) = τ + S (α, 1, τ) = τ + S p (α, 1)(1 − e−kτ )
ипродифференцируем ее:
ϕ′(τ) = 1 + S p (α, 1) ke−kτ .
Изображением, полученным с помощью (2.6), будет ϕ'= 1 + S p s
Теперь полученные изображения можно подставить в изображение закона управления мощностью:
(2.6)
k
s + k
.
(2.3) и найти
36
|
|
|
|
|
|
|
|
|
|
|
|
θп |
|
|
|
|
|
|
|
|
|
|
|
|
s + k |
|
|
|
|
|
|
|
s + k |
|
|
|
|
|
s + k |
|
|
||||||||||
Ki = |
|
|
|
|
|
|
|
|
|
|
|
|
|
= θп |
|
|
|
|
= |
|
|
|
|
|
= θп |
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k + (1 + kS p )s |
|
. |
|||||||||||||||||||||||||||||
|
|
1 |
+ S p |
k |
|
|
|
s + k + kS ps |
k + ϕ′(0) s |
||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
s |
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
s |
|
|
|
|
s + k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
′ |
|
|
|
|
|
+ kS p . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Здесь ϕ (0) = 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
Чтобы перейти от изображения к оригиналу, необходимо из дроби выде- |
|||||||||||||||||||||||||||||||||||||||||||||||||||
лить целую часть и разложить остаток на простые дроби: |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
|
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
s + k + |
|
|
|
|
|
|
|
− |
|
|
|
|
|
|
|
|
|
|
|
|
s |
+ |
|
|
|
+ k 1 |
− |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ϕ ' 0 |
ϕ ' 0 |
|
|
||||||||||||||||
|
|
|
|
|
|
|
θ |
|
|
|
|
ϕ |
|
'( |
0) |
|
ϕ '( |
0) |
|
|
θ |
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
п |
|
|
|
|
|
|
|
|
|
|
|
|
п |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
Ki = |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
|
|
|
|
|
|
|
|
( ) |
|
|
|
( ) |
= |
|
|||||||||||||||||||
|
|
ϕ '( 0) |
|
|
s + |
|
|
|
k |
|
|
|
|
|
|
|
|
ϕ '( 0) |
|
|
|
|
|
s + |
|
k |
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ϕ '( 0) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ϕ '( 0) |
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 − |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
θП |
|
|
|
ϕ ' |
0 |
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
|
|
|
|
1 + k |
|
|
|
|
|
( |
) |
|
|
. |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
s + |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ϕ '( 0) |
|
|
|
|
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ϕ '( |
0) |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Оригиналом этого изображения будет |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
θп |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Ki (τ) |
= |
|
|
|
|
|
δ0 |
(τ) + k 1 − |
|
|
|
e−kτ/ϕ′(0) . |
|
|
(2.7) |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
ϕ′(0) |
ϕ′(0) |
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Как видно из (2.7), искомый закон управления мощностью имеет импульсную составляющую бесконечной амплитуды, действующую нулевое время, благодаря которой мгновенно устанавливается заданная температура
на поверхности. Другая составляющая (см. рис. 2.4), |
|||||||
|
θ |
|
|
|
1 |
|
|
Ki(0+ + τ) = |
|
п |
k 1 |
− |
|
|
e−kτ ϕ′(0) , |
′ |
(0) |
′ |
|
||||
|
ϕ |
|
|
ϕ |
(0) |
||
убывающая по экспоненте, обеспечивает постоянство температуры на поверхности.
Для того чтобы найти время нагрева τн до заданного перепада темпера-
туры T между температурой поверхности Tп |
и температурой центра Tц , |
|||||
вводятся понятия средней температуры Tср = 0,5(Tп + Tц ) и отношения сред- |
||||||
ней температуры к температуре поверхности |
θср |
= |
Tср − Tнач |
, где T |
- ис- |
|
|
|
|||||
|
θп |
|
|
|
нач |
|
|
|
|
Tп − Tнач |
|
||
ходная (перед нагревом) температура, обычно 293 К. Закон роста средней
37
температуры был определен ранее как
|
τ |
|
|
|
|
|
|
|
θ |
п |
|
|
|
τ |
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
||||||
θ*ср = |
∫ Ki (τ)d τ = |
|
|
|
|
|
1 |
+ ∫ k 1 − |
|
|
|
e− kτ |
ϕ′(0)d τ |
= |
||||||||||||||||||||||
|
|
′ |
|
|
|
|
|
|
′ |
|
||||||||||||||||||||||||||
|
0 |
|
|
|
|
|
|
ϕ |
(0) |
|
0 |
|
|
|
|
|
|
ϕ (0) |
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
′ |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
ϕ (0) |
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
θ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
′ |
|
|
|
1 − e− kτ ϕ (0) |
|
|
||||||||||||
= |
|
|
п |
|
|
1 + k |
|
|
1 |
− |
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
ϕ′(0) |
|
k |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
ϕ′(0) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
θ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
′ |
(0) |
|
|
|
= |
|
п |
|
|
1 |
+ (ϕ′(0) −1) − (ϕ′(0) −1)e |
|
|
|
|
= |
|
||||||||||||||||||||||||
|
ϕ′(0) |
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
θ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− kτ |
|
′ |
|
|
|
|
|
|
|
|
|
= |
|
|
п |
|
ϕ′(0) − (ϕ′( |
0) −1)e |
|
|
|
|
|
= |
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
ϕ′(0) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(0) |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
′ |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
= θ |
|
1 − |
|
1 − |
|
|
|
|
|
|
e− kτ ϕ |
|
|
. |
|
|
|
(2.8) |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
п |
|
|
|
|
|
|
ϕ′( |
0) |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
В выражении (2.8) средняя температура θ еще не определена, известно |
||||||||||||||||||||||||||||||||||||
только, что θср |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ср |
|
|
|
|
|
|
|
|
|
θп ≤ 1. Чтобы найти время нагрева до заданного перепада |
|||||||||||||||||||||||||||||||||||
температуры между поверхностью и центром нагреваемого сплошного цилиндра τн , надо задать соответствующее перепаду конкретное значение θср и
решить уравнение θ |
|
|
= θ |
|
|
|
− 1 − |
|
|
|
1 |
|
e−kτн |
ϕ′(0) |
|
|
относительно τ |
|
|
||||||||||||||||||||||||
|
|
|
|
1 |
|
|
|
|
|
|
: |
||||||||||||||||||||||||||||||||
ср |
|
|
|
|
|
|
|
|
|
|
|
н |
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
п |
|
|
ϕ′(0) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
θср |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
θср |
|
|
||||
|
|
= 1 |
− 1 − |
|
|
|
|
|
|
e−kτ ϕ′(0) ; 1 − 1 − |
|
|
|
|
|
e−kτ ϕ′(0) = |
1 − |
|
; |
|
|||||||||||||||||||||||
|
θп |
|
′ |
(0) |
|
|
|
′ |
θп |
|
|||||||||||||||||||||||||||||||||
|
|
|
|
ϕ |
|
|
|
( )) |
|
|
|
|
ϕ (0) |
|
|
|
( |
|
ср |
|
|
|
п ) |
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
( |
|
|
|
|
|
|
|
|
|
|
|
1 − |
θ |
|
θ |
|
|
|
|
|||||||||||||
|
|
ekτн ϕ′(0) |
= |
1 − 1 ϕ′ 0 |
|
|
|
; |
e−kτн |
ϕ′(0) = |
|
|
|
|
|
|
|
; |
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
1 |
− ( |
θср θп ) |
1 − 1 ϕ′ |
( |
0 |
)) |
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
( |
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
′ |
(0) |
) |
|
|
|
|
|
′ |
(0) |
|
|
1 − |
|
|
|
′ |
(0) |
) |
|
|
|
|
|||||||||
|
|
|
kτн |
|
|
|
1 − 1 ϕ |
|
|
|
|
|
|
|
|
1 ϕ |
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
= ln |
|
|
|
( |
|
|
|
; |
τн = |
ϕ |
ln |
|
|
|
|
|
( |
|
|
|
|
|
|
|
. |
|
|
|
|||||||||||
|
|
|
ϕ′(0) |
|
|
|
|
|
|
|
|
|
|
|
1 − (θср θп ) |
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
1 − (θср θп ) |
|
|
|
|
|
|
k |
|
|
|
|
||||||||||||||||||||||||||
2.2.2. Управление мощностью с начальной ступенькой
Рассматрим режим управления мощностью более близкий к реальности по сравнению с предыдущим режимом. Управление нагревом в нем осуществляется в два этапа. На первом этапе мощность задается в виде ступеньки. Уровень ее определяется технико-экономическими соображениями. На втором этапе управление мощностью осуществляется по закону, обеспечивающему поддержание заданной температуры на поверхности.
38
Ставится задача определения длительности τ0 первого этапа нагрева, в течение которого при заданном уровне мощности Ki0 ступеньки температура поверхности достигнет заданного значения θп . На втором этапе устанавливается закон изменения мощности Ki (τ − τ0 ) , обеспечивающий поддержа-
ние достигнутой температуры, и определяется его длительность до момента |
|||||
установления заданного перепада темпера- |
Ki, θ |
|
|
||
туры между поверхностью и центром нагре- |
|
|
|||
|
|
|
|
|
|
ваемого тела (рис. 2.5). |
Ki0 |
|
|
||
|
Как и прежде, воспользуемся критери- |
|
|
||
|
θп |
|
|
||
альной формой |
уравнения теплопроводно- |
|
|
||
|
|
Ki(τ – τ0) |
|||
|
∂θ - Ñ2q = Ki (t) w(a, b) . При Ki (τ) = |
|
|
||
сти |
|
|
|
||
|
¶t |
|
|
τ0 |
τ |
= Ki0 = const и |
τ ≤ τ0 решением этого урав- |
|
|||
|
Рис. 2.5 |
|
|||
нения будет |
|
|
|
||
|
|
|
|
||
|
|
θ (β, τ) = Ki0ϕ (α, β, τ) . |
|
|
|
|
|
|
(2.9) |
|
|
|
При β = 1 к концу первого этапа ( τ = τ0 ) решение (2.9) имеет вид |
||||
|
|
θ (β = 1, τ0 ) = θп = Ki0ϕ1 (α, 1, τ0 ) . |
|
(2.10) |
|
Здесь ϕ1 (α, β, τ) – переходная характеристика первого этапа. В дальнейшем, для краткости она будет обозначаться ϕ1 (τ) .
Если же Ki (τ) – любая гладкая функция, содержащая в том числе ступеньку, то решение уравнения можно найти с помощью интеграла наложения
q(b, t) = ∫τ Ki (t)j¢(t - x)dx, который при τ > τ0 можно разделить на два инте-
0
грала:
|
|
τ0 |
|
′ |
|
|
|
|
|
τ |
|
|
′ |
|
θ (τ) = ∫ |
|
|
|
|
|
|
|
|
|
|
||||
Ki0ϕ1 (τ − ξ) d ξ + ∫ Ki (ξ) |
ϕ2 (τ − ξ) d ξ = |
|
||||||||||||
|
|
0 |
|
|
|
|
|
|
τ |
τ0 |
|
|
|
|
|
|
|
′ |
|
|
ξ=τ0 |
|
|
Ki ( |
′ |
(τ − ξ) d ξ = |
|
||
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|||||
= −Ki0ϕ1 |
(τ − ξ) |
|
ξ=0 |
+ ∫ |
ξ)ϕ2 |
|
||||||||
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
τ0 |
|
τ |
|
|
|
= Ki |
|
ϕ′ |
( |
τ) − ϕ |
(τ − τ |
|
) |
+ |
Ki (ξ)ϕ′ (τ − ξ) d ξ . |
|
||||
0 |
0 |
∫ |
(2.11) |
|||||||||||
|
1 |
|
1 |
|
|
|
|
|
|
2 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
τ0 |
|
|
|
|
|
|
|
|
|
|
|
39 |
|
|
|
|
|
|
Здесь ϕ′ (τ ) – переходная характеристика на втором этапе нагрева.
2
Введем новую координатную систему τ′ с началом координат на оси τ в точке τ = τ0 . Таким образом, τ′ = τ − τ0 , τ = τ0 + τ′ . Теперь при τ′ > 0 и β = 1 с учетом (2.10) получаем:
′ |
|
= Ki |
ϕ (τ |
) = Ki |
ϕ (τ |
|
′ |
|
+ |
τ |
′ |
′ |
|
|
|
|
|
||||||||||
θ(τ ) = θ |
0 |
+ τ ) |
∫ |
Ki (ξ)ϕ |
2 |
(τ − ξ)d ξ . (2.12) |
|||||||
|
п |
|
0 1 0 |
|
0 1 |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
Перевод (2.12) в изображения дает
θп = |
Ki0ϕ1 (τ0 ) |
= Ki |
|
ϕ |
(τ |
|
) − ϕ |
+ |
|
sϕ |
|
|
|
|
0 |
0 |
Ki |
2 |
. |
(2.13) |
|||||||||
|
||||||||||||||
s |
s |
1 |
|
1 |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
||||
Далее формула (2.13) приводится к виду, удобному для обратного преобразования Ki в оригинал:
|
|
|
|
|
|
|
|
|
|
ϕ (τ |
) |
+ ϕ1 − ϕ1 (τ0 ) |
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
1 0 |
|
= |
ϕ1 (τ0 ) + sϕ1 − sϕ1 (τ0 ) |
. |
|
|
|
||||||||||
|
|
|
|
|
|
= |
|
|
|
|
s |
|
|
|
|
|||||||||||||
|
|
|
|
|
Ki |
|
(2.14) |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
sϕ2 |
|
|
|
|
|
|
|
s2ϕ2 |
|
|
|
|
|
|
Результат преобразуется с целью получения оригинала Ki (τ′) . |
|
|
|
|||||||||||||||||||||||
|
|
Воспользуемся |
теоремой |
о дифференцировании изображения |
( |
|
)′ = |
|||||||||||||||||||||
|
|
f |
||||||||||||||||||||||||||
= |
|
|
|
− f (0) , |
|
|
|
= |
|
′ + f (0) |
и преобразуем переходные характеристики и |
|||||||||||||||||
f |
′ = sf |
sf |
|
f |
||||||||||||||||||||||||
их изображения. |
Напомним, что sϕ |
i |
= ϕ′ |
и ϕ |
i |
(0) = 0 , i = 1, 2. Преобразуем |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
||
слагаемые числителя и знаменателя (2.14) с целью понижения порядка относительно s :
sϕ |
(τ |
|
) = ϕ′ |
(τ |
|
) + ϕ |
(τ |
|
) ; sϕ = ϕ′ |
; s |
2 |
ϕ |
|
= s (sϕ |
|
′ |
. |
||
0 |
0 |
0 |
|
2 |
2 |
) = sϕ |
2 |
||||||||||||
1 |
|
1 |
|
1 |
|
1 |
1 |
|
|
|
|
|
|
||||||
После подстановки (2.15) в (2.14) получим:
|
|
|
′ |
′ |
(τ0 ) |
|
|
|
|
||||
Ki = Ki |
|
ϕ1 |
− ϕ1 |
. |
||
0 |
|
sϕ′ |
|
|||
|
|
|
|
|
||
|
|
|
|
2 |
|
|
(2.15)
(2.16)
Чтобы найти изображения производных переходных характеристик, используют их экспоненциальные аппроксимации. Оригиналом изображения
(2.16) будет
|
k ϕ (τ ) − τ |
|
−( |
k |
|
′ |
( |
0 |
)) |
τ′ |
|
|
|
|
|
1 1 0 |
0 |
2 |
|
−k τ′ |
|||||||||
Ki (τ′) = Ki0 |
|
|
Ae |
|
|
|
|
|
|
|
|
+ Be |
1 |
. |
′ |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ϕ2 (0) |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
40 |
|
|
|
|
|
|
|
|
|
|
|
|
Таким образом, установлен закон изменения мощности после ступеньки. Для определения длительности процесса нагрева до заданного перепада температуры между центром и поверхностью можно использовать закон роста средней температуры, который формулируется следующим образом:
τ′
θср = Ki0τ0 + ∫ Ki (τ′)d τ′ .
0
При τ′ = τ′н средняя температура достигает заданного значения. После
интегрирования получаем:
|
|
|
|
|
|
|
|
|
|
θср = Ki0τ0 + Ki0 × |
|
|
|
|||||||||
|
k |
ϕ |
(τ |
0 |
) − τ |
|
|
ϕ′ |
(0) |
|
|
−( |
k |
2 |
′ |
( |
0 |
′ |
|
|
1 |
|
× |
1 |
1 |
|
0 |
|
2 |
|
|
− e |
|
ϕ2 |
|
))τн |
+ |
|
|||||||
|
|
|
|
|
|
|
|
|
|
A 1 |
|
|
|
|
|
|
|
|
|
|||
|
|
′ |
|
|
|
k2 |
|
|
|
|
|
|
|
|
k1 |
|||||||
|
|
|
ϕ2 (0) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
B (1 − e−k1τ′н ) . (2.17)
Поскольку процесс адиабатический, при τ′н → ∞ средняя температура θср → θп . Действительно, если учесть, что в этом случае экспоненты устре-
мятся к единице и что Ki0ϕ1 (τ0 ) = θп , |
то после подстановки в правую часть |
||||||||||||||||||||||||||||||||||||||
(2.17) значений A и B и некоторых преобразований получаем: |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
k |
2 |
|
|
|
k1 |
ϕ1 (τ0 ) − τ0 |
ϕ′ (0) k − k |
2 |
|
|
|
|
1 k |
2 |
− k |
|
|
|||||||||||||||||||||
Ki0τ0 + Ki0 |
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
0 |
|
|
|
+ |
|
|
|
|
|
|
|
1 |
|
= |
||||||||||
′ |
(0) |
|
|
k2 |
|
|
|
|
k2 k0 − k1 |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
ϕ2 |
|
|
|
|
|
|
|
|
|
|
k1 k0 − k1 |
|
||||||||||||||||||||||||||
|
|
|
|
k0k1 ϕ1 (τ0 ) − τ0 |
|
1 |
|
k − k |
2 |
|
|
|
1 k |
2 |
− k |
|
|
|
|
|
|||||||||||||||||||
= Ki0τ0 + Ki0 |
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
+ |
|
|
|
|
|
|
|
|
1 |
|
= |
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
k |
2 |
|
|
|
|
k |
0 |
|
k − k |
k |
k |
0 |
− k |
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
1 |
|
|
|
1 |
|
|
|
|
|
1 |
|
|
|
|
|||||||
|
|
|
k0k1 |
ϕ1 (τ0 ) − τ0 |
|
k k |
− k k |
+ k |
|
k − k k |
|
|
|
|
|
||||||||||||||||||||||||
= Ki0τ0 + Ki0 |
|
|
|
|
|
|
|
|
|
|
|
1 0 |
1 2 |
|
|
|
2 0 |
|
|
|
1 0 |
|
= |
|
|||||||||||||||
|
|
|
|
k2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k1k0 |
(k0 − k1 ) |
|
|
|
|
|
|
|
|
|
||||||||||||||
= Ki0τ0 + Ki0ϕ1 (τ0 ) − Ki0τ0 = θп .
Аналитического решения трансцендентного уравнения (2.17) относительно τ′н не существует. Можно найти приближенное решение в аналитическом виде, приняв B = 0, т. е. A = 1 и
k = k |
|
= Ki |
|
τ |
|
+ Ki |
ϕ |
|
(τ |
|
|
) |
− Ki |
|
|
τ |
|
− Ki |
|
ϕ |
(τ |
|
) |
− τ |
|
e |
− |
k |
ϕ′ |
( |
0) τ′ |
|
|
|||||||||||||||||
2 |
0 |
0 |
|
0 |
0 |
0 |
0 |
0 |
0 |
|
( |
|
2 2 |
|
) н = |
|
||||||||||||||||||||||||||||||||||
1 |
|
|
|
|
|
|
|
|
0 1 |
|
|
|
|
|
|
|
|
|
|
|
1 |
(0) |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
= θ − Ki |
|
|
ϕ (τ |
|
|
) − τ |
|
e |
− |
k |
ϕ′ |
τ′ |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
0 |
0 |
0 |
|
( |
|
2 2 |
|
|
) |
н ; |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
п |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
θп − θср |
|
|
|
|
|
|
|
|
|
|
(0) |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− |
k |
ϕ′ |
(0) |
τ′ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− |
k |
ϕ′ |
τ′ |
||||||||||||
θ − θ = Ki |
|
ϕ |
|
(τ |
|
) − τ |
|
e |
|
( |
|
2 2 |
|
|
|
) |
н ; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= e |
( |
|
2 2 |
|
) н ; |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ki |
|
ϕ |
(τ |
|
) − τ |
|
|
|
|
|
|
|||||||||||||||||||||||||
п |
ср |
|
|
0 |
|
1 |
|
0 |
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
0 |
0 |
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
41 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
