Qué se reproduce
Esta primera versión implementa el equilibrio directo de la sección 2.1, Fig. 2 y Ec. (7). La red tiene un nudo libre y cuatro apoyos fijos. Su topología sigue el artículo; sus dimensiones y cargas se eligen aquí explícitamente, sin extraerlas de la figura. El artículo desarrolla diseño inverso diferenciable; esa contribución queda fuera de esta versión.
El método directo
La densidad de fuerza es la tracción axial dividida por la longitud. Densidades positivas dan una matriz de equilibrio definida positiva a una red anclada. El crate de investigación ensambla la matriz y usa Cholesky disperso del motor de producción para las tres coordenadas. Las fuerzas resultan de f = qℓ.
x* = (Σ qⱼ sⱼ + p) / Σ qⱼ
Los apoyos están en (±2, ±2, 0) m. La carga base hacia abajo es 4 kN y cada densidad es 1 kN/m. La Ec. (7) da z = −1 m, longitud = 3 m y tracción = 3 kN. Son valores derivados para nuestros datos, no observaciones numéricas publicadas.
Evidencia y tolerancias
Cinco casos comparan cada coordenada y fuerza con jax_fdm sin modificar en float64. Las tolerancias absolutas son 1e-12 m para coordenadas, 1e-12 kN para fuerzas de referencia y residuos, y 1e-9 kN para fuerzas de la cercha de producción. El archivo de validación registra cada error observado. La reproducción recorre cargas con el mismo solver Rust; no representa iteraciones de convergencia.
Error máximo de posición respecto a la referencia: 2,2 × 10−16 m; Error máximo de fuerza respecto a la referencia: 0 kN.
Referencia: jax_fdm, MIT, revisión 9a17fd379071fb3511ca914f46e8d1e8c286460e. El manifiesto registra SHA-256 del archivo fuente, hashes de resultados y entorno instalado. No se ha confirmado que esta revisión corresponda a la publicación.
Qué verifica el motor de producción
El editor recibe la geometría, las traslaciones fijas de los apoyos y la misma carga vertical como cercha elástica lineal 3D. Con rigideces iguales, la simetría hace que sus fuerzas axiales coincidan con FDM. Los tests también verifican el desplazamiento vertical con la solución elástica exacta. Es un análisis separado de barras inicialmente sin tensión: no se transfieren pretensión ni rigidez geométrica.
Aún no reproducido
No incluye adjuntos, gradientes, optimización, cáscaras, tracción/compresión mixta, cargas dependientes de la forma, pandeo ni afirmaciones de velocidad. No se reproducen formas optimizadas ni valores tabulados del artículo. La vista previa es una proyección indicada; el editor contiene el modelo espacial completo.
Investigación, no una herramienta de diseño. La concordancia con otro código no es validación experimental ni certificación de ingeniería.
Preguntas para los autores
¿Qué commit y entorno generaron el artículo? ¿Están disponibles los datos exactos de la Fig. 2? ¿Qué red pequeña y coordenadas, fuerzas y gradientes de referencia recomiendan para una posterior reproducción del diseño inverso?
Citar esto
Stabileo Research, “Cuatro cables, un equilibrio”, reproducción de Rafael Pastrana, Deniz Oktay, Kai-Uwe Bletzinger, Ryan P. Adams, Sigrid Adriaenssens (2026), doi:10.1016/j.cma.2026.118783. Implementación: research/form-finding/pastrana-2026 en el commit 8068cd16ea2c. Motor: el motor de Stabileo en el mismo commit. Se prevén enlaces permanentes a esta compilación exacta.