Quais os métodos mais adequados para calcular a rigidez de materiais que não obedecem a Lei de Hooke?

Passo 1

O que o enunciado pede: identificar os métodos mais adequados para calcular a rigidez (ou stiffness) de materiais que não obedecem à Lei de Hooke.
Informações dadas (implícitas pela pergunta): não temos um caso específico — o material pode ser não linear elástico, viscoelástico, plástico, hiperelástico, dependente de taxa, composto, etc. A resposta deve cobrir métodos práticos (experimentais) e teóricos/computacionais (modelos e técnicas de cálculo).

Passo 2

O que significa “rigidez” em regime não linear — conceitos básicos.
  • Em materiais lineares a rigidez é simplesmente a constante E de Hooke (módulo elástico): \sigma = E\,\varepsilon.
  • Para materiais não lineares precisamos distinguir duas noções úteis: o módulo secante e o módulo tangente.
    • Módulo secante: é a razão entre tensão e deformação em um ponto, ou seja \displaystyle E_s=\frac{\sigma}{\varepsilon}. Ele dá uma “rigidez média” entre a origem e o ponto considerado.
    • Módulo tangente: é a derivada local da curva tensão-deformação, ou seja \displaystyle E_t=\frac{d\sigma}{d\varepsilon}. Ele é a rigidez instantaneous que aparece ao linearizar o comportamento num pequeno incremento.
      Explicação intuitiva: pense na curva \sigma(\varepsilon). O secante é a inclinação da linha do zero até o ponto; o tangente é a inclinação da reta que toca a curva naquele ponto. Na prova, quando pedirem “rigidez” de material não linear quase sempre querem o módulo tangente (para análises incrementais) ou explicitar se usam o secante (em modelos simplificados).

Passo 3

Métodos experimentais simples para obter rigidez não linear (o que medir e como extrair).
  • Teste monotônico (tração/compresión) registrando par (\varepsilon,\sigma). A curva experimental fornece E_s e E_t.
    • Para E_s: calcule \displaystyle E_s(\varepsilon)=\frac{\sigma(\varepsilon)}{\varepsilon}.
    • Para E_t: calcule numericamente a derivada; por exemplo usando diferença finita centrada:
    E_t(\varepsilon)\approx\frac{\sigma(\varepsilon+\Delta\varepsilon)-\sigma(\varepsilon-\Delta\varepsilon)}{2\Delta\varepsilon}.Use \Delta\varepsilon pequeno, compatível com ruído experimental.
  • Ensaios dinâmicos (DMA — Dynamic Mechanical Analysis) para materiais dependentes de frequência: mede módulo complexo \displaystyle E^*(\omega)=E'(\omega)+\mathrm{i}E''(\omega). Aqui E'(\omega) é o módulo de armazenamento (rigidez efetiva), E''(\omega) é a perda (amortecimento).
  • Nanoindentação para regiões locais (obtém E efetivo por análise do carregamento/descarga).

Passo 4

Modelagem constitutiva: classes de modelos e como deles extrair rigidez. (aula prática sobre cada uma)
  • Modelos hiperelásticos (materiais elásticos grandes deformações, ex.: borrachas). São definidos por uma energia de deformação W(\mathbf{F}) onde \mathbf{F} é o gradiente de deformação. A tensão é obtida por derivada: \mathbf{S}=\dfrac{\partial W}{\partial \mathbf{E}} (forma de tensões depende da medida escolhida). A rigidez tangente consiste na segunda derivada da energia: \mathbb{C}=\dfrac{\partial^2 W}{\partial \mathbf{E}\partial \mathbf{E}}. Exemplo simples (unidimensional, para didática): modelo de Lei de potência/elastômero simples \displaystyle \sigma=K\varepsilon^n. Então
E_t(\varepsilon)=\frac{d\sigma}{d\varepsilon}=K\,n\,\varepsilon^{\,n-1}, \qquad E_s(\varepsilon)=K\,\varepsilon^{\,n-1}.
Mostra como calcular direto a partir da lei constitutiva.
  • Modelos viscoelásticos (dependência temporal): descrevem comportamento com operadores diferenciais ou integrais de convolução. Forma canônica em tempo:
