Enunciado

Pilares de concreto reforçados com aço são usados em grandes construções. Falhas estruturais podem ocorrer em altas temperaturas em função de um incêndio por causa do enfraquecimento do núcleo metálico. Considere um pilar composto com 200   m m de espessura constituído por um núcleo central de aço (espessura de 50   m m ) entre duas paredes de concreto com 75   m m de espessura. O pilar encontra-se em uma temperatura inicial uniforme de Ti = 27°C e é subitamente exposto a produtos de combustão a T∞ = 900°C e h   =   40   W / ( m 2   ·   K ) nas duas superfícies expostas. A temperatura da vizinhança também é de 900 ° C .

(a) Usando um método implícito de diferenças finitas, com Δ x   =   10   m m e Δ t   =   100   s , determine as temperaturas na superfície exposta do concreto e no centro da placa de aço em t   =   10.000   s . As propriedades do aço são: k a ç o   =   55   W / ( m   ·   K ) ,   ρ a ç o   =   7850   k g / m 3 ,   e   c a ç o   =   450   J / ( k g   ·   K ) . As propriedades do concreto são: k c o n   =   1,4   W / ( m   ·   K ) ,   ρ c o n   =   2300   k g / m 3 ,   c c o n   =   880   J / ( k g   ·   K ) ,   e   ε   =   0,90 . Represente graficamente as temperaturas máxima e mínima do concreto juntamente com as temperaturas máxima e mínima do aço ao longo do p e r í o d o   0   ≤   t   ≤   10.000   s .

(b) Repita a parte (a), mas agora leve em conta uma resistência térmica de contato de R t ,   c ' ' =   0,20   m 2 · K / W na interface concreto-aço.

(c) Em t = 10.000 s, o fogo é extinto e as temperaturas da vizinhança e do ambiente retornam a T ∞   =   T v i z   =   27 ° C . Usando o mesmo coeficiente de transferência de calor e a mesma efetividade das partes (a) e (b), determine a temperatura máxima do aço e o tempo crítico, no qual ocorre a temperatura máxima no aço para os casos com e sem a resistência térmica de contato. Represente graficamente a temperatura da superfície do concreto, a temperatura do concreto adjacente ao aço e as temperaturas do aço no período 10.000   ≤   t   ≤   20.000   s .

Passo 1

Vamos lá! Nosso problema pede os gráficos de temperatura, mas antes não vamos esquecer do nosso velho esboço pra facilitar nosso problema.

Agora precisamos fazer nossas famosas considerações. Hehe (u.u)

I - C o n d u ç ã o   u n i d i m e n s i o n a l

I I - P r o p r i e d a d e s   c o n s t a n t e s

Passo 2

(a)

Vamos lá! Primeiro que tal fazermos o balanço de energia para todos os nós e seus respectivos volumes de controle:

1 :

ρ c o n c c o n ⋅ Δ x 2 ⋅ T 1 p + 1 - T 1 p Δ t = h + h r T ∞ - T 1 p + 1 + k c o n T 2 p + 1 - T 1 p + 1 Δ x

Onde h r = ε σ T s - T 1 p + 1 T s 2 - T 1 p + 1 2 e T s = T ∞ :

Para os nós 2 - 7 :

ρ c o n c c o n Δ x ⋅ T n p + 1 - T n p Δ t = k c o n T n + 1 p + 1 - T n p + 1 Δ x + k c o n T n - 1 p + 1 - T n p + 1 Δ x

Para o nó 8 :

ρ c o n c c o n Δ x ⋅ T 8 p + 1 - T 8 p Δ t = k c o n T 9 p + 1 - T 8 p + 1 R t o t ' ' + k c o n T 7 p + 1 - T 8 p + 1 Δ x

Onde R t o t ' ' = Δ x 2 k c o n + Δ x 2 k a ç o + R t , c ' ' , com R t ,   c ' ' desprezível:

Para o nó 9 :

ρ a ç o c a ç o Δ x ⋅ T 9 p + 1 - T 9 p Δ t = k a ç o T 8 p + 1 - T 9 p + 1 R t o t ' ' + k a ç o T 10 p + 1 - T 9 p + 1 Δ x

Para o nó 10 :

