⚡ electrical engineering

A Physically Constrained Energy-Based Excess Pore Pressure Model for Saturated Sands with an Ordinary Differential Equation Interpretation

本研究は、飽和したきれいな砂に対するエネルギーに基づく過剰間隙水圧モデルを、コーシー・オイラー常微分方程式の解として解釈することで再定式化し、それによって物理的な有界性を確保し、校正を単一のパラメータへと削減し、かつ様々な相対密度にわたって高い精度を実現するものである。

原著者： Ferruh Gündoğan, Ayfer Erken

公開日 2026-06-30

📖 1 分で読めます☕ さくっと読める

CC BY 4.0

原著者： Ferruh Gündoğan, Ayfer Erken

原論文は CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/) でライセンスされています。 ✨ これは以下の論文のAI生成解説です。著者が執筆または承認したものではありません。技術的な正確性については原論文を参照してください。免責事項の全文を読む

大きな全体像：なぜ砂の城は崩れるのか？

ビーチで砂の城を作っているところを想像してみてください。ただその上に立っているだけなら、形は保たれます。しかし、リズムよくその上でジャンプし始めると、砂の粒がガタガタと動き出します。水に浸かった砂の城（飽和状態の砂）の場合、この揺れによって、粒の間に閉じ込められた水が、砂を押し返そうとする力がどんどん強まっていきます。

やがて、水圧が非常に強くなり、砂の粒を完全に押し広げてしまいます。砂は固体としての性質を失い、ドロドロした液体（ミルクシェイクのようなもの）のように振る舞い始めます。これが「液状化」と呼ばれる現象であり、地震の際に建物が沈んだり傾いたりする原因となります。

エンジニアには、これが「いつ」起こるのかを正確に予測する方法が必要です。彼らは「散逸エネルギー（dissipated energy）」という概念を使います。これは、砂が自らを再配置するために費やす「ガタガタと揺れる」仕事量、あるいは総量だと考えてください。より激しく揺らすほど、より多くのエネルギーを吸収し、水圧は高まっていきます。

古い問題点：壊れた定規

数十年にわたり、エンジニアは水圧を予測するために、有名な公式（1985年にベリルとデイビスによって作成されたもの）を使用してきました。それは、「投入されるエネルギーが増えるほど、水圧も高くなる」と示す、単純な定規のようなものです。

しかし、この古い定規には2つの大きな欠陥がありました。

天井がない： この公式は永遠に上昇し続けます。現実の世界では、水圧は上にある土の重さを超えることはできません（100%を超えることはありません）。古い公式は時として、150%という不可能な数値を予測してしまい、物理法則を無視してしまいます。
絡まったつまみ： この公式には、本来は独立しているはずの2つの設定（パラメータ）がありますが、それらが実は互いにくっついてしまっています。一方を変えると、もう一方も変わってしまうのです。これにより、異なる種類の砂に対してモデルを調整（キャリブレーション）することが非常に困難になっています。

新しい解決策：スマートで境界のあるモデル

この論文の著者たち（GündوğanとErken）は、この「定規」を直すことに決めました。彼らは、2種類の砂を用いた実際のラボテストのデータを確認しました。一つはOttawa F65（多くの研究で使用される標準的なきれいな砂）、もう一つはAkpınar sand（イスタンブールの地元の砂）です。

単にグラフ上に直線を描くのではなく、彼らは、揺れによるエネルギーと水圧の関係は、もっと**「3段階の旅」**のようなものであることに気づきました。

開始時： ゆっくりとした緩やかな上昇。
中間期： 急激で急峻な上昇。
終了時： 限界値（100%の水圧）でちょうど平坦になる滑らかな曲線。

これを捉えるために、彼らは数学を一つの**「レシピ」**として扱いました。彼らは公式を2つの部分に分解しました。

成長エンジン： 急激な上昇を処理する基本的な項。
調停者（メディエーター）： 限界に近づくにつれて成長を抑制し、決して100%の天井を突破しないように制御する「交通整理役」。

「合成的」な秘密：数学的な物語

ここが彼らの手法の巧妙な点です。彼らは、自分たちの新しい曲線の形が、**常微分方程式（ODE）**と呼ばれる特定の物理方程式の解と全く同じ形をしていることに気づきました。

常微分方程式を、システムの「時間の経過に伴う変化についての物語」だと考えてください。

