Ir para o conteúdo

Integração numérica

Os vetores de forças internas do elemento \(\boldsymbol{q}^e, \boldsymbol{Q}^e\) e a matriz de rigidez do elemento \(\boldsymbol{K}^e\) obtidos em Discretização do trabalho virtual interno e Trabalho virtual externo e montagem da equação global assumem a forma de integrais sobre o domínio do elemento \(\Omega^e\) ou \(\Omega^e_0\). O FrontISTR avalia essas integrais numericamente usando quadratura de Gauss.

Quadratura de Gauss e mudança de variáveis

A quadratura de Gauss aproxima uma integral sobre o domínio de referência \(\Xi\) por uma combinação linear dos valores do integrando nos pontos de integração \(\boldsymbol{\xi}_i \in \Xi\) e dos pesos \(w_i\). Sua aplicação ao domínio do elemento \(\Omega^e\) envolve uma mudança de variáveis por meio do mapeamento \(\boldsymbol{x}: \Xi \to \Omega^e\), resultando em

\[ \int_{\Omega^e} f(\boldsymbol{x})\, dv \approx \sum_{i=1}^{n_q} w_i\, f(\boldsymbol{x}(\boldsymbol{\xi}_i))\, J_{\xi_i}, \qquad J_{\xi_i} = \left.\det\!\left(\frac{\partial \boldsymbol{x}}{\partial \boldsymbol{\xi}}\right)\right|_{\boldsymbol{\xi}_i} \]

onde \(n_q\) é o número de pontos de integração e \(J_{\xi_i}\) é o determinante do Jacobiano da transformação. O domínio de referência \(\Xi\) é definido para cada tipo de elemento (para um hexaedro, \([-1,1]^3\); para triângulos, tetraedros e prismas triangulares, as respectivas formas de referência), e os pontos de integração \(\boldsymbol{\xi}_i\) e os pesos \(w_i\) são fornecidos em tabelas numéricas. Integrais de superfície são tratadas da mesma forma, mapeando-se uma face do elemento a partir de um domínio de referência bidimensional.

A seguir são apresentados números representativos de pontos de integração usados pelo FrontISTR (para a correspondência com os tipos de elemento, consulte Sistema de numeração de elementos e biblioteca de funções de forma).

Tipo de elemento Regra de integração Número de pontos de integração
Tetraedro de 4 nós (tet4n) Regra de 1 ponto 1
Tetraedro de 10 nós (tet10n) Regra de 4 pontos 4
Prisma triangular de 6 nós (prism6n) Regra de 2 pontos 2
Prisma triangular de 15 nós (prism15n) Regra de 9 pontos 9
Hexaedro de 8 nós (hex8n) Gauss-Legendre 2×2×2 8
Hexaedro de 20 nós (hex20n) Gauss-Legendre 3×3×3 27
Quadrilátero de 4 nós (quad4n) Gauss-Legendre 2×2 4
Quadrilátero de 8 nós (quad8n) Gauss-Legendre 3×3 9
Triângulo de 3 nós (tri3n) Regra de 1 ponto 1
Triângulo de 6 nós (tri6n) Regra de 3 pontos 3

Para elementos hexaédricos, quadrilaterais e de linha, são usados produtos tensoriais da regra de Gauss-Legendre em cada direção coordenada. Triângulos, tetraedros e prismas triangulares usam regras específicas adequadas a geometrias do tipo simplex (disposições de pontos que integram polinômios exatamente sobre o próprio triângulo).

Aplicação à integração do elemento

Na integração do elemento, as coordenadas naturais \(\boldsymbol{r}\) são usadas como coordenadas do domínio de referência (\(\boldsymbol{r} = \boldsymbol{\xi}\)), e o mapeamento para as coordenadas físicas é dado pela interpolação das coordenadas nodais usando as funções de forma. Dependendo da escolha da configuração de referência (Estrutura da análise incremental), usam-se as formulações a seguir.

Método de Lagrange Total (integração sobre a configuração de referência \(\Omega^e_0\)): o mapeamento e o Jacobiano são

\[ \boldsymbol{X}^e(\boldsymbol{r}) = \sum_{\alpha} N_\alpha^e(\boldsymbol{r})\, \boldsymbol{X}^e_\alpha, \qquad J_{r_i} = \left.\det\!\left(\frac{\partial \boldsymbol{X}}{\partial \boldsymbol{r}}\right)\right|_{\boldsymbol{r}_i} \]

e a força interna do elemento e a matriz de rigidez são aproximadas por

\[ \boldsymbol{Q}^e \approx \sum_{i=1}^{n_q} w_i\, (\boldsymbol{B}_L + \boldsymbol{B}_{NL})^T\, \boldsymbol{S}\, J_{r_i}, \qquad \boldsymbol{K}^e \approx \sum_{i=1}^{n_q} w_i\, \boldsymbol{K}^e_{x}\, J_{r_i} \]

Todos os termos \(\boldsymbol{B}_L, \boldsymbol{B}_{NL}, \boldsymbol{S}, \boldsymbol{K}^e_{x}\) são avaliados no ponto de integração \(\boldsymbol{r}_i\).

Método de Lagrange Atualizado (integração sobre a configuração atual \(\Omega^e\)): o mapeamento e o Jacobiano são

\[ \boldsymbol{x}^e(\boldsymbol{r}) = \sum_{\alpha} N_\alpha^e(\boldsymbol{r})\, \boldsymbol{x}^e_\alpha = \sum_{\alpha} N_\alpha^e(\boldsymbol{r})\, (\boldsymbol{X}^e_\alpha + \boldsymbol{u}^e_\alpha), \qquad J_{r_i} = \left.\det\!\left(\frac{\partial \boldsymbol{x}}{\partial \boldsymbol{r}}\right)\right|_{\boldsymbol{r}_i} \]

e

\[ \boldsymbol{q}^e \approx \sum_{i=1}^{n_q} w_i\, \boldsymbol{B}_L^T\, \boldsymbol{\sigma}\, J_{r_i}, \qquad \boldsymbol{K}^e \approx \sum_{i=1}^{n_q} w_i\, \boldsymbol{K}^e_{x}\, J_{r_i} \]

são usados como aproximações.

A única diferença entre as duas formulações é se as coordenadas nodais fornecidas ao mapeamento são \(\boldsymbol{X}^e_\alpha\) ou \(\boldsymbol{x}^e_\alpha\); os pontos de integração, os pesos e a estrutura do loop sobre os pontos de integração são comuns.

Integração completa e integração reduzida

A integração que utiliza pontos de integração suficientes para integrar exatamente o grau polinomial do integrando é chamada de integração completa, enquanto a integração que utiliza um nível a menos de pontos de integração é chamada de integração reduzida. A integração reduzida é usada para aliviar o locking de cisalhamento e volumétrico, mas exige tratamento de modos de deformação espúrios, como modos hourglass. O número de pontos de integração de cada tipo de elemento e a escolha entre integração completa e reduzida são tratados em Sistema de numeração de elementos e biblioteca de funções de forma e nas seções subsequentes, bem como em Formulações avançadas de elementos.

Tópicos relacionados