1. Einführung
Eine genaue Bewertung der Windeinwirkungen auf Konstruktionen hängt nicht nur vom verwendeten Lastmodell ab, sondern vor allem davon, wie die einströmende atmosphärische Grenzschicht (ABL) dargestellt wird. In der praktischen Bauplanung ist Wind ein äußerst instabiles, turbulentes und räumlich korreliertes Phänomen. Die meisten normbasierten Ansätze, darunter auch das Standard-Windmodell des Eurocodes, reduzieren diese Komplexität jedoch auf äquivalente statische Größen, um Einfachheit, Robustheit und konservative Sicherheitsmargen zu gewährleisten.
Bei Windsimulationen auf Basis der numerischen Strömungsmechanik (CFD) ist die Definition der Anströmrandbedingungen ein entscheidender Faktor für die Genauigkeit der vorhergesagten Oberflächendrücke, Spitzenlasten und dynamischen Effekte.
Die hier beschriebene Methode implementiert einen zufälligen, räumlich korrelierten turbulenten Einströmungsgenerator auf einer vertikalen Anströmebene (YZ-Ebene) und ist vorgesehen für:
- Instationäre CFD-Windsimulationen
- Bewertung von lastrelevanten Druck- und Kraftwerten
- Validierung unter windkanalähnlichen Anströmbedingungen
- Fortgeschrittene Anwendungen, die über den konservativen Anwendungsbereich des Eurocodes hinausgehen
Dieser Ansatz geht deutlich über das Standard-Windmodell des Eurocodes (EC) hinaus, das in erster Linie für die Ableitung statischer Bemessungslasten und nicht für die zeitaufgelöste Strömungsphysik konzipiert ist.
2. Überblick über die implementierte Methode
Die Einström-Generierung besteht aus fünf gekoppelten Komponenten:
2.1 Räumlich diskretisierte Anströmebene
- Anströmung definiert auf einer YZ-Ebene bei konstantem 𝑥=𝑥0
- Regelmäßiges Raster: 𝑛𝑦 × 𝑛𝑧 Punkte (z.B. 20 × 100 → 2.000 Proben)
- Jeder Rasterpunkt repräsentiert ein zeitabhängiges Geschwindigkeitssignal
- Es wird keine Mittelwertbildung über die Spannweite angewandt (wichtig für die Erhaltung von Fluktuationen)
Dadurch lassen sich vertikale und seitliche Kohärenzeffekte physikalisch darstellen.
2.2 Mittleres Windprofil (Potenzgesetz, physikalisch erzwungen)
Das mittlere Windprofil wird anhand einer Potenzgesetzformulierung definiert, einer empirischen ingenieurtechnischen Näherung, die den Anstieg der Windgeschwindigkeit mit der Höhe innerhalb der atmosphärischen Grenzschicht abbildet. Das Modell setzt die Geschwindigkeit in Beziehung zu einer Referenzwindbedingung und einem geländespezifischen Exponenten. Obwohl der Potenzgesetzansatz einfacher ist als die aus der Ähnlichkeitstheorie abgeleitete logarithmische Grenzschichtformulierung, bietet er eine praktische und weit verbreitete Definition der Einströmung für ingenieurtechnische Anwendungen. Die Formulierung gewährleistet eine konsistente vertikale Geschwindigkeitsverteilung, bezogen auf eine festgelegte Höhe.
|
U(z) |
Mittlere Windgeschwindigkeit in Höhe z über Geländeoberkante |
|
Uref |
Referenzwindgeschwindigkeit |
|
z |
Höhe über Grund |
|
𝑧ref |
Bezugshöhe |
|
α |
Potenzgesetz-Exponent (Gelände-Exponent) |
2.3 Turbulenzmodell: Kaimal-Spektrum
Die Turbulenz an jedem Anströmpunkt wird unter Verwendung des Kaimal-Spektrums generiert, das das Verhalten der atmosphärischen Grenzschicht abbildet und festlegt, wie sich die Turbulenzenergie über die Frequenzen verteilt. Dieser Ansatz erfasst realistische frequenzabhängige Eigenschaften von Windschwankungen, anstatt von einem konstanten Verhalten auszugehen. Das Modell berücksichtigt die Turbulenzintensität und charakteristische Längenskalen, um die Struktur der Geschwindigkeitsschwankungen zu beschreiben. Dadurch enthält die generierte Turbulenz einen physikalisch konsistenten Energiegehalt, der mit der Frequenz variiert. An jedem Anströmpunkt wird die Turbulenz unter Verwendung des Kaimal-Autospektrums generiert, in Übereinstimmung mit der Theorie der atmosphärischen Grenzschicht [1]:
|
Su(f) |
Spektrale Leistungsdichte der Längsgeschwindigkeit bei Frequenz |
|
σu |
Standardabweichung der Längsgeschwindigkeitsschwankungen |
|
𝐼𝑢 |
Turbulenzintensität |
|
Lu |
Integrale Längenskala in Längsrichtung |
Dadurch ergibt sich ein frequenzabhängiger Energiegehalt, der im statischen Ansatz des Eurocodes nicht berücksichtigt wird.
2.4 Räumliche Kohärenz und Kreuzspektralmatrix
Es wird ein räumliches Kohärenzmodell angewendet, um eine physikalisch realistische Korrelation zwischen den Anströmpunkten sicherzustellen, indem beschrieben wird, wie die Turbulenzähnlichkeit mit der Entfernung und der Frequenz abnimmt. Diese Kohärenzfunktion dient dazu, die Beziehung zwischen Signalen an verschiedenen Orten zu definieren, wodurch räumlich korrelierte Windschwankungen widergespiegelt werden. Auf der Grundlage dieser Kohärenz wird eine vollständige Kreuzspektraldichtematrix konstruiert, die einzelne Spektren mit ihrer räumlichen Korrelation kombiniert. Die resultierende Matrix repräsentiert sowohl die Energieverteilung als auch die räumliche Kopplung von Geschwindigkeitsschwankungen über den Einlass hinweg. Um eine physikalisch realistische Korrelation zwischen den Anströmpunkten zu gewährleisten, wird ein entfernungsabhängiges Kohärenzmodell angewendet:
|
γij(f) |
Frequency-dependent coherence function between inlet points i , j and |
|
f |
Frequenz |
|
dij |
Euklidischer Abstand zwischen Anströmpunkten |
|
𝐶coh |
Kohärenzabnahmekoeffizient |
Daraus ergibt sich eine vollständige Kreuzspektraldichtematrix:
|
Sij(f) |
Kreuzspektraldichte zwischen Geschwindigkeitssignalen an Anströmpunkten i und j |
|
Sii(f), Sjj(f) |
Autospektraldichten an Punkten i und j |
Wesentliche Merkmale:
- Vollständig räumlich korrelierte Turbulenz
- Frequenzabhängiger Korrelationsabfall
- Erfasst realistische Böenstrukturen über dem Einlass
2.5 Robuste Cholesky-Zerlegung (numerische Stabilität)
Eine robuste Cholesky-Zerlegung wird angewendet, um die numerische Stabilität zu gewährleisten, wenn die Kreuzspektralmatrix schlecht konditioniert wird, insbesondere in Bodennähe. Die Matrix wird zunächst symmetrisiert, um Konsistenz zu gewährleisten, gefolgt von einer adaptiven diagonalen Belastung zur Verbesserung der Stabilität. Bei Bedarf wird eine automatische Wiederholungsstrategie mit zunehmender Regularisierung verwendet, um die positive Definierbarkeit sicherzustellen. Dieser Ansatz ermöglicht die stabile und zuverlässige Generierung korrelierter Zeitreihendaten. Da die Kreuzspektralmatrix (insbesondere in Bodennähe) schlecht konditioniert werden kann, wird eine robuste Cholesky-Faktorisierung verwendet:
- Hermitische Durchsetzung:
- Adaptive diagonale Belastung
- Automatische Wiederholungsstrategie mit zunehmender Regularisierung
Dies gewährleistet positive Definitheit und eine stabile Generierung von Zeitreihen.
2.6 Konstruieren von Zeitreihen (inverse FFT)
Es werden zufällige Phasen eingefügt, um Turbulenzen im Signal darzustellen, während das Spektrum skaliert und strukturiert wird, um ein physikalisch realistisches, reellwertiges Ergebnis zu gewährleisten. Anschließend wird eine inverse FFT angewendet, um die Informationen aus dem Frequenzbereich in ein zeitabhängiges Geschwindigkeitssignal umzuwandeln. Das Endergebnis ist ein Geschwindigkeitsfeld, das sich aus einer mittleren Strömung in Kombination mit fluktuierenden Komponenten zusammensetzt.
- Pro Frequenz werden Zufallsphasen angewendet
- Einseitiges Spektrum, skaliert mit sqrt(2Δ𝑓)
- Hermitische Symmetrie erzwungen
- Inverse FFT → zeitaufgelöstes Geschwindigkeitssignal an jedem Anströmpunkt
Ergebnis:
3. Visualisierung und Diagnostik
Die Methode ermöglicht:
- Sofortige vertikale Profile 𝑈(𝑧,𝑡o)
- Vergleiche mehrerer Probenlinien (ohne Mittelungsartefakte)
- Vollständige Animationen der Konturen in der Anströmebene 𝑈(𝑦,𝑧,𝑡)
Dadurch lässt sich die Strömungsphysik im Gegensatz zu statischen Profilen auf EC-Basis direkt beobachten.
4. Vergleich mit dem Standard-Windmodell des Eurocodes
Tabelle 1 fasst die grundlegenden konzeptionellen und methodischen Unterschiede zwischen dem Standard-Windmodell des Eurocodes (EC) und der derzeitigen zufälligen CFD-Einstromgenerierungsmethode zusammen. Obwohl beide Ansätze darauf abzielen, Windeinwirkungen auf Konstruktionen darzustellen, basieren sie auf unterschiedlichen Modellierungsphilosophien und dienen unterschiedlichen ingenieurtechnischen Zwecken.
Tabelle 1: Vergleich mit dem Standard-Eurocode-Windmodell
| Aspekt | Eurocode (EC) | Aktuelle CFD-Einströmungsmethode |
|---|---|---|
| Natur | Statisch/quasistatisch | Vollständig instationär |
| Turbulenz | Implizit über Faktoren | Explizite Zeitreihen |
| Räumliche Korrelation | Nicht aufgelöst | Vollständig aufgelöst |
| Frequenzinhalt | Nicht vorhanden | Kaimal-Spektrum |
| Böenstruktur | Äquivalente statische Böe | Physikalisch entwickelnd |
| Druckspitzen | Empirische Faktoren | Aus Strömung entstehend |
| Lastpfade | Konservative Umhüllende | Physikbasiert |
| Eignung | Normkonformität | Erweiterte Analyse & Validierung |
5. Hauptunterschiede in der Interpretation
Der Eurocode-Ansatz zielt in erster Linie darauf ab, sichere Bemessungslasten zu gewährleisten, indem Turbulenzeffekte durch Teilsicherheitsbeiwerte berücksichtigt werden; er liefert jedoch keine detaillierten Erkenntnisse über den zeitlichen Ablauf oder die Mechanismen der Lastentwicklung. Im Gegensatz dazu ermöglicht die vorliegende Methode ein tieferes physikalisches Verständnis, indem sie aufzeigt, wann, wo und warum Spitzendrücke auftreten, einschließlich Windböen, durch Ablösung verursachter Druckspitzen und räumlich korrelierter Belastungseffekte. Dies macht sie besonders wertvoll für Leichtbaukonstruktionen, Fassaden- und Verkleidungselemente, Membransysteme, auf Dächern montierte Anlagen sowie für Untersuchungen zur dynamischen Strukturantwort:
📘Eurocode
- Ausgelegt für die Erzeugung sicherer Bemessungslasten
- Turbulenzeffekte sind in den Teilsicherheitsbeiwerten enthalten
- Keine Angaben dazu, wann oder wie Lasten auftreten
🌬️Aktuelle Methode
- Ermittelt, wann, wo und warum Druckspitzen entstehen
- Böendurchdringung
- Durch Ablösung verursachte Druckspitzen
- Korrelierte Belastung über große Flächen
⭐Besonders wichtig für:
- Leichtbaukonstruktionen
- Verkleidung und Fassadenelemente
- Membranen und Dachaufbauten
- Untersuchungen zur dynamischen Reaktion
6. Technische Implikationen
Die vorgeschlagene Methode zielt nicht darauf ab, den Eurocode für die normgerechte Bemessung zu ersetzen; vielmehr ist sie als ergänzender Ansatz zu verstehen, der zusätzliche physikalische Einblicke liefert, die über die typischerweise in den Eurocode-Verfahren enthaltenen Annahmen hinausgehen. Sie eignet sich besonders für Sensitivitätsanalysen, Forschungs- und Validierungszwecke, Optimierungsprozesse sowie zur Erklärung von Abweichungen, die zwischen Eurocode-basierten Ergebnissen und CFD-Simulationen auftreten können. In der praktischen Ingenieuranwendung hilft dieser Ansatz oft dabei, zu klären, warum sich CFD-basierte Druckwerte von den Eurocode-Vorhersagen unterscheiden, insbesondere hinsichtlich lokaler Spitzendrücke und Effekte von instationärer Strömung, ohne dass dies zwangsläufig auf einen Fehler in einer der beiden Methoden hindeuten muss.
7. Zusammenfassung
Die vorgestellte Methode zur Generierung von Einströmungen kombiniert eine physikalisch konsistente Turbulenzmodellierung mit einer räumlich kohärenten und zeitaufgelösten Darstellung des Windfeldes, gestützt durch eine robuste numerische Umsetzung. Im Vergleich zum Standardansatz des Eurocodes verlagert die Methode den analytischen Fokus von konservativen Lastumhüllenden hin zu einem physikalisch fundierten Verständnis des Strömungsverhaltens und ermöglicht so tiefere Einblicke und fortgeschrittenere windtechnische Analysen.