Redes de resortes auxéticas

Qué muestra el artículo

Muchos materiales, desde redes de metamateriales hasta membranas celulares, se modelan como redes de resortes de Hooke. Grossman y Boudaoud deducen el continuo elástico equivalente usando solo la geometría y la topología de la red, y calculan los desplazamientos no afines : las desviaciones locales de un estiramiento uniforme que hacen que las redes desordenadas sean más blandas que la estimación de deformación uniforme. Contrastan la teoría con simulaciones de una tira estirada para tres familias: redes triangulares ordenadas, un panal que se vuelve reentrante y una «espuma» formada al desplazar aleatoriamente los nodos de una red triangular, que se vuelve auxética (coeficiente de Poisson negativo: se ensancha al estirarla) al aumentar el desorden.

Probalo: red perfecta, η = 0 · η = 0.3 · η = 0.5, auxética

El modelo que resolvemos

Resortes e unen nodos; cada uno tiene rigidez k_e, longitud de reposo ℓ̄_e y longitud actual ℓ_e, y la energía es E = ½ Σ_e k_e (ℓ_e − ℓ̄_e)² (Ec. 3 del artículo). Una red de resortes es un reticulado articulado, así que el panel derecho la muestra en el visor de Stabileo: cada barra coloreada por su fuerza axial, sobre la forma deformada. Calculamos dos cosas.

• El ensayo de tira del artículo. Una tira de 13 × 52 (unos 800 nodos), con las filas superior e inferior sujetas y la superior desplazada según la deformación impuesta; el coeficiente de Poisson se obtiene del cambio de ancho entre los dos nodos de borde más cercanos a media altura, como en la Sec. V A del artículo. Se resuelve como un reticulado lineal con la factorización dispersa de Cholesky del motor de Stabileo; el mismo modelo resuelto por el análisis de reticulados 2D del motor de producción coincide hasta 10⁻¹⁴.
• Módulos efectivos exactos. Una celda periódica bajo deformación macroscópica impuesta, con todos los nodos libres para relajarse: la rigidez restante es el tensor elástico C de la red, incluido el ablandamiento no afín, sin efectos de borde. El tensor de deformación uniforme (afín) es la misma suma sin relajación.
• Pretensión (un agregado, ausente en los cálculos del artículo): longitudes de reposo acortadas una fracción, la red relajada hasta el equilibrio por el método de Newton y la respuesta lineal tomada alrededor de ese estado, incluido el aumento de rigidez por tracción.

Comparación con el artículo

Espuma: la transición auxética

El artículo informa que ν disminuye con el desorden η, cruza cero en η ≃ 0.46 y alcanza −0.1 cuando η → 0.5. Nuestros módulos periódicos exactos (10 muestras de 26 × 26 celdas por cada η) cruzan cero en η = 0.446 y dan ν = −0.148 en η = 0.5. Es una comparación entre protocolos: nuestra media de celda periódica queda dentro de todas las barras de error de las simulaciones de tiras sujetas de 13 × 52 de la Fig. 5 para η ≤ 0.45; no valida nuestro protocolo de tira. El marcador del artículo en η = 0.49 es una flecha debajo del gráfico (aproximadamente −0.14); nuestro valor allí es −0.117.

Nuestra implementación del protocolo de tira (trazo discontinuo) da 0.371 para la red perfecta, donde el valor exacto es 1/3 y la Fig. 5 del artículo muestra 0.331. El exceso es un efecto de borde libre de orden 1/ancho: las tiras de anchos 26 y 52 dan 0.354 y 0.345, extrapolando a 0.336. La tira del artículo no lo muestra, y su texto no detalla lo suficiente la medición para saber por qué.

Panal: Ec. (23)

Un panal de aristas unitarias cuyos vértices opuestos están separados 2ℓ es regular en ℓ = 1 y reentrante por debajo de ℓ = ½. El artículo triangula cada hexágono con un nodo central y resortes de rigidez 10⁻³ respecto de las aristas, y compara con la fórmula ν(ℓ) = −(2 − 4ℓ)(ℓ + ½)/(3 + 4ℓ − 4ℓ²). Nuestra red periódica coincide hasta 3.4 × 10⁻³ con esos resortes y hasta 3.4 × 10⁻⁶ con resortes de rigidez relativa 10⁻⁶: la diferencia es la rigidez propia de los resortes blandos. Mostrar ℓ = 0.25 (reentrante) · ℓ = 0.9