方程式の左辺は、砂の内部的な挙動（どのように抵抗するか）を表します。
右辺は、外部からの揺れ（地震）を表します。

彼らのモデルをこの方程式として書き直すことで、厳格な「交通ルール（物理的制約）」を適用することができました。

ルール1： 圧力は常に上昇しなければならない（砂を「逆方向に揺らす」ことはできないため）。
ルール2： 圧力は100%に達した瞬間に、必ず上昇を停止しなければならない（液状化）。
ルール3： 100%に達した瞬間、上昇の速度はゼロでなければならない（鋭いスパイクではなく、滑らかに水平にならなければならない）。

数学にこれらのルールを守らせることで、「絡まったつまみ」の問題を排除することができました。彼らは複雑なモデル全体を、見つけ出すべきたった一つの数字（bと呼ばれるもの）へと集約させたのです。残りの数学部分は、物理法則によって固定されています。

彼らが発見したこと

彼らはこの新しい「ワンナンバー・モデル（一つの数字によるモデル）」をラボのデータと比較検証しました。

精度： 実験結果と驚くほど一致しました（精度99%）。
安全性： 不可能な数値を予測することはありませんでした（0%から100%の間に収まります）。
簡潔さ： 調整すべき設定が一つだけなので、古いモデルよりもはるかに使いやすくなっています。

「アハー！」の瞬間：失敗の2つのステージ

数学は、砂がどのように破壊されるかについて興味深い事実を明らかにしました。

初期段階： 砂自身の内部構造（粒がいかに密に詰まっているか）が、圧力の上昇をコントロールしています。これが方程式の「同次（Homogeneous）」の部分です。
後期段階： 砂が液状化に非常に近づくと（圧力の75%付近）、砂は「記憶」を失います。内部構造が崩壊し、圧力の上昇はほぼ完全に外部からの揺れエネルギーによって駆動されるようになります。これが「特殊（Particular）」の部分です。

まとめ

著者たちは、砂の液状化を予測するための、乱雑で上限がなく、複雑な公式を取り壊し、再構築しました。彼らは、モデルが現実の物理学のように振る舞う（滑らかに上昇し、限界でピタリと止まる）ように、数学的な「物語（方程式）」を用いることで、モデルを構築しました。その結果、地震の際に砂がいつ液体に変わるかを理解するための、よりシンプルで、より安全で、より正確なツールが完成しました。

技術要約：飽和砂における物理的制約を設けたエネルギーに基づく過剰間隙水圧モデル

問題提起
土壌液状化は極めて重要な地盤工学的リスクであり、従来、応力ベース、ひずみベース、間隙水圧ベース、およびエネルギーベースの手法を用いて評価されてきた。その中でも、エネルギーベースの手法（EBM）は、累積散逸エネルギーと過剰間隙水圧の発生を関連付けることで、等価サイクルへの変換に伴う不確実性を排除し、不規則な地震動への対応が可能であるという明確な利点を持つ。最も影響力のあるEBMの定式化は、ベリル・デイビス（Berrill–Davis, 1985）の冪関数（ $r_u = \alpha W^\beta$ ）であるが、本研究ではこのモデルにおける2つの根本的な限界を特定している：

非有界な応答： 冪関数は数学的に非有界であり、過剰間隙水圧比（ $r_u$ ）が1を超えることを許容してしまう。これは、液状化において $r_u$ が1.0に漸近すべき飽和砂においては、物理的に不可能である。
パラメータの相互依存性： 校正パラメータである $\alpha$ と $\beta$ はしばしば独立したものとして扱われるが、経験的研究によれば、両者は互いに相関している。この依存性は、校正を複雑にし、物理的な解釈性を低下させ、外挿を困難にする。

手法
著者らは、これらの制限に対処するため、均一で清浄な飽和砂を対象とした再定式化モデルを提案する。そのプロセスは以下の多段階の導出過程を経る：

