Ir para o conteúdo

Tensões e leis de balanço

Este capítulo organiza as definições e as relações entre tensores de tensão, juntamente com as equações de equilíbrio e a simetria das tensões derivadas da conservação da massa, da quantidade de movimento linear e do momento angular. Os três tensores de tensão — tensão de Cauchy \(\boldsymbol{\sigma}\), primeira tensão de Piola-Kirchhoff \(\boldsymbol{P}\) e segunda tensão de Piola-Kirchhoff \(\boldsymbol{S}\) — são introduzidos conforme a configuração usada para a superfície de referência e o vetor de força. Consulte Movimento, deformação e medida de deformação para as definições de configuração, movimento e gradiente de deformação, e Lista de símbolos de grandezas físicas para a lista de símbolos.

Tensor de tensão de Cauchy

Considere uma superfície infinitesimal no ponto \(\boldsymbol{x}\) da configuração atual (área \(d\Gamma\), normal unitária exterior \(\boldsymbol{n}\)) sobre a qual atua a força \(d\boldsymbol{f}\), e defina o vetor de tração \(\boldsymbol{t}\) como a força por unidade de área, \(d\boldsymbol{f} = \boldsymbol{t}\, d\Gamma\). \(\boldsymbol{t}\) depende tanto da posição \(\boldsymbol{x}\) quanto da normal \(\boldsymbol{n}\). A aplicação da conservação da quantidade de movimento a um tetraedro infinitesimal mostra que \(\boldsymbol{t}\) é linear em \(\boldsymbol{n}\); portanto, existe um tensor de segunda ordem \(\boldsymbol{\sigma}(\boldsymbol{x}, t)\) tal que

\[ \boldsymbol{t} = \boldsymbol{\sigma} \boldsymbol{n}, \qquad t_i = \sigma_{ij} n_j \]

para qualquer \(\boldsymbol{n}\) (teorema da tensão de Cauchy). Esse tensor \(\boldsymbol{\sigma}\) é chamado de tensor de tensão de Cauchy e também de tensão verdadeira, pois representa a força por unidade de área medida em relação à geometria da configuração atual. O componente \(\sigma_{ij}\) representa o componente \(x_i\) da força por unidade de área que atua sobre uma superfície infinitesimal normal ao eixo coordenado \(x_j\) na configuração atual.

Como mostrado adiante a partir da conservação do momento angular, a tensão de Cauchy é um tensor simétrico:

\[ \boldsymbol{\sigma} = \boldsymbol{\sigma}^T, \qquad \sigma_{ij} = \sigma_{ji} \]

Em três dimensões, portanto, possui seis componentes independentes (\(\sigma_{11}, \sigma_{22}, \sigma_{33}, \sigma_{12}, \sigma_{23}, \sigma_{31}\)). A vetorização da tensão na notação de Voigt pressupõe essa simetria (consulte Notação tensorial e fundamentos matemáticos). A tensão de Cauchy é usada na saída de tensões da análise geometricamente não linear do FrontISTR com o método Updated Lagrange e na análise de pequenas deformações.

Primeira e segunda tensões de Piola-Kirchhoff

Em deformações finitas, muitas vezes é conveniente expressar a tensão usando superfícies e vetores da configuração de referência; por isso são introduzidos dois tensores de tensão de Piola-Kirchhoff. Seja \(d\Gamma_0\) a área de uma superfície infinitesimal na configuração de referência e \(\boldsymbol{N}\) sua normal unitária exterior.

A primeira tensão de Piola-Kirchhoff (primeira tensão PK, tensão nominal) \(\boldsymbol{P}\) é definida como o tensor de tensão obtido quando a força \(d\boldsymbol{f}\) na configuração atual atua sobre uma superfície infinitesimal na configuração de referência:

\[ d\boldsymbol{f} = \boldsymbol{P} \boldsymbol{N}\, d\Gamma_0 \]

Usando a fórmula de Nanson \(\boldsymbol{n}\, d\Gamma = J \boldsymbol{F}^{-T} \boldsymbol{N}\, d\Gamma_0\) juntamente com o teorema da tensão de Cauchy, obtém-se a relação com a tensão de Cauchy

\[ \boldsymbol{P} = J\, \boldsymbol{\sigma} \boldsymbol{F}^{-T}, \qquad \boldsymbol{\sigma} = \frac{1}{J} \boldsymbol{P} \boldsymbol{F}^T \]

A primeira tensão PK é, em geral, não simétrica.

A segunda tensão de Piola-Kirchhoff (segunda tensão PK) \(\boldsymbol{S}\) é definida como o tensor de tensão obtido ao trazer a força da configuração atual \(d\boldsymbol{f}\) de volta à configuração de referência por \(\boldsymbol{F}^{-1}\) e aplicá-la a uma superfície infinitesimal na configuração de referência:

\[ \boldsymbol{F}^{-1}\, d\boldsymbol{f} = \boldsymbol{S} \boldsymbol{N}\, d\Gamma_0 \]

Tanto o vetor de força quanto a superfície sobre a qual ele atua são expressos por grandezas da configuração de referência, tornando-a invariante sob rotação de corpo rígido e simétrica. Sua relação com a primeira tensão PK e a transformação a partir da tensão de Cauchy são

\[ \boldsymbol{P} = \boldsymbol{F} \boldsymbol{S}, \qquad \boldsymbol{\sigma} = \frac{1}{J} \boldsymbol{F} \boldsymbol{S} \boldsymbol{F}^T, \qquad \boldsymbol{S} = J\, \boldsymbol{F}^{-1} \boldsymbol{\sigma} \boldsymbol{F}^{-T} \]

