Exercícios Resolvidos de Modelos térmicos

Ver Teoria

Enunciado

Considere o reservatório aquecido mostrado na figura abaixo. Um fluxo de processo está sendo aquecido no reservatório por um aquecedor elétrico. A proporção da transferência calorífica, q ( t ) para o fluido de processo está relacionado ao sinal, m ( t ) , por:

q t = a m ( t )

Você pode supor que o reservatório de aquecimento está bem isolado, que o fluido está bem misturado no reservatório e que a capacidade calorífica e a densidade do fluido são constantes. Desenvolva o modelo matemático que descreva como a temperatura de entrada, T i ( t ) , o fluxo do processo, f ( t ) , e o sinal, m ( t ) , afetam a temperatura de saída T ( t ) . Depois determine as funções de transferência e esboce o diagrma de blocos.

Passo 1

Nesse exercício o enunciado foi bonzinho e já fez todas as considerações necessárias. Então vamos direto pro balanço de energia porque a gente tá sem tempo, irmão.

A c ú m u l o   d e   e n e r g i a = E n t r a - S a i + G e r a d a - { C o n s u m i d a }

Aqui, além das correntes de entrada e saída do reservatório, também teremos uma taxa de geração de energia promovido pelo aquecedor elétrico.

A c ú m u l o   d e   e n e r g i a = E n t r a - S a i + { G e r a d o }

Vamos escrever as expressões para cada um dos termos acima:

A c ú m u l o   d e   e n e r g i a = d U ( t ) d t = m c p d T ( t ) d t = ρ V c p d T ( t ) d t

E n t r a = c p m ˙ t T i t = c p ρ f ( t ) T i t

s a i = c p m ˙ ( t ) T t = c p ρ f ( t ) T ( t )

E o enunciado nos deu a expressã da energia do gerador elétrico:

G e r a d o = a m ( t )

Juntando tudo:

ρ V c p d T ( t ) d t = c p ρ f ( t ) T i t - c p ρ f ( t ) T ( t ) - a m ( t )

ρ V c p d T ( t ) d t = c p ρ f t T i t - f t T t - a m ( t )

Passo 2

E olha só. De novo caímos em uma EDO não linear. Isso acontece porque temos dois termos onde temos duas funções de tempo se multiplicando. Precisamos resolver isso.

Para fazer essa linearização, usaremos a linearização de Taylor. Quando temos um caso como esse de produto de funções, a linearização é feita da seguinte forma:

f x = h x . g ( x )

f x ≅ f x - + ∂ f ∂ h x - . h - h - 1 ! + ∂ f ∂ g x - . g - g - 1 !

Primeiro a gente escreve a função em torno do ponto onde estamos linearizando, depois nós vamos derivar a função duas vezes, uma vez em relação a cada função que compõe a função principal.

Sendo assim:

f t T i t = f - T i - + f - T i t - T i - + T i - [ f t - f - ]

E analogamente:

f t T t = f - T - + f - T t - T - + T - [ f t - f - ]

Substituindo no balanço

ρ V c p d T ( t ) d t = c p ρ f - T i - + f - T i t - T i - + T i - f t - f - - f - T - - f - T t - T - - T - f t - f - + a m t               ( 1 )

Passo 3

O próximo passo é colocar a nossa expressão em variáveis desvio. Para isso, vamos escrever a expressão ( 1 ) com todos os termos em seus estados iniciais.

ρ V c p d T - d t = c p ρ f - T i - - f - T - + a m -               ( 2 )

Olha só, a maioria dos termos se anulou porque as diferenças ficaram igual a 0 . Fazendo 1 - ( 2 ) :

ρ V c p d [ T t - T - ] d t = c p ρ f - T i t - T i - + T i - f t - f - - f - T t - T - - T - f t - f - + a m t - m -             ( 3 )

Quando fazemos 1 - ( 2 ) os termos que são puramente constantes se anulam e sobram apenas os que são função do tempo.

Vamos definir nossas variáveis desvio:

Γ t = T t - T -

Γ i t = T i t - T i -

F t = f t - f -

M t = m t - m -

Substituindo em ( 3 ) :

ρ V c p d Γ ( t ) d t = c p ρ f - Γ i t + T i - F t - f - Γ t - T - F t + a M ( t )

Rearranjando:

ρ V c p d Γ ( t ) d t = c p ρ f - Γ i t - f - Γ t + T - i - T - F t + a M ( t )

Passo 4

Agora, vamos passar tudo que depende de Γ ( t ) para a esquerda:

ρ V c p d Γ ( t ) d t + c p ρ f - Γ t = c p ρ f - Γ i t + c p ρ T - i - T - F t + a M ( t )

Show! Agora já podemos aplicar Laplace:

ρ V c p s Γ s + c p ρ f - Γ s = c p ρ f - Γ i s + c p ρ T - i - T - F s + a M ( s )

Colocando Γ ( s ) em evidência:

( ρ V c p s + c p ρ f - ) Γ s = c p ρ f - Γ i s + c p ρ T - i - T - F s + a M ( s )

Por fim:

Γ s = c p ρ f - ρ V c p s + c p ρ f - Γ i s + c p ρ T - i - T - ρ V c p s + c p ρ f - F s + a ρ V c p s + c p ρ f - M ( s )

Passo 5

Agora só nos resta esboçar o diagrama de blocos.

Como a saída Γ ( s ) depende de três entradas que estão se somando:

Resposta

ρ V c p d T ( t ) d t = c p ρ f t T i t - f t T t - a m ( t )

Γ s = c p ρ f - ρ V c p s + c p ρ f - Γ i s + c p ρ T - i - T - ρ V c p s + c p ρ f - F s + a ρ V c p s + c p ρ f - M ( s )