\sigma(t)=\int_0^t G(t-\tau)\,\dot\varepsilon(\tau)\,d\tau,
onde G(t) é o módulo relaxação. No domínio de frequência a rigidez é o módulo complexo E^*(\omega). Para uso incremental numérico, usa-se séries de Prony ou modelos de Maxwell/Kelvin–Voigt; a rigidez instantânea é G(0^+) e o tangente depende do tempo e da taxa de carregamento. Para provas: aprenda a transformar um modelo Prony em operadores diferenciais e extrair a resposta na rampa/degrau.
  • Modelos elastoplásticos: combinação elástica + critério de escoamento + regra de fluxo + endurecimento. Procedimento típico para obter rigidez incremental: usar linearização algorítmica (return-mapping) e derivar o operador tangente consistente (o chamado “algorithmic tangent modulus” C_\text{alg}). Em uniaxial, o módulo elastoplástico tangente pode ser calculado como
E_t=\frac{E}{1+\dfrac{E}{H_\text{tan}}},
onde E é a rigidez elástica e H_\text{tan} é um parâmetro de endurecimento (isso é um caso típico com endurecimento linear isotrópico; se H_\text{tan}\to\infty recupera elástico, se H_\text{tan}=0 dá plasticidade perfeita com E_t\to0 na região plástica). Na prova, o procedimento é: (i) calcular tensão elástica “trial”, (ii) verificar condição de yield, (iii) projetar de volta ao domínio admissível (return mapping), (iv) derivar a atualização para obter E_t.
  • Modelos constitutivos empíricos (curve-fitting): ajustar função matemática a curva \sigma(\varepsilon) (ex.: polinômios, funções exponenciais, lawas de tipo Ramberg–Osgood) e então derivar para obter E_t. Exemplo Ramberg–Osgood (simplificado):
\varepsilon=\frac{\sigma}{E}+K\left(\frac{\sigma}{\sigma_0}\right)^n.
Dado isso, pode-se implicitamente derivar \,d\sigma/d\varepsilon para achar E_t.

Passo 5

Cálculo prático do módulo tangente a partir de uma lei dada — exemplo passo a passo (didático).
Exemplo: suponha lei uniaxial \displaystyle \sigma(\varepsilon)=K\,\varepsilon^n. Queremos E_t e E_s em \varepsilon=\varepsilon_0.
  • Passo A: escrever a função (já dada).
  • Passo B: derivar para obter o tangente:
E_t(\varepsilon)=\frac{d\sigma}{d\varepsilon}=K\,n\,\varepsilon^{\,n-1}.
  • Passo C: avaliar em \varepsilon_0:
E_t(\varepsilon_0)=K\,n\,\varepsilon_0^{\,n-1}.
  • Passo D: secante em \varepsilon_0:
E_s(\varepsilon_0)=\frac{\sigma(\varepsilon_0)}{\varepsilon_0}=K\,\varepsilon_0^{\,n-1}.
Interpretação: E_t=n\,E_s para esse modelo. Isso mostra como, numa prova, você pode transformar rapidamente uma lei dada em rigidez instantânea.

Passo 6

Como proceder numericamente (FEM e análises incrementais) — o que o examinador geralmente espera.
  • Análise incremental–iterativa: para materiais não lineares faz-se pequenos incrementos de carga/deformação e resolve-se de forma iterativa (Newton–Raphson). A matriz de rigidez tangente global é
\mathbf{K}_\text{t}=\int_\Omega \mathbf{B}^\top\,\mathbb{C}_t\,\mathbf{B}\,dV,
onde \mathbf{B} é a matriz de forma e \mathbb{C}_t é o tensor de rigidez tangente (derivada da tensão em relação à deformação). Logo, para conseguir convergência eficiente você precisa do “consistent tangent” (linearização exata da atualização constitutiva). Em prova: explique o conceito e, se pedirem, escreva a expressão integral acima e diga como obter \mathbb{C}_t derivando a lei constitutiva.
  • Se não houver tempo para derivar o consistente, pode usar aproximação numérica (diferença finita) para \mathbb{C}_t, mas isso encarece e pode prejudicar convergência.

Passo 7

Métodos para materiais dependentes de escala ou heterogêneos (compósitos): homogenização.
  • Regras de mistura simples: Voigt (deformação uniforme) e Reuss (tensão uniforme) para obter módulos efetivos E_\text{Voigt},E_\text{Reuss}. Por exemplo:
