📄 other

Multiscale Calculation of Equivalent Elastic Modulus for Rock Mass Based on Perturbation Theory

本論文は、岩盤の等価弾性係数を正確に決定するために、初期損傷とメソスケールの不均質性を組み込んだ摂動論に基づくマルチスケール計算手法を提案し、室内試験および現場試験の結果を用いた検証を通じて、従来の手法と比較して優れた精度を実証するものである。

原著者： Lei Wen, Wu Yang, Huayun Yang

公開日 2026-07-24

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

CC BY 4.0

原著者： Lei Wen, Wu Yang, Huayun Yang

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

巨大で崩れかけの石造りの城が、圧力を受けたときにどのように耐えるかを予測しようとしている場面を想像してみてください。しかし、ここにはひねりがあります。その城は、滑らかなものもあれば、ギザギザしたものもある岩の乱雑な積み重ねであり、それらはセメントのパッチによって繋ぎ止められ、中にはクモの巣のように隠れた亀裂が走っています。これが、「岩盤（ロックマス）」工学の現実の世界です。エンジニアがトンネル、鉱山、あるいはダムを建設する際、研究所で完璧で小さな岩の立方体をテストし、その結果をもって山全体が同じように振る舞うと仮定することはできません。それは、ビル全体が地震に耐えられるかを確認するために、レンガ一個をテストするようなものです。問題は、岩石が乱雑であり、その曲がり方や壊れ方に影響を与える微細な欠陥や亀裂に満ちていることです。科学者たちは、微視的な微細な亀裂の世界と、巨大な岩石形成の世界との間の溝を埋めようと長い間試みてきましたが、彼らの古い地図はしばしば誤った方向へと導きました。彼らには、木を見て森も同時に見ることができる、新しい方法が必要でした。

「摂動論に基づく岩盤の等価弾性係数のマルチスケール計算（Multiscale Calculation of Equivalent Elastic Modulus for Rock Mass Based on Perturbation Theory）」と題されたこの論文は、エンジニアが、凹凸のある亀裂の入った岩がどのように押しつぶされ、伸びるのかを正確に把握するための、新しい、非常にスマートな眼鏡のようなものです。著者である Lei Wen、Wu Yang、そして Huayun Yang は、「マルチスケール計算」という巧妙な数学的トリックを開発しました。これは、岩石の初期損傷（そこにある微細な亀裂）を起点として扱うものです。彼らは岩石が完璧であると仮定する代わりに、「摂動論」と呼ばれる手法を用いました。これは、基本的には「完璧なシステムをわずかに揺らして、どのように揺らぐかを見る」という高度な方法です。これを用いて、微細な亀裂が岩の剛性にどのように影響を与えるかを算出します。彼らは、これらの微細な亀裂によって生じる歪み（押しつぶされること）を、微視的な世界と巨視的な世界を繋ぐ「架け橋」と呼んでいます。

チームはただコンピュータに向き合っていたわけではありません。彼らは中国南部の深い金属鉱山へと赴き、実際の岩石を採取しました。彼らは、これらの岩石を異なる強度の力（具体的には、完璧な岩石の強度の0.25、0.4、0.5、0.65、および0.75倍）で押しつぶすことで、意図的に損傷を与え、「初期損傷」を作り出しました。その後、彼らはこれらの損傷した岩石に対して、2つのテストを実施しました。一つは「ブラジル分割試験」（岩の両側から押しつぶしてオレンジのように割れるまで圧縮するテスト）、もう一つは標準的な「一軸圧縮試験」（上から下へと押しつぶすテスト）です。さらに、彼らは音波（P波速度）とCTスキャンを使用して、岩石内部の目に見えない亀裂を観察しました。

