Ir para o conteúdo

Elasticidade linear

Este capítulo aborda a lei constitutiva (lei de Hooke) para materiais elásticos lineares. Para detalhes sobre a seleção do modelo e as especificações de entrada, consulte a secção de funções 03_material.

Fenómeno (resposta elástica linear)

Um corpo elástico linear é um material que, num intervalo em que a deformação sob carga é suficientemente pequena, apresenta uma relação linear entre tensão e deformação e uma resposta reversível e independente do histórico, na qual toda a deformação regressa ao estado original após a descarga. Ao contrário da plasticidade ou da fluência, não possui variáveis internas que representem o estado (como deformação plástica ou deformação viscosa), e a tensão é determinada unicamente pela deformação atual.

No intervalo de pequenas deformações, utilizam-se a tensão de Cauchy \(\boldsymbol{\sigma}\) e a deformação infinitesimal \(\boldsymbol{\varepsilon}\); no intervalo de deformação finita (problemas que incluem grandes rotações, mas em que a própria deformação é pequena), utilizam-se a segunda tensão de Piola-Kirchhoff \(\boldsymbol{S}\) e a deformação de Green-Lagrange \(\boldsymbol{E}\) (material de St.Venant-Kirchhoff).

Lei constitutiva

Material elástico linear isotrópico

No intervalo de pequenas deformações, a lei de Hooke isotrópica é expressa utilizando as constantes de Lamé \(\lambda, \mu\) como

\[ \boldsymbol{\sigma} = \lambda\, \mathrm{tr}(\boldsymbol{\varepsilon})\, \boldsymbol{I} + 2\mu\, \boldsymbol{\varepsilon} \]

Na forma de componentes,

\[ C_{ijkl} = \lambda\, \delta_{ij}\delta_{kl} + \mu\, \bigl( \delta_{ik}\delta_{jl} + \delta_{il}\delta_{jk} \bigr) \]

e a dupla contração com a deformação dá \(\sigma_{ij} = C_{ijkl}\,\varepsilon_{kl}\). As constantes de Lamé relacionam-se com o módulo de Young \(E\) e o coeficiente de Poisson \(\nu\) da seguinte forma:

\[ \lambda = \frac{E\nu}{(1+\nu)(1-2\nu)}, \qquad \mu = \frac{E}{2(1+\nu)}. \]

Na notação de Voigt, a relação entre o vetor de deformações \(\hat{\varepsilon}\) e o vetor de tensões \(\hat{\sigma}\) é escrita como \(\hat{\sigma} = D\, \hat{\varepsilon}\), e a matriz de material tridimensional \(D\) é

\[ D = \begin{bmatrix} \lambda + 2\mu & \lambda & \lambda & 0 & 0 & 0 \\ \lambda & \lambda + 2\mu & \lambda & 0 & 0 & 0 \\ \lambda & \lambda & \lambda + 2\mu & 0 & 0 & 0 \\ 0 & 0 & 0 & \mu & 0 & 0 \\ 0 & 0 & 0 & 0 & \mu & 0 \\ 0 & 0 & 0 & 0 & 0 & \mu \end{bmatrix} \]

Esta representação matricial é utilizada na implementação.

Material elástico linear ortotrópico

Para materiais cujas constantes elásticas diferem consoante a direção dos eixos principais, utiliza-se um modelo ortotrópico com nove constantes independentes: os módulos de Young nas direções dos eixos principais \(E_1, E_2, E_3\), os coeficientes de Poisson entre os eixos principais \(\nu_{12}, \nu_{23}, \nu_{31}\) e os módulos de corte \(G_{12}, G_{23}, G_{31}\). Para a forma explícita da matriz \(D\), consulte referências normalizadas.

Intervalo de deformação finita: material de St.Venant-Kirchhoff

Para problemas que incluem grandes rotações, mas em que a própria deformação é pequena, utiliza-se o material de St.Venant-Kirchhoff como lei elástica linear que relaciona a segunda tensão de Piola-Kirchhoff \(\boldsymbol{S}\) e a deformação de Green-Lagrange \(\boldsymbol{E}\):

\[ \boldsymbol{S} = \boldsymbol{C} : \boldsymbol{E}, \qquad C_{ijkl} = \lambda\, \delta_{ij}\delta_{kl} + \mu\, \bigl( \delta_{ik}\delta_{jl} + \delta_{il}\delta_{jk} \bigr). \]

As definições das constantes de Lamé \(\lambda, \mu\) são as mesmas que no caso de pequenas deformações. Note-se, contudo, que, como as definições do par tensão-deformação são diferentes, esta é uma lei constitutiva diferente da lei de Hooke para pequenas deformações.

