Redes de molas auxéticas

O que o artigo mostra

Muitos materiais, de redes de metamateriais a membranas celulares, são modelados como redes de molas de Hooke. Grossman e Boudaoud deduzem o contínuo elástico equivalente usando apenas a geometria e a topologia da rede e calculam os deslocamentos não afins : os desvios locais de um alongamento uniforme que tornam as redes desordenadas mais flexíveis que a estimativa de deformação uniforme. Comparam a teoria com simulações de uma faixa tracionada para três famílias: redes triangulares ordenadas, uma colmeia que se torna reentrante e uma «espuma» formada ao deslocar aleatoriamente os nós de uma rede triangular, que se torna auxética (coeficiente de Poisson negativo: alarga ao ser tracionada) com o aumento da desordem.

Experimente: rede perfeita, η = 0 · η = 0.3 · η = 0.5, auxética

O modelo que resolvemos

Molas e unem nós; cada uma tem rigidez k_e, comprimento de repouso ℓ̄_e e comprimento atual ℓ_e, e a energia é E = ½ Σ_e k_e (ℓ_e − ℓ̄_e)² (Eq. 3 do artigo). Uma rede de molas é uma treliça articulada, então o painel direito a mostra no visualizador do Stabileo: cada barra colorida pela força axial, sobre a forma deformada. Calculamos duas coisas.

• O ensaio de faixa do artigo. Uma faixa de 13 × 52 (cerca de 800 nós), com as fileiras superior e inferior presas e a superior deslocada conforme a deformação imposta; o coeficiente de Poisson é obtido da mudança de largura entre os dois nós de borda mais próximos da meia altura, como na Seç. V A do artigo. Resolvida como treliça linear pela fatoração esparsa de Cholesky do motor do Stabileo; o mesmo modelo resolvido pela análise de treliças 2D do motor de produção concorda até 10⁻¹⁴.
• Módulos efetivos exatos. Uma célula periódica sob deformação macroscópica imposta, com todos os nós livres para relaxar: a rigidez restante é o tensor elástico C da rede, incluindo o amolecimento não afim, sem efeitos de borda. O tensor de deformação uniforme (afim) é a mesma soma sem relaxamento.
• Pré-tensão (um acréscimo, ausente nos cálculos do artigo): comprimentos de repouso encurtados por uma fração, a rede relaxada até o equilíbrio pelo método de Newton e a resposta linear tomada nesse estado, incluindo o enrijecimento por tração.

Comparação com o artigo

Espuma: a transição auxética

O artigo relata ν decrescendo com a desordem η, cruzando zero em η ≃ 0.46 e chegando a −0.1 quando η → 0.5. Nossos módulos periódicos exatos (10 amostras de 26 × 26 células para cada η) cruzam zero em η = 0.446 e dão ν = −0.148 em η = 0.5. É uma comparação entre protocolos: nossa média da célula periódica fica dentro de todas as barras de erro das simulações de faixas presas de 13 × 52 da Fig. 5 para η ≤ 0.45; isso não valida nosso protocolo de faixa. O marcador do artigo em η = 0.49 é uma seta abaixo do gráfico (cerca de −0.14); nosso valor ali é −0.117.

Nossa implementação do protocolo de faixa (tracejada) dá 0.371 para a rede perfeita, onde o valor exato é 1/3 e a Fig. 5 do artigo mostra 0.331. O excesso é um efeito de borda livre de ordem 1/largura: faixas de larguras 26 e 52 dão 0.354 e 0.345, extrapolando para 0.336. A faixa do artigo não o mostra, e o texto não detalha suficientemente a medição para sabermos por quê.

Colmeia: Eq. (23)

Uma colmeia de arestas unitárias cujos vértices opostos estão separados por 2ℓ é regular em ℓ = 1 e reentrante abaixo de ℓ = ½. O artigo triangula cada hexágono com um nó central e molas de rigidez 10⁻³ em relação às arestas, e compara com a fórmula ν(ℓ) = −(2 − 4ℓ)(ℓ + ½)/(3 + 4ℓ − 4ℓ²). Nossa rede periódica concorda até 3.4 × 10⁻³ com essas molas e até 3.4 × 10⁻⁶ com molas de rigidez relativa 10⁻⁶: a diferença é a própria rigidez das molas flexíveis. Mostrar ℓ = 0.25 (reentrante) · ℓ = 0.9