ρ a ç o c a ç o Δ x ⋅ T 10 p + 1 - T 10 p Δ t = k a ç o T 9 p + 1 - T 10 p + 1 Δ x + k a ç o T 1 1 p + 1 - T 10 p + 1 Δ x

Para o nó 10 :

ρ a ç o c a ç o ⋅ Δ x 2 ⋅ T 11 p + 1 - T 11 p Δ t = k a ç o T 10 p + 1 - T 11 p + 1 Δ x

Passo 3

Agora ele pede o histórico de temperatura para o núcleo de metal e superfície de concreto, para o nó 1 1 e 1 , normalmente podemos resolver isso em forma de tabela, mas devido à natureza da questão fica melhor fazer com gráfico.

Assim usando o método das diferenças finitas temos.

P.S: Pode desenvolver em qualquer software de sua preferência.

Devido à alta condutividade térmica do metal a temperatura é quase uniforme, entretanto para o concreto devido à baixa condutividade térmica temos uma discrepância entre a temperatura de superfície do concreto T m a x e a temperatura interna T m i n .

Assim temos que as temperaturas são após 10.000   s :

T 1 = 882 ° C ,   T 11 = 594   ° C

Passo 4

(b) Adicionando a resistência térmica de contato agora para o nosso problema temos que nosso gráfico fica assim:

A explicação segue a mesma lógica da letra (a), a discrepância entre a temperatura de superfície e a temperatura interna do concreto, e a temperatura praticamente uniforme do aço.

Entretanto a resistência térmica de contato reduz bastante a quantidade de calor que atravessa o concreto até o aço, assim a temperatura do aço cai drasticamente.

Assim temos que as temperaturas são após 10.000   s :

T 1 = 88 9 ° C ,   T 11 = 231   ° C

Passo 5

(c)

Repetindo o passo 1 entretanto agora as temperaturas nodais para T 10.000 s , com T ∞ = 300 K , considerando o tempo inicial os gráficos vão ficar assim:

S e m   r e s i s t ê n c i a   d e   c o n t a t o

C o m   r e s i s t ê n c i a   d e   c o n t a t o

Agora temos que:

Sem resistência de contato

T a ç o s e m ,   R ' ' t , c = 617,5   ° C

t c r i t s e m ,   R ' ' t , c = 1.100   s

Com resistência de contato

T a ç o c o m ,   R ' ' t , c = 305,7   ° C

t c r i t c o m ,   R ' ' t , c = 6.000   s

Isso é tudo pessoal! \o/

Resposta

(a) T 1 = 882 ° C , T 11 = 594   ° C

(b)   T 1 = 88 9 ° C , T 11 = 231   ° C

(c)

Sem resistência de contato

T a ç o s e m ,   R ' ' t , c = 617,5   ° C

t c r i t s e m ,   R ' ' t , c = 1.100   s

Com resistência de contato

T a ç o c o m ,   R ' ' t , c = 305,7   ° C

t c r i t c o m ,   R ' ' t , c = 6 . 00 0   s

Ver Outros Exercícios desse livro

Exercícios de Livros Relacionados

Formule, por diferenças finitas, sistemas de equações cuja solução aproxime a solução dos seguintes problemas de contorno: y ' ' = y s e n ( y ) + t y y ( 0 ) = 1 y ( 1 ) = 5 y ' ' = 2 y + y 3 - t y (
Considere a rede para um sistema bidimensional sem geração volumétrica interna com temperaturas nodais mostradas abaixo. Se o espaçamento da grade for de 100m m e a condutividade térmica do material f
Considere a ferramenta de corte do Problema 4.40, mas com um ângulo de ponta de 60 ° . Derive a equação de diferença finita nodal para o nó m ,n .
Em livros-texto sobre CFD [5, 23-27], frequentemente se substitui os lados esquerdos da equação (8.119b e c) pelas duas expressões seguintes, respectivamente: ∂ ∂ x u 2 + ∂ ∂ y ( v u ) ∂ ∂ x u v + ∂ ∂
A equação de Laplace em coordenadas polares planas, equação 8.11, é complicada pelo raio variável. Considere a malha de diferenças finitas na figura abaixo, com nós (i,j) igualmente espaçados de Δ θ e