Redes ordenadas: una observación sobre las constantes elásticas

Para redes triangulares ordenadas, el artículo da una fórmula, la Ec. suplementaria (9), para ν en función de la forma de la celda (φ, ψ) y la orientación θ, graficada en la Fig. 4. Una red de Bravais no tiene desplazamientos no afines, así que la respuesta exacta es la de deformación uniforme, que calculamos independientemente. La Ec. suplementaria (9) y la línea de cálculo de la Fig. 7 (módulo de Young de la red (3, 2), leído de los datos vectoriales del PDF y multiplicado por el área de un triángulo) coinciden con ese resultado afín para constantes elásticas k_e ∝ 1/ℓ̄_e², hasta 10⁻¹⁴ y 6 × 10⁻⁴. Con constantes elásticas iguales obtenemos ν = 1/6 en (φ, ψ, θ) = (1, 0.5, 0), donde la Ec. (9) da 8/27. Esto podría reflejar una convención que no identificamos; está en nuestra lista de preguntas para los autores. Mostrar (1, 0.5) · (2, 0.5) con θ = 0.7, peso igual

Elegí una red ordenada (sin dilución ni pretensión) para comparar la Ec. (9) con ambas leyes de resortes.

No pudimos reproducir los puntos de simulación de la Fig. 7: el artículo no indica cómo mide el módulo de Young su simulación de tira, y nuestra lectura (fuerza vertical a media altura dividida por el ancho y la deformación) no reproduce esos puntos.

Rigidez y conectividad (un control independiente)

Fuera del artículo: quitar enlaces aleatoriamente de la red triangular reduce la coordinación media z desde 6. El conteo de Maxwell (2 grados de libertad por nodo frente a z/2 enlaces por nodo) indica que una red genérica no puede ser rígida por debajo de z = 4. Nuestros módulos se anulan hasta el redondeo en z = 3 y permanecen finitos por encima de aproximadamente z = 3.9, el umbral de percolación de rigidez (por ejemplo, Jacobs y Thorpe, PRL 75, 4051, 1995). Diluir hasta z ≈ 3.7 · …y ponerla en tracción

Todos los controles

Limitaciones

• No hay código de referencia público. Las comparaciones se basan en el texto del artículo, sus fórmulas, los datos vectoriales de las Figs. 7–9 y una digitalización de la Fig. 5 rasterizada (±0.002).
• No implementamos el método semianalítico propio del artículo (los tensores no afines W de su Ec. 20 en la aproximación de campo medio); comparamos con sus resultados.
• Solo respuesta lineal (pequeña deformación). El artículo minimizó la energía no lineal a una deformación de 1%; a esa deformación la diferencia es de orden 1%.
• Las resoluciones en vivo y el barrido automático se ejecutan actualmente de forma síncrona en el hilo principal. Las redes grandes o pretensadas pueden bloquear la navegación y Pausa hasta que termine el cálculo. Quedan pendientes los trabajadores interrumpibles y la reproducción precalculada.
• Los bordes libres de la tira sesgan el ν medido en O(1/ancho); los módulos periódicos no tienen ese sesgo.
• Nuestro conjunto (10 muestras por η) estima la media con aproximadamente ±0.005 en η ≤ 0.3 y ±0.01 cerca de 0.5.
• Coincidir con un artículo no es validar contra experimentos, y nada aquí es un resultado de diseño de ingeniería.

Código y referencias

• Implementación, pruebas, valores de referencia y su procedencia: research/metamaterials/grossman-2026 (AGPL-3.0).
• Artículo y material suplementario: doi:10.1103/17sg-gzhl, CC BY 4.0. Las figuras de esta página provienen de nuestras propias ejecuciones; los puntos del artículo son sus datos publicados.
• Construcción de espuma a partir de J. Liu, Y. Nie, H. Tong y N. Xu, Phys. Rev. Materials 3, 055607 (2019), Ref. 27 del artículo.

Citar esto

Stabileo Research, “Redes de resortes auxéticas”, reproducción de Doron Grossman, Arezki Boudaoud (2026), doi:10.1103/17sg-gzhl. Implementación: research/metamaterials/grossman-2026 en el commit 41d1ab4f8a8a. Motor: el motor de Stabileo en el mismo commit. Se prevén enlaces permanentes a esta compilación exacta.