データ分析と分解： オタワF65砂（LEAP-2017）およびアクピナル（Akpınar）砂の繰り返し三軸試験データを用い、正規化された間隙水圧比（ $r_u$ ）と正規化された散逸エネルギー（ $W_N$ ）の関係を対数–対数スケールで分析する。著者らは、データが単一の直線ではなく、少なくとも3つの線形セグメントで構成されていることを観察した。
多項式近似： データは定数項のない6次多項式によって近似される。この関数は、主要な成長項（ $W_N^b$ ）と媒介関数 $g(W_N)$ に分解される。
テイラー展開と簡略化： 媒介関数をさらに分析する。指数成分のテイラー展開により、3次多項式がデータの傾向を適切に表現できることが示された。これにより、候補となるモデル形式 $r_u = W_N^b (a_0 + a_1 W_N + a_2 W_N^2 + a_3 W_N^3)$ が得られる。
常微分方程式（ODE）による解釈： 数学的構造を認識し、モデルを非斉次コーシー・オイラー常微分方程式の解として再解釈する。左辺は系の内部応答を表し、右辺は外部荷重条件を表す。
物理的制約： 非物理的な変動とパラメータの相互依存性を排除するため、モデルには液状化点（ $W_N = 1, r_u = 1$ $W_{N} = 1, r_{u} = 1$ ）における特定の境界条件が課される：
- $W_N = 1$ において $r_u = 1$ 。
- $W_N = 1$ において水平接線（ $r_u' = 0$ ）を持ち、曲線が限界に達する際に平坦化することを保証する。
- 荷重プロセス全体を通じて単調増加（ $r_u' > 0$ ）すること。
- 液状化時に非斉次項が消失することから導かれる、係数間の特定の関係（ $a_2 = -3a_3$ ）。

これらの制約により、媒介多項式の係数が単一のパラメータ $b$ の関数として解析的に決定され、モデルは $b$ のみに支配される形式へと簡素化される。

主な貢献

物理的に制約された定式化： 提案されたモデルは、 $0 \le r_u \le 1$ を本質的に満たすように構築されており、ベリル・デイビス・モデルに見られるような過剰間隙水圧の非物理的な過大予測を防ぐ。
単一の校正パラメータ： 媒介係数を成長指数 $b$ に解析的に結びつけることで、複数のパラメータを独立して校正する必要性を排除し、従来の定式化で見られた相互依存性の問題を解決した。
ODEフレームワーク： 本研究は、経験的な間隙水圧の関係を、非斉次コーシー・オイラー常微分方程式の解として捉える新しい解釈を提供する。このフレームワークは、応答を、初期段階で支配的な内部材料特性を司る斉次部分と、外部荷重条件を司る非斉次部分へと分離する。
検証： モデルは、様々な相対密度（緩い状態から非常に密な状態まで）におけるオタワF65砂およびアクピナル砂の高精度な繰り返し三軸試験データに対して検証され、 $R^2 \ge 0.99$ を達成した。

結果

簡略化されたモデル（式20）は、オタワF65砂およびアクピナル砂の測定された間隙水圧履歴を正確に再現した。
パラメータ $b$ は相対密度によって変化することが判明した（例：相対密度51%では $b \approx 0.36$ 、80–91%では $b \approx 0.43$ ）。
モデルは、液状化に近づくにつれての初期の急速な上昇と、その後のプラトー（停滞）を含む、非線形な間隙水圧の進化を捉えている。
元のベリル・デイビス・モデルとの比較により、提案された定式化は、特に完全な液状化という終端条件付近において、元のモデルが大きなばらつきや非有界な成長を示すのに対し、より一貫した間隙水圧進化の記述を提供することが示された。
残りのLEAP-2017データセットへの適用により、局所的な極大値を避けるための条件 $b \ge 0.27$ を満たせば、様々な相対密度および繰返しせん断応力比にわたって間隙水圧発生をシミュレートできることが確認された。

意義
本論文は、本研究が、均一で清浄な砂を用いた液状化研究において、古典的なベリル・デイビスの冪関数に対する、実行可能でコンパクト、かつ物理的に一貫した代替案を提供すると主張している。応答を制限し、ODEベースの導出を通じてパラメータの相互依存性を排除することで、高い精度を維持しながら校正を簡素化している。著者らは、この定式化は経験的なものであるものの、厳格な物理的制約の下での回帰から導出されているため、液状化の根底にあるメカニズム（内部の粒状特性に支配される初期相（斉次解）と、材料が機械的な記憶を失い、完全な液状化に至るために多大なエネルギーを必要とする相（非斉次解）の連続）を捉えることができると強調している。なお、適用範囲は均一で清浄な飽和砂に限定されている。