Material hipoelástico

No método de Lagrange atualizado, utiliza-se um material hipoelástico que adota a lei elástica linear entre a taxa de Jaumann do tensor de tensão relativa de Kirchhoff \(\hat{\boldsymbol{\sigma}}^{\nabla J}\) e o tensor taxa de deformação \(\boldsymbol{D}\):

\[ \hat{\boldsymbol{\sigma}}^{\nabla J} = \boldsymbol{C} : \boldsymbol{D} \]

Também aqui, a definição de \(\boldsymbol{C}\) em termos das constantes de Lamé é a mesma que acima, mas, como o par tensão-deformação é diferente, trata-se de uma lei constitutiva distinta. Para detalhes sobre a forma como \(\boldsymbol{\sigma}_{t_{n+1}}\) é atualizado através do procedimento de integração temporal (integração de Euler progressiva e aproximação por diferenças centrais), consulte tmptexdocs/11b_continuum_mechanics.tex.

Redução para análise bidimensional e elementos de casca

Estado plano de tensão (\(\sigma_{33} = 0\))

Para estado plano de tensão, que representa o carregamento de placas finas e cascas, assume-se \(\sigma_{33} = \sigma_{13} = \sigma_{23} = 0\) e utiliza-se \(D\) com o componente de deformação fora do plano \(\varepsilon_{33}\) eliminado algebricamente:

\[ D = \frac{E}{1 - \nu^2} \begin{bmatrix} 1 & \nu & 0 \\ \nu & 1 & 0 \\ 0 & 0 & \dfrac{1-\nu}{2} \end{bmatrix}. \]

Estado plano de deformação (\(\varepsilon_{33} = 0\))

Para estado plano de deformação, que representa a secção transversal de uma estrutura longa restringida na direção da espessura, assume-se \(\varepsilon_{33} = \varepsilon_{13} = \varepsilon_{23} = 0\):

\[ D = \frac{E}{(1+\nu)(1-2\nu)} \begin{bmatrix} 1-\nu & \nu & 0 \\ \nu & 1-\nu & 0 \\ 0 & 0 & \dfrac{1-2\nu}{2} \end{bmatrix}. \]

Axissimétrico

Na análise axissimétrica, que trata problemas uniformes na direção \(\theta\) em coordenadas cilíndricas \((r, \theta, z)\), são mantidos os quatro componentes \(\boldsymbol{\sigma} = (\sigma_{rr}, \sigma_{\theta\theta}, \sigma_{zz}, \sigma_{rz})^T\) e \(\boldsymbol{\varepsilon} = (\varepsilon_{rr}, \varepsilon_{\theta\theta}, \varepsilon_{zz}, 2\varepsilon_{rz})^T\), obtendo-se

\[ D = \frac{E}{(1+\nu)(1-2\nu)} \begin{bmatrix} 1-\nu & \nu & \nu & 0 \\ \nu & 1-\nu & \nu & 0 \\ \nu & \nu & 1-\nu & 0 \\ 0 & 0 & 0 & \dfrac{1-2\nu}{2} \end{bmatrix} \]

.

Elementos de casca (estado plano de tensão + corte transversal)

Para elementos de casca, a tensão no plano é tratada como estado plano de tensão e, para os componentes de deformação de corte transversal \((2\varepsilon_{13}, 2\varepsilon_{23})\), utiliza-se \(\kappa\, G\) com o fator de correção de corte \(\kappa\) (geralmente \(\kappa = 5/6\)). Para a formulação detalhada de cada elemento, consulte Formulações avançadas de elementos. A especificação na secção de funções 03_material de que «os elementos de casca suportam apenas elasticidade linear» baseia-se nesta simplificação de tratar o estado no plano como estado plano de tensão.

Dependência da temperatura

O módulo de Young \(E\), o coeficiente de Poisson \(\nu\) e o coeficiente de dilatação linear \(\alpha\) são geralmente especificados como funções da temperatura \(T\). No FrontISTR, introduzem-se valores para uma sequência finita de pontos de temperatura \(T_1 < T_2 < \cdots < T_n\), e a temperatura \(T\) durante o cálculo é avaliada por interpolação linear entre pontos adjacentes. Para temperaturas fora do intervalo (\(T < T_1\) ou \(T > T_n\)), utiliza-se para extrapolação o valor no extremo mais próximo, \(T_1\) ou \(T_n\) (extrapolação constante).

Para a dependência da temperatura das propriedades térmicas e de condução de calor (densidade, calor específico e condutividade térmica), consulte Propriedades térmicas. As regras de interpolação são comuns a estas propriedades.

Tópicos relacionados