E_\text{Voigt}=\sum_i v_i E_i,\qquad \frac{1}{E_\text{Reuss}}=\sum_i \frac{v_i}{E_i},
onde v_i são frações volúmicas.
  • Homogenização mais avançada: método auto-consistente, Mori–Tanaka, ou FE ^2 (solve problemas de microestrutura por elementos finitos e extrair tensor efetivo). Esses métodos fornecem uma rigidez efetiva não-linear que pode depender de \varepsilon, taxa, etc.

Passo 8

Viscoelasticidade e análise no domínio da frequência — como obter “rigidez” útil.
  • Se o material é dependente de tempo/velocidade, rigidez instantânea e de longo-prazo são diferentes: E_\text{inst}=G(0^+), E_\infty=G(\infty) (no caso de módulo de relaxação).
  • Em regime harmônico, medir E'(\omega) (módulo armazenamento) e usar isso como rigidez efetiva para análises dinâmicas. Fórmula (definição):
E^*(\omega)=E'(\omega)+\mathrm{i}E''(\omega)=\frac{\sigma_0}{\varepsilon_0}\mathrm{e}^{\mathrm{i}\varphi},
onde \varphi é defasagem entre tensão e deformação. Para provas, explique que rigidez para vibrações é E'(\omega).

Passo 9

Receita prática, passo-a-passo, que você pode seguir numa prova quando pedirem “calcular a rigidez” de um material não linear:
  1. Identifique o regime de comportamento (não linear elástico, viscoelástico, plastificante, grande deformação, composto, dependente de taxa).
  2. Se tiver uma lei constitutiva analítica, derive explicitamente E_t=\dfrac{d\sigma}{d\varepsilon}. Faça simplificações caso haja simetria/unidimensionalidade. (Ex.: veja Passo 5.)
  3. Se tiver só curva experimental, calcule E_s=\dfrac{\sigma}{\varepsilon} e E_t por diferença finita. Use suavização se houver ruído.
  4. Se for uso em FEM/iteração, derive o operador tangente \mathbb{C}_t (segunda derivada de energia ou linearização da lei constitutiva) e monte \mathbf{K}_t=\int \mathbf{B}^\top\mathbb{C}_t\mathbf{B}dV. Explique no enunciado que precisa do consistent tangent.
  5. Se o material for dependente de tempo, diga qual módulo usar (instantâneo, longo-prazo, ou módulo de armazenamento E'(\omega)) dependendo do problema.
  6. Se o material for heterogêneo, escolha método de homogenização apropriado (Voigt/Reuss para estimativas simples; Mori–Tanaka ou FE ^2 para maior precisão).

Passo 10

Dicas rápidas e armadilhas do exame (empatia e prática):
  • Não confunda módulo secante com tangente; sempre declare qual está usando.
  • Em problemas de pequenas variações num FEM, o examinador provavelmente quer o módulo tangente (linearização local). Se não tiver certeza, escreva ambos e justifique.
  • Se o enunciado disser “rigidez efetiva para vibração”, lembre-se de falar de E'(\omega) ou usar ensaio dinâmico.
  • Quando o enunciado der uma lei em termos de energia (hiperelástico), calcule tensões por derivada e rigidez por segunda derivada (mostre isso passo a passo).
  • Se faltarem dados experimentais, proponha um método experimental simples (ensaio monotônico ou DMA) e escreva como extrair E_t e E_s.

Resposta

Módulo tangente (calcular E_t=\dfrac{d\sigma}{d\varepsilon}) obtido por derivação da lei constitutiva.
Módulo secante (calcular E_s=\dfrac{\sigma}{\varepsilon}) extraído da curva tensão–deformação.
Modelagem hiperelástica (energia W\Rightarrow tensão pela derivada e rigidez pelo tensor \mathbb{C}=\partial^2 W/\partial\mathbf{E}^2).
Modelos viscoelásticos (Prony/Maxwell/Kelvin–Voigt) e uso de módulo complexo E^*(\omega)=E'(\omega)+\mathrm{i}E''(\omega).
Elastoplástico com return-mapping e cálculo do algoritmo de rigidez tangente consistente (algorithmic tangent).
Homogenização para materiais heterogêneos (Voigt/Reuss, Mori–Tanaka, FE ^2).
Ensaios experimentais (tração compressão para E_s,E_t; DMA para E'(\omega); nanoindentação para regiões locais).

Exercícios relacionados