建築物の風荷重解析
理論と物理
概要
先生、建物まわりの風の解析って何を求めるんですか？
大きく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で使われる乱流モデルを整理しよう。
建築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と風洞実験を組み合わせた検証を経て採用されたこの工夫が、世界一高い自立式電波塔の安全性を支えています。
各項の物理的意味
- 時間項 $\partial(\rho\phi)/\partial t$：蛇口をひねった瞬間を思い浮かべてください。最初は水がバタバタと不安定に出て、しばらくすると安定した流れになりますよね？ この「変化している最中」を記述するのが時間項です。心臓の拍動で血流が脈打つのも、エンジンのバルブが開閉するたびに流れが変動するのも、すべて非定常現象。では定常解析とは？ 「十分時間が経って流れが落ち着いた後」だけを見る——つまりこの項をゼロにする。計算コストが大幅に下がるため、まず定常で解いてみるのがCFDの基本戦略です。
- 対流項 $\nabla \cdot (\rho \mathbf{u} \phi)$：川に落ち葉を落としたらどうなりますか？ 流れに乗って下流に運ばれますよね。これが「対流」——流体の動きが物を運ぶ効果です。暖房の温風が部屋の端まで届くのも、空気という「運び屋」が熱を対流で輸送しているから。ここが面白いところ——この項は「速度×速度」を含むため非線形です。つまり、流れが速くなるとこの項が急激に強くなり、制御が難しくなる。これが乱流の根本原因です。よくある勘違い：「対流と伝導は同じようなもの」→ 全然違います！ 対流は流れが運ぶ、伝導は分子が伝える。桁違いの効率差があります。
- 拡散項 $\nabla \cdot (\Gamma \nabla \phi)$：コーヒーにミルクを入れて放置したことはありますか？ かき混ぜなくても、しばらく経つと自然に混ざりますよね。あれが分子拡散です。では次の質問——ハチミツとお水、どちらが流しやすいですか？ 当然お水ですよね。ハチミツは粘性（$\mu$）が高いから流れにくい。粘性が大きいと拡散項が強くなり、流体は「もったりした」動きになります。レイノルズ数が小さい流れ（ゆっくり、ドロドロ）では拡散が支配的。逆にRe数が大きい流れでは対流が圧倒し、拡散は脇役になります。
- 圧力項 $-\nabla p$：注射器のピストンを押すと、液体が針先から勢いよく出ますよね？ なぜでしょう？ ピストン側が高圧、針先が低圧——この圧力差が流体を押す力になるからです。ダムの放水も同じ原理。天気図で等圧線がギュッと密になっている場所では？ そう、強風が吹きます。「圧力差があるところに流れが生まれる」——これがナビエ-ストークス方程式の圧力項の物理的意味。ここでの勘違いポイント：CFDの「圧力」は絶対圧ではなくゲージ圧のことが多い。圧縮性解析に切り替えたとたんに結果がおかしくなる場合、絶対圧/ゲージ圧の混同が原因かもしれません。
- ソース項 $S_\phi$：暖められた空気が上に昇る——なぜでしょう？ 周囲より軽く（密度が低く）なったから、浮力で押し上げられるのです。この浮力はソース項として方程式に追加されます。他にも、ガスコンロの炎で化学反応熱が発生する、工場の電磁ポンプで金属溶湯にローレンツ力がかかる…これらはすべて「外部から流体にエネルギーや力を注入する」作用であり、ソース項で表現します。ソース項を忘れるとどうなるか？ 自然対流の解析で浮力を入れ忘れると、流体は一切動かない——冬の部屋で暖房をつけたのに暖かい空気が上に行かない、という物理的にありえない結果になります。
仮定条件と適用限界
- 連続体仮定：クヌッセン数 Kn < 0.01（分子平均自由行程 ≪ 代表長さ）で成立
- ニュートン流体仮定：せん断応力と歪み速度が線形関係（非ニュートン流体では粘度モデルが必要）
- 非圧縮性仮定（Ma < 0.3の場合）：密度を一定として扱う。マッハ数0.3以上では圧縮性効果を考慮
- ブシネスク近似（自然対流）：密度変化を浮力項のみで考慮し、他の項では一定密度を使用
- 適用外ケース：希薄気体（Kn > 0.1）、超音速・極超音速流れ（衝撃波捕捉が必要）、自由表面流れ（VOF/Level Set等が必要）
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 速度 $u$ | m/s | 入口条件で体積流量から換算する際、断面積の単位に注意 |
| 圧力 $p$ | Pa | ゲージ圧と絶対圧の区別。圧縮性解析では絶対圧を使用 |
| 密度 $\rho$ | kg/m³ | 空気: 約1.225 kg/m³@20°C、水: 約998 kg/m³@20°C |
| 粘性係数 $\mu$ | Pa·s | 動粘性係数 $\nu = \mu/\rho$ [m²/s] との混同に注意 |
| レイノルズ数 $Re$ | 無次元 | $Re = \rho u L / \mu$。層流/乱流遷移の判定指標 |
| CFL数 | 無次元 | $CFL = u \Delta t / \Delta x$。時間刻みの安定性に直結 |
数値解法と実装
計算領域とメッシュ
建物まわりの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倍になるが、局所的なピーク負圧を正しく予測できるメリットは大きいよ。