Espectros de Respuesta de Diseño
Esta gráfica compara la aceleración espectral máxima (adimensional $A_d$, en g) a la que se diseñaría la estructura bajo ambas normas según las condiciones geológicas elegidas. Las líneas verticales punteadas indican el período fundamental ($T_{fund}$) inicial estimado para los edificios.
Aceleración del Terreno ($a_g$)
Desplazamiento del Techo ($x_N$)
Curva de Histéresis del Primer Piso (Cortante de Piso vs Deriva de Piso)
Memoria de Cálculo (Tiempo Real)
Esta sección muestra las memorias de cálculo detalladas en tiempo real. Aquí puede auditar las inercias bruta y agrietada transformada, relaciones modulares, masas tributarias, coeficientes de rigidez y amplificaciones torsionales del modelo estructural dinámico. Al finalizar la simulación, se anexará automáticamente al final de esta sección un resumen detallado con el estado de las rótulas plásticas y derivas máximas de cada columna individual.
Sismos Recientes en la Región de Venezuela
Monitoreo en tiempo real de los últimos eventos sísmicos registrados en el territorio de Venezuela y sus alrededores. Los datos provienen del catálogo público global del USGS filtrados geográficamente (cobertura de un radio de 1000 km en torno al centro del país).
Deslinde de Responsabilidad (Disclaimer)
Esta aplicación web es una herramienta puramente educativa, interactiva y de divulgación académica. Los cálculos y simulaciones presentados se basan en formulaciones simplificadas y no consideran la complejidad geotécnica, topográfica o estructural detallada de un proyecto real. Bajo ninguna circunstancia esta herramienta sustituye un estudio de ingeniería sísmica, un análisis de vulnerabilidad ni un diseño estructural profesional realizado por profesionales licenciados bajo las normativas nacionales vigentes.
Simplificaciones y Consideraciones de Cálculo
Modelo de Edificio de Corte (Shear Building)
El edificio se modela como un sistema de múltiples grados de libertad (MDOF) donde toda la masa está concentrada en los niveles de losa y las vigas son infinitamente rígidas a flexión comparadas con las columnas. En el caso de "Secciones Físicas", se incorpora un coeficiente de flexibilidad de vigas (η) basado en la relación de rigidez viga-columna para modelar de manera simplificada la flexibilidad del nodo, omitiendo deformaciones axiales y rotacionales acopladas.
Análisis Dinámico de Tiempo-Historia No Lineal
Para resolver las ecuaciones de movimiento, se utiliza el método de integración numérica de Newmark (Beta = 0.25, Gamma = 0.50). El comportamiento no lineal de los entrepisos se simula mediante una ley constitutiva histerética bilineal (elastoplástica con endurecimiento por deformación), en lugar de un análisis de elementos finitos de plasticidad distribuida o modelos de fibras tridimensionales.
Espectro de Diseño vs. Aceleración del Terreno
El simulador opera bajo una doble metodología: la sacudida dinámica de la estructura y el cálculo de desplazamientos paso a paso están basados directamente en la señal del acelerograma del terreno (gráfico de aceleración). Sin embargo, la resistencia lateral ante la fluencia ($V_y$) y la rigidez inicial de diseño de los entrepisos se definen y calibran a partir del espectro de diseño de cada norma para el período fundamental calculado del edificio ($T_1$), acoplando el análisis espectral con la integración en el dominio del tiempo.
Generación Estocástica de la Onda Sísmica (Tiempo-Historia)
En cada simulación o cambio de configuración, el simulador genera un nuevo acelerograma sintético estocástico único. Para esto se emplea la Transformada de Box-Muller para producir ruido blanco gaussiano, el cual luego pasa por un Filtro de Kanai-Tajimi para moldear su frecuencia dominante según el período espectral del suelo seleccionado ($T_{soil}$). Finalmente, se modula en amplitud usando una envolvente de Jennings para simular el inicio, fase fuerte y atenuación del sismo. Esto implica que la señal en el tiempo varía en cada ejecución, manteniendo constantes únicamente su contenido de frecuencia promedio y aceleración pico (PGA).
Relación entre Magnitud ($M_w$ 7.1 vs. 7.5) e Intensidad (PGA)
El simulador representa la diferencia de magnitud entre el Sismo 1 ($M_w$ 7.1) y el Sismo 2 ($M_w$ 7.5) a través del escalamiento de la aceleración pico del terreno (PGA) y el daño acumulado. En la física sísmica, un sismo de mayor magnitud libera significativamente más energía y genera aceleraciones de pico más altas en el sitio. Esto se modela aplicando un factor de amplificación de 1.5 en el PGA del segundo evento ($\text{PGA}_2 = 1.5 \times \text{PGA}_1$). Adicionalmente, el sismo sucesivo aprovecha el ablandamiento y daño estructural (índice de Park-Ang) acumulado durante el primer evento, simulando de manera muy realista cómo la pérdida de rigidez inicial hace al edificio mucho más vulnerable al embate de una sacudida de mayor magnitud.
Práctica Profesional: Suites de Acelerogramas (ASCE 7)
Debido a la naturaleza caótica y aleatoria de la fase de las ondas, los códigos de diseño internacionales modernos (como el ASCE 7 Capítulo 16 o ACI 318) exigen que el análisis tiempo-historia no lineal se realice utilizando una suite de al menos 7 a 11 acelerogramas independientes. Si se corre un grupo de 7 o más sismos, la norma permite diseñar la estructura para la respuesta promedio obtenida. Si se corren menos (por ejemplo, 3), se debe diseñar obligatoriamente para la respuesta máxima absoluta, lo que resulta sumamente conservador.
¿Cómo generar una onda nueva en el simulador? Cada vez que presionas el botón "Reiniciar" (Reset) o cambias cualquier parámetro de la configuración, el sistema ejecuta de forma transparente el generador estocástico y produce un acelerograma completamente nuevo e independiente de la corrida anterior. Para evaluar el edificio de forma rigurosa, se aconseja realizar al menos 7 corridas consecutivas usando el botón Restablecer y promediar los resultados de deriva y daño.
Modelo de Daño de Park-Ang por Elemento y Degradación
El daño estructural se evalúa de manera continua en cada piso utilizando el Índice de Daño de Park-Ang. Para mayor precisión, el simulador rastrea individualmente la deriva en las columnas de esquina (las más solicitadas por torsión) para calcular el daño máximo local (Dmax). Este índice gobierna la formación de rótulas plásticas y el colapso del entrepiso cuando se alcanza un daño del 99% en cualquier columna individual, reflejando fallas locales críticas.
Inercias de Secciones Físicas y Masa de Diseño
Al activar "Secciones Físicas", la inercia de las columnas se calcula usando la sección transformada total (concreto + acero simétrico), y la de las vigas empleando la sección agrietada transformada en flexión simple. La masa de diseño se escala de forma realista según el área tributaria de la losa usando la relación m = mref × (0.3 + 0.7 × Área / 25), donde 0.3 representa el peso estructural fijo y 0.7 la losa y cargas variables asociadas.
Filosofía Columna Fuerte - Viga Débil (SCWB)
Cuando se activan las secciones personalizadas, el sistema calcula los momentos nominales de vigas (Mn,vig) y columnas (Mn,col) en cada nudo. Si se cumple el requisito de diseño dúctil (ΣMn,col ≥ 1.2 ΣMn,vig), las rótulas plásticas se forman preferentemente en las vigas, protegiendo las columnas y reduciendo la degradación de rigidez de entrepiso en un 60%. Si no se cumple, el daño y la fluencia ocurren directamente en las columnas, induciendo un mecanismo de piso blando frágil.
Efecto de Torsión en Planta y Excentricidad
Se simula la torsión acoplando el giro de la losa (θz) al desplazamiento traslacional en función de la excentricidad torsional elegida (ea). Este giro amplifica geométricamente la deriva en las columnas perimetrales más alejadas del centroide, lo que a su vez se traduce en un incremento proporcional del daño y la aceleración del colapso en el resolvedor físico.
Independencia de Ondas en Sismos Sucesivos
Al simular sismos sucesivos (doble sismo), por facilidad de cálculo se asume que ambos eventos son independientes. No se considera la superposición o interferencia destructiva/constructiva de las ondas sísmicas en el suelo de fundación; el segundo sismo actúa directamente sobre la estructura en su estado de rigidez y resistencia remanentes degradadas por el primer sismo.
Simulación de Evacuación de Personas y Vulnerabilidad
El simulador incluye un agente que modela el comportamiento y escape de 4 personas que inician en el techo (último piso) de cada edificio. El sistema evalúa dos fases: un retardo inicial de toma de decisión de 10.0 segundos ante el sismo activo (omitido en el modo 'después del sismo', el cual inicia a los 30.0 segundos cuando cesa el sismo 1), y un descenso constante a razón de 8.0 segundos por nivel. Las personas se acoplan cinemáticamente al balanceo lateral e inercial de la losa correspondiente a su ubicación. Si el resolvedor detecta el colapso del edificio, los ocupantes por encima de la losa fallada quedan atrapados y caen físicamente junto con los escombros de la estructura.
Vector de Fuerza Cortante Basal (Corte Basal 3D)
Representa la fuerza cortante total en la base del edificio (V_0 = sum F_i) mediante un vector 3D dinámico. La longitud de la flecha se escala en tiempo real en proporción al peso total del edificio (una longitud de 3.5 metros representa un corte basal equivalente al 30% del peso reactivo de la estructura), y su sentido indica la dirección en que empuja la inercia del sismo en cada instante del resolvedor físico.
Formación de Rótulas Plásticas
Visualiza las articulaciones plásticas (zonas de fluencia por flexión) en los extremos superiores e inferiores de cada columna de forma individual. Si la deriva local de la columna supera el límite elástico (1.5%), se dibuja un anillo amarillo. Si la deformación progresa y supera la zona de seguridad estructural (1.8%), el anillo cambia a color rojo de peligro, mostrando cómo se concentra la ductilidad y el daño inelástico a lo largo del sismo.
Recomendaciones para el Cálculo de la Excentricidad Torsional
La excentricidad torsional es la separación entre el Centro de Masas (CM) y el Centro de Rigidez (CR) de un entrepiso. Para ingenieros en formación o sin suficiente experiencia práctica, se sugieren los siguientes criterios de estimación y cálculo bajo normativas sismorresistentes (como la COVENIN 1756):
1. Excentricidad Estática o Teórica ($e_s$)
Es la distancia física real entre el baricentro de las cargas de la losa (CM, gobernado por la distribución de tabiques, losas, equipos) y el centro de cortante de los elementos verticales de soporte (CR, gobernado por la ubicación y sección de columnas y muros). Se calcula en cada nivel de entrepiso mediante:
x_CM = Σ(m_i · x_i) / Σm_i |
x_CR = Σ(K_y,j · x_j) / ΣK_y,j
e_s = |x_CM - x_CR|
2. Excentricidad Accidental ($e_a$) obligatoria
Debido a la incertidumbre en la posición real de las cargas vivas, variación del módulo de elasticidad del concreto, y efectos torsionales dinámicos no modelados, las normas exigen considerar una excentricidad accidental mínima. En la norma venezolana **COVENIN 1756**, se fija típicamente en:
e_a = 0.06 · L (Norma 2001) | e_a = 0.05 · L (Norma 2019)
Donde L es la dimensión de la planta del edificio medida en dirección perpendicular a la acción sísmica aplicada. Para un edificio simétrico ($e_s = 0$), el valor mínimo de diseño a introducir es esta excentricidad accidental (ej. 5% a 6%).
3. Excentricidad de Diseño ($e_d$) para Análisis Estático
Al realizar análisis estáticos equivalentes, la excentricidad teórica se amplifica dinámicamente para proteger los elementos extremos. Se calculan dos hipótesis de diseño en cada eje:
e_d1 = 1.5 · e_s + e_a (Amplifica el torque y castiga
el extremo alejado)
e_d2 = 0.5 · e_s - e_a (Reduce el torque para revisar
el lado menos solicitado)
4. Recomendaciones e Intuición Práctica
Evite irregularidad torsional: Si la excentricidad calculada supera el 10% del ancho del edificio (0.10), el edificio se comportará de manera muy inestable en planta, concentrando rotación e induciendo derivas severas en las columnas de esquina (efecto látigo torsional). Lo ideal es redistribuir las columnas o añadir pantallas (muros de cortante) para acercar el CR al CM y reducir la excentricidad a menos del 5%.
Limitación del Modelo frente a una Estructura Real
Es crucial comprender que el simulador es generalmente más optimista (menos estricto) que el comportamiento de un edificio real ante el mismo terremoto. Esto se debe a que el modelo simplificado omite modos de falla críticos que reducen la capacidad de supervivencia de la estructura:
Oclusión de Fallas Frágiles (Cortante y Adherencia)
El simulador asume una ductilidad ideal continua donde los elementos fluyen y se deforman cíclicamente hasta derivas altas (3.5% o 4.5%). En la realidad, si los estribos transversales no tienen el detallado moderno o están muy espaciados, se producirá un pandeo del acero o falla por cortante frágil a derivas mucho menores (entre 1.5% y 2.0%), colapsando la estructura de forma súbita.
Efectos Bidireccionales y Sismo Vertical
El simulador modela las fuerzas sísmicas en un único eje horizontal unidireccional. Un terremoto real actúa en tres dimensiones ortogonales simultáneas, induciendo flexión esviada bidireccional y oscilaciones dinámicas de la carga axial por sismo vertical, disminuyendo drásticamente la capacidad resistente real de las columnas.
Torsión en Planta e Irregularidades
Aunque el simulador ahora incorpora una torsión estática-dinámica acoplada (que rota las losas y castiga las columnas extremas), en la realidad la torsión se acopla dinámicamente de forma tridimensional para cada modo de vibración, amplificando el daño de manera no lineal y asimétrica de acuerdo a las irregularidades geométricas específicas de la estructura.
Diferencias y Consideraciones por Normativa
COVENIN 1756:2001 (Norma Tradicional)
- Espectro de Diseño: Basado en factores de aceleración de zona tradicional ($A_0 = 0.30$ para Caracas, Zona 5) y formas espectrales rígidas (suelos S1, S2, S3 o S4).
- Límite de Deriva de Colapso: Establecido de forma simplificada en un 3.5% para el análisis dinámico no lineal en el simulador.
- Severidad de Degradación: Mayor susceptibilidad a la pérdida de rigidez y resistencia ante deformaciones plásticas cíclicas, emulando los detalles de confinamiento menos estrictos del siglo pasado.
COVENIN 1756:2019 (Actualización y Propuestas)
- Espectro de Diseño: Emplea una formulación moderna basada en mapas de microzonificación e intensidades de aceleración espectral para períodos cortos ($A_0$) y de 1 segundo ($A_1$), permitiendo una mayor definición del tipo de sitio.
- Límite de Deriva de Colapso: Ampliado a un 4.5%, debido a los mayores requerimientos de ductilidad exigidos en las normativas modernas.
- Severidad de Degradación: Reducida a la mitad en la simulación no lineal, representando la excelente capacidad de confinamiento y retención de resistencia de los marcos de concreto armado con detallado dúctil moderno (ND3).
¿Tienes sugerencias o comentarios?
Este simulador es de código abierto. Si deseas proponer mejoras técnicas, reportar alguna inconsistencia en las fórmulas o sugerir nuevas funciones para las normas de sismo de Venezuela, por favor abre una sugerencia en mi repositorio de GitHub.
Registrar sugerencia en GitHub Issues