建築物の風荷重解析
建築物の風荷重の理論基礎
概要
先生、建物まわりの風の解析って何を求めるんですか？
大きく3つの目的がある。(1)構造設計のための風荷重算定、(2)歩行者レベルの風環境評価（ペデストリアンコンフォート）、(3)自然換気計画だ。
超高層ビルでは風荷重が構造設計の支配的要因になる。建築基準法では風力係数を使って設計風圧を算定するが、複雑な建物形状や周辺建物との干渉がある場合はCFD解析が求められるんだ。
風洞試験の代わりにCFDを使うケースが増えてるんですよね。
その通り。ただし建築分野のCFDは風洞試験の完全代替ではなく、相互補完の関係だ。日本建築学会の「建築物荷重指針」でもCFD解析のガイドラインが整備されている。
支配方程式
建物まわりの風を記述する方程式はどうなりますか？
非圧縮性Navier-Stokes方程式が基本だ。建物まわりの風速はM < 0.3なので圧縮性は無視できる。
風圧係数はこう定義される。
ここで$p$は局所圧力、$p_\infty$は基準圧力、$V_H$は建物高さでの基準風速だ。
大気境界層の風速プロファイルはべき法則で表されることが多い。
ここで$\alpha$は地表面粗度に依存するべき指数だ。市街地で$\alpha \approx 0.25$--$0.35$、海上で$\alpha \approx 0.10$--$0.15$程度になる。
なるほど。入口境界条件として大気境界層プロファイルを与えるんですね。
乱流モデル
建築CFDで使われる乱流モデルを整理しよう。
| モデル | 特徴 | 建築風解析での適性 |
|---|---|---|
| 標準k-epsilon | 等方性乱流。計算コスト低 | 鈍頭物体の剥離を過小予測する傾向 |
| RNG k-epsilon | 渦度依存の粘性。剥離改善 | 角柱まわりの流れに改善効果 |
| SST k-omega | 壁面近傍の精度良好 | 建物表面の風圧分布に推奨 |
| LES (Smagorinsky) | 大スケール渦を直接解く | ピーク風圧・変動風圧に必須 |
| DES/DDES | RANS+LESハイブリッド | 実務的なコストで変動風圧を予測 |
k-epsilonだと建物の剥離を正しく予測できないんですか？
標準k-epsilonは鈍頭物体（角柱や直方体）の後方渦を過小予測する傾向がある。建物屋上角部の再付着長さが実験と合わないケースが多い。RNG k-epsilonや実現可能（Realizable）k-epsilonで改善されるが、ピーク風圧の予測にはLESが望ましいんだ。
歩行者レベル風環境
ビル風の評価基準ってどうなっていますか？
日本建築学会では風環境評価尺度を定めている。歩行者の高さ（地上1.5m）での風速を対象にする。
| ランク | 年間累積超過確率 | 環境の目安 |
|---|---|---|
| 1（良好） | 10m/s超が1%未満 | 住宅地・公園 |
| 2（許容） | 10m/s超が5%未満 | 一般市街地 |
| 3（やや不良） | 10m/s超が10%未満 | 商業地域 |
| 4（不良） | 10m/s超が10%以上 | 対策が必要 |
風環境の評価には年間の風向頻度分布も考慮するんですよね。
その通り。16風向（22.5度刻み）のCFDを実施し、対象地点のアメダスデータの風向頻度と組み合わせて年間の超過確率を算出するのが標準的な手法だ。
東京スカイツリーが「三角形を高さで回転」させた空力的理由
東京スカイツリーの断面は下部で正三角形ですが、高さとともに断面が少しずつ回転し、上部では円形に近づく設計になっています。これはただのデザインではなく、カルマン渦による共振（ビル風での揺れ）を抑制するための空力設計です。円柱形や単純な三角形断面は特定の風速でカルマン渦が同期し大きく振動しますが、高さで断面を変えることで渦が全高で同期しにくくなる。CFDと風洞実験を組み合わせた検証を経て採用されたこの工夫が、世界一高い自立式電波塔の安全性を支えています。
建築物の風荷重の数値計算手法
計算領域とメッシュ
建物まわりのCFDで計算領域はどのくらいの大きさにすればいいですか？
AIJ（日本建築学会）ガイドラインに基づく推奨値がある。
| パラメータ | 推奨値 | 備考 |
|---|---|---|
| 入口から建物まで | 5H以上 | Hは建物高さ |
| 建物から出口まで | 15H以上 | 後流発達のため |
| 側面まで | 5H以上 | 閉塞率5%以下 |
| 上面まで | 5H以上 | 閉塞率5%以下 |
| 閉塞率 | 3%以下推奨 | 建物断面/計算領域断面 |
閉塞率を低くする必要があるんですね。
そうだ。閉塞率が高いと人工的な加速効果が生じて風圧を過大評価する。3%以下が理想で、最大でも5%を超えないようにすべきだ。
入口境界条件
大気境界層の入口条件はどう設定するんですか？
べき法則またはログ法則のプロファイルを与える。乱流諸量も同時に指定する必要があるんだ。
速度プロファイル（ログ法則）:
乱流エネルギー:
散逸率:
ここで$u_*$は摩擦速度、$\kappa = 0.41$はカルマン定数、$z_0$は粗度長、$C_\mu = 0.09$だ。
$z_0$（粗度長）はどう決めるんですか？
地表面の粗度カテゴリに対応する値を使う。
| 地表面区分 | $z_0$ [m] | べき指数$\alpha$ | 例 |
|---|---|---|---|
| I（海上） | 0.0002--0.005 | 0.10 | 海岸、空港 |
| II（田園） | 0.01--0.05 | 0.15 | 田畑、低層住宅 |
| III（郊外） | 0.1--0.5 | 0.20 | 中層市街地 |
| IV（市街地） | 0.5--2.0 | 0.27 | 高層ビル群 |
メッシュ戦略
建物まわりのメッシュ生成のポイントを整理しよう。
- 建物面上: 最低10分割/辺（角部は細分化）
- 地表面近傍: $y^+ < 1$（壁面せん断応力の精度確保）
- 建物周辺リファインメント: 建物高さの2倍の範囲を細分化
- 後流域: 建物後方10Hまでは粗くしすぎない
- セル成長率: 1.2以下
SnappyHexMeshで建物まわりのメッシュを作ることが多いですよね。
OpenFOAMのsnappyHexMeshは建築CFDで広く使われている。STL形式で建物形状を読み込み、自動的に局所細分化とプリズム層追加ができる。STAR-CCM+のトリムメッシュも同様のアプローチで効率的だよ。
定常RANSとLESの使い分け
どんな場合にLESが必要になりますか？
使い分けの目安はこうだ。
| 目的 | 推奨手法 | 理由 |
|---|---|---|
| 平均風圧分布 | 定常RANS | 実務で十分な精度 |
| ピーク風圧 | LES/DES | 変動成分の予測が必要 |
| 歩行者風環境（平均） | 定常RANS | 16風向の効率的計算 |
| 渦励振評価 | LES | 渦放出周波数の予測 |
| 自然換気 | 非定常RANS/LES | 開口部の変動風圧が重要 |
平均的な風圧分布だけならRANSで十分なんですね。
そうだ。ただしRANSのピーク風圧予測は過小になる傾向がある。外装材の設計ではピーク風圧が必要になるから、その場合はLESかDESが推奨される。計算コストはRANSの50--100倍になるが、局所的なピーク負圧を正しく予測できるメリットは大きいよ。
建築CFDの「入口乱流プロファイル」はどこから来るのか
建物周りのCFDで最も地味だけど重要なのが「入口の大気境界層プロファイル設定」です。現実の風は地表摩擦で高さとともに増速するプロファイルを持ちますが、これを正しく設定しないと計算途中で崩壊してしまいます（インホモジニティ問題）。Richards-Hoxeyの整合条件式を使い、$k$-$\varepsilon$モデルの壁関数と矛盾しないプロファイルを導出する必要がある。この「乱流入口条件の設定方法」はAIJや国際的なガイドライン（COST Action C14など）に詳細が記載されており、建築CFDに固有の重要な知識です。
建築物の風荷重の実務適用
解析フロー
建築CFDの実務的な手順を教えてください。
歩行者風環境評価を例に説明しよう。
1. 対象建物と周辺建物のモデル化: 半径500m程度の周辺建物をGISデータから3Dモデル化
2. 計算領域設定: AIJガイドラインに準拠（閉塞率3%以下）
3. メッシュ生成: 500万--2000万セル（RANS）。建物周辺は2--5m、歩行者高さ1.5mに最低3層
4. 境界条件: 入口にべき法則プロファイル、出口に0Pa固定、上面・側面にスリップ壁
5. 16風向の計算: 22.5度刻みで全風向をカバー
6. 風環境評価: アメダスデータの風向頻度と組み合わせて超過確率を算出
7. 報告書作成: 風速分布図、超過確率マップを作成
16風向全部計算するのは大変ですね。
定常RANSなら1風向あたり数時間で計算できる。16風向で2--3日あれば完了する。LESだと1風向あたり数日かかるから、16風向では非現実的だ。そのため通常は卓越風向の数方向だけLESを実施する。
周辺建物のモデル化
周辺の建物はどこまで詳細にモデル化すべきですか？
モデル化の範囲と詳細度には目安がある。
| 距離 | モデル化方針 | 理由 |
|---|---|---|
| 0--100m | 詳細モデル（外壁凹凸含む） | 対象建物への直接的な影響 |
| 100--300m | 簡略モデル（直方体近似） | ビル風の増速効果 |
| 300--500m | さらに簡略化 | 大気境界層の形成 |
| 500m以上 | 地表面粗度で代替 | 直接的な影響は小さい |
ビル風の増速効果って具体的にはどういう現象ですか？
2棟の高層ビルの間の隙間で風速が増幅される「ビルドワインド効果」が代表的だ。自由流風速の1.5--2.0倍に達することもある。対象建物の風上に高層ビルがある場合は、その下降気流（ダウンウォッシュ）も歩行者レベルに大きな影響を与えるんだ。
検証と妥当性確認
建築CFDの結果をどう検証すればいいですか？
AIJが公開しているベンチマークケースで検証するのが標準的だ。
- 単体角柱: 2:1:1角柱まわりの風圧分布（AIJデータベース）
- 2棟並列配置: 2棟間の増速効果の検証
- 4:4:1角柱: 高層建物を想定したベンチマーク
- 新宿副都心モデル: 実市街地の風洞実験データ
これらのベンチマークでCFDの設定を検証してから、本番の建物解析に移行するのがベストプラクティスだ。
よくある失敗と対策
| 症状 | 原因 | 対策 |
|---|---|---|
| 入口プロファイルが維持されない | 壁関数と入口条件の不整合 | Richards-Hoxey型の整合条件を使用 |
| 風圧係数が過大 | 閉塞率が高い | 計算領域を拡大（閉塞率3%以下） |
| 後流の非対称解 | メッシュの非対称、不十分な収束 | メッシュ対称化、収束基準を強化 |
| 地表面近傍の風速が過小 | プリズム層の不足 | 地面にも$y^+=1$相当のプリズム層を配置 |
入口プロファイルが維持されないって、よく聞く問題ですよね。
k-epsilonモデルでは入口の乱流プロファイルが建物に到達する前に崩れることがある。これはRichards-Hoxeyの整合条件を満たしていないことが原因だ。入口の$k$と$\varepsilon$のプロファイルが壁関数の仮定と矛盾しないように設定する必要があるんだよ。
歩行者風環境基準がCFDで審査されるようになった経緯
日本では2000年代から超高層建物の建築確認申請にCFDによる風環境評価が実質的に必要になってきました。「ビルの隅で突風が吹いて歩行者が倒れる」問題への対応です。AIJ（日本建築学会）が風環境の評価基準を整備し、16風向×複数風速のCFD解析と気象統計を組み合わせた「年間超過頻度」による評価手法が普及。今では「CFD無しで超高層の確認申請は通らない」とも言える状況で、建築確認行政にCFDが組み込まれた特異な分野になっています。
建築物の風荷重のソフトウェア比較
主要ツールの比較
建築風解析に使えるCFDソフトにはどんなものがありますか？
| ツール | 開発元 | 特徴 | 建築風解析での強み |
|---|---|---|---|
| Ansys Fluent | Ansys | 汎用CFD | RANS/LES両対応、豊富な乱流モデル |
| STAR-CCM+ | Siemens | 自動メッシュ | トリムメッシュで大規模都市モデルに強い |
| OpenFOAM | OSS | 無償、カスタマイズ可 | 研究・中小事務所で普及 |
| WindPerfectDX | 環境シミュレーション | 建築専用 | 日本の基準に準拠、GUI操作 |
| STREAM | Cradle (MSC) | 熱流体統合 | 建築設備設計との連携 |
| Autodesk CFD | Autodesk | BIM連携 | Revitからの直接インポート |
WindPerfectDXって日本の建築業界では結構使われていますよね。
そうだ。WindPerfectDXはAIJガイドラインに準拠した大気境界層入口条件がプリセットされていて、建築技術者がCFDの専門知識なしでもある程度使えるように設計されている。ただしカスタマイズ性は汎用ソルバーに劣るよ。
OpenFOAMの活用
OpenFOAMで建築風解析をやるケースは多いですか？
大学・研究機関を中心に非常に多い。ライセンスコストがかからないのが最大の利点だ。
建築CFDでよく使われるOpenFOAMのソルバー:
- simpleFoam: 定常RANS（歩行者風環境）
- pimpleFoam: 非定常RANS/LES
- buoyantSimpleFoam: 浮力効果を含む場合（温熱環境）
大気境界層の入口条件にはatmBoundaryLayerInletVelocity とatmBoundaryLayerInletEpsilon が標準で用意されている。
BIM連携
BIMモデルからCFD解析への連携はどうなっていますか？
BIM（Revit等）からCFDへの連携は発展途上だ。現状の主なワークフローは以下の通り。
| ワークフロー | 特徴 | 課題 |
|---|---|---|
| Revit → STL → CFD | 最も一般的 | 形状簡略化が手動 |
| Revit → Autodesk CFD | 直接連携 | 大規模都市モデルには不向き |
| GIS → STL → CFD | 周辺建物の効率的モデル化 | LODの調整が必要 |
| CityGML → CFD | 都市スケール | データ整備状況に依存 |
建物の形状をどこまで簡略化していいかの判断が難しそうですね。
バルコニーや庇は歩行者レベルの風環境に影響するから残すべきだ。一方、窓枠や手すりなど翼弦の1/100以下の微細形状は省略しても影響は小さい。判断に迷ったら感度分析を行うのが確実だよ。
ツール選定の指針
| 用途 | 推奨ツール | 理由 |
|---|---|---|
| 歩行者風環境（実務） | WindPerfectDX, Fluent | 基準準拠、実績豊富 |
| 大規模都市CFD | STAR-CCM+ | 自動メッシュ、大規模対応 |
| 研究・教育 | OpenFOAM | 無償、カスタマイズ性 |
| 設備設計連携 | STREAM | 温熱・換気の統合解析 |
| BIM連携 | Autodesk CFD | Revitとの直接連携 |
BIM連携で建築CFDが「設計サイクル内」に入ってきた
かつて建築風環境CFDは「設計ほぼ確定後に専門家に外注する」ものでした。でも最近は建築設計ソフトRevitのBIMモデルから直接CFDモデルを生成できるAutodesk CFDやSimScaleが普及し、設計初期から建築家自身がCFDを回す時代になりつつあります。「方角を5度変えると風環境がどう変わるか」を打合せ中にリアルタイムで確認できるようになった。専門家が数週間かけていた解析が数時間に短縮されることで、設計の意思決定フローが根本から変わりつつあります。
建築物の風荷重の先端研究
都市スケールCFD
都市全体を対象としたCFDって可能なんですか？
計算機の発展により、数km四方の都市スケールCFDが実現しつつある。数億セルのLESも可能になってきた。
主な研究テーマ:
- マイクロクライメート予測: ヒートアイランド現象の解析。建物表面温度と大気の熱交換をモデル化
- 汚染物質拡散: 交通排気ガス、工場排煙の都市内拡散予測
- パンデミック対策: 建物換気口からのエアロゾル拡散シミュレーション
- 再生可能エネルギー: 都市内での小型風力発電ポテンシャル評価
都市スケールだとメッシュ数がとんでもないことになりそうですね。
1km四方の都市で建物を解像するには最低でも1000万セル、LESなら1億セル以上が必要だ。GPUソルバー（ProLB、XFlowなど格子ボルツマン法ベース）が計算効率の面で注目されているよ。
格子ボルツマン法
格子ボルツマン法（LBM）は建築分野でどう使われていますか？
LBMはN-S方程式を直接解くのではなく、分子の衝突・移流をメゾスケールでモデル化する手法だ。建築風解析には大きなメリットがある。
代表的なLBMソルバー:
- XFlow (Dassault): 商用。建築風解析の実績豊富
- ProLB: 商用。GPUネイティブで超高速
- Palabos: OSS。研究用途
風力-構造連成
超高層ビルの風応答解析はどうやりますか？
CFDで算出した変動風力を構造モデルに入力して、風応答（揺れ）を予測する。
アプローチは2通りある。
1. ワンウェイ連成: CFDの変動風力→FEMの動的解析。建物変形が風力に影響しない場合
2. 双方向FSI: CFDとFEMを時間ステップごとに連成。渦励振など変形がフィードバックする場合
渦励振が問題になるのはどんな建物ですか？
建築物の風荷重のトラブル対応
よくあるトラブルと対策
建築風CFDでよくある問題を教えてください。
1. 入口境界層プロファイルが維持されない
症状: 入口で設定した風速プロファイルが建物到達前に変化してしまう
原因: k-epsilonモデルの壁関数と入口条件の不整合、または地面粗度の不適切な設定
対策:
- Richards-Hoxeyの整合条件を使用: $k$, $\varepsilon$ のプロファイルが平衡状態を満たすように設定
- 地面の壁関数に対応する粗度長$z_0$をKs（等価砂粗さ）に変換: $K_s = \frac{9.793 z_0}{C_s}$（$C_s$はFluent粗度定数、デフォルト0.5）
- 入口から建物まで十分な距離（5H以上）を確保
粗度長からKsへの変換は結構ハマりやすいですよね。
そうだ。FluentとOpenFOAMで粗度パラメータの定義が異なるから注意が必要だ。OpenFOAMではatmBoundaryLayer系の境界条件を使えば$z_0$を直接指定できる。
2. 歩行者レベルの風速が実測と合わない
症状: 地上1.5mの風速がCFDで過大または過小
対策:
- 地面と建物壁面のプリズム層を適切に設定（$y^+ \approx 30$--300で壁関数使用、または$y^+ < 1$で壁面解像）
- 歩行者高さに最低3層のセルを確保
- 街路樹や低層構造物のモデル化を検討（多孔体モデルで代替可能）
- 風向の微小な変化で結果が大きく変わる場合があるため、前後2風向の感度を確認
3. 建物後方の渦が非対称
対称な建物なのに後流が非対称になることがあるんですが。
症状: 定常RANSで対称建物の後方渦が片側に偏る
原因: メッシュの微小な非対称性がトリガーとなり、双安定状態の一方に収束
対策:
- 対称面を利用した半分モデルで計算（ただし非対称渦放出は捉えられない）
- 非定常解析に切り替えて時間平均を取る
- メッシュを厳密に対称化する
4. LESで非物理的な結果
症状: LESで壁面近傍にスパイク状の速度/圧力が発生
対策:
- 壁面近傍の格子がLESの解像要件を満たしているか確認（$\Delta x^+ < 50$, $\Delta z^+ < 20$）
- サブグリッドスケールモデルの選択を見直す（動的Smagorinskyが推奨）
- 時間刻みがCFL < 1を満たしているか確認
- 入口の乱流生成法を確認（Vortex Method、Digital Filter等）
品質保証チェックリスト
報告書を出す前にチェックすべき項目を教えてください。
建築CFDの品質保証チェックリスト:
- 計算領域の閉塞率は3%以下か
- 入口プロファイルは建物位置で維持されているか（空走テスト実施）
- メッシュ収束性は3水準で確認したか
- AIJベンチマーク（単体角柱等）で検証済みか
- 結果は物理的に妥当か（ビル谷間の増速率1.5--2.0倍程度）
- 16風向の計算は完了しているか
- アメダスデータの統計期間は10年以上か
空走テストって何ですか？
建物を置かずに空の計算領域で計算し、入口から出口まで境界層プロファイルが維持されるかを確認するテストだ。プロファイルの変化が5%以内であれば入口条件の設定は適切と判断できる。建築CFDでは必須の検証ステップだよ。
「角地の風」が予測より強い理由——コリドー効果の罠
建物が密集した市街地で、2棟のビルの隙間に走る強風（コリドー効果またはビルの風のベンチュリ効果）は、CFDで繰り返し過小評価されることで知られています。原因は2つあります。一つは隣接建物の高さや形状の簡略化で流路断面積が過大になること。もう一つは定常RANSが非定常の突風ピークを予測できないこと。竣工後に「想定外の強風」でクレームが入るケースのかなりの割合がこのコリドー効果の予測不足です。LESと実測風速計の設置の組み合わせで事後対策するケースが後を絶ちません。
関連トピック
なった
詳しく
報告