A segunda tensão PK é derivada de uma função de energia de deformação hiperelástica \(W(\boldsymbol{C})\) como \(\boldsymbol{S} = 2\,\partial W / \partial \boldsymbol{C}\) e, no método Total Lagrange, é usada em par de trabalho conjugado com a deformação de Green-Lagrange \(\boldsymbol{E}\), isto é, \((\boldsymbol{S}, \boldsymbol{E})\).

As configurações de referência e as propriedades de simetria são resumidas abaixo. No limite de pequenas deformações (\(\boldsymbol{F} \to \boldsymbol{I}\), \(J \to 1\)), as três tensões coincidem.

Tensor de tensão Superfície de referência Vetor de força Simetria Aplicação
Tensão de Cauchy \(\boldsymbol{\sigma}\) Configuração atual \(d\Gamma, \boldsymbol{n}\) Configuração atual \(d\boldsymbol{f}\) Simétrica Método Updated Lagrange / análise de pequenas deformações
Primeira tensão PK \(\boldsymbol{P}\) Configuração de referência \(d\Gamma_0, \boldsymbol{N}\) Configuração atual \(d\boldsymbol{f}\) Geralmente não simétrica Equações de equilíbrio na configuração de referência
Segunda tensão PK \(\boldsymbol{S}\) Configuração de referência \(d\Gamma_0, \boldsymbol{N}\) Configuração de referência \(\boldsymbol{F}^{-1} d\boldsymbol{f}\) Simétrica Método Total Lagrange / hiperelasticidade

Conservação da massa e da quantidade de movimento e equações de equilíbrio

Seja \(\rho\) a densidade de massa na configuração atual e \(\rho_0\) a densidade de massa na configuração de referência. A lei de conservação da massa \(\int_{\Omega} \rho\, dv = \int_{\Omega_0} \rho_0\, dV\), juntamente com a transformação do elemento de volume \(dv = J\, dV\), reduz-se à forma local

\[ \rho_0 = J \rho \]

Esta é a relação local de conservação da massa.

Sejam \(\boldsymbol{g}\) a força de corpo por unidade de massa e \(\boldsymbol{a}\) a aceleração. Aplicando o teorema da tensão de Cauchy e o teorema da divergência de Gauss à conservação da quantidade de movimento linear (primeira lei de movimento de Euler), obtém-se a equação de equilíbrio (equação do movimento) local na configuração atual:

\[ \nabla_x \cdot \boldsymbol{\sigma} + \rho \boldsymbol{g} = \rho \boldsymbol{a}, \qquad \frac{\partial \sigma_{ij}}{\partial x_j} + \rho g_i = \rho a_i \]

Reescrevendo a integral na configuração de referência usando a fórmula de Nanson e \(\boldsymbol{P} = J \boldsymbol{\sigma} \boldsymbol{F}^{-T}\), obtém-se a forma local na configuração de referência:

\[ \nabla_X \cdot \boldsymbol{P} + \rho_0 \boldsymbol{g} = \rho_0 \boldsymbol{a}, \qquad \frac{\partial P_{ij}}{\partial X_j} + \rho_0 g_i = \rho_0 a_i \]

As duas formas são equivalentes pela regra de transformação de tensões e por \(\rho_0 = J\rho\). Desprezando o termo de inércia, obtém-se

\[ \nabla_x \cdot \boldsymbol{\sigma} + \rho \boldsymbol{g} = \boldsymbol{0}, \qquad \nabla_X \cdot \boldsymbol{P} + \rho_0 \boldsymbol{g} = \boldsymbol{0} \]

que é a equação de equilíbrio estático e constitui o ponto de partida da análise estática do FrontISTR (linear e não linear).

Conservação do momento angular e simetria das tensões

Combinando a conservação do momento angular (segunda lei de movimento de Euler) com a equação de equilíbrio obtida da conservação da quantidade de movimento linear, obtém-se a simetria do tensor de tensão de Cauchy:

\[ \boldsymbol{\sigma} = \boldsymbol{\sigma}^T, \qquad \sigma_{ij} = \sigma_{ji} \]

Assim, ele possui seis componentes independentes em três dimensões. Tomando a transposta de ambos os lados da transformação \(\boldsymbol{\sigma} = J^{-1} \boldsymbol{F} \boldsymbol{S} \boldsymbol{F}^T\) e usando a não singularidade de \(\boldsymbol{F}\), conclui-se que a segunda tensão PK também é simétrica:

\[ \boldsymbol{S} = \boldsymbol{S}^T \]

e, portanto, também possui seis componentes independentes. Em contraste, a primeira tensão PK \(\boldsymbol{P} = \boldsymbol{F} \boldsymbol{S}\) geralmente não é simétrica; vale apenas a relação \(\boldsymbol{P} \boldsymbol{F}^T = \boldsymbol{F} \boldsymbol{P}^T\) (isto é, \(\boldsymbol{P} \boldsymbol{F}^T\) é simétrico), e ela possui nove componentes independentes.

Essas simetrias constituem a base para expressar a tensão como um vetor de seis componentes na notação de Voigt. Consulte Notação tensorial e fundamentos matemáticos e Elasticidade linear para as convenções da notação de Voigt e a construção das matrizes de material.

Tópicos relacionados