彼らが発見したのは、彼らの新しい数学的手法が、従来の標準的なコンピュータ・シミュレーションよりも現実を予測する上ではるかに優れているということです。彼らが結果を実際のラボテストと比較したところ、彼らの「マルチスケール」アプローチは、実際の岩石とほぼ完璧に一致する曲線を描きました。従来のメソッドは、ピーク強度（壊れるまでにどれほど強く押せるか）を推測することには成功していましたが、曲線の形状を間違えていました。新しい手法は、引張荷重－変位曲線（分割試験によるもの）が、特定の「厳密に低い凸（strictly lower convex）」の形で曲がることを示しており、これは実生活と一致していますが、古いシミュレーションはこのニュアンスを見落としていました。対照的に、圧縮応力－ひずみ曲線（押しつぶしテストによるもの）は、異なる「凸（convex）」の形状に従います。彼らはさらに、この手法を、崩落が発生した鉛亜鉛鉱山の実際の建設現場へと適用しました。彼らは、「注入固結体（grouting solidified bodies）」、つまり、セメントのスラリーで再び接着された壊れた岩石の強度を計算するために、この手法を用いました。数学は、その接着された塊が破壊点に達する前の挙動を正確に予測しました。

しかし、著者たちは、自分たちの手法が魔法ではないことも注意深く述べています。彼らのモデルは、岩石が壊れる前の挙動については完璧に捉えましたが、岩石が完全に崩壊し始めて「ひずみ軟化（strain-softening）」が起こる現象を正確に予測することには苦戦しました。これは、彼らのモデルが、岩石が完全に破壊される際に起こる乱雑なプロセスをまだ考慮していないためです。しかし、最も重要な部分、すなわち、損傷した岩盤がいかに硬く、いかに強いかを知るための、微視的世界と巨視的世界の間のこの新しい「架け橋」は、従来のツールよりも優れた成果を上げています。これは、岩石の初期の傷跡を認め、この摂動数学を用いることで、エンジニアが岩盤が持ちこたえるか崩壊するかについて、より明確なイメージを得られることを示唆しており、私たちのトンネルや鉱山をより安全なものにしています。

技術要約：摂動論に基づく岩盤の等価弾性係数のマルチスケール計算

問題提起
等価弾性係数の正確な決定は、岩盤工学の安全性および安定性評価において極めて重要である。現在の手法には重大な限界が存在する。 intact（無傷）な岩石試料を用いた室内試験では、工学的岩盤に固有のマイクロおよびマクロなジョイントや亀裂を考慮できていない。一方で、現場試験は空間的範囲が限定されているため、代表性に欠けることが多い。さらに、室内試験の結果を現場条件へとスケールアップするために用いられる従来の低減法は、信頼性の低い低減係数に依存しており、岩盤の破壊過程における力学的応答プロセスを捉えることができない。したがって、intactな岩石力学と、損傷した岩盤の不均質挙動との間の溝を埋める、体系的なマルチスケール結合アプローチが求められている。

手法
著者らは、岩盤における初期損傷を考慮するために設計された、摂動論に基づくマルチスケール計算手法を提案している。この手法は、内部亀裂のマルチスケール展開プロセスを反映する小さなパラメータ（ $\delta$ ）を用いて、岩盤を初期損傷を持つ摂動系として扱う。

主な手法のステップは以下の通りである：

理論的枠組み： メソ（中間）スケールの不均質性に起因する歪みテンソルを、マクロおよびミクロのスケールを接続する「架け橋」として利用する。変位場に対して漸近展開を行い、変数をマクロ座標（ $x$ ）とミクロ近傍（ $y$ ）に分離する。
支配方程式： 微小弾性損傷因子（ $D^\delta$ ）を組み込んだ構成方程式と平衡微分方程式を結合することにより、支配方程式を導出する。初期損傷は、初期弾性損傷テンソルによって表現される。
均質化： 異なる物理的次元スケール下での運動方程式を解くことにより、均質化されたマクロ等価弾性係数（ $E_{mnkl}^\delta$ ）を導出する。これには、変分原理に基づく有限要素法を用いて、テンソル固有関数（ $\chi$ ）を解くプロセスが含まれる。
数値実装： 損傷係数の反復計算にはMATLABを、有限要素モデリングにはCOMSOLを用い、本アプローチを実装した。モデルは、進化する損傷状態を反映するために、各ステップで荷重（ $Q$ ）と変位（ $u$ ）を更新する。
検証： 本アプローチは、主に2つのケースを通じて検証された：
1. 室内検証： 中国南部の深部金属鉱山から採取された石灰岩試料を用いて、一軸圧縮試験およびブラジル型割裂引張試験を実施した。試料には様々な程度（0.25 $\sigma_{cs}$ から0.75 $\sigma_{cs}$ ）の事前損傷を与え、P波速度およびCTスキャンによって特性評価を行った。
2. 現場検証： 鉛・亜鉛鉱山の坑内の充填固結体に本手法を適用した。ここで、充填材は不規則なジョイント裂罅を持つ岩盤として扱われた。

