1. Introducción
La evaluación precisa de las acciones del viento sobre estructuras depende no solo del modelo de carga aplicado, sino fundamentalmente de cómo se representa la capa límite atmosférica (ABL) entrante. En el diseño de ingeniería práctico, el viento es un fenómeno altamente inestable, turbulento y correlacionado espacialmente. Sin embargo, la mayoría de los enfoques basados en códigos, incluido el modelo estándar de viento del Eurocódigo, reducen esta complejidad a cantidades estáticas equivalentes para garantizar simplicidad, robustez y márgenes de seguridad conservadores.
En las simulaciones de viento basadas en Dinámica de Fluidos Computacional (CFD), la definición de la condición de contorno de entrada es un factor decisivo para la exactitud de las presiones de superficie predichas, cargas máximas y efectos dinámicos.
El método documentado aquí implementa un generador de inflow turbulento aleatorio y correlacionado espacialmente en un plano de entrada vertical (plano YZ) y está destinado para:
- Simulaciones transitorias de viento por CFD
- Evaluación de presiones y fuerzas relevantes para cargas
- Validación frente a condiciones de inflow similares a túneles de viento
- Aplicaciones avanzadas más allá del alcance conservador del Eurocódigo
Este enfoque va significativamente más allá del modelo de viento del Eurocódigo (EC) estándar, que está diseñado principalmente para la derivación de cargas de diseño estáticas, no para la física del flujo resuelto en el tiempo.
2. Resumen del Método Implementado
La generación de inflow consta de cinco componentes acoplados:
2.1 Plano de Entrada Espacialmente Discretizado
- Entrada definida en un plano YZ a constante 𝑥=𝑥0
- Malla regular: puntos 𝑛𝑦 × 𝑛𝑧 (ej. 20 × 100 → 2,000 sondas)
- Cada punto de la malla representa una señal de velocidad dependiente del tiempo
- No se aplica promediación en el sentido del vano (importante para preservar fluctuaciones)
Esto permite que los efectos de coherencia vertical y lateral se representen físicamente.
2.2 Perfil de Viento Medio (Ley de Potencia, Impuesta Físicamente)
El perfil de viento medio se define utilizando una formulación de ley de potencia, una aproximación de ingeniería empírica que representa el aumento de la velocidad del viento con la altura dentro de la capa límite atmosférica. El modelo relaciona la velocidad con una condición de viento de referencia y un exponente dependiente del terreno. Aunque es más simple que la formulación logarítmica de la capa límite derivada de la teoría de similitud, el enfoque de ley de potencia proporciona una definición de inflow práctica y ampliamente utilizada para aplicaciones de ingeniería. La formulación asegura una distribución de velocidad vertical coherente referenciada a una altura especificada.
|
U(z) |
Velocidad media del viento a la altura z sobre el nivel del suelo |
|
Uref |
Velocidad básica del viento |
|
z |
Altura sobre el suelo |
|
𝑧ref |
Altura de referencia |
|
α |
Exponente de la ley de potencias (exponente del terreno) |
2.3 Modelo de Turbulencia: Espectro de Kaimal
La turbulencia en cada punto de entrada se genera utilizando el espectro de Kaimal, que representa el comportamiento de la capa límite atmosférica y define cómo se distribuye la energía turbulenta a través de las frecuencias. Este enfoque captura características realistas dependientes de frecuencia de las fluctuaciones del viento en lugar de asumir un comportamiento constante. El modelo refleja la intensidad de turbulencia y escalas de longitud características para describir la estructura de variaciones de velocidad. Como resultado, la turbulencia generada incluye contenido energético físicamente consistente que varía con la frecuencia. En cada punto de entrada, la turbulencia se genera usando el autospectro de Kaimal, consistente con la teoría de la capa límite atmosférica [1]:
|
Su(f) |
Densidad espectral de potencia de velocidad longitudinal a frecuencia |
|
σu |
Desviación estándar de las fluctuaciones de velocidad longitudinal |
|
𝐼𝑢 |
Intensidad de turbulencia |
|
Lu |
Escala de longitud integral |
Esto produce contenido energético dependiente de frecuencia, que está ausente en el enfoque estático del Eurocódigo.
2.4 Coherencia Espacial y Matriz Espectral de Cruce
Se aplica un modelo de coherencia espacial para asegurar una correlación físicamente realista entre los puntos de entrada describiendo cómo disminuye la similitud de la turbulencia con la distancia y la frecuencia. Esta función de coherencia se utiliza para definir la relación entre señales en diferentes ubicaciones, reflejando fluctuaciones de viento espacialmente correlacionadas. Basado en esta coherencia, se construye una matriz completa de densidad espectral de cruce, combinando espectros individuales con su correlación espacial. La matriz resultante representa tanto la distribución de energía como el acoplamiento espacial de las fluctuaciones de velocidad a través de la entrada. Para asegurar una correlación físicamente realista entre los puntos de entrada, se aplica un modelo de coherencia dependiente de la distancia:
|
γij(f) |
Función de coherencia dependiente de la frecuencia entre los puntos de entrada i, j y |
|
f |
Frecuencia |
|
dij |
Distancia euclidiana entre puntos de entrada |
|
𝐶coh |
Coeficiente de decaimiento de coherencia |
Esto resulta en una matriz completa de densidad espectral de cruce:
|
Sij(f) |
Densidad de cruce del espectro entre señales de velocidad en puntos de entrada i y j |
|
Sii(f), Sjj(f) |
Densidades auto-espectrales en los puntos i y j |
Características clave:
- Turbulencia completamente correlacionada espacialmente
- Decaimiento de la correlación dependiente de la frecuencia
- Captura estructuras de ráfagas realistas a través de la entrada
2.5 Descomposición de Cholesky Robusta (Estabilidad Numérica)
Se aplica una descomposición de Cholesky robusta para mantener la estabilidad numérica cuando la matriz espectral de cruce se vuelve mal condicionada, especialmente cerca del suelo. La matriz se simetriza primero para asegurar consistencia, seguida de una carga diagonal adaptativa para mejorar la estabilidad. Se utiliza una estrategia de reintento automático con aumento de regularización si es necesario, asegurando positividad definida. Este enfoque permite una generación estable y confiable de datos de series de tiempo correlacionados. Debido a que la matriz espectral de cruce puede volverse mal condicionada (especialmente cerca del suelo), se utiliza una factorización de Cholesky robusta:
- Cumplimiento Hermitiano:
- Carga diagonal adaptativa
- Estrategia de reintento automático con regulación creciente
Esto garantiza positividad definida y generación estable de series de tiempo.
2.6 Construcción de Series de Tiempo (FFT Inversa)
Se introducen fases aleatorias para representar la turbulencia dentro de la señal, mientras que el espectro se escala y estructura para asegurar un resultado realista físicamente real. Luego se aplica una FFT inversa para transformar la información en el dominio de frecuencia en una señal de velocidad dependiente del tiempo. El resultado final es un campo de velocidad compuesto por un flujo medio combinado con componentes fluctuantes.
- Fases aleatorias aplicadas por frecuencia
- Espectro unilateral escalado por sqrt(2Δ𝑓)
- Simetría Hermitiana aplicada
- FFT inversa → señal de velocidad resuelta en el tiempo en cada punto de entrada
Resultado:
3. Visualización y Diagnóstico
El método permite:
- Perfiles verticales instantáneos 𝑈(𝑧,𝑡o)
- Comparaciones de múltiples líneas de sondas (sin artefactos de promediación)
- Animaciones de contorno en todo el plano de entrada 𝑈(𝑦,𝑧,𝑡)
Esto hace que la física del flujo sea directamente observable, a diferencia de los perfiles estáticos basados en el EC.
4. Comparación con el Modelo Estándar de Viento del Eurocódigo
La Tabla 1 resume las diferencias conceptuales y metodológicas fundamentales entre el modelo estándar de viento del Eurocódigo (EC) y el método actual de generación de inflow aleatorio por CFD. Mientras que ambos enfoques buscan representar los efectos del viento sobre estructuras, se basan en filosofías de modelado diferentes y sirven a propósitos de ingeniería distintos.
Tabla 1: Comparación con el Modelo Estándar de Viento del Eurocódigo
| Aspecto | Eurocódigo (EC) | Método Actual de Inflow por CFD |
|---|---|---|
| Naturaleza | Estática / casi estática | Totalmente transitoria |
| Turbulencia | Implícita a través de factores | Series de tiempo explícitas |
| Correlación espacial | No resuelta | Totalmente resuelta |
| Contenido de frecuencia | No presente | Espectro de Kaimal |
| Estructura de ráfaga | Ráfaga estática equivalente | Evolución física |
| Picos de presión | Factores empíricos | Emergente del flujo |
| Rutas de carga | Envolvente conservadora | Basada en la física |
| Adecuación | Cumplimiento del código | Análisis avanzado y validación |
5. Diferencias Clave en la Interpretación
El enfoque del Eurocódigo está diseñado principalmente para asegurar cargas de diseño seguras incorporando efectos de turbulencia a través de factores de seguridad parciales, pero no proporciona un conocimiento detallado del momento o mecanismos de desarrollo de las cargas. En contraste, el método actual permite una comprensión física más profunda al resolver cuándo, dónde y por qué ocurren los picos de presión, incluyendo penetración de ráfagas, picos de presión inducidos por separación, y efectos de carga correlacionados espacialmente. Esto lo hace particularmente valioso para estructuras ligeras, elementos de la fachada y revestimiento, sistemas de membranas, equipos montados en techos y estudios que involucran respuesta estructural dinámica:
📘Eurocódigo
- Diseñado para producir cargas de diseño seguras
- Efectos de turbulencia incorporados en factores de seguridad parciales
- Sin información sobre cuándo o cómo ocurren las cargas
🌬️Método Actual
- Resuelve cuándo, dónde y por qué se forman los picos de presión
- Penetración de ráfagas
- Picos de presión inducidos por separación
- Carga correlacionada sobre áreas grandes
⭐Particularmente importante para:
- Estructuras ligeras
- Elementos de fachada y revestimiento
- Membranas y equipos en techos
- Estudios de respuesta dinámica
6. Implicaciones de Ingeniería
El método propuesto no pretende reemplazar al Eurocódigo para diseño reglamentario; en cambio, debe entenderse como un enfoque complementario que proporciona un conocimiento físico adicional más allá de las suposiciones típicamente incorporadas en los procedimientos del Eurocódigo. Es particularmente adecuado para estudios de sensibilidad, propósitos de investigación y validación, procesos de optimización, y para explicar discrepancias que pueden surgir entre resultados basados en el Eurocódigo y simulaciones por CFD. En aplicaciones de ingeniería práctica, este enfoque a menudo ayuda a aclarar por qué los valores de presión derivados por CFD difieren de las predicciones del Eurocódigo, especialmente en cuanto a picos locales de presión y efectos transitorios del flujo, sin necesariamente indicar un error en ninguna de las metodologías.
7. Resumen
El método de generación de inflow presentado introduce un modelado de turbulencia físicamente consistente combinado con coherencia espacial y representación de campo de viento resuelto en el tiempo, respaldado por una implementación numérica robusta. En comparación con el enfoque estándar del Eurocódigo, el método cambia la perspectiva analítica de los envolventes de carga conservadora hacia una comprensión basada en la física del comportamiento del flujo, permitiendo así un conocimiento más profundo y análisis de ingeniería del viento más avanzados.