Redes ordenadas: uma observação sobre as constantes elásticas

Para redes triangulares ordenadas, o artigo fornece uma fórmula, a Eq. suplementar (9), para ν em função da forma da célula (φ, ψ) e da orientação θ, mostrada na Fig. 4. Uma rede de Bravais não tem deslocamentos não afins, portanto a resposta exata é a de deformação uniforme, que calculamos independentemente. A Eq. suplementar (9) e a linha de cálculo da Fig. 7 (módulo de Young da rede (3, 2), lido dos dados vetoriais do PDF e multiplicado pela área de um triângulo) coincidem com esse resultado afim para constantes elásticas k_e ∝ 1/ℓ̄_e², até 10⁻¹⁴ e 6 × 10⁻⁴. Com constantes elásticas iguais obtemos ν = 1/6 em (φ, ψ, θ) = (1, 0.5, 0), onde a Eq. (9) dá 8/27. Isso pode refletir uma convenção que não identificamos; está na nossa lista de perguntas aos autores. Mostrar (1, 0.5) · (2, 0.5) com θ = 0.7, peso igual

Selecione uma rede ordenada (sem diluição ou pré-tensão) para comparar a Eq. (9) com ambas as leis de molas.

Não conseguimos reproduzir os pontos de simulação da Fig. 7: o artigo não informa como sua simulação de faixa mede o módulo de Young, e nossa leitura (força vertical à meia altura dividida pela largura e pela deformação) não reproduz esses pontos.

Rigidez e conectividade (uma verificação independente)

Fora do artigo: remover ligações aleatoriamente da rede triangular reduz a coordenação média z a partir de 6. A contagem de Maxwell (2 graus de liberdade por nó contra z/2 ligações por nó) indica que uma rede genérica não pode ser rígida abaixo de z = 4. Nossos módulos se anulam até o arredondamento em z = 3 e permanecem finitos acima de cerca de z = 3.9, o limiar de percolação de rigidez (por exemplo, Jacobs e Thorpe, PRL 75, 4051, 1995). Diluir até z ≈ 3.7 · …e colocá-la em tração

Todas as verificações

Limitações

• Não há código de referência público. As comparações se baseiam no texto do artigo, suas fórmulas, os dados vetoriais das Figs. 7–9 e uma digitalização da Fig. 5 rasterizada (±0.002).
• Não implementamos o método semianalítico próprio do artigo (os tensores não afins W da sua Eq. 20 na aproximação de campo médio); comparamos com seus resultados.
• Apenas resposta linear (pequena deformação). O artigo minimizou a energia não linear a uma deformação de 1%; nessa deformação a diferença é de ordem 1%.
• As resoluções ao vivo e a varredura automática são executadas atualmente de forma síncrona na thread principal. Redes grandes ou pré-tensionadas podem bloquear a navegação e Pausar até o fim do cálculo. Workers interrompíveis e reprodução pré-calculada estão pendentes.
• As bordas livres da faixa enviesam o ν medido em O(1/largura); os módulos periódicos não têm esse viés.
• Nosso conjunto (10 amostras por η) estima a média com cerca de ±0.005 em η ≤ 0.3 e ±0.01 perto de 0.5.
• Concordar com um artigo não é validar contra experimentos, e nada aqui é um resultado de projeto de engenharia.

Código e referências

• Implementação, testes, valores de referência e sua origem: research/metamaterials/grossman-2026 (AGPL-3.0).
• Artigo e material suplementar: doi:10.1103/17sg-gzhl, CC BY 4.0. As figuras desta página vêm das nossas próprias execuções; os pontos do artigo são seus dados publicados.
• Construção de espuma a partir de J. Liu, Y. Nie, H. Tong e N. Xu, Phys. Rev. Materials 3, 055607 (2019), Ref. 27 do artigo.

Citar isto

Stabileo Research, “Redes de molas auxéticas”, reprodução de Doron Grossman, Arezki Boudaoud (2026), doi:10.1103/17sg-gzhl. Implementação: research/metamaterials/grossman-2026 no commit 41d1ab4f8a8a. Motor: o motor do Stabileo no mesmo commit. Estão previstos links permanentes para esta compilação exata.