0
1
Teoría avanzada de estructuras
estáticamente determinadas
Escalante Espinoza, Nelver Javier; Calderón
Rodríguez, Luis Carlos; Tume Ramirez,
Hipolito; Lindo Santos, Prudencio Filotero;
Quillos Ruiz, Serapio Agapito; Del Castillo
Villacorta, Henry Joseph; Sparrow Alamo,
Edgar Gustavo
© Escalante Espinoza, Nelver Javier;
Calderón Rodríguez, Luis Carlos; Tume
Ramirez, Hipolito; Lindo Santos, Prudencio
Filotero; Quillos Ruiz, Serapio Agapito; Del
Castillo Villacorta, Henry Joseph; Sparrow
Alamo, Edgar Gustavo, 2026
Primera edición (1.ª ed.): abril, 2026
Editado por:
Editorial Mar Caribe ®
www.editorialmarcaribe.es
Av. Gral. Flores 547, 70000 Col. del
Sacramento, Departamento de Colonia,
Uruguay.
Diseño de carátula e ilustraciones: Editorial
Mar Caribe
Libro electrónico disponible en:
https://editorialmarcaribe.es/ark:/10951/i
sbn.9789915698960
Formato: Electrónico
ISBN: 978-9915-698-96-0
ARK: ark:/10951/isbn.9789915698960
Editorial Mar Caribe (OASPA): Como
miembro de la Open Access Scholarly
Publishing Association, apoyamos el
acceso abierto de acuerdo con el código de
conducta, la transparencia y las mejores
prácticas de OASPA para la publicación de
libros académicos y de investigación.
Estamos comprometidos con los más altos
estándares editoriales en ética y
deontología, bajo la premisa de «Ciencia
Abierta en América Latina y el Caribe»
Editorial Mar Caribe, rmante Nº 795 de
12.08.2024 de la Declaración de Berlín
... Nos sentimos obligados a abordar los
retos de Internet como un medio funcional
emergente para la distribución del
conocimiento. Obviamente, estos avances
pueden modicar signicativamente la
naturaleza de la publicación cientíca, así
como el actual sistema de garantía de
calidad... (Max Planck Society, ed. 2003, pp.
152-153).
CC BY-NC 4.0
Los autores pueden autorizar al público en
general a reutilizar sus obras únicamente
con nes no lucrativos, los lectores pueden
utilizar una obra para generar otra,
siempre que se dé crédito a la investigación,
y conceden al editor el derecho a publicar
primero su ensayo bajo los términos de la
licencia CC BY-NC 4.0.
Editorial Mar Caribe se adhiere a la
Recomendación relativa a la preservación
del patrimonio documental, comprendido el
patrimonio digital, y el acceso al mismo de
la UNESCO y a la Norma Internacional de
referencia para un sistema abierto de
información archivística (OAIS-ISO 14721).
Este libro está preservado digitalmente por
datasegura.info
DATA SEGURA
2
Editorial Mar Caribe
Teoría avanzada de estructuras
estáticamente determinadas
Colonia, Uruguay
2026
4
Índice
Introducción ........................................................................................... 7
Capítulo 1 ............................................................................................. 11
Teoría avanzada de estructuras estáticamente determinadas:
fundamentos, rigor matricial y aplicaciones complejas .................. 11
Formulación matricial dual y los cuatro subespacios
fundamentales .................................................................................. 13
Rigor matricial mediante la descomposición en valores
singulares (SVD) ............................................................................... 16
Aplicaciones complejas: reticulados periódicos, tensegridad y
eciencia computacional ................................................................. 19
Capítulo 2 ............................................................................................. 23
Fundamentos topológicos: determinación estática y cinemática
avanzada ............................................................................................... 23
Dualidad estático-cinemática y operadores estructurales
fundamentales .................................................................................. 24
Los cuatro subespacios vectoriales de Pellegrino y Calladine .... 25
Formulación grafo-teórica y bases topológicas en estructuras
reticuladas ........................................................................................ 27
Rigidez de primer orden, estructuras autotensadas y
tensegridad ....................................................................................... 29
Comparativa de marcos formulativos topológicos y matriciales 31
Capítulo 3 ............................................................................................. 34
Estabilidad y condicionamiento de sistemas isostáticos: un análisis
matricial, numérico y topológico ....................................................... 34
Fundamentos matriciales y criterio de Maxwell-Calladine .......... 35
Descomposición en valores singulares y espacios vectoriales ... 36
Condicionamiento numérico e inestabilidad geométrica ............ 38
5
Estabilidad por ................................................................................. 39
Avances topológicos y aplicaciones en metamateriales .............. 42
Capítulo 4 ............................................................................................. 46
Métodos variacionales y energéticos en la determinación de
deformaciones en sistemas isostáticos ............................................. 46
El principio de los trabajos virtuales y el método de la carga
unitaria .............................................................................................. 49
Principios variacionales y teoremas de Castigliano ..................... 52
Equivalencia formal entre Castigliano II y el principio de los
trabajos virtuales ............................................................................. 55
Técnica de cargas cticias o modales ............................................ 55
Inuencia de la deformación por cortante y factores de forma . 56
Simplicación del cálculo mediante integración gráca: La Regla
de Vereschaguin ............................................................................... 59
Análisis comparativo de métodos energéticos y variacionales ... 63
Capítulo 5 ............................................................................................. 67
Formulación matricial y vectorial del equilibrio en sistemas
discretos ............................................................................................... 67
Introducción y discretización de sistemas estructurales ............ 68
Dualidad cinemática-estática y el principio de los trabajos
virtuales ............................................................................................ 69
Transformación de coordenadas y mecanismos de ensamble
directo ............................................................................................... 71
Capítulo 6 ............................................................................................. 75
El marco canónico de Gilbert Strang y la caracterización
subespacial de Pellegrino-Calladine ................................................. 75
Los cuatro subespacios fundamentales de Pellegrino y Calladine
............................................................................................................ 76
6
Subespacio nulo o núcleo de la matriz de equilibrio null .......... 77
Rigidez ............................................................................................... 80
Estabilidad y condicionamiento de sistemas isostáticos: un
análisis matricial, numérico y topológico...................................... 82
Fundamentos matriciales y criterio de Maxwell-Calladine .......... 83
Descomposición en valores singulares y espacios vectoriales ... 84
Condicionamiento numérico e inestabilidad geométrica ............ 86
Estabilidad por preesfuerzo y rigidez geométrica........................ 87
Avances ............................................................................................. 90
Conclusión ............................................................................................ 94
Bibliografía ........................................................................................... 97
7
Introducción
El análisis y diseño de estructuras constituyen la columna
vertebral de la ingeniería civil y mecánica. Tradicionalmente, la
enseñanza y la literatura técnica han categorizado el comportamiento
estructural en dos grandes dominios: las estructuras estáticamente
determinadas (isostáticas) y las estáticamente indeterminadas
(hiperestáticas). Mientras que estas últimas suelen acaparar el interés
cientíco reciente debido a su redundancia estructural y a la difusión
masiva del Método de Elementos Finitos (MEF), las estructuras
isostáticas son a menudo relegadas a un nivel introductorio en la
formación de grado.
Sin embargo, esta simplicación oculta una realidad fundamental:
las estructuras estáticamente determinadas representan el estado
límite primario del equilibrio, la base fundamental para el análisis de
comportamiento no lineal, la trayectoria de cargas limpia y la referencia
insustituible para la validación de modelos numéricos complejos.
Esta obra, Teoría avanzada de estructuras estáticamente
determinadas, busca redenir la frontera del conocimiento en torno a los
sistemas isostáticos. A través de un enfoque riguroso que integra la
mecánica del continuo, la formulación matricial avanzada, la
optimización topológica y el análisis no lineal geométrico y material,
este libro demuestra que la isostaticidad no es un tema básico superado,
sino un campo vivo de investigación con aplicaciones críticas en la
ingeniería de vanguardia.
8
La necesidad de desarrollar un tratado cientíco y técnico
enfocado en la teoría avanzada de estructuras isostáticas responde a
tres factores determinantes:
i. Brecha teórico-práctica en la modelación numérica: el uso
indiscriminado de software comercial de cálculo ha generado
una dependencia metodológica donde los ingenieros confían
en resultados matriciales sin una vericación analítica rigurosa.
Las estructuras estáticamente determinadas ofrecen la única
plataforma analítica donde la trayectoria de fuerzas internas
es independiente de la rigidez del material y de las
deformaciones (en primer orden), lo que las convierte en la
herramienta superior para auditar y vericar algoritmos
complejos de cálculo.
ii. Optimización y estructuras de peso mínimo: en la ingeniería
aeroespacial, biomédica y en el diseño de cubiertas de grandes
luces, las estructuras isostáticas (como las cerchas o
reticulados óptimos tipo Michell) representan las soluciones de
peso mínimo teórico para sistemas de cargas dados.
Profundizar en su comportamiento avanzado permite diseñar
sistemas más sostenibles y ecientes.
iii. Comportamiento no lineal y dinámico: cuando las estructuras
isostáticas experimentan grandes desplazamientos (no
linealidad geométrica) o entran en el rango inelástico, el mito
de que sus reacciones no dependen de las deformaciones se
rompe. El análisis de cómo la variación geométrica altera las
9
ecuaciones de equilibrio en sistemas isostáticos ofrece una
visión profunda sobre fenómenos de inestabilidad, pandeo y
bifurcación que raras veces se abordan en la literatura
convencional.
En síntesis, revisitar la isostaticidad con las herramientas
analíticas y computacionales del siglo XXI no es un ejercicio regresivo,
sino una necesidad cientíca para fundamentar el diseño estructural
avanzado. El alcance de este trabajo abarca la investigación teórica y
analítica de sistemas estáticamente determinados bajo condiciones no
convencionales:
i. No linealidad geométrica y material: análisis del
comportamiento isostático bajo grandes deformaciones,
cambio de geometría en la trayectoria de carga y plasticidad
localizada.
ii. Estructuras articuladas y sistemas cinemáticamente exibles:
estudio dinámico y estático de mecanismos de cinemática
nita y estructuras desplegables.
iii. Optimización topológica isostática: algoritmos para la
generación de morfologías óptimas basadas en redes
isostáticas.
Como eje articulador de la investigación que fundamenta esta
obra, se formulan las siguientes interrogantes centrales: ¿De qué
manera los métodos de análisis avanzado (no linealidad geométrica,
optimización topológica y formulaciones matriciales de orden superior)
10
transforman la comprensión tradicional de las estructuras
estáticamente determinadas y potencian su aplicación en el diseño
estructural óptimo y la validación numérica de ingeniería?, ¿Cómo afecta
la gran deformación a la independencia entre la rigidez del elemento y
la distribución de esfuerzos en un sistema puramente isostático?.
Este libro unica y profundiza conceptos de mecánica del sólido
continuo, geometría diferencial aplicada a estructuras y métodos
energéticos avanzados sobre sistemas isoestáticos. Permite establecer
las fronteras teóricas donde termina el equilibrio estático puro y donde
la geometría de la estructura determina inequívocamente su respuesta.
Y sirve como un puente conceptual entre la estática elemental y el
análisis computacional no lineal, proporcionando a estudiantes de
posgrado e investigadores un marco que desarrolla intuición física y
rigor matemático.
11
Capítulo 1
Teoría avanzada de estructuras estáticamente
determinadas: fundamentos, rigor matricial y
aplicaciones complejas
El análisis de la estabilidad y la determinabilidad en la ingeniería
de estructuras ha experimentado una profunda transformación analítica
desde las formulaciones clásicas basadas en recuentos escalares hasta
los enfoques contemporáneos basados en el álgebra lineal.
Históricamente, la determinación del carácter isostático o hiperestático
de un sistema reticulado dependía del criterio enunciado por Maxwell,
el cual establece una relación aritmética simple entre el número de
barras , el número de nudos y el número de restricciones cinemáticas
impuestas por los apoyos . En un espacio euclidiano de dimensión
(donde para estructuras planas y para estructuras
espaciales), la condición necesaria para la determinabilidad se expresa
según la ecuación (Morales et al., 2024):
Sin embargo, la mecánica estructural moderna ha demostrado de
manera concluyente que la ecuación clásica de Maxwell constituye una
condición puramente necesaria pero insuciente para garantizar la
estabilidad geométrica y la determinación estática de una estructura.
Un sistema reticulado puede satisfacer de forma estricta la igualdad
escalar global y, no obstante, albergar simultáneamente redundancias
12
estáticas locales (estados de autotensión) e inestabilidades geométricas
regionales o globales (mecanismos inextensionales).
La superación de estas deciencias conceptuales comenzó con las
contribuciones de Calladine, quien extendió la regla de Maxwell al
considerar el rango de la matriz de equilibrio, y se consolidó mediante
la teoría desarrollada por Pellegrino y Calladine. Desde la perspectiva
de la mecánica matricial, la determinabilidad estática y la
determinabilidad cinemática no son propiedades independientes, sino
dos facetas duales que se derivan de las propiedades espectrales y
vectoriales de la matriz de equilibrio del sistema .
En términos rigurosos, una estructura es estáticamente
determinada cuando las fuerzas internas de todos sus miembros
(tensiones en las barras) pueden calcularse de manera única y biunívoca
a partir de las ecuaciones de equilibrio estático frente a cualquier
sistema de cargas externas aplicadas en los nudos, sin necesidad de
recurrir a las ecuaciones constitutivas del material ni a las condiciones
de compatibilidad de deformación (Zienkiewicz, 1977).
Complementariamente, un ensamblaje es cinemáticamente
determinado cuando la posición de cada uno de sus nudos está jada
de manera unívoca por la longitud de sus barras, garantizando que no
existan modos de deformación innitesimal o de mecanismo
inextensional.
Mediante la introducción del rango algebraico de la matriz de
equilibrio , el número de estados independientes de autotensión y el
13
número de mecanismos inextensionales independientes quedan
formalmente denidos por las relaciones:
donde representa el número total de grados de libertad nodales
del sistema. La sustracción algebraica de ambas ecuaciones conduce a
la generalización exacta de la regla de Maxwell-Calladine:
Esta relación fundamental establece que un ensamblaje
estructural es simultánea y estrictamente estática y cinemáticamente
determinado si y solo si y , lo que exige numéricamente que
el rango algebraico de la matriz de equilibrio sea igual tanto al número
de barras como al número de grados de libertad nodales (
), convirtiendo la matriz de equilibrio en un operador cuadrado y no
singular (Escalante et al., 2026).
Formulación matricial dual y los cuatro subespacios
fundamentales
La descripción rigurosa del comportamiento mecánico de un
ensamblaje de barras y uniones articuladas exige la construcción de un
marco matricial dual que conecte las variables de fuerza con las
variables de desplazamiento. Considérese una estructura general en un
espacio -dimensional con nudos, barras y restricciones de apoyo,
resultando en grados de libertad nodales (véase la tabla 1).
14
El estado de fuerzas internas se dene mediante el vector de
tensiones de las barras , mientras que la acción exterior aplicada
sobre los nudos se representa a través del vector de cargas nodales
. La condición de equilibrio estático en cada nudo del sistema da
lugar al sistema de ecuaciones lineales (Medina y Zavadskij, 2004):
donde es la matriz de equilibrio global de la estructura. Los
coecientes de la matriz están constituidos por los cosenos directores
de las barras referidos al sistema de ejes cartesianos globales.
De forma estrictamente dual, el comportamiento cinemático se
describe mediante el vector de desplazamientos nodales y el
vector de elongaciones de las barras . Bajo la hipótesis de
desplazamientos innitesimales, la compatibilidad geométrica entre los
desplazamientos de los nudos y las deformaciones del elemento se
expresa mediante la relación:
donde
siendo la matriz de compatibilidad cinemática, la cual es la
transpuesta exacta de la matriz de equilibrio . Esta dualidad estático-
cinemática permite categorizar exhaustivamente el comportamiento
mecánico del sistema mediante el análisis de los cuatro subespacios
vectoriales fundamentales del operador matricial .
15
Tabla 1: Interpretación mecánica y estructural de los cuatro subespacios
fundamentales
Subespacio
Vectorial
Dimensión
Algebraica
Interpretación Mecánica y Estructural
Espacio Nulo
Nul
Estados de Autotensión ( ): Dene las
distribuciones de tensiones en las barras
que satisfacen . Representa conjuntos
de fuerzas internas en equilibrio con cargas
externas nulas.
Espacio
Columna
Col
Cargas Admisibles de Equilibrio: Representa
el espacio de vectores de fuerza externa
que pueden ser equilibrados directamente
por tensiones internas en las barras ().
Espacio Nulo
Izquierdo
Nul
Mecanismos Inextensionales (): Dene los
vectores de desplazamiento nodal que no
producen elongación en las barras ().
Incluye movimientos rígidos y mecanismos
internos.
Espacio Fila
Row
Elongaciones Compatibles: Representa el
conjunto de vectores de deformación
derivados de desplazamientos nodales
continuos que no constituyen modos de
mecanismo.
Nota: Interpretado desde Voutetaki (2024).
La caracterización cinemática exige además segregar los
modos inextensionales en dos categorías con signicado físico
diferenciado: movimientos de cuerpo rígido y mecanismos internos. Si
la estructura se encuentra insucientemente vinculada a la cimentación,
hasta 6 modos (en tres dimensiones) o 3 modos (en dos dimensiones)
corresponden a desplazamientos de cuerpo rígido . Los modos
16
restantes representan mecanismos internos , los cuales
constituyen deformaciones geométricas en las que las barras rotan en
torno a sus articulaciones sin experimentar cambios de longitud a
primer orden (Ramalingam, 2009).
Rigor matricial mediante la descomposición en valores
singulares (SVD)
Aunque métodos como la eliminación gaussiana o la factorización
QR permiten calcular la dimensión de los subespacios vectoriales, la
descomposición en valores singulares (SVD) introducida formalmente en
este campo por Pellegrino ofrece la herramienta computacional y teórica
más sólida para el análisis cuantitativo de la matriz de equilibrio (Fallas
et al., 2021).
Para cualquier matriz de equilibrio de rango , la
descomposición SVD se expresa como:
donde y son matrices ortogonales cuyas columnas
constituyen bases ortonormales para los espacios de desplazamientos
nodales y de tensiones en las barras, respectivamente. La matriz
es una matriz diagonal cuyos primeros elementos son los
valores singulares estrictamente positivos ,
estando el resto poblado por ceros.
La estructura ortogonal de la SVD efectúa una partición directa
de los cuatro subespacios fundamentales dividiendo las matrices y
en bloques ordenados:
17
En esta formulación, el subbloque contiene vectores
columna que generan una base ortonormal para el espacio nulo Nul,
aislando los estados de autotensión. El subbloque alberga
vectores columna que denen la base ortonormal del espacio nulo
izquierdo Nul , revelando la totalidad de los mecanismos
inextensionales del sistema. Por su parte, las columnas de
corresponden al espacio de cargas equilibrables Col, mientras que
las columnas de delimitan el espacio de elongaciones
compatibles Row.
Esta descomposición facilita además el estudio de la interacción
no lineal entre los estados de autotensión y los mecanismos
inextensionales. Cuando un ensamblaje cinemáticamente
indeterminado activa un modo de mecanismo Nul, las barras
experimentan rotaciones que inducen elongaciones de segundo orden.
Si la estructura posee un estado de autotensión prévio Nul, la
interacción entre estas fuerzas internas y las variaciones geométricas
genera el fenómeno denominado fuerzas producto (product-forces).
El trabajo de segundo orden realizado por el estado de
autotensión a lo largo del desplazamiento del mecanismo se evalúa
mediante la matriz de rigidez geométrica :
Cuando para cualquier desplazamiento no nulo asociado
al modo , el estado de autotensión imparte rigidez de primer orden
18
al mecanismo, rigidizando cinemáticamente el ensamblaje
(comportamiento característico de sistemas de tensegridad, redes de
cables y domos preesforzados). Si , el mecanismo corresponde a
una inestabilidad geométrica de orden superior que permite
distorsiones nitas sin resistencia elástica inicial.
A partir del análisis de las dimensiones de estos subespacios
(véase la tabla 2), Affan y Calladine (1989) clasicaron la totalidad de
los ensamblajes estructurales en cuatro tipos fundamentales:
Tabla 2: Ensamblajes estructurales tipo
Tipo
Estructural
Estado de
Autotensión
(s)
Mecanismos
(m)
Condición de
Determinabilidad y
Comportamiento Mecánico
Tipo I
Estáticamente Determinado
y Cinemáticamente
Determinado: Sistema
isostático puro con respuesta
elástica única. La matriz es
cuadrada y no singular.
Tipo II
Estáticamente
Indeterminado y
Cinemáticamente
Determinado: Estructura
hiperestática convencional.
Posee capacidad de sobre-
restricción pero carece de
modos inestables.
Tipo III
Estáticamente Determinado
y Cinemáticamente
Indeterminado: Mecanismo
puro no estabilizable. Incapaz
19
Tipo
Estructural
Estado de
Autotensión
(s)
Mecanismos
(m)
Condición de
Determinabilidad y
Comportamiento Mecánico
de admitir estados de
pretensado inicial.
Tipo IV
Estáticamente
Indeterminado y
Cinemáticamente
Indeterminado: Ensamblaje
complejo (tensegridad). Los
mecanismos pueden
rigidizarse mediante la
aplicación de autotensión.
Aplicaciones complejas: reticulados periódicos,
tensegridad y eciencia computacional
El análisis algebraico de la matriz de equilibrio adquiere un papel
insustituible en el diseño de estructuras complejas contemporáneas,
como los metamateriales mecánicos, los reticulados periódicos y los
sistemas desplegables. Un teorema fundamental formulado por Guest y
Hutchinson demuestra que no existen estructuras reticulares innitas,
periódicas y articuladas en dos o tres dimensiones que sean
simultáneamente estática y cinemáticamente determinadas. Todo
reticulado innito periódico que satisfaga localmente la regla de
Maxwell albergará invariablemente modos blandos o estados de
autotensión no localizados (Vargas, 2025).
A pesar de esta restricción en el dominio innito, es posible
diseñar estructuras reticulares nitas e isostáticas mediante la
20
modicación rigurosa de las condiciones de contorno o la introducción
de perturbaciones geométricas. En la red plana de Kagome, un
ensamblaje periódico innito satisface numéricamente la ecuación de
Maxwell pero contiene líneas de mecanismos innitesimales. Al recortar
un dominio nito de dicha red, la pérdida de restricciones de frontera
genera múltiples mecanismos inestables. Para restablecer la
determinabilidad estática y cinemática estricta (), se deben
añadir elementos perimetrales denominados barras de parche (patch
bars), los cuales eliminan los modos exibles sin inducir hiperestaticidad.
En contraposición, en reticulados tridimensionales de alta
conectividad como la Malla Octogonal (Octet Truss), la elevada
conectividad nodal ( ) genera un grado masivo de redundancia
estática ( ). Para transformar esta topología en un sistema
estáticamente determinado manteniendo sus prestaciones resistentes,
se introduce una ruptura de simetría mediante la alteración de la
posición de los nudos centrales, derivando en la Malla Octogonal
Modicada (MOT). Esta perturbación desrumpe los planos de simetría,
eliminando las redundancias estáticas internas sin activar mecanismos
de colapso.
Esta cualidad de determinabilidad cinemática en mallas de
Kagome modicadas resulta indispensable para la creación de
estructuras adaptativas y sistemas de actuación. Como el espacio nulo
de la matriz de compatibilidad está libre de tensiones elásticas internas
durante la activación de ciertos grados de libertad, las barras pueden
utilizarse como actuadores lineales capaces de alterar la forma global
21
de la estructura sin almacenar energía elástica indeseada por distorsión
del reticulado.
En el ámbito del cálculo computacional y la optimización
topológica, la descomposición SVD de la matriz de equilibrio actúa como
un preprocesador geométrico de alta eciencia. En los ujos
tradicionales de diseño asistido por ordenador, la estabilidad geométrica
se evalúa ensamblando la matriz de rigidez global , asignando
propiedades mecánicas y resolviendo el sistema de ecuaciones elásticas.
Si la estructura es inestable, la matriz resulta singular,
interrumpiendo el proceso numérico.
El análisis directo de los subespacios de la matriz de equilibrio
permite detectar inestabilidades geométricas antes de denir la sección
transversal de los elementos ( ) o el módulo de elasticidad del
material ( ). La evaluación aislada de elimina la necesidad de
multiplicar , reduciendo sustancialmente la carga de cómputo en
las fases iniciales de exploración morfológica y garantizando la
viabilidad estructural previa al cálculo de tensiones y deformaciones.
El diagnóstico de la estabilidad geométrica y la determinabilidad
estática en estructuras articuladas complejas exige superar los criterios
escalares de recuento numérico para adoptar el marco del álgebra lineal
matricial. La matriz de equilibrio y su transpuesta, la matriz de
compatibilidad , constituyen los operadores fundamentales cuya
estructura de rangos determina unívocamente el comportamiento
estático y cinemático del sistema.
22
La aplicación de la Descomposición en Valores Singulares (SVD)
sobre la matriz de equilibrio proporciona una metodología rigurosa para
aislar los subespacios vectoriales asociados a los estados de
autotensión y a los mecanismos inextensionales. Esta caracterización
espectral no solo permite clasicar cualquier ensamblaje en los cuatro
tipos mecánicos de Pellegrino-Calladine, sino que posibilita evaluar la
capacidad de los estados de pretensado para impartir rigidez de primer
orden a mecanismos innitesimales mediante interacciones de rigidez
geométrica.
Así, el análisis de la matriz de equilibrio demuestra su valor
operativo en el diseño de metamateriales, reticulados periódicos nitos
y algoritmos de optimización estructural. Al permitir la detección
prematura de mecanismos locales sin necesidad de construir la matriz
de rigidez elástica global, el rigor matricial se consolida como una
herramienta analítica indispensable para el desarrollo de estructuras
ligeras, adaptativas e isostáticas de alta eciencia (Parras et al., 1985).
23
Capítulo 2
Fundamentos topológicos: determinación estática
y cinemática avanzada
La determinación estática y cinemática constituye uno de los ejes
fundamentales de la mecánica de estructuras discretas y la topología
computacional. Históricamente, la evaluación de la estabilidad e
hiperestaticidad estructural se abordaba mediante conteos escalares
simplicados formulados por James Clerk Maxwell. Sin embargo, la
aparición de sistemas reticulados complejos, ensamblajes autotensados
de tensegridad y estructuras espacialmente sobredeterminadas
demostró la insuciencia de los métodos puramente numéricos en
escenarios donde la geometría nodal y la conectividad topológica
interactúan no linealmente (Medina y Zavadskij, 2004).
El estudio avanzado de la respuesta estática y cinemática exige
un marco conceptual unicado basado en el álgebra lineal de espacios
vectoriales, la teoría de grafos topológicos y la teoría de homología. La
formulación rigurosa de la matriz de equilibrio, la matriz de
compatibilidad cinemática y la descomposición matricial mediante
descomposición en valores singulares (SVD), junto con la teoría de bases
cíclicas y cocíclicas, permite caracterizar sistemáticamente los estados
de autotensión, los mecanismos inextensionales e innitesimales, y el
acondicionamiento numérico de los operadores estructurales globales
(Monleón y Lázaro, 2011).
24
Dualidad estático-cinemática y operadores estructurales
fundamentales
La respuesta mecánica de cualquier ensamblaje estructural
formado por barras conectadas en nudos se rige por la relación dual
entre el equilibrio de fuerzas y la compatibilidad de deformaciones. En
una estructura articulada tridimensional con nudos activos y
elementos unidimensionales, las ecuaciones del equilibrio estático
linealizado en la conguración no deformada se expresan mediante el
sistema matricial (Monleón y Lázaro, 2011):
donde (con grados de libertad nodales y
esfuerzos internos en las barras) representa la matriz de equilibrio
global. El vector contiene los esfuerzos axiales generalizados en
los elementos, mientras que denota el vector de fuerzas externas
aplicadas sobre los nudos.
De forma simétrica, las condiciones cinemáticas dictan que el
vector de elongaciones axiales de los elementos se deriva del
vector de desplazamientos nodales mediante la matriz
cinemática o de compatibilidad :
El Teorema de los Trabajos Virtuales establece que el trabajo
realizado por las fuerzas internas debe ser igual al trabajo realizado por
las cargas externas para cualquier campo de desplazamientos y fuerzas
virtualmente admisibles. Esta equivalencia energética impone que la
25
matriz de compatibilidad sea exactamente la transpuesta de la matriz
de equilibrio :
La dualidad estático-cinemática queda sellada conceptualmente
mediante el par adjunto . Toda la información referente a la
estabilidad estática, la indeterminación cinemática, la rigidez latente y
los modos de deformación sin extensión de la estructura se encuentra
codicada en la estructura topológica y geométrica de la matriz .
Los cuatro subespacios vectoriales de Pellegrino y
Calladine
La aplicación del Teorema Fundamental del Álgebra Lineal a la
matriz de equilibrio transforma la comprensión del comportamiento
estático y cinemático. Para cualquier matriz de equilibrio de
dimensiones y rango rank, el espacio vectorial del sistema
se particiona en cuatro subespacios ortogonales fundamentales. Desde
la perspectiva de Farias et al. (2024) y Medina y Zavadskij (2004):
El primer subespacio es el Espacio Nulo (null), de dimensión
, el cual abarca todos los vectores de esfuerzo interno que
satisfacen la ecuación en ausencia de cargas externas ( ).
Físicamente, este subespacio dene los estados independientes de
autotensión o pretensado del sistema, donde el parámetro representa
el grado de indeterminación estática o hiperestaticidad.
El segundo subespacio es el Espacio Nulo Izquierdo (null ),
cuya dimensión es , donde representa el número de
26
restricciones cinemáticas impuestas por los apoyos de la cimentación.
Este subespacio está integrado por los vectores de desplazamiento
nodal que producen una elongación nula en las barras, es decir,
. Físicamente, estos vectores representan los modos de deformación
inextensional o mecanismos.
Los dos subespacios restantes corresponden al Espacio Fila
(range) y al Espacio Columna (range), ambos de dimensión . El
Espacio Fila contiene el conjunto de elongaciones de barra
cinemáticamente compatibles con un campo de desplazamientos
nodales continuos, mientras que el Espacio Columna está integrado por
las fuerzas externas nodales que pueden ser equilibradas por
combinaciones lineales de los esfuerzos axiales de los elementos.
Al vincular numéricamente las dimensiones de estos subespacios con
los recuentos topológicos, se obtiene la formulación generalizada de la
relación de Maxwell y Calladine:
Esta formulación demuestra que el número de estados de
autotensión y el número de mecanismos no pueden ser
determinados individualmente mediante recuentos de componentes
aislados. Dos ensamblajes con idéntico número de elementos ( ) y
nudos () pueden presentar valores radicalmente opuestos de y si
existen simetrías geométricas o alineaciones nodales que alteren el
rango de la matriz de equilibrio.
27
Para aislar numéricamente estos subespacios y caracterizar su
comportamiento cinemático, se recurre a la Descomposición en Valores
Singulares (SVD) de la matriz de equilibrio:
donde y son matrices ortogonales compuestas por
los subbloques y , y la matriz diagonal
contiene los valores singulares dispuestos en orden decreciente.
Las columnas de constituyen una base ortonormal del
espacio de autotensiones, mientras que las columnas de
forman la base ortonormal del espacio de mecanismos inextensionales.
Asimismo, este tratamiento matricial permite segregar los
mecanismos inextensionales en movimientos globales de cuerpo rígido
allowed por la cimentación () y mecanismos internos verdaderos
(), tales que .
Formulación grafo-teórica y bases topológicas en
estructuras reticuladas
Mientras que la formulación matricial directa opera sobre
coordenadas cartesianas nodales, la teoría de grafos estructurada por
Ali Kaveh analiza la topología de la conectividad para optimizar la
formulación matricial, la esparsidad y el acondicionamiento numérico
de los métodos clásicos de análisis estructural.
Una estructura discreta se modela topológicamente mediante un
grafo , donde el conjunto de vértices representa los nudos
o articulaciones, y el conjunto de aristas representa los elementos
28
unidimensionales. El primer número de Betti o rango del circuito
dene el número de ciclos independientes del grafo:
donde denota el número de componentes conexas del grafo. Para
el caso de pórticos tridimensionales continuos (estructuras con uniones
rígidas que transmiten momentos), el Grado de Indeterminación
Estática ( o ) depende de la topología cíclica según la
formulación de Kaveh:
En el caso de reticulados espaciales con articulaciones esféricas
perfeccionadas (cerchas tridimensionales), la ecuación adapta las
restricciones nodales según la expresión:
En la aplicación del Método de las Fuerzas, la matriz de exibilidad
del sistema se obtiene mediante la proyección de los elementos
exibles sobre una Base Cíclica Generalizada (GCB). La esparsidad y la
estructura de banda de la matriz son equivalentes al patrón de la
matriz de adyacencia de ciclos del grafo subyacente.
Cuando una estructura presenta variaciones marcadas en las
rigideces relativas de sus elementos, las elegidas arbitrariamente como
bases cíclicas inducen matrices de exibilidad con un elevado número
de condición , lo que genera errores de redondeo en el proceso de
inversión. Para resolver este problema, los algoritmos desarrollados por
Kaveh seleccionan bases cíclicas submínimas y subóptimas. En lugar de
29
minimizar únicamente la longitud de los ciclos (como en el algoritmo
estándar de Greedy), estos algoritmos minimizan las superposiciones o
solapamientos (overlaps) entre los ciclos de la base, lo que garantiza la
dominancia diagonal de la matriz de exibilidad y optimiza el
acondicionamiento numérico de la solución.
De forma dual, en el Método de los Desplazamientos, la matriz de
rigidez global se construye mediante una Base Cococíclica o conjunto
de cortadura (Cut-Set Basis). La generación de conjuntos de cortadura
asociados a los vértices aisla las ecuaciones de equilibrio nodal de
manera sistemática. La reordenación nodal óptima mediante la
construcción de grafos -totales permite reducir el ancho de banda y el
perl de la matriz de rigidez, maximizando la eciencia computacional
en la resolución de sistemas algebraicos lineales mediante eliminación
gaussiana.
Rigidez de primer orden, estructuras autotensadas y
tensegridad
La presencia de mecanismos inextensionales () no implica
indefectiblemente el colapso cinemático de una estructura. La existencia
simultánea de estados de autotensión ( ) introduce un efecto de
rigidez geométrica que estabiliza los modos de deformación
innitesimales (Pérez, 2015).
La matriz de rigidez tangencial de un ensamblaje pretensado
se compone de la suma de la matriz de rigidez elástica de primer orden
y la matriz de rigidez geométrica o de esfuerzo inicial :
30
donde la componente elástica responde a la formulación ,
siendo la matriz de rigidez constitutiva del material. Para un modo
inextensional denido por un vector de desplazamiento null, la
contribución elástica se anula ( ). En tales circunstancias, la
estabilidad de primer orden del mecanismo depende del escalar
resultante de la proyección de la rigidez geométrica:
Un valor indica que la autotensión conere rigidez de
primer orden al mecanismo innitesimal, transformándolo en un modo
cinemático rígido frente a acciones externas. Si , el modo
permanece como un mecanismo nito de orden superior. Por el contrario,
si , la autotensión induce una rigidez geométrica negativa que
desestabiliza la respuesta cinemática.
Esta rigidez inducida por pretensado constituye el principio de
funcionamiento de las estructuras de tensegridad—sistemas formados
por componentes de compresión discontinuos (puntales) suspendidos
dentro de una red continua de componentes en tracción (cables)—
popularizadas por Buckminster Fuller y Kenneth Snelson.
El diseño de estos sistemas requiere resolver el problema de
búsqueda de forma (form-nding), consistente en identicar posiciones
nodales y distribuciones de fuerza que garanticen el equilibrio estático
autotensado. La formulación de densidad de fuerza introducida por
Schek recongura la matriz de equilibrio deniendo las variables
31
, donde es el esfuerzo y la longitud del elemento . Esto
transforma el sistema en un operador lineal gobernado por la matriz de
densidad de fuerza . Para garantizar la superestabilidad de un sistema
de tensegridad, la matriz debe ser semidenida positiva con un rango
(donde es el número de nudos y la dimensión
espacial), asegurando la estabilidad global independientemente de las
propiedades elásticas de los cables y puntales.
En el ámbito del análisis estructural avanzado, investigaciones
recientes basadas en álgebras homológicas y cohaces celulares (cellular
cosheaves) han extendido estos principios topológicos a marcos rígidos
continuos que transmiten momentos ectores y esfuerzos cortantes.
Mediante el establecimiento de secuencias exactas largas de homología,
se ha logrado unicar la contabilidad clásica de Maxwell-Calladine
diseñada para cerchas articuladas con el primer número de Betti ()
representativo de los ciclos en marcos con momento. Este enfoque
resuelve las singularidades geométricas presentes en las cerchas
articuladas (como la conguración de Desargues) y absorbe las
inestabilidades innitesimales mediante campos de deformación por
exión local.
Comparativa de marcos formulativos topológicos y
matriciales
En la tabla 3 se presentan las características distintivas de los
cuatro marcos computacionales más relevantes empleados en la
determinación estática y cinemática avanzada, desde la interpretación
de los autores, con base en los fundamentos de Petkov (2026):
32
Tabla 3: Comparación de marcos formulativos topológicos y matriciales según Petkov
(2026)
Marco Teórico
Operador
Fundament
al
Representaci
ón Espacial /
Topológica
Entidad
Objetivo
Calibrada
Ventaja
Computacio
nal y Física
SVD de
Pellegrino-
Calladine
Matriz de
Equilibrio
Geometría
euclidiana
cartesiana
nodal
directriz
Subespaci
os null
y null
()
Caracteriza
mecanismos
de cuerpo
rígido,
mecanismos
internos y
estados de
autotensión.
Método de las
Fuerzas por
Bases Cíclicas
(Kaveh)
Matriz de
Flexibilidad
Global
Grafo del
sistema
e
hipergrafos
Bases
cíclicas
submínim
as y
subóptima
s (GCB)
Minimiza la
superposició
n (overlaps)
entre ciclos,
optimizando
el número de
condición.
Método de
Desplazamien
tos por Bases
Cococíclicas
(Kaveh)
Matriz de
Rigidez
Global
Grafo dual y
particiones de
vértices (Cut-
Sets)
Reducción
de ancho
de banda y
perl
Optimiza el
ordenamient
o nodal
mediante
grafos -
totales para
soluciones
matriciales
directas.
Cohaces
Celulares y
Homología
Estructural
Secuencias
Exactas de
Homología
Complejos
celulares
topológicos
-
dimensionale
s
Unicació
n de la
regla
Maxwell-
Calladine
con Betti
()
Integra la
exión y el
cortante;
elimina
singularidad
es
33
Marco Teórico
Operador
Fundament
al
Representaci
ón Espacial /
Topológica
Entidad
Objetivo
Calibrada
Ventaja
Computacio
nal y Física
geométricas
locales.
La determinación estática y cinemática en sistemas estructurales
ha evolucionado desde los esquemas de recuento escalar hacia un
marco matemático integrado en el álgebra lineal de espacios vectoriales,
la teoría de grafos y la homología topológica.
La descomposición ortogonal de la matriz de equilibrio mediante
la descomposición en valores singulares (SVD) constituye una
metodología rigurosa para la caracterización de sistemas discretos,
permitiendo diferenciar la movilidad cinemática de los estados de
autotensión y evaluando el papel de la rigidez geométrica en la
estabilización de mecanismos innitesimales (Fallas et al., 2021).
Asimismo, los métodos topológicos basados en grafos
demuestran que la eciencia computacional y la estabilidad numérica
no dependen exclusivamente de los algoritmos de resolución algebraica,
sino de la selección de bases cíclicas y cocíclicas que minimicen la
superposición topológica de los elementos.
34
Capítulo 3
Estabilidad y condicionamiento de sistemas
isostáticos: un análisis matricial, numérico y
topológico
La mecánica de estructuras y la física de materia condensada
convergen de manera fundamental en el estudio de los sistemas
isostáticos. Denidos en la estática clásica como aquellos sistemas
articulados donde el número de restricciones cinemáticas se equipara
de forma exacta con el número de grados de libertad independientes,
los sistemas isostáticos se sitúan en el límite crítico entre la exibilidad
mecánica y la rigidez estructural.
El análisis moderno de estas estructuras supera la regla de conteo
escalar tradicional, empleando formulaciones rigurosas del álgebra
lineal computacional como la descomposición en valores singulares
(SVD) de la matriz de equilibrio, la teoría de operadores diferenciales
discretos y el estudio de la topología de bandas en metamateriales
elásticos.
El comportamiento mecánico de una estructura isostática es
intrínsecamente sensible a su geometría exacta, a la conguración
espacial de sus conexiones y al condicionamiento numérico de sus
matrices gobernantes. Una pequeña perturbación geométrica en la
posición de sus nudos puede transformar un sistema nominalmente
35
estable en una estructura singularmente inestable o susceptible a
mecanismos innitesimales de energía cero (Lizarraga, 1990).
Por este motivo, la comprensión detallada de la estabilidad
geométrica, la rigidez por preesfuerzo y el condicionamiento matricial
resulta imprescindible para el diseño de estructuras desplegables, redes
de tensegridad de alto rendimiento y metamateriales mecánicos
programables.
Fundamentos matriciales y criterio de Maxwell-Calladine
El equilibrio estático y la compatibilidad cinemática de un sistema
articulado compuesto por barras y nudos en un espacio euclídeo de
dimensión se formulan mediante ecuaciones matriciales lineales que
vinculan los desplazamientos nodales , las elongaciones axiales
de las barras , las fuerzas externas aplicadas y las
tensiones o compresiones internas . La relación de equilibrio en
la conguración no deformada viene dada por la ecuación matricial
(Rocks y Mehta, 2024):
donde representa la matriz de equilibrio de la estructura. Por
aplicación directa del principio de los trabajos virtuales, la
compatibilidad cinemática dual que expresa el alargamiento o
acortamiento de las barras en función de los desplazamientos nodales
queda denida por la matriz de compatibilidad ,
cumpliéndose rigurosamente que .
36
El criterio histórico planteado por Maxwell establece que un
sistema espacial sin restricciones de apoyo es isostático cuando el
número de grados de libertad no restringidos satisface
, siendo el número de movimientos rígidos en el
espacio -dimensional. No obstante, la reformulación generalizada
introducida por Calladine demostró que la existencia simultánea de
estados de autoesfuerzo y mecanismos innitesimales invalida la regla
contable simple, imponiendo la relación vectorial:
donde denota el número de estados de autoesfuerzo independientes
y representa la dimensión total del espacio nulo cinemático, el cual
incluye tanto los movimientos de cuerpo rígido como los
mecanismos internos de deformación inextensional ( ).
Una estructura se dene como rigurosamente isostática en sentido
estático y cinemático cuando y , lo que conduce a un número
medio de coordinación o valencia nodal igual al valor crítico bajo
condiciones de frontera abiertas o periódicas.
Descomposición en valores singulares y espacios
vectoriales
La caracterización completa del comportamiento cinemático y
estático de un ensamblaje articulado se obtiene mediante la
descomposición en valores singulares (SVD) de la matriz de equilibrio
, de acuerdo con el marco teórico consolidado por Pellegrino y
Calladine. La matriz se factoriza en la forma (Alessandrini et al.,
2021):
37
donde y son matrices ortogonales cuya estructura
dene bases orquestadas para los espacios nodales y de barras,
respectivamente, mientras que contiene los valores
singulares , con rango (véase la tabla 4).
Tabla 4: Subespacio matricial ortogonal
Subespacio
Fundamental
Matricial
Denición
Matemática
Dimensión
Signicado Estructural y
Mecánico
Núcleo de
Null
Estados de autoesfuerzo
autosostenidos () en ausencia de
cargas externas.
Imagen de
Range
rango
Conjunto de fuerzas nodales
externas garantizadamente
equilibrables.
Núcleo de
Null
Mecanismos de deformación
inextensional ( ) y movimientos
rígidos.
Imagen de
Range
rango
Deformaciones elásticas de barras
cinemáticamente compatibles.
En una estructura isostática perfecta de rango completo, el rango
equivale al número total de barras , lo que anula la dimensión del
núcleo de ( ) y restringe la dimensión del núcleo de
estrictamente a los desplazamientos rígidos globales (). La
descomposición ortogonal permite aislar con precisión matemática si un
38
mecanismo detectado corresponde a una exibilidad nita o
innitesimal, aportando una herramienta de cálculo unicada para
resolver problemas estáticos e indeterminaciones cinemáticas sin
recurrir a suposiciones heurísticas.
Condicionamiento numérico e inestabilidad geométrica
La viabilidad práctica y la estabilidad física de un sistema
isostático no dependen únicamente del valor discreto de su rango
matricial, sino también de la sensibilidad de la matriz de equilibrio
frente a perturbaciones geométricas en la localización de sus nudos.
El condicionamiento numérico de la estructura se cuantica
formalmente mediante el número de condición de , el cual se expresa
como la razón entre su valor singular máximo y su valor singular mínimo
no nulo (Rodríguez et al., 2017):
Un sistema geométricamente bien congurado presenta un
número de condición moderado, lo que asegura que las fuerzas internas
calculadas varíen de forma suave y proporcional frente a desviaciones
menores en las cargas aplicadas. No obstante, cuando la disposición
geométrica de las barras se aproxima a una alineación crítica —por
ejemplo, cuando tres o más articulaciones se vuelven casi colineales o
cuando múltiples barras convergentes se orientan de forma cuasi-
coplanar— el valor singular mínimo tiende asintóticamente hacia
cero. En estas condiciones, aunque la matriz mantenga formalmente
39
su rango pleno desde un punto de vista puramente simbólico, el sistema
se vuelve numéricamente mal condicionado o singularmente
inestable.
El mal condicionamiento matricial desencadena fenómenos
mecánicos severos que comprometen la integridad de la estructura. En
primer lugar, se produce una amplicación extrema de las fuerzas
internas; cargas externas de magnitud modesta aplicadas en la
dirección del vector singular asociado a demandan tensiones y
compresiones de escala , superando con frecuencia los límites
de uencia o induciendo el colapso localizado.
En segundo lugar, el sistema demuestra una extrema sensibilidad
a las tolerancias de fabricación e imperfecciones del montaje, donde
desviaciones milimétricas en la posición de los nudos alteran
radicalmente la distribución interna de esfuerzos. Finalmente, la
estructura experimenta una pérdida drástica de rigidez efectiva frente a
ciertos patrones de carga, exhibiendo una respuesta cinemática
excesivamente blanda que simula mecánicamente la presencia de un
mecanismo innitesimal no restringido.
Estabilidad por preesfuerzo y rigidez geométrica
En sistemas que exhiben mecanismos innitesimales no
restringidos de primer orden (), la rigidez elástica convencional
resulta insuciente para resistir perturbaciones en las direcciones de
deformación de energía cero (véase la tabla 5). Sin embargo, la
estructura puede recuperar la estabilidad estática mediante la
40
introducción deliberada de un estado de autoesfuerzo previa aplicación
de cargas externas, principio fundamental que rige el comportamiento
de las estructuras de tensegridad y los reticulados hiperestáticos
pretensados (Agredo, 2016).
El comportamiento elástico linealizado de un sistema articulado
sometido a desplazamientos nodales viene gobernado por la matriz de
rigidez tangente , la cual se descompone aditivamente en dos
términos diferenciados:
donde es la matriz de rigidez elástica del material y representa la
matriz de rigidez geométrica o de esfuerzo inicial. La componente
elástica se ensambla mediante la matriz de compatibilidad y la matriz
diagonal de rigideces axiales diagsegún la formulación
. Dado que para cualquier modo de deformación
innitesimal Nullse satisface , la contribución elástica
a la energía de deformación es nula (
), dejando la
estabilidad cinemática del mecanismo supeditada de manera exclusiva
al término geométrico .
La matriz de rigidez geométrica se genera a partir del
vector de autoesfuerzo Nullo mediante la matriz de densidad de
fuerzas , adoptando la estructura tensorial , donde
representa la matriz identidad -dimensional.
41
Tabla 5: Estado de estabilidad, respuesta mecánica y dinámica
Estado de
Estabilidad
Condición Matemática de la
Matriz
Respuesta Mecánica y Dinámica
Rigidez Elástica
Primaria
Estabilidad por
Preesfuerzo
Superestabilidad
Inestabilidad por
Pandeo
De acuerdo con el teorema de Calladine y Pellegrino, la condición
necesaria y suciente para que un sistema sea estable por preesfuerzo
es que exista al menos un estado de autoesfuerzo de
tracción/compresión que haga que la forma cuadrática
sea
estrictamente denida positiva para todo mecanismo innitesimal no
nulo.
No obstante, la elevación indiscriminada de los niveles de
preesfuerzo introduce efectos no lineales contraintuitivos. Mientras que
incrementar la tensión en los miembros de tracción aumenta la rigidez
geométrica inicial, induce simultáneamente un incremento proporcional
en las fuerzas de compresión presentes en las barras rígidas. Cuando la
42
carga axial de compresión se aproxima a los límites de estabilidad
elástica o de pandeo euleriano de los elementos individualizados, los
autovalores negativos introducidos por la componente de compresión
en neutralizan los autovalores positivos de , conduciendo
inevitablemente al colapso por pandeo crítico del ensamblaje.
Avances topológicos y aplicaciones en metamateriales
El estudio de la isostaticidad y el condicionamiento matricial ha
experimentado una profunda transformación con la introducción de los
conceptos de la topología de bandas en la ciencia de materiales y los
metamateriales mecánicos. Los reticulados periódicos isostáticos,
conocidos ampliamente como redes de Kane-Lubensky, constituyen
estructuras innitas o repetitivas donde el número medio de
coordinación por nudo satisface exactamente la condición crítica de
Maxwell .
El análisis elastodinámico de estas redes se realiza transformando
la matriz de compatibilidad al espacio de momentos mediante el
teorema de Bloch, permitiendo obtener la matriz de compatibilidad de
onda en función del vector de onda dentro de la zona de Brillouin.
Cuando un reticulado isostático posee una estructura gapeada, su
espectro de fonones presenta una brecha acústica completa a
frecuencias no nulas, reservando los modos de energía cero únicamente
para las translaciones globales en el centro de la zona de Brillouin (
). La estructura global del cristal se caracteriza mediante un invariante
topológico no trivial denominado vector de polarización topológica ,
expresado como:
43
donde son los vectores primitivos de la red cristalina y representan
los índices de winding topológicos calculados a partir de las fases del
determinante de .
La existencia de una polarización topológica no nula altera
radicalmente la respuesta mecánica del metamaterial en sus fronteras
físicas y discontinuidades internas. En los bordes libres de la estructura,
la polarización topológica rompe la simetría traslacional de la rigidez,
provocando la acumulación localizada de mecanismos cinemáticos
exibles (oppy modes) en una supercie especíca, mientras que el
borde opuesto queda desprovisto de mecanismos y concentra los
estados de autoesfuerzo de frontera.
Esta polarización adquiere una relevancia crítica en el diseño de
interfaces y paredes de dominio. Al ensamblar dos dominios de
metamaterial isostático caracterizados por vectores de polarización
opuestos ( ), la pared de contacto alberga estados de
autoesfuerzo topológicamente protegidos de forma intrínseca. Dichos
estados canalizan y concentran las fuerzas externas aplicadas
exclusivamente a lo largo de la línea de interfaz, desencadenando
procesos de pandeo localizado altamente controlables cuando la carga
supera el umbral crítico de compresión.
Esta capacidad de guiar la deformación elástica de manera no
local permite fabricar dispositivos de protección contra impactos,
absorbedores de energía programables y estructuras metamateriales
44
adaptativas que conservan su funcionalidad aun en presencia de
desorden estructural o pequeñas variaciones geométricas.
En el campo de las estructuras desplegables y la robótica espacial,
los fundamentos de la isostaticidad modicada se aplican en el análisis
de sistemas de tensegridad asimétricos equipados con actuadores
agrupados o continuos. Al integrar cables continuos que pasan a través
de poleas o articulaciones lisas accionadas por un número reducido de
motores, se reconguran las ecuaciones de equilibrio estático de los
cables, reduciendo sustancialmente los grados de libertad de control y
manteniendo la rigidez geométrica global mediante la preservación
estricta de la condición denida positiva de la matriz en el subespacio
cinemático del sistema.
El análisis integrado de la estabilidad y el condicionamiento de los
sistemas isostáticos revela que la rigidez de un ensamblaje articulado
no constituye un atributo meramente cuantitativo derivado de la
densidad de barras, sino una propiedad algebraica y geométrica
altamente condicionada. El empleo de criterios escalares simples como
la regla de Maxwell resulta maniestamente insuciente para
diagnosticar la viabilidad de estructuras complejas en las proximidades
de conguraciones singulares o bajo perturbaciones nodales (Vergara,
2023).
La descomposición en valores singulares (SVD) de la matriz de
equilibrio se consolida como la herramienta analítica fundamental para
evaluar de manera unicada la estaticidad y la cinemática de cualquier
45
reticulado. El análisis sistemático de los cuatro subespacios vectoriales
fundamentales y del número de condición permite identicar no
solo la presencia de mecanismos innitesimales o redundancias
estáticas, sino también prever la susceptibilidad del sistema a
concentraciones destructivas de tensión e inestabilidades numéricas
provocadas por pequeñas desviaciones de montaje.
Asimismo, la introducción de preesfuerzo mediante la matriz de
rigidez geométrica proporciona un mecanismo efectivo para
estabilizar modos cinemáticos no restringidos, aunque exige un control
riguroso de las fuerzas de compresión resultantes para no inducir el
colapso por pandeo crítico. Finalmente, la convergencia entre la
mecánica matricial clásica y la topología de bandas abarca una de las
fronteras más prometedoras de la ingeniería de estructuras,
posibilitando el desarrollo de metamateriales capaces de guiar la
deformación, localizar estados de autoesfuerzo y ofrecer respuestas
mecánicas programables con una robustez intrínseca ante el desorden
geométrico.
46
Capítulo 4
Métodos variacionales y energéticos en la
determinación de deformaciones en sistemas
isostáticos
El análisis del comportamiento elástico de las estructuras
isostáticas bajo la acción de cargas externas se fundamenta en los
principios de conservación de la energía y en las leyes de la
termodinámica de sólidos deformables. Cuando un sistema isostático se
somete a un conjunto de fuerzas o momentos externos aplicados de
manera estática y gradual, las fuerzas realizan un trabajo externo ()
sobre la estructura.
Como respuesta a este estímulo, las bras del material
experimentan deformaciones internas, acumulando energía en forma de
energía elástica de deformación (). En ausencia de pérdidas disipativas
signicativas —tales como efectos térmicos no reversibles o uencia—,
el principio de conservación de la energía establece la equivalencia
estricta entre el trabajo externo ejecutado y la energía interna
almacenada ().
El trabajo externo realizado por una carga puntual aplicada
gradualmente desde un valor nulo hasta su magnitud nal,
incrementando linealmente el desplazamiento , se expresa
formalmente como (Pan et al., 2020; Zhang et al., 2015):
47
Para entender la respuesta estructural completa es preciso
diferenciar entre la energía de deformación ( ) y la energía
complementaria (, a menudo denotada como ). Grácamente, en un
diagrama carga-desplazamiento, la energía de deformación
representa el área situada bajo la curva, calculada como la integral de
la fuerza con respecto al desplazamiento:
Por su parte, la energía complementaria representa el área
comprendida por encima de la curva de carga-desplazamiento, limitada
por el eje de las cargas:
En materiales con comportamiento elástico lineal que obedecen
la Ley de Hooke, la relación entre fuerza y deformación es estrictamente
proporcional. Bajo esta condición cinemática y constitutiva, se cumple
que la energía de deformación y la energía complementaria son
numéricamente equivalentes ( ). Esta identidad es la piedra
angular que permite la aplicación intercambiable de teoremas
energéticos basados en desplazamientos o en fuerzas.
A nivel microscópico y seccional, la energía elástica de
deformación acumulada dentro del volumen del sólido deformable ()
48
se calcula integrando la densidad de energía elástica , la cual depende
del tensor de tensiones () y del tensor de deformaciones ():
En el contexto de la teoría de barras para sistemas isostáticos
unidimensionales, como vigas, pórticos y cerchas, la integración
tridimensional de la densidad energética se simplica mediante la
integración sobre la longitud de la barra () de los esfuerzos internos
resultantes: esfuerzo axil (), momento ector (), esfuerzo cortante ()
y momento de torsión () (véase la tabla 6).
Tabla 6: Modelo diferencial e integral para estimar la energía elástica de deformación
Solicitación
Interna
Expresión Diferencial de
Energía de Deformación
(dU)
Formulación Integral a lo
Largo del Elemento (U)
Esfuerzo Axil
()
Momento
Flector ()
Esfuerzo
Cortante ()
Momento
Torsor ()
En las formulaciones descritas, representa el módulo de
elasticidad longitudinal, es el módulo de elasticidad transversal o de
49
cortante, es el área de la sección transversal, es el momento
de inercia de la sección, denota la constante de rigidez a la torsión,
y es el factor de forma por cortante.
Dado que la energía de deformación depende del cuadrado de los
esfuerzos internos, el principio de superposición no puede aplicarse
directamente a la energía total. La energía producida por la acción
simultánea de dos sistemas de cargas no equivale a la suma algebraica
de las energías producidas por cada sistema actuando de forma aislada,
debido a la presencia de términos cruzados de trabajo mutuo. No
obstante, la linealidad constitutiva sí permite superponer los esfuerzos
internos y los desplazamientos resultantes, sobre los cuales se
fundamentan las formulaciones variacionales.
El principio de los trabajos virtuales y el método de la carga
unitaria
El Principio de los Trabajos Virtuales (PTV) constituye una de las
herramientas más potentes y versátiles de la mecánica de sólidos
deformables, pues no exige la existencia de una función de energía
potencial estricta y es aplicable tanto a análisis lineales como a
conguraciones no lineales. Para la determinación de desplazamientos
y giros en sistemas isostáticos, el PTV se operacionaliza mediante el
denominado Método de la Carga Unitaria o sistema de fuerzas virtuales
(Chowdhury y Summerscales, 2024).
50
El fundamento conceptual del PTV para el cálculo de deformaciones
se basa en la imposición de dos estados completamente independientes
pero compatibles entre sí (Wang et al., 2021):
• Estado Real (Sistema de Deformación): Representa la estructura
isostática real sometida al sistema de cargas externas reales,
variaciones térmicas o asentamientos, de donde se derivan las
deformaciones reales del material (, , ) y los desplazamientos
reales (, ).
• Estado Virtual (Sistema de Fuerzas): Consiste en la misma
estructura isostática en equilibrio, pero despojada de las cargas
reales y sometida exclusivamente a una carga virtual unitaria sin
dimensión () aplicada en el nudo o sección transversal en la
que se desea determinar el desplazamiento, y actuando en la
dirección especíca de dicho desplazamiento. Si se busca calcular
un giro angular (), la acción virtual corresponde a un momento
par unitario virtual ().
Al aplicar el principio de conservación del trabajo virtual, el
trabajo virtual externo realizado por la carga unitaria sobre el
desplazamiento real (
) se iguala al trabajo virtual interno realizado
por las fuerzas internas virtuales sobre las deformaciones reales del
material (
):
51
Donde
,
, y representan las leyes de esfuerzos
internos generadas únicamente por la carga o momento unitario virtual,
mientras que las deformaciones reales se expresan en función de las
cargas reales y la rigidez de la estructura:
Sustituyendo estas relaciones constitutivas en la formulación del
trabajo virtual, se obtiene la Ecuación General de Deformación por Carga
Unitaria para vigas y pórticos isostáticos:
En el caso especíco de cerchas o armaduras articuladas
isostáticas, donde los elementos estructurales solo experimentan
esfuerzos axiales puros y sus secciones transversales permanecen
constantes a lo largo de cada barra, la integral se reduce a una suma
discreta sobre el total de barras () de la estructura:
Una ventaja determinante del PTV frente a métodos geométricos
directos estriba en su capacidad para abordar efectos no mecánicos
52
como gradientes térmicos, errores de fabricación e imperfecciones
geométricas iniciales.
Cuando una barra experimenta una variación térmica homogénea
( ), su deformación longitudinal real viene dada por ,
donde es el coeciente de dilatación térmica. Si además existe un
gradiente térmico lineal entre la bra superior e inferior de una viga
( ), se induce una curvatura real
,
donde es el canto de la sección. Del mismo modo, si una barra de una
cercha presenta una imperfección longitudinal de fabricación (), esta
actúa directamente como una deformación impuesta.
La ecuación del trabajo virtual que integra estos efectos en un
sistema isostático adopta la siguiente forma ampliada:
Puesto que los sistemas isostáticos no poseen restrinjencias
redundantes, los cambios térmicos uniformes y los errores de
fabricación no generan esfuerzos internos adicionales (, ). Sin
embargo, producen deformaciones geométricas puras que el PTV evalúa
mediante la proyección de dichos desplazamientos reales sobre la red
de esfuerzos internos virtuales (
,
) dictados por la geometría estática
de la estructura.
Principios variacionales y teoremas de Castigliano
Los métodos variacionales formalizan la búsqueda de
conguraciones de equilibrio a través de la minimización o
53
estacionariedad de funcionales energéticos. El Principio de la Energía
Potencial Total ( ) establece que, de todas las
conguraciones cinemáticamente admisibles que un sólido deformable
puede tomar, la conguración que satisface el equilibrio estático es
aquella que minimiza el funcional de energía potencial total.
A partir del desarrollo del cálculo variacional aplicado a
estructuras elásticas, el ingeniero italiano Alberto Castigliano postuló
dos teoremas fundamentales que relacionan las derivadas parciales de
las funciones de energía con las fuerzas externas y los desplazamientos.
Primer teorema de Castigliano
Establece que si la energía elástica de deformación de una
estructura elástica se expresa en función de los desplazamientos
nodales , la derivada parcial de la función de energía de deformación
con respecto a un desplazamiento determinado es igual a la fuerza
externa actuante en ese punto y en la dirección de dicho
desplazamiento:
Análogamente, la derivada parcial con respecto a una rotación
angular produce el momento ector o par aplicado en dicho punto:
El Primer Teorema es de carácter general y resulta aplicable tanto
a sistemas elásticos lineales como a sistemas elásticos no lineales. Su
54
uso se orienta primariamente a la formulación de ecuaciones de rigidez
en el método de los elementos nitos y en el análisis de sistemas no
lineales, siendo menos directo para hallar desplazamientos puntuales
en estructuras isostáticas calculadas manualmente.
Segundo teorema de Castigliano
El Segundo Teorema de Castigliano postula que si la energía
complementaria de deformación de una estructura elástica se expresa
en función de las fuerzas aplicadas , la derivada parcial de la energía
complementaria respecto a una fuerza puntual externa es igual al
desplazamiento del punto de aplicación de dicha fuerza en la
dirección de su línea de acción:
Para el caso particular de materiales con respuesta elástico-lineal,
donde la energía de deformación es idéntica a la energía
complementaria ( ), el teorema permite reemplazar la energía
complementaria por la energía elástica convencional:
Sustituyendo la integral de energía de deformación por exión en
la formulación del Segundo Teorema, y aplicando la regla de la cadena
para la derivación bajo el signo integral, se obtiene:
55
Equivalencia formal entre Castigliano II y el principio de los
trabajos virtuales
La relación matemática entre el Segundo Teorema de Castigliano
y el Método de la Carga Unitaria muestra una coincidencia de fondo. Al
diferenciar el momento ector total producido por un conjunto de
cargas reales con respecto a una fuerza especíca , el término
analítico resultante corresponde exactamente a la ley de momentos
ectores generada por una carga unitaria aplicada en la posición de :
Por consiguiente, ambas expresiones convergen analíticamente
en la integral
. La diferencia radica en el procedimiento de
cálculo: mientras que el PTV requiere plantear dos sistemas estáticos
independientes (real y virtual), el Segundo Teorema de Castigliano opera
derivando algebraicamente la función parametrizada del momento
ector con respecto a las magnitudes de carga actora.
Técnica de cargas cticias o modales
El Segundo Teorema de Castigliano exige la presencia de una
fuerza o momento concentrado real en la posición exacta y en la
dirección del desplazamiento o giro que se desea calcular. Para evaluar
desplazamientos en secciones no cargadas o determinar giros en puntos
intermedios, se aplica la técnica de la carga cticia o auxiliar.
Este procedimiento consiste en introducir una carga puntual
cticia (o un momento cticio ) en la posición y dirección del
56
desplazamiento de interés. Posteriormente, se determinan las leyes de
esfuerzos internos incorporando la variable en las
ecuaciones de equilibrio, y se evalúa la derivada parcial de dichas leyes
respecto a la variable auxiliar (
). Finalmente, se ejecuta la integración
sustituyendo inmediatamente después de operar la
diferenciación algebraicamente, anulando el impacto de la carga cticia
sobre el estado de tensiones reales:
Esta metodología resulta equivalente en rigor al método de la
carga unitaria pero mantiene una representación matemática continua
mediante derivadas parciales.
Inuencia de la deformación por cortante y factores de
forma
En la teoría clásica de vigas de Euler-Bernoulli, se asume la
hipótesis cinemática de que las secciones transversales originalmente
planas y perpendiculares al eje neutro permanecen planas y
perpendiculares a dicho eje tras la deformación. Esta simplicación
ignora la distorsión angular provocada por las tensiones tangenciales
de cortante (). No obstante, en la teoría avanzada de vigas de
Timoshenko, la distorsión por cortante se incorpora explícitamente,
separando la pendiente total de la deformada (
) en dos componentes:
la rotación de la sección por exión pura () y la distorsión angular por
cortante () (David et al., 2020):
57
Puesto que las tensiones tangenciales ( ) no se distribuyen
uniformemente sobre el área de la sección transversal —variando de
acuerdo con la fórmula de Collignon-Jourawski
— , la bra
neutra y las bras extremas experimentan distintos niveles de
deformación por cortante. Para corregir esta disparidad al calcular la
energía de deformación interna mediante una formulación
unidimensional, se introduce el factor de forma por cortante ( ),
también denominado coeciente de corrección por cortante o factor de
distorsión:
Donde es el área total de la sección transversal, es el
momento de inercia centroidal, es el primer momento de área de
la sección cortada por encima del nivel , y es el ancho de la sección
transversal en dicha posición.
La integración de este factor en la formulación variacional permite
expresar la energía de deformación por cortante como
.
En secciones rectangulares macizas, las tensiones cortantes
siguen una distribución parabólica simétrica a lo largo del canto. El
esfuerzo alcanza un valor nulo en las bras extremas superior e inferior,
mientras que en la bra neutra centroidal registra su valor máximo,
equivalente a veces el esfuerzo cortante promedio ().
58
Al integrar la distribución parabólica del primer momento de área en la
fórmula del factor de forma, se obtiene exactamente .
En secciones circulares macizas, las tensiones tangenciales
presentan una distribución parabólica modicada con un esfuerzo
máximo en el centroide de veces el valor medio (
),
arrojando un factor de corrección de (véase la tabla 7). En
perles laminados en I o H (como los perles IPE o HEB), el esfuerzo
cortante es absorbido casi en su totalidad por el alma, por lo que el
factor de forma se aproxima por la relación entre el área total de la
sección y el área efectiva del alma (Manoli et al., 2022).
Tabla 7: Impacto de la deformación por tensión cortante y los factores de forma
Geometría de la
Sección
Transversal
Relación de
Tensión Cortante
Máxima (τmax
/τprom)
Factor de
Forma por
Cortante (α)
Área Efectiva
de Cortante
(Av=A/α)
Rectangular
Maciza
()
()
Circular Maciza
()
()
Tubo Circular de
Pared Delgada
Perl I (IPE /
HEB) - Flexión Eje
Mayor
59
La importancia relativa de la deformación por cortante respecto a
la exión depende de la esbeltez del elemento. Para una viga
simplemente apoyada con una carga puntual central , longitud ,
sección rectangular de canto y módulo de Poisson , la razón
entre la deexión por cortante y la deexión por exión viene dada por:
Esta relación matemática evidencia el impacto de la geometría:
• Vigas Esbeltas (): La deexión por cortante representa
menos del de la deformación total por exión, validando la
omisión del término de cortante según la hipótesis de Euler-
Bernoulli.
• Vigas de Gran Canto ( ): La deformación por cortante
supera el de la exión, pudiendo ser dominante en vigas muy
cortas (). En estas condiciones, las integrales energéticas
deben incluir el término de cortante con su correspondiente factor
de forma para no subestimar la deformación real.
Simplicación del cálculo mediante integración gráca: La
Regla de Vereschaguin
La determinación de desplazamientos mediante el Principio de los
Trabajos Virtuales requiere evaluar integrales del producto de dos
funciones de esfuerzo a lo largo de la estructura: la función de
momentos reales y la función de momentos virtuales
(Manglik, 2013):
60
En sistemas isostáticos, las fuerzas concentradas y distribuidas
producen diagramas de momentos reales descritos por funciones
polinómicas de grado o . Puesto que el sistema virtual está cargado
únicamente por una fuerza o momento unitario concentrado, el
diagrama de momentos virtuales
es lineal (polinomio de primer
grado ) dentro de cada tramo.
Demostración analítica del Teorema de Vereschaguin
Considérese la integral del producto de una función
(diagrama de momentos reales) y una función lineal
(diagrama de momentos virtuales) en el intervalo (Manglik,
2013):
Por denición geométrica, la integral
representa el área
total () del diagrama :
Asimismo, la posición del centro de gravedad () del diagrama
se dene mediante su momento estático de primer orden:
61
Sustituyendo estas relaciones geométricas en la expresión de la
integral, se obtiene:
El término corresponde a la función lineal
evaluada exactamente en la posición del centro de gravedad del
primer diagrama. Denotando esta ordenada virtual como , la
Regla de Vereschaguin se formula como:
Esta formulación establece que la integral del producto de dos
diagramas, siendo al menos uno lineal, es igual al área del primer
diagrama () multiplicada por la ordenada del diagrama lineal ()
medida en la vertical que pasa por el centro de gravedad del primer
diagrama.
Requisitos y condiciones formales de aplicabilidad
Para aplicar la Regla de Vereschaguin, deben cumplirse las
siguientes condiciones:
• La rigidez a la exión ( ) debe ser constante en el tramo de
integración. Si varía, el tramo debe subdividirse en segmentos de
rigidez uniforme.
• Al menos uno de los dos diagramas integrados debe ser
estrictamente lineal en el tramo analizado.
62
• El diagrama cuyo centro de gravedad se evalúa debe ser continuo.
Si presenta cambios de signo o discontinuidades, se fracciona en
guras geométricas elementales (rectángulos, triángulos,
parábolas) (véase la tabla 8).
Tabla 8: Aplicación de la regla de Vereschaguin
Geometría
del
Diagrama
Real M(x)
Geometría
del Diagrama
Virtual Mˉ(x)
Área del
Diagrama
Real (ΩM)
Ubicación
Centroidal
(xC desde
origen)
Valor de la
Integral
(EIMMˉdx)
Rectángulo
(base ,
altura )
Rectángulo
(base ,
altura
)
Rectángulo
(base ,
altura )
Triángulo
(base ,
altura máx
)
Triángulo
(base ,
altura máx
)
Triángulo
(vértices
alineados)
Triángulo
(base ,
altura máx
)
Triángulo
(vértices
opuestos)
Parábola de
2º Grado
(Carga C. D.)
Rectángulo
(base ,
altura
)
Parábola de
2º Grado
(Carga C. D.)
Triángulo
(base ,
altura máx
)
63
Análisis comparativo de métodos energéticos y
variacionales
La elección del método de cálculo para evaluar deformaciones en
estructuras isostáticas depende de la tipología estructural (vigas,
cerchas, pórticos), del tipo de solicitación aplicada y de si se requiere
una función matemática continua de la deformada o el desplazamiento
puntual en una sección especíca (véase la tabla 9).
Tabla 9: Análisis comparativo de métodos energéticos y variacionales
Criterio de
Comparació
n
Principio
de los
Trabajos
Virtuales
(PTV)
Segundo
Teorema de
Castigliano
Regla de
Vereschagui
n
Métodos
Geométricos
Directos
(Doble
Integración)
Fundamento
Mecánico
Variacional
/ Equilibrio
estático y
compatibili
dad
cinemática.
Variacional /
Derivada de
la energía
complement
aria elástica.
Integración
gráca del
PTV basada
en
propiedades
geométricas.
Integración
diferencial
directa de la
ecuación de
Euler-
Bernoulli.
Campo de
Aplicación
Principal
Cerchas,
pórticos
complejos y
vigas bajo
cargas
arbitrarias.
Vigas, arcos
y pórticos
con cargas
puntuales
parametriza
das.
Vigas e
hiperestática
s resueltas
mediante el
método de la
fuerza.
Vigas
sencillas con
pocas
discontinuida
des de carga
o sección.
Manejo de
Cambios
Térmicos e
Directo:
Incorpora
deformacio
nes
Complejo:
Exige
expresar la
energía
Moderado:
Requiere
transformar
la curvatura
Laborioso:
Modica la
ecuación
diferencial
64
Criterio de
Comparació
n
Principio
de los
Trabajos
Virtuales
(PTV)
Segundo
Teorema de
Castigliano
Regla de
Vereschagui
n
Métodos
Geométricos
Directos
(Doble
Integración)
Imperfeccio
nes
impuestas
directamen
te en el
trabajo
interno.
elástica en
función de
potenciales
térmicos.
térmica en
momentos
equivalentes.
imponiendo
condiciones
de borde
adicionales.
Determinaci
ón de
Deformadas
Continuas
Ineciente
(requiere
aplicar una
carga
unitaria
para cada
punto).
Ineciente
(exige
parametrizar
la fuerza
cticia como
variable de
posición).
Ineciente
(diseñado
para
desplazamie
ntos
discretos).
Directo:
Obtiene la
ecuación
algebraica
continua
.
Complejidad
Algebraica
Media
(resolución
de
integrales
denidas
simples).
Media-Alta
(derivación
parcial bajo
el signo
integral).
Baja:
Reemplaza la
integración
por
productos de
áreas y
ordenadas.
Alta
(resolución
de
constantes
de
integración
en múltiples
tramos).
Este cuadro comparativo reeja la complementariedad operativa
de los métodos. Para la formulación computacional y el desarrollo del
método de elementos nitos, las aproximaciones variacionales
(Castigliano I y Trabajo Virtual) resultan superiores debido a su
integración con el álgebra matricial. En contraste, para el cálculo
manual de estados límite de servicio (echas y giros) en la práctica del
65
diseño estructural, el Principio de los Trabajos Virtuales combinado con
la Regla de Vereschaguin proporciona el método más expedito y menos
propenso a errores de integración matemática (Luévanos 2013).
El análisis de los métodos variacionales y energéticos aplicados a
estructuras isostáticas conduce a las siguientes conclusiones analíticas:
Los métodos variacionales proporcionan un marco conceptual unicado
para el análisis de deformaciones provocadas por cargas mecánicas,
gradientes térmicos y errores de fabricación. La energía elástica de
deformación actúa como la función funcional generatriz de la cual
derivan tanto las condiciones de equilibrio estático como las relaciones
de compatibilidad cinemática.
Existe una equivalencia matemática directa entre el Segundo
Teorema de Castigliano y el Método de la Carga Unitaria. La derivada
parcial del momento ector respecto a una carga (
) es idéntica a la
distribución de momentos virtuales (
), conrmando que ambas
formulaciones representan expresiones alternativas del mismo principio
energético en el régimen elástico lineal.
La Regla de Vereschaguin simplica el cálculo analítico al
transformar integrales de productos de funciones en multiplicaciones
de áreas de momentos reales por ordenadas virtuales evaluadas en la
posición de sus centroides (). Esta metodología constituye una
herramienta eciente para el cálculo manual de deformaciones y sienta
las bases para la resolución de sistemas hiperestáticos mediante el
método de las exibilidades (Manglik, 2013).
66
Por último, el dominio de los principios variacionales y energéticos
permite comprender la rigidez global de las estructuras y fundamenta
los algoritmos utilizados por los programas modernos de análisis
estructural por elementos nitos.
67
Capítulo 5
Formulación matricial y vectorial del equilibrio en
sistemas discretos
El análisis del equilibrio en sistemas discretos constituye una de
las columnas vertebrales de la mecánica computacional y la ingeniería
estructural moderna. A diferencia de los medios continuos, cuya
respuesta mecánica se rige por ecuaciones diferenciales parciales
denidas sobre un dominio innito-dimensional, los sistemas discretos
—tales como armaduras articuladas, redes de cables, reticulados
espaciales y estructuras tensegríticas— reducen el comportamiento
global a un número nito de grados de libertad nodales mediante la
concentración de masas, fuerzas y rigideces en elementos
discretizados.
La formulación del equilibrio en estos sistemas puede abordarse
desde dos esquemas conceptuales complementarios: el enfoque
vectorial clásico, basado en el balance de fuerzas y momentos según las
leyes de Newton-Euler, y la formulación matricial algebraico-lineal, que
unica el equilibrio, la compatibilidad cinemática y las relaciones
constitutivas. La integración de ambos métodos mediante el Principio de
los Trabajos Virtuales revela una dualidad matemática profunda entre
la matriz de equilibrio y la matriz de compatibilidad cinemática,
estableciendo la base para la resolución sistemática de problemas
estática e hiperestáticamente indeterminados (Yilmaz, 2026).
68
Introducción y discretización de sistemas estructurales
La discretización de un sistema mecánico implica la reducción de
un problema continuo a una representación algebraica equivalente en
un número nito de puntos coordenados denominados nodos. En una
estructura discreta compuesta por elementos barras unidas en nodos
en un espacio de dimensión (), el comportamiento geométrico
global del sistema queda gobernado por el vector de desplazamientos
nodales , donde representa el número total de grados
de libertad no restringidos y corresponde al número de restricciones
cinemáticas impuestas por la cimentación o apoyos.
El balance estático en cada nodo del sistema exige que la suma
vectorial de las fuerzas internas ejercidas por los elementos
concurrentes y las fuerzas externas aplicadas sea idénticamente nula.
Para un nodo genérico , la ecuación vectorial de equilibrio estático se
representa como:
donde es el vector de fuerza interna ejercido por el elemento que
conecta el nodo con el nodo , representa el conjunto de nodos
adyacentes al nodo , y es el vector de carga externa aplicada
directamente en dicho nodo.
Cuando los elementos estructurales corresponden a barras
articuladas, capaces de responder únicamente a solicitaciones axiales
69
de tracción o compresión, la fuerza interna se descompone en el
producto de la magnitud de su tensión escalar y el vector unitario
director orientado a lo largo del eje geométrico del elemento :
siendo los vectores de posición espacial de los nodos y en la
conguración no deformada. Esta expresión vectorial local es la base
para la construcción global de los operadores algebraicos del sistema.
Dualidad cinemática-estática y el principio de los trabajos
virtuales
Al ensamblar escalarmente las ecuaciones vectoriales de
equilibrio para la totalidad de los grados de libertad del sistema, se
obtiene un sistema lineal de ecuaciones algebraicas que relaciona las
fuerzas internas de los elementos con el vector de fuerzas externas
aplicadas (Garcés et al., 2002):
donde es la matriz de equilibrio del sistema discreto,
es el vector de solicitaciones internas de los elementos,
y representa el vector global de cargas nodales.
La matriz de equilibrio es un operador puramente geométrico
cuyos componentes dependen de los cosenos directores de la
orientación de los elementos en el espacio. Simultáneamente, el
70
comportamiento cinemático del sistema discretizado describe cómo los
desplazamientos nodales producen elongaciones o
deformaciones en los elementos, compiladas en el vector de
deformaciones . La relación de compatibilidad
cinemática se expresa linealmente como:
donde representa la matriz de compatibilidad cinemática del
sistema.
La relación entre la matriz de equilibrio y la matriz de
compatibilidad cinemática está acoplada a través del Principio de los
Trabajos Virtuales. Considerando un estado de desplazamientos nodales
virtuales compatible con las restricciones del sistema, las
deformaciones virtuales inducidas en los elementos corresponden a
.
El trabajo virtual externo realizado por las fuerzas nodales
aplicadas a lo largo de los desplazamientos virtuales se expresa
como:
Por otro lado, el trabajo virtual interno absorbido por los
elementos sometidos a tensiones internas bajo las elongaciones
virtuales se dene como:
71
De acuerdo con el Principio de los Trabajos Virtuales, para
cualquier sistema en equilibrio estático, el trabajo virtual externo debe
igualar al trabajo virtual interno para todo desplazamiento virtual
admisible:
Sustituyendo la ecuación de compatibilidad cinemática y
la ecuación de equilibrio global en la condición de equivalencia
energética:
Dado que esta identidad debe satisfacerse de forma unánime para
cualquier campo de desplazamientos virtuales arbitrario y
cualquier vector de tensiones internas , se deduce la condición de
dualidad cinemática-estática fundamental:
Esta dualidad demuestra que la matriz de compatibilidad
cinemática de un sistema discreto es la traspuesta de la matriz de
equilibrio estático. Así, la relación de compatibilidad cinemática adopta
la forma canónica:
Transformación de coordenadas y mecanismos de ensamble
directo
Las ecuaciones mecánicas de un elemento discreto individual se
formulan de manera natural en su sistema de coordenadas locales,
donde el eje longitudinal concuerda con la directriz del elemento. Sin
72
embargo, el equilibrio global del sistema requiere que todos los vectores
de fuerza y desplazamiento de los elementos se expresen en un marco
de referencia ortogonal común denominado sistema de coordenadas
globales (Esquivel et al., 2014):
Matrices de transformación de coordenadas
Para un elemento unidimensional en un espacio bidimensional,
con una inclinación angular respecto al eje global , la matriz de
rotación local a global para un nodo individual se dene mediante
transformaciones ortogonales:
Para un elemento articulado con dos nodos, con cuatro grados de
libertad globales y desplazamientos locales
, la matriz de transformación de coordenadas
expandida se compone en bloques diagonales:
Debido a la naturaleza ortogonal de las transformaciones
geométricas rotacionales, la matriz cumple rigurosamente con la
propiedad:
y
73
Esta propiedad permite proyectar vectores de fuerza y
desplazamiento inversamente entre ambos sistemas sin alterar la
magnitud norma del vector ni escalar la energía del sistema. La
conversión cinemática local se evalúa como , mientras que la
proyección de fuerzas nodales locales al sistema global viene dada por
.
Ensamble directo de la matriz de rigidez global (Blanco et al., 2012)
La rigidez intrínseca de un elemento barra individual en su marco
local viene descrita por su matriz de rigidez elemental local :
donde es el módulo de elasticidad del material, es el área de la
sección transversal y es la longitud inicial del elemento. Aplicando las
ecuaciones de transformación local-global, la matriz de rigidez del
elemento referida al sistema de coordenadas globales se
obtiene mediante la ecuación de congruencia energética:
Sustituyendo los componentes de la matriz de rotación, la rigidez
elemental global de una barra bidimensional adopta la expresión
explícita:
74
El ensamble del sistema estructural global consiste en la
superposición topológica de las matrices dentro de la matriz de
rigidez global del sistema atendiendo a los mapas de
conectividad nodales de cada elemento:
donde representa la matriz de asignación booleana que relaciona los
grados de libertad locales del elemento con el vector de grados de
libertad globales del sistema. El comportamiento global queda denido
nalmente mediante la ecuación matricial fundamental de rigidez:
75
Capítulo 6
El marco canónico de Gilbert Strang y la
caracterización subespacial de Pellegrino-Calladine
En la mecánica de medios discretos, Gilbert Strang formalizó la
interacción entre las tres ecuaciones fundamentales del
comportamiento estructural (estática, cinemática y constitutiva) en una
estructura operativa de operadores matriciales conocida como el marco
.
El ujo operacional de este marco vincula secuencialmente los
espacios vectoriales del sistema. Los desplazamientos nodales se
transforman en elongaciones de elementos mediante la matriz de
compatibilidad . Dichas elongaciones generan tensiones internas a
través de la matriz constitutiva . Finalmente, las tensiones internas se
proyectan en fuerzas nodales externas mediante la matriz de equilibrio
. Las tres ecuaciones constitutivas del sistema se expresan como
(Brenner, 1967):
1. Ecuación Cinemática:
2. Ecuación Constitutiva:
3. Ecuación de Equilibrio:
donde es la matriz constitutiva del sistema, una matriz diagonal
constituida por la rigidez axial inherente a cada barra individual
diag
.
76
Componiendo de forma directa las tres etapas operacionales se
obtiene la relación entre fuerzas y desplazamientos nodales:
Se sintetiza de manera directa que la matriz de rigidez global del
sistema discreto admite una factorización canónica del tipo:
Esta formulación compacta garantiza la simetría y la
semidenición positiva de la matriz de rigidez , ya que para cualquier
matriz constitutiva denida positiva y cualquier operador real , la
forma cuadrática .
Los cuatro subespacios fundamentales de Pellegrino y
Calladine
El análisis algebraico de la matriz de equilibrio mediante
el Teorema Fundamental del Álgebra Lineal fue introducido por Sergio
Pellegrino y Christopher Calladine en 1986. Este enfoque descompone
los espacios de fuerzas y desplazamientos en los cuatro subespacios
vectoriales fundamentales de , permitiendo clasicar y cuanticar el
grado de indeterminación estática y cinemática de cualquier ensamblaje
de barras.
Si rank representa el rango algebraico de la matriz de
equilibrio , la estructura espacial del sistema se desglosa en los
subespacios descritos a continuación (Dugas et al., 1997; Brenner y
Butler, 1965):
77
Subespacio nulo o núcleo de la matriz de equilibrio null
Este subespacio vectorial reside en el espacio de tensiones
internas y posee una dimensión igual a:
null
Las bases de null representan los estados de autotensión
(states of self-stress) del sistema. Un estado de autotensión es una
combinación lineal de fuerzas internas nulltal que equilibran un
vector de carga externa nulo ( ). Estos estados permiten la
introducción de pretensado en la estructura sin alterar el balance de
cargas externas, siendo un requerimiento operativo en estructuras de
tenso-estructura y tensegridad (véase la tabla 10).
Subespacio Nulo a la Izquierda o Núcleo de la Traspuesta null
Este subespacio vectorial reside en el espacio de desplazamientos
nodales y su dimensión se calcula como:
null
Las bases de null representan los mecanismos o modos
inextensionales del sistema. Si un vector de desplazamiento
null, la elongación producida en la totalidad de las barras es nula
( ). El parámetro mide el número de grados de libertad
cinemáticos o el grado de indeterminación cinemática del sistema.
El espacio de mecanismos se segrega en dos componentes
físicas independientes:
78
donde es el número de movimientos de cuerpo rígido globales
permitidos por las restricciones de apoyo (hasta 3 en el plano y 6 en el
espacio tridimensional) y representa el número de mecanismos
internos cinemáticos de la estructura.
Subespacio Imagen o Espacio Columna range
Subespacio perteneciente a con dimensión igual a .
Representa el conjunto de todos los vectores de cargas externas que
pueden ser equilibrados directamente por el trabajo elástico interno de
las barras de la estructura.
Subespacio Fila o Espacio Imagen de la Traspuesta range
Subespacio perteneciente a con dimensión igual a .
Representa el conjunto de elongaciones de las barras que son
compatibles con un campo de desplazamientos nodales continuos
(Brenne y Butler, 1965).
Tabla 10: Subespacio vectorial. Interpretación mecánica y estructural
Subespacio
Vectorial
Espacio Dominio
Dimensión
Algebraica
Interpretación Mecánica y
Estructural
Núcleo de
(null)
(Tensiones)
Estados de autotensión o
pretensado estáticamente
admisible sin cargas externas.
Núcleo de
(null)
(Desplazamientos)
Modos de deformación
inextensional (Mecanismos
internos y cuerpo rígido ).
79
Subespacio
Vectorial
Espacio Dominio
Dimensión
Algebraica
Interpretación Mecánica y
Estructural
Imagen de
(range)
(Fuerzas nodales)
rank
Vectores de carga externa
resolubles mediante resistencia
elástica primaria.
Imagen de
(range)
(Elongaciones)
rank
Estados de deformabilidad
interna compatibles
cinemáticamente.
A partir de estas relaciones subespaciales se deduce la regla
combinatoria fundamental para la caracterización topológica
estructural de reticulados:
Esta ecuación sustituye la fórmula clásica de la estática basada
en el recuento simple de barras y nudos ( ), generalizando el
análisis para estructuras geométricamente complejas o con
singularidades locales donde el recuento lineal simple resulta
insuciente (véase la tabla 11) (Brenner y Butler, 1965).
Tabla 11: Operador matricial de relación para estructuras geométricamente complejas
o con singularidades locales
Dominio /
Nivel
Dominio Espacial /
Coordenadas
Vector de Estado /
Variable
Operador Matricial de
Relación
Cinemática
Local
Coordenadas
Locales de
Elemento
Desplazamientos
locales ( ) /
Elongaciones ()
Matriz de Compatibilidad
Local ()
80
Dominio /
Nivel
Dominio Espacial /
Coordenadas
Vector de Estado /
Variable
Operador Matricial de
Relación
Cinemática
Global
Coordenadas
Globales de
Estructura
Desplazamientos
nodales globales ()
Matriz de Transformación
Ortogonal ()
Respuesta
Constitutiva
Material / Sección
Transversal
Tensiones internas del
elemento ()
Matriz Constitutiva
Elasticidad ()
Equilibrio
Estático
Sistema Global
Ensamblado
Vector de Fuerzas
Exteriores ()
Matriz de Equilibrio (
) / Rigidez Global (
)
Rigidez Geométrica y Mecanismos de Primer Orden
Un aspecto derivado del análisis subespacial es el
comportamiento de estructuras donde existen mecanismos
inextensionales ( ). Un mecanismo inextensional es de primer
orden si no produce elongación de primer orden en las barras bajo un
desplazamiento innitesimal . Sin embargo, si la estructura posee
simultáneamente un estado de autotensión null, la interacción
entre la tracción preexistente de los elementos y los cambios
geométricos cuadráticos generados por el mecanismo origina una
matriz de rigidez geométrica o matriz de rigidez por tensión inicial .
La matriz de rigidez tangente global del sistema discreto queda
denida por la suma explícita (Pierre, 2022):
81
donde es la matriz de rigidez elástica primaria y es la
matriz de rigidez geométrica dependiente del nivel de esfuerzo del
estado de autotensión . Si la matriz es denida positiva sobre el
subespacio de mecanismos null, el estado de autotensión rigidiza
mecánicamente el mecanismo, dotando a la estructura de estabilidad
de segundo orden sin necesidad de restricciones de apoyo adicionales.
Este fenómeno constituye la base fundamental para el diseño estático y
geométrico de cables pretensados, mallas tensadas y estructuras
tensegríticas.
La formulación matricial y vectorial del equilibrio en sistemas
discretos proporciona un marco matemático riguroso que unica el
análisis de estructuras convencionales con el de sistemas avanzados. La
integración de la dualidad cinemático-estática mediante el Principio de
los Trabajos Virtuales demuestra formalmente que la compatibilidad
cinemática () y el equilibrio estático () son aspectos duales de la
misma geometría estructural.
El marco canónico propuesto por Strang, acoplado con la
descomposición en cuatro subespacios vectoriales de Pellegrino y
Calladine, transforma la evaluación de la estabilidad e hiperestatidad en
un problema computacional de álgebra lineal. La determinación del
rango y de los subespacios nulos de la matriz de equilibrio permite
identicar con precisión no solo la hiperestatidad () y la existencia de
mecanismos ( ), sino también la capacidad de estabilización de un
sistema mediante estados de pretensado ().
82
En el ámbito de la ingeniería computacional moderna, este
enfoque sistemático es la base para el desarrollo de algoritmos de
optimización topológica, simulaciones de grandes desplazamientos
mediante formulaciones corrotacionales y la resolución del equilibrio no
lineal en estructuras reticuladas exibles. La versatilidad de este
paradigma garantiza su aplicación continuada en el análisis estructural
y el diseño de nuevos sistemas mecánicos complejos.
Estabilidad y condicionamiento de sistemas isostáticos: un
análisis matricial, numérico y topológico
La mecánica de estructuras y la física de materia condensada
convergen de manera fundamental en el estudio de los sistemas
isostáticos. Denidos en la estática clásica como aquellos sistemas
articulados donde el número de restricciones cinemáticas se equipara
de forma exacta con el número de grados de libertad independientes,
los sistemas isostáticos se sitúan en el límite crítico entre la exibilidad
mecánica y la rigidez estructural (Sánchez, 2020).
El análisis moderno de estas estructuras supera la regla de conteo
escalar tradicional, empleando formulaciones rigurosas del álgebra
lineal computacional como la descomposición en valores singulares
(SVD) de la matriz de equilibrio, la teoría de operadores diferenciales
discretos y el estudio de la topología de bandas en metamateriales
elásticos.
El comportamiento mecánico de una estructura isostática es
intrínsecamente sensible a su geometría exacta, a la conguración
83
espacial de sus conexiones y al condicionamiento numérico de sus
matrices gobernantes. Una pequeña perturbación geométrica en la
posición de sus nudos puede transformar un sistema nominalmente
estable en una estructura singularmente inestable o susceptible a
mecanismos innitesimales de energía cero. Por este motivo, la
comprensión detallada de la estabilidad geométrica, la rigidez por
preesfuerzo y el condicionamiento matricial resulta imprescindible para
el diseño de estructuras desplegables, redes de tensegridad de alto
rendimiento y metamateriales mecánicos programables.
Fundamentos matriciales y criterio de Maxwell-Calladine
El equilibrio estático y la compatibilidad cinemática de un sistema
articulado compuesto por barras y nudos en un espacio euclídeo
de dimensión se formulan mediante ecuaciones matriciales lineales
que vinculan los desplazamientos nodales , las elongaciones
axiales de las barras , las fuerzas externas aplicadas y
las tensiones o compresiones internas . La relación de equilibrio
en la conguración no deformada viene dada por la ecuación matricial
(Rocks y Mehta, 2024):
donde representa la matriz de equilibrio de la estructura.
Por aplicación directa del principio de los trabajos virtuales, la
compatibilidad cinemática dual que expresa el alargamiento o
acortamiento de las barras en función de los desplazamientos nodales
queda denida por la matriz de compatibilidad ,
cumpliéndose rigurosamente que .
84
El criterio histórico planteado por Maxwell establece que un
sistema espacial sin restricciones de apoyo es isostático cuando el
número de grados de libertad no restringidos satisface
, siendo el número de movimientos rígidos en el
espacio -dimensional. No obstante, la reformulación generalizada
introducida por Calladine demostró que la existencia simultánea de
estados de autoesfuerzo y mecanismos innitesimales invalida la regla
contable simple, imponiendo la relación vectorial:
donde denota el número de estados de autoesfuerzo independientes
y representa la dimensión total del espacio nulo cinemático, el cual
incluye tanto los movimientos de cuerpo rígido como los
mecanismos internos de deformación inextensional (). Una
estructura se dene como rigurosamente isostática en sentido estático
y cinemático cuando y , lo que conduce a un número medio
de coordinación o valencia nodal igual al valor crítico bajo
condiciones de frontera abiertas o periódicas.
Descomposición en valores singulares y espacios
vectoriales
La caracterización completa del comportamiento cinemático y
estático de un ensamblaje articulado se obtiene mediante la
descomposición en valores singulares (SVD) de la matriz de equilibrio
, de acuerdo con el marco teórico consolidado por Pellegrino y
Calladine. La matriz se factoriza en la forma (Demmel, 1997):
85
donde y son matrices ortogonales cuya
estructura dene bases orquestadas para los espacios nodales y de
barras, respectivamente, mientras que contiene los valores
singulares , con rango (véase la tabla 12).
Tabla 12: Descomposición en valores singulares y espacios vectoriales
Subespacio
Fundamental
Matricial
Denición
Matemática
Dimensión
Signicado Estructural y
Mecánico
Núcleo de
Null
Estados de autoesfuerzo
autosostenidos ( ) en
ausencia de cargas
externas.
Imagen de
Range
rango
Conjunto de fuerzas
nodales externas
garantizadamente
equilibrables.
Núcleo de
Null
Mecanismos de
deformación inextensional
( ) y movimientos
rígidos.
Imagen de
Range
rango
Deformaciones elásticas
de barras
cinemáticamente
compatibles.
En una estructura isostática perfecta de rango completo, el rango
equivale al número total de barras , lo que anula la dimensión del
86
núcleo de ( ) y restringe la dimensión del núcleo de
estrictamente a los desplazamientos rígidos globales ( ). La
descomposición ortogonal permite aislar con precisión matemática si un
mecanismo detectado corresponde a una exibilidad nita o
innitesimal, aportando una herramienta de cálculo unicada para
resolver problemas estáticos e indeterminaciones cinemáticas sin
recurrir a suposiciones heurísticas.
Condicionamiento numérico e inestabilidad geométrica
La viabilidad práctica y la estabilidad física de un sistema
isostático no dependen únicamente del valor discreto de su rango
matricial, sino también de la sensibilidad de la matriz de equilibrio
frente a perturbaciones geométricas en la localización de sus nudos
(Fabré et al., 2007). El condicionamiento numérico de la estructura se
cuantica formalmente mediante el número de condición de , el cual
se expresa como la razón entre su valor singular máximo y su valor
singular mínimo no nulo:
Un sistema geométricamente bien congurado presenta un
número de condición moderado, lo que asegura que las fuerzas internas
calculadas varíen de forma suave y proporcional frente a desviaciones
menores en las cargas aplicadas. No obstante, cuando la disposición
geométrica de las barras se aproxima a una alineación crítica —por
ejemplo, cuando tres o más articulaciones se vuelven casi colineales o
cuando múltiples barras convergentes se orientan de forma cuasi-
87
coplanar— el valor singular mínimo tiende asintóticamente
hacia cero. En estas condiciones, aunque la matriz mantenga
formalmente su rango pleno desde un punto de vista puramente
simbólico, el sistema se vuelve numéricamente mal condicionado o
singularmente inestable.
El mal condicionamiento matricial desencadena fenómenos
mecánicos severos que comprometen la integridad de la estructura. En
primer lugar, se produce una amplicación extrema de las fuerzas
internas; cargas externas de magnitud modesta aplicadas en la
dirección del vector singular asociado a demandan tensiones y
compresiones de escala , superando con frecuencia los límites
de uencia o induciendo el colapso localizado. En segundo lugar, el
sistema demuestra una extrema sensibilidad a las tolerancias de
fabricación e imperfecciones del montaje, donde desviaciones
milimétricas en la posición de los nudos alteran radicalmente la
distribución interna de esfuerzos (Fabré et al., 2007).
Ahora bien, la estructura experimenta una pérdida drástica de
rigidez efectiva frente a ciertos patrones de carga, exhibiendo una
respuesta cinemática excesivamente blanda que simula mecánicamente
la presencia de un mecanismo innitesimal no restringido.
Estabilidad por preesfuerzo y rigidez geométrica
En sistemas que exhiben mecanismos innitesimales no
restringidos de primer orden (), la rigidez elástica convencional
resulta insuciente para resistir perturbaciones en las direcciones de
88
deformación de energía cero. Sin embargo, la estructura puede
recuperar la estabilidad estática mediante la introducción deliberada de
un estado de autoesfuerzo previa aplicación de cargas externas,
principio fundamental que rige el comportamiento de las estructuras de
tensegridad y los reticulados hiperestáticos pretensados (véase la tabla
13).
El comportamiento elástico linealizado de un sistema articulado
sometido a desplazamientos nodales viene gobernado por la matriz de
rigidez tangente , la cual se descompone aditivamente en dos
términos diferenciados:
donde es la matriz de rigidez elástica del material y representa la
matriz de rigidez geométrica o de esfuerzo inicial. La componente
elástica se ensambla mediante la matriz de compatibilidad y la matriz
diagonal de rigideces axiales diag según la formulación
. Dado que para cualquier modo de deformación
innitesimal Null se satisface , la contribución
elástica a la energía de deformación es nula (
), dejando la
estabilidad cinemática del mecanismo supeditada de manera exclusiva
al término geométrico .
La matriz de rigidez geométrica se genera a partir del
vector de autoesfuerzo Null o mediante la matriz de densidad de
fuerzas , adoptando la estructura tensorial , donde
representa la matriz identidad -dimensional.
89
Tabla 13: Estabilidad por preesfuerzo y rigidez geométrica
Estado de
Estabilidad
Condición Matemática
de la Matriz KT
Respuesta Mecánica y
Dinámica
Rigidez Elástica
Primaria
es denida positiva
sobre Null
La geometría imparte rigidez
directa sin requerir
preesfuerzo previo.
Estabilidad por
Preesfuerzo
Null
El autoesfuerzo aporta
rigidez de segundo orden a
los mecanismos
innitesimales.
Superestabilidad
es semi-denida
positiva globalmente y
no requiere elasticidad
Estabilidad incondicional e
independiente del módulo de
Young de los elementos.
Inestabilidad por
Pandeo
Las fuerzas de compresión
superan los umbrales de
estabilidad elástica global.
De acuerdo con el teorema de Calladine y Pellegrino, la condición
necesaria y suciente para que un sistema sea estable por preesfuerzo
es que exista al menos un estado de autoesfuerzo de
tracción/compresión que haga que la forma cuadrática
sea
estrictamente denida positiva para todo mecanismo innitesimal no
nulo.
No obstante, la elevación indiscriminada de los niveles de
preesfuerzo introduce efectos no lineales contraintuitivos. Mientras que
incrementar la tensión en los miembros de tracción aumenta la rigidez
geométrica inicial, induce simultáneamente un incremento proporcional
90
en las fuerzas de compresión presentes en las barras rígidas. Cuando la
carga axial de compresión se aproxima a los límites de estabilidad
elástica o de pandeo euleriano de los elementos individualizados, los
autovalores negativos introducidos por la componente de compresión
en neutralizan los autovalores positivos de , conduciendo
inevitablemente al colapso por pandeo crítico del ensamblaje.
Avances topológicos y aplicaciones en metamateriales
El estudio de la isostaticidad y el condicionamiento matricial ha
experimentado una profunda transformación con la introducción de los
conceptos de la topología de bandas en la ciencia de materiales y los
metamateriales mecánicos. Los reticulados periódicos isostáticos,
conocidos ampliamente como redes de Kane-Lubensky, constituyen
estructuras innitas o repetitivas donde el número medio de
coordinación por nudo satisface exactamente la condición crítica de
Maxwell (Zhoufei et al., 2025).
El análisis elastodinámico de estas redes se realiza transformando
la matriz de compatibilidad al espacio de momentos mediante el
teorema de Bloch, permitiendo obtener la matriz de compatibilidad de
onda en función del vector de onda dentro de la zona de Brillouin.
Cuando un reticulado isostático posee una estructura gapeada, su
espectro de fonones presenta una brecha acústica completa a
frecuencias no nulas, reservando los modos de energía cero únicamente
para las translaciones globales en el centro de la zona de Brillouin (
). La estructura global del cristal se caracteriza mediante un invariante
91
topológico no trivial denominado vector de polarización topológica ,
expresado como (Qian et al., 2025):
donde son los vectores primitivos de la red cristalina y representan
los índices de winding topológicos calculados a partir de las fases del
determinante de .
La existencia de una polarización topológica no nula altera
radicalmente la respuesta mecánica del metamaterial en sus fronteras
físicas y discontinuidades internas. En los bordes libres de la estructura,
la polarización topológica rompe la simetría traslacional de la rigidez,
provocando la acumulación localizada de mecanismos cinemáticos
exibles (oppy modes) en una supercie especíca, mientras que el
borde opuesto queda desprovisto de mecanismos y concentra los
estados de autoesfuerzo de frontera.
Esta polarización adquiere una relevancia crítica en el diseño de
interfaces y paredes de dominio. Al ensamblar dos dominios de
metamaterial isostático caracterizados por vectores de polarización
opuestos ( ), la pared de contacto alberga estados de
autoesfuerzo topológicamente protegidos de forma intrínseca. Dichos
estados canalizan y concentran las fuerzas externas aplicadas
exclusivamente a lo largo de la línea de interfaz, desencadenando
procesos de pandeo localizado altamente controlables cuando la carga
supera el umbral crítico de compresión.
92
Esta capacidad de guiar la deformación elástica de manera no
local permite fabricar dispositivos de protección contra impactos,
absorbedores de energía programables y estructuras metamateriales
adaptativas que conservan su funcionalidad aun en presencia de
desorden estructural o pequeñas variaciones geométricas.
En el campo de las estructuras desplegables y la robótica espacial,
los fundamentos de la isostaticidad modicada se aplican en el análisis
de sistemas de tensegridad asimétricos equipados con actuadores
agrupados o continuos. Al integrar cables continuos que pasan a través
de poleas o articulaciones lisas accionadas por un número reducido de
motores, se reconguran las ecuaciones de equilibrio estático de los
cables, reduciendo sustancialmente los grados de libertad de control y
manteniendo la rigidez geométrica global mediante la preservación
estricta de la condición denida positiva de la matriz en el subespacio
cinemático del sistema.
El análisis integrado de la estabilidad y el condicionamiento de los
sistemas isostáticos revela que la rigidez de un ensamblaje articulado
no constituye un atributo meramente cuantitativo derivado de la
densidad de barras, sino una propiedad algebraica y geométrica
altamente condicionada. El empleo de criterios escalares simples como
la regla de Maxwell resulta maniestamente insuciente para
diagnosticar la viabilidad de estructuras complejas en las proximidades
de conguraciones singulares o bajo perturbaciones nodales (Suresh et
al., 2021).
93
La descomposición en valores singulares (SVD) de la matriz de
equilibrio se consolida como la herramienta analítica fundamental para
evaluar de manera unicada la estaticidad y la cinemática de cualquier
reticulado. El análisis sistemático de los cuatro subespacios vectoriales
fundamentales y del número de condición permite identicar no
solo la presencia de mecanismos innitesimales o redundancias
estáticas, sino también prever la susceptibilidad del sistema a
concentraciones destructivas de tensión e inestabilidades numéricas
provocadas por pequeñas desviaciones de montaje.
Asimismo, la introducción de preesfuerzo mediante la matriz de
rigidez geométrica proporciona un mecanismo efectivo para
estabilizar modos cinemáticos no restringidos, aunque exige un control
riguroso de las fuerzas de compresión resultantes para no inducir el
colapso por pandeo crítico. Finalmente, la convergencia entre la
mecánica matricial clásica y la topología de bandas abarca una de las
fronteras más prometedoras de la ingeniería de estructuras,
posibilitando el desarrollo de metamateriales capaces de guiar la
deformación, localizar estados de autoesfuerzo y ofrecer respuestas
mecánicas programables con una robustez intrínseca ante el desorden
geométrico.
94
Conclusión
Teoría avanzada de estructuras estáticamente determinadas,
consolidó una premisa fundamental para la ingeniería estructural
moderna: la determinidad estática no constituye únicamente una fase
introductoria de la mecánica de materiales, sino que representa el
núcleo riguroso donde el equilibrio, la eciencia geométrica y la física
del comportamiento de los cuerpos se integran con absoluta precisión.
De la evidencia, se ha trascendido la resolución trivial de sistemas
bidimensionales para abordar la formulación matemática y
metodológica de sistemas espaciales complejos, métodos energéticos
avanzados y la respuesta frente a acciones variables. A diferencia de los
sistemas hiperestáticos, donde la rigidez y las deformaciones
condicionan la distribución de esfuerzos, las estructuras isostáticas
dependen exclusivamente de la geometría y de las leyes fundamentales
de Newton:
y
Esta independencia respecto a las propiedades del material o
variaciones térmicas convierte al análisis determinado en la prueba de
fuego denitiva para comprender los caminos de carga (load paths) en
cualquier tipología estructural.
El estudio avanzado ha demostrado que la simple condición
numérica entre grados de libertad y restricciones ( o ) es
insuciente. La investigación ratica la importancia crítica de evitar los
95
mecanismos innitesimales y la inestabilidad cinemática implícita,
aspectos determinantes en el diseño seguro de reticulados espaciales,
arcos y estructuras articuladas.
Además, el desarrollo del principio de los trabajos virtuales (PTV)
y de los teoremas de energía aplicados a sistemas determinados
favorece el cálculo analítico de desplazamientos y sienta las bases para
el análisis matricial y computacional. Asimismo, en términos rigurosos,
una estructura es estáticamente determinada cuando las fuerzas
internas de todos sus miembros (tensiones en las barras) pueden
calcularse de manera única y biunívoca a partir de las ecuaciones de
equilibrio estático frente a cualquier sistema de cargas externas
aplicadas en los nudos, sin necesidad de recurrir a las ecuaciones
constitutivas del material ni a las condiciones de compatibilidad de
deformación.
En la práctica profesional, el dominio del comportamiento
isostático es el mecanismo intelectual con el que el proyectista evalúa
con rapidez la viabilidad de una forma estructural antes de someterla a
modelos computacionales complejos. Por ende, las estructuras
estáticamente determinadas poseen la propiedad intrínseca de no
desarrollar esfuerzos secundarios ante asentamientos diferenciales en
los apoyos o dilataciones térmicas, convirtiéndolas en soluciones
estratégicas para contextos geotécnicos complejos o juntas de dilatación.
En conclusión, el análisis de estructuras estáticamente
determinadas conserva una vigencia incalculable en la era del modelado
96
numérico por elementos nitos. Más allá de calcular reacciones e
hiperestáticas internas, este libro rearma que la intuición estructural
no surge del software, sino de la comprensión profunda de cómo la
forma, el espacio y las fuerzas dialogan a través del equilibrio estricto.
La teoría avanzada aquí expuesta entrega las herramientas formales
para que el investigador y el diseñador dominen ese diálogo con
precisión matemática y juicio crítico.
97
Bibliografía
Affan, A., and Calladine, C. R. (1989). Initial Bar Tensions in Pin-
Jointed Assemblies. International Journal of Space Structures. 4(1), 1-
16. https://doi.org/10.1177/026635118900400101
Agredo Chávez, A. M., Sarmiento Nova, S. J., & Viviescas Jaime, A.
(2016). Evaluación de la rigidez a exión de puentes de viga y losa en
concreto presforzado a partir de pruebas de carga. Caso de estudio:
Puente La Parroquia vía La Renta - San Vicente de Chucurí. Revista UIS
Ingenierías, 15(2), 145-159.
https://www.redalyc.org/pdf/5537/553776590010.pdf
Alessandrini, M., Biagetti, G., Crippa, P., Falaschetti, L., Manoni, L., &
Turchetti, C. (2021). Singular Value Decomposition in Embedded
Systems Based on ARM Cortex-M Architecture. Electronics, 10(1), 34.
https://doi.org/10.3390/electronics10010034
Blanco Claraco, J. L., González Herrera, A., y García-Manrique, J. M.
(2012). Análisis estático de estructuras por el método matricial. Málaga:
Universidad de Málaga.
https://ingmec.ual.es/~jlblanco/papers/blanco2012calculo_matricial_e
structuras.pdf
Brenner, S. (1967). Endomorphism Algebras of Vector Spaces with
Distinguished Sets of Subspaces. Journal of Algebra, 6, 100-114.
https://acortar.link/0bG5T1
98
Brenner, S. and Butler, M.C.R. (1965). Endomorphism Rings of Vector
Spaces and Torsion Free Abelian Groups. Journal of the London
Mathematical Society, 40, 183-187. https://doi.org/10.1112/jlms/s1-
40.1.183
Chowdhury, I. R., & Summerscales, J. (2024). Tejidos para el refuerzo
de materiales compuestos: una revisión. Journal of Composites
Science, 8 (7), 280. https://doi.org/10.3390/jcs8070280
David Müzel, S., Bonhin, E. P., Guimarães, N. M., & Guidi, E. S. (2020).
Application of the Finite Element Method in the Analysis of Composite
Materials: A Review. Polymers, 12(4), 818.
https://doi.org/10.3390/polym12040818
Demmel, J. W. (1997). Other Titles in Applied Mathematics Applied
Numerical Linear Algebra. Philadelphia: Society for Industrial and
Applied Mathematics (SIAM).
https://doi.org/10.1137/1.9781611971446
Dugas, M., and Göbel, R. (1997). Endomorphism rings ofB 2-groups of
innite rank. Isr. J. Math., 101, 141–156.
https://doi.org/10.1007/BF02760926
Escalante Espinoza, N. J., Osorio Rojas, E. A., Quillos Ruiz, S. A., Tume
Ramirez, H., Sparrow Alamo, E. G., López Carranza, A. R., & Calderón
Rodríguez, L. C. (2026). Análisis estructural estáticamente
determinadas. Colonia del Sacramento: Editorial Mar Caribe.
https://doi.org/10.17613/k5chg-z9q46
99
Esquivel Cárdenas, J. A., Ruiz Tello, J., & Rincón Maltos, G. (2014).
Metodología para encontrar Matrices de Transformación Homogénea
para la Cinemática de Mecanismos con Matlab y VRML. PROSPECTIVA,
12(1), 64-69. https://www.redalyc.org/pdf/4962/496250639008.pdf
Fabré Sentile, C., Sánchez Noa, M., González Fernández, V., Wellesley-
Bourke Funcasta, J., Veliz Marrero, L. M., García Domínguez, M. E., &
Rodríguez Delgado, T. (2007). Análisis de pérdida de estabilidad no
lineal de elementos estructurales sometidos a compresión. Ingeniería
Mecánica, 10(1), 45-50.
https://www.redalyc.org/articulo.oa?id=225117649006
Fallas-Monge, J. J., Chavarría-Molina, J., & Soto-Quiros,
P. (2021). Descomposición en valores singulares de una matriz: Un
repaso por los fundamentos teóricos y sus aplicaciones en el
procesamiento de imágenes. Investigacion Operacional, 42(2), 148-173.
https://orion.tec.ac.cr/es/publications/descomposici%C3%B3n-en-
valores-singulares-de-una-matriz-un-repaso-por-/
Farias, J. G., De Pieri, E., & Martins, D. (2024). A Review on the
Applications of Dual Quaternions. Machines, 12(6), 402.
https://doi.org/10.3390/machines12060402
Garcés, F., Genatios, C., Lafuente, M., & Mebarki, A.. (2002). Reajuste
de matrices de rigidez y exibilidad para sistemas aporticados. Boletín
Técnico , 40(3), 19-26.
http://ve.scielo.org/scielo.php?script=sci_arttext&pid=S0376-
723X2002000300002&lng=en&tlng=es.
100
Lizarraga, I. (1990). Estructuras isostáticas. Ciudad de México: McGraw-
Hill
Luévanos Rojas, A. (2013). Método de aproximaciones sucesivas para
vigas estáticamente indeterminadas incluyendo una nueva variable.
Revista de Arquitectura e Ingeniería, 7(1), 1-14.
https://www.redalyc.org/articulo.oa?id=193927502002
Manglik, V. K. (2013). Elements of Mechanical Engineering. New Delhi:
PHI Learning Private Limited
Manoli, C. K., Papatzani, S., & Mouzakis, D. E. (2022). Exploring the
Limits of Euler–Bernoulli Theory in Micromechanics. Axioms, 11(3), 142.
https://doi.org/10.3390/axioms11030142
Medina, G., & Zavadskij, A.G. (2004). The four subspace problem: An
elementary solution. Linear Algebra and its Applications, 392, 11-23.
https://doi.org/10.1016/j.laa.2004.05.010
Monleón, S., y Lázaro, C. (2011). El concepto de hiperviga como
generalización de los modelos unidimensionales en la estática de piezas
alargadas. Revista Internacional de Métodos Numéricos para Cálculo y
Diseño en Ingeniería, 27(4), 294-303. https://www.elsevier.es/es-revista-
revista-internacional-metodos-numericos-calculo-338-pdf-download-
S0213131511000162 18
Morales Robles, A. J., Ortiz Dominguez, M., Cardoso Legorreta, E., &
Cruz Avilés, A. (2024). Comparación de dos métodos para el análisis
estructural de un eje estáticamente indeterminado. Ingenio Y Conciencia
101
Boletín Cientíco De La Escuela Superior Ciudad Sahagún, 11(22), 53–68.
https://doi.org/10.29057/escs.v11i22.12589
Pan, C., Han, Y., & Lu, J. (2020). Diseño y optimización de estructuras
reticulares: una revisión. Applied Sciences , 10 (18), 6374.
https://doi.org/10.3390/app10186374
Parras Galán, L., Montes Tubío, M., García Guzmán, A., Entrenas
Angulo, J. A., & de Dios Palomares, R. (1985). Un método de
optimización de estructuras articuladas estáticamente
determinadas. Informes De La Construcción, 37(372), 65–73.
https://doi.org/10.3989/ic.1985.v37.i372.1857
Pérez Montiel, H. (2015). Física general. Ciudad de México: Grupo
Editorial Patria
Pierre Sanon, S. (2022). Endomorphism Rings of Ordinary Abelian
Varieties. [Tesis de doctorado]. Technische Universität Kaiserslautern.
https://d-nb.info/1255183004/34
Qian, C., Kaminer, I. y Chen, H. (2025). Una guía para metamateriales
inteligentes e inteligencia de metamateriales. Nat Commun, 16, 1154.
https://doi.org/10.1038/s41467-025-56122-3
Ramalingam, K. K. (2009). Handbook of Mechanical Engineering Terms.
New Delhi: New Age International.
Rocks, J.W., and Mehta, P. (2024). Integrating local energetics into
Maxwell-Calladine constraint counting to design mechanical
102
metamaterials. Phys. Rev. E., 110, 025002.
https://doi.org/10.1103/PhysRevE.110.025002
Rodríguez Calderón, W., Rojas Aguero, R., Yépez Aguirre, J., &
Pallares - Muñoz, M. (2017). Análisis de sensibilidad y de estabilidad
numérica en el cálculo de factores de intensidad de tensiones en un caso
de mecánica de fractura. Revista UIS Ingenierías, 16(2), 151–160.
https://doi.org/10.18273/revuin.v16n2-2017014
Sánchez Beitia, S. (2021). Estabilidad e isostaticidad en la arquitectura.
Madrid: Ediciones Díaz de Santos.
https://www.editdiazdesantos.com/wwwdat/pdf/9788490522752.pdf
Suresh Kumar, N., Naidu, K. C. B., Banerjee, P., Anil Babu, T., &
Venkata Shiva Reddy, B. (2021). A Review on Metamaterials for Device
Applications. Crystals, 11(5), 518.
https://doi.org/10.3390/cryst11050518
Vargas Silva, G. (2025). Teoría de Estructuras en Ingeniería Mecánica
o cómo acercar la industria al mundo académico. Anales De Ingeniería
Mecánica, 1(24). https://doi.org/10.63450/aim.1.195.2025
Vergara Arredondo, F. (2023). Topological mechanics of hyperstatic
metamaterials. [Tesis de maestría]. Universidad de Chile.
https://repositorio.uchile.cl/handle/2250/194920
Voutetaki, M. E. (2024). Application and Assessment of an Experiential
Deformation Approach as a Didactive Tool of Truss Structures in
103
Architectural Engineering. Education Sciences, 14(4), 354.
https://doi.org/10.3390/educsci14040354
Wang, H., Wu, S., & Shao, Z. (2021). Analytical Compliance Equations
of Generalized Elliptical-Arc-Beam Spherical Flexure Hinges for 3D
Elliptical Vibration-Assisted Cutting Mechanisms. Materials, 14(20),
5928. https://doi.org/10.3390/ma14205928
Yılmaz, Y. (2026). On Quasilinear Algebra of Linear Interval Equations
and Interval Cramer’s Rule. Mathematics, 14(6), 1018.
https://doi.org/10.3390/math14061018
Zhang, P., Toman, J., Yu, Y., Biyikli, E., Kirca, M., Chmielus, M., and To,
A. C. (2015). Efcient Design-Optimization of Variable-Density
Hexagonal Cellular Structure by Additive Manufacturing: Theory and
Validation. ASME. J. Manuf. Sci. Eng., 137(2),
021004. https://doi.org/10.1115/1.4028724
Zhoufei, L., Peng Jin, M. L., Chengmeng, W., Pengfei Zhuang, P.
T., Jian-Hua Jiang, F. M., Jiping, H. (2025). Topology in Thermal,
Particle, and Plasma Diffusion Metamaterials. Chem. Rev., 24, 125(18),
8655–8730. https://doi.org/10.1021/acs.chemrev.4c00912
Zienkiewicz, O. C. (1977). The nite element method. New York: McGraw
Hill
104
De esta edición de Teoría avanzada de estructuras estáticamente
determinadas, se terminó de editar en la ciudad de Colonia del
Sacramento en la República Oriental del Uruguay el 15 de abril
de 2026
105
Editorial Mar Caribe
Teoría avanzada de estructuras estáticamente
determinadas
ISBN: 978-9915-698-96-0
