Ir para o conteúdo

Análise transiente de condução de calor

Esta seção apresenta a discretização temporal e o método iterativo de solução para a análise de condução de calor em sólidos pelo método dos elementos finitos (MEF). Para as equações governantes e as condições de contorno do meio contínuo, consulte Equação de condução de calor.

Equação discretizada (ponto de partida)

Ao discretizar a equação de condução de calor (equação de condução de calor (gov_he_main)) pelo método de Galerkin, obtém-se

\[\begin{equation} K T + M \frac{\partial T}{\partial t} = F \label{eq:2.4.8} \end{equation}\]

onde

\[\begin{equation} K = \int\left( k_x \frac{\partial N^T}{\partial x}\frac{\partial N}{\partial x} + k_y \frac{\partial N^T}{\partial y}\frac{\partial N}{\partial y} + k_z \frac{\partial N^T}{\partial z}\frac{\partial N}{\partial z} \right) dV + \int hc N^T N ds + \int hr N^T N ds \label{eq:2.4.9} \end{equation}\]
\[\begin{equation} M = \int \rho c N^T N dV \label{eq:2.4.10} \end{equation}\]
\[\begin{equation} F = \int Q N^T dV - \int q_s N^T dS + \int{hc} T c N^T dS + \int{hcTr} ({T+Tr}) ({T^2 + T r^2}) N^T dS \label{eq:2.4.11} \end{equation}\]
\[\begin{equation} N = (N^1, N^2, \ldots, Ni) \label{eq:2.4.12} \end{equation}\]

Aqui, \(K\), \(M\), \(F\) e \(N\) são, respectivamente, a matriz de condução de calor (incluindo as contribuições de convecção e radiação na fronteira), a matriz de massa, o vetor de carga térmica e a matriz de funções de forma. As definições dos símbolos das propriedades do material (\(\rho\), \(c\), \(k_x, k_y, k_z\), \(Q\), \(hc\), \(hr\) etc.) seguem Equação de condução de calor.

Discretização temporal e solução iterativa

A equação \(\eqref{eq:2.4.8}\) é uma equação não linear e transiente. Usando o método de Euler implícito para a discretização temporal, quando a temperatura no instante \(t=t_0\) é conhecida, a temperatura no instante \(t=t_0+\Delta t\) é calculada usando a seguinte equação.

\[\begin{equation} K_{t=t_0+\Delta t} T_{t=t_0+\Delta t} + M_{t=t_0+\Delta t} \frac{T_{t=t_0+\Delta t} - T_{t=t_0}}{\Delta t} = F_{t=t_0+\Delta t} \label{eq:2.4.13} \end{equation}\]

Considere melhorar o vetor de temperaturas \(T_{t=t_0+\Delta t}^{(i)}\) que satisfaz aproximadamente a Eq.\(\eqref{eq:2.4.13}\), a fim de obter a solução mais precisa \(T_{t=t_0+\Delta t}^{(i)+1}\).

Para isso, primeiro expresse o vetor de temperaturas da seguinte forma.

\[\begin{equation} T_{t=t_0+\Delta t}= T_{t=t_0+\Delta t}^{(i)} + \Delta T_{t=t_0+\Delta t}^{(i)} \label{eq:2.4.14} \end{equation}\]

O produto da matriz de condução de calor pelo vetor de temperaturas, a matriz de massa e os termos relacionados são representados aproximadamente da seguinte forma.

\[\begin{equation} K_{t=t_0+\Delta t} T_{t=t_0+\Delta t} = K_{t=t_0+\Delta t}^{(i)} T_{t=t_0+\Delta t}^{(i)} + \frac{\partial \big(K_{t=t_0+\Delta t}^{(i)} T_{t=t_0+\Delta t}^{(i)}\big) } {\partial T_{t=t_0+\Delta t}^{(i)} } \Delta T_{t=t_0+\Delta t}^{(i)} \label{eq:2.4.15} \end{equation}\]
\[\begin{equation} M_{t=t_0+\Delta t} = M_{t=t_0+\Delta t}^{(i)} + \frac{\partial M_{t=t_0+\Delta t}^{(i)}}{\partial T_{t=t_0+\Delta t}^{(i)}} \Delta T_{t=t_0+\Delta t}^{(i)} \label{eq:2.4.16} \end{equation}\]