主な貢献

マルチスケール摂動モデル： 本論文は、メソスケールの不均質性に起因する歪みテンソルがミクロとマクロのスケールを接続するメカニズムとして機能する理論的枠組みを確立し、初期損傷を考慮した等価弾性係数の計算を可能にした。
動的な損傷進化： 弾性係数を定数として扱うことが多い従来のモデルとは異なり、本アプローチは、荷重過程における等価弾性係数の動的な発達、およびマクロおよびミクロの両方の損傷の進化を反映している。
精度の向上： 本研究は、提案されたマルチスケール計算が、従来の有限要素シミュレーション手法と比較して、実際の室内試験および現場試験のデータにより近い結果をもたらすことを示している。

結果

室内試験の結果：
- P波速度（ $V_p$ ）は初期損傷が増加するにつれて減少することが判明し、損傷指標としての有用性が確認された。
- CTスキャンにより、高い初期損傷は空隙率および亀裂密度の増加と相関していることが明らかになった。
- ブラジル型割裂引張試験： 荷重－変位曲線は厳密に下凸の関数を示した。マルチスケール計算の結果は、特にピーク前領域において、従来のシミュレーションよりも実験曲線によく一致した。
- 一軸圧縮試験： 応力－歪み曲線は単調増加する凹関数であった。マルチスケールアプローチは等価弾性係数を正確に予測し、その値は従来のシミュレーションよりも試験結果に近いものであった。損傷の進化は、初期には緩やかに増加し、その後急速に発達して安定化する様子を示した。
現場試験の結果：
- 充填固体の現場圧縮試験において、マルチスケールアプローチは最大一軸圧縮応力に至るまでの応力－歪み曲線を正確に予測した。
- 充填固体の算出された等価弾性係数は現場の試験データと一致しており、複雑な工学的シナリオ（砕石やセメントスラリーの集合体）に対する本手法の適用可能性を実証した。
観察された限界： 本研究では、ピーク前の結果は正確であったものの、充填固体のポストピーク（ピーク後）の応力－歪み曲線が現場の結果から乖離していることが指摘された。この不一致は、深部の脆性岩石とは異なる挙動を示す材料のポストフェイラー（破壊後）挙動を捉えるために必要な、歪み軟化パラメータがモデルに欠けていることに起因する。

意義および主張
本論文は、提案されたマルチスケール計算アプローチが、岩石工学における複雑な問題を解決するための堅牢な理論的手法を提供すると主張している。初期損傷を考慮し、摂動論を利用することで、本手法はミクロ構造の不均質性とマクロな力学的挙動の間の溝を克服することに成功した。著者らは、このアプローチが、岩盤の等価弾性係数を決定するための従来のシミュレーション手法に代わる、より正確な選択肢を提供すると断言している。また、マクロとミクロのスケール間の適切な架け橋さえ特定できれば、この手法を他の不均質材料の力学的パラメータ決定へと拡張できる可能性があると示唆している。本研究は、特に初期損傷やスケール効果が重要な要因となるエンジニアリングプロジェクトにおいて、岩盤の安定性評価を改善するためのリファレンスとして提示されている。

関連論文