Substituindo a Eq.\(\eqref{eq:2.4.14}\), a Eq.\(\eqref{eq:2.4.15}\) e a Eq.\(\eqref{eq:2.4.16}\) na Eq.\(\eqref{eq:2.4.13}\) e omitindo os termos de segunda ordem e superiores, obtém-se a seguinte equação.

\[\begin{equation} \bigg(\frac{M_{t=t_0+\Delta t}^{(i)}}{\Delta t} + \frac {\partial M_{t=t_0+\Delta t}^{(i)} } { \partial T_{t=t_0+\Delta t}^{(i)} } \frac{T_{t=t_0+\Delta t}^{(i)} - T_{t=t_0}}{\Delta t} + \frac{\partial \big(K_{t=t_0+\Delta t}^{(i)} T_{t=t_0+\Delta t}^{(i)}\big)} {\partial T_{t=t_0+\Delta t}^{(i)}} \bigg) \Delta T_{t=t_0+\Delta t}^{(i)} \\\ = F_{t=t_0+\Delta t} - M_{t=t_0+\Delta t}^{(i)} \frac{T_{t=t_0+\Delta t}^{(i)} - T_{t=t_0}}{\Delta t} - K_{t=t_0+\Delta t}^{(i)} T_{t=t_0+\Delta t}^{(i)} \label{eq:2.4.17} \end{equation}\]

A matriz de coeficientes do lado esquerdo é então avaliada aproximadamente usando a seguinte equação.

\[\begin{equation} K^{(i)} = \frac{M_{t=t_0+\Delta t}^{(i)}}{\Delta t} + \frac{\partial \big( K_{t=t_0+\Delta t}^{(i)} T_{t=t_0+\Delta t}^{(i)} \big)}{\partial T^{(i)}_{t=t_0+\Delta t}} = \frac{M_{t=t_0+\Delta t}^{(i)}}{\Delta t} + K_{T_{t=t_0+\Delta t}}^{(i)} \label{eq:2.4.18} \end{equation}\]

Aqui, \(K_{T_{t=t_0+\Delta t}}^{(i)}\) é a matriz de rigidez tangente.

Por fim, a temperatura no instante \(t=t_0+\Delta t\) pode ser calculada executando cálculos iterativos com a seguinte equação.

\[\begin{equation} K^{(i)} \Delta T_{t=t_0+\Delta t}^{(i)} = F_{t=t_0+\Delta t} - M_{t=t_0+\Delta t}^{(i)} \frac{T_{t=t_0+\Delta t}^{(i)} - T_{t=t_0}}{\Delta t} - K^{(i)} T_{t=t_0+\Delta t}^{(i)} \label{eq:2.4.19} \end{equation}\]

Em particular, na análise em regime permanente, os cálculos iterativos são realizados usando a seguinte equação.

\[ K_T^{(i)} \Delta T_{t=\infty}^{(i)} = F_{t=\infty} - K_T^{(i)} \Delta T_{t=\infty}^{(i)} \]
\[\begin{equation} T_{t=\infty}^{(i+1)} = T_{t=\infty}^{(i)} + \Delta{T}_{t=\infty}^{(i)} \label{eq:2.4.20} \end{equation}\]

Na análise transiente, como é usado um método implícito para a discretização temporal, a escolha do incremento de tempo \(\Delta t\) geralmente não está sujeita a uma restrição quanto à sua magnitude. Entretanto, se o incremento de tempo \(\Delta t\) for grande demais, aumentará o número de iterações necessárias para a convergência. Em geral, um incremento de tempo \(\Delta t\) excessivamente grande aumenta o número de iterações. Na implementação, o controle automático do incremento monitora a magnitude do vetor residual, reduzindo \(\Delta t\) quando a convergência é lenta e aumentando \(\Delta t\) quando o número de iterações é pequeno (→ consulte Controle de passo para obter detalhes).

Tópicos relacionados