技術資料
Feel&Think
第15回 めざせ三軸試験のシミュレーション
今回からは、三軸圧縮試験結果を表現できる構成方程式(応力とひずみの関係を表現する理論式)を組み立ててみます。また、この理論式を用いて三軸圧縮試験をシミュレートし、実際の結果と比較します。この過程で、これまで行ってきた岩石の変形メカニズムの推定の妥当性が検証されます。さらには、岩石の挙動から地下の岩盤の挙動を想像する手掛かりを得ることができると考えます。
最初に行うのは、岩石試料に軸差応力を加えて行った際に見られた、特徴的な体積ひずみの表現です。軸差応力(圧縮応力)を加えていくと、最初は体積が減少しますが途中から膨張に転じ、有効拘束圧が大きくなると膨張が発生しにくくなることは、これまで示してきたとおりです(図-1)。
![図-1 軸差応力と体積ひずみの関係の一例](/wp/wp-content/uploads/2023/07/rm15_z01.png)
この体積膨張が、破壊した初期クラック端より発生した二次クラックの、開口を伴った伸長に起因すると考え、破壊力学をもとに二次クラックの伸長方向を検討しました。この結果、軸差応力が有効拘束圧に比べ充分に大きい場合には、開口した二次クラックは最大圧縮応力方向に伸長すると近似でき、このときの二次クラックの相対長は次式で与えられました(図-2参照)。
![数式1](/wp/wp-content/uploads/2023/07/rm15_s01.png)
すなわち、開口した二次クラックが最大圧縮応力方向に伸長すると近似すれば、二次クラックの長さは軸差応力に比例し、有効拘束圧に反比例することになります。
![図-2 最大圧縮応力方向へ伸長する二次クラック](/wp/wp-content/uploads/2023/07/rm15_z02.png)
さて、最大圧縮応力方向に二次クラックが伸長する場合には、岩石の体積ひずみはどうなるのでしょうか。体積ひずみを求めるためには、二次クラックの変形から体積変化を求め、さらに二次クラックの個数分体積変化を加え合わせ、全体の体積で除すといった手続きが必要です。これを順に行っていきましょう。
まず、二次クラックの開口に伴う体積変化を求めます。図-2(b)のモデルでは、二次クラック面に作用するそれぞれの力Pn、σ2、φによる二次クラック面の変形が、弾性論より次のように表されます。
![数式2](/wp/wp-content/uploads/2023/07/rm15_s02.png)
![数式3](/wp/wp-content/uploads/2023/07/rm15_s03.png)
![数式4](/wp/wp-content/uploads/2023/07/rm15_s04.png)
ここに、
![数式5](/wp/wp-content/uploads/2023/07/rm15_s05.png)
であり、E*とν*は、それぞれクラック以外の部分のヤング率とポアソン比です。ただし、これらの値は計測できません。
この変形を二次クラック面に沿って積分すれば、体積変化となります。
![数式6](/wp/wp-content/uploads/2023/07/rm15_s06.png)
それぞれについて体積変化を求めると、次のようになります。
![数式7](/wp/wp-content/uploads/2023/07/rm15_s07.png)
![数式8](/wp/wp-content/uploads/2023/07/rm15_s08.png)
![数式9](/wp/wp-content/uploads/2023/07/rm15_s09.png)
これらの和が、二次クラックがもたらす体積変化です。式(1)も考慮して、
![数式10](/wp/wp-content/uploads/2023/07/rm15_s10.png)
![数式11](/wp/wp-content/uploads/2023/07/rm15_s11.png)
が得られます。ここで、
![数式12](/wp/wp-content/uploads/2023/07/rm15_s12.png)
であり、γが初期クラックの長さaとその角度θの関数であることから、次のように書くことができます。
![数式13](/wp/wp-content/uploads/2023/07/rm15_s13.png)
ここに、 であり、三軸圧縮試験では軸差応力に相当します。
二次クラックの開口による体積ひずみは、二次クラックの個数分体積変化を加え合わせ、これを全体の体積Vで除すことで求まります。
![数式14](/wp/wp-content/uploads/2023/07/rm15_s14.png)
γの総和は、γの平均値 を二次クラックの個数N分加え合わせたもので近似できるものとします。
![数式15](/wp/wp-content/uploads/2023/07/rm15_s15.png)
二次クラックの個数Nですが、これは破壊した初期クラックの数に他なりません。そこで、ある軸差応力で破壊する初期クラックの数を表し、全体の総和を1とする破壊密度関数Dを導入します。これに、破壊する可能性のある単位体積あたりの初期クラックの数を表すρと全体の体積を乗ずると、ある軸差応力で破壊する初期クラックの個数となります。ある軸差応力に達するまでに破壊した初期クラックの総数は、破壊密度関数をその軸差応力まで積分して求められます。
![数式16](/wp/wp-content/uploads/2023/07/rm15_s16.png)
最終的に、二次クラックの開口を伴った伸長による体積ひずみは次式となります。すなわち、二次クラックによる体積ひずみは有効拘束圧に反比例します。
![数式17](/wp/wp-content/uploads/2023/07/rm15_s17.png)
ここまでで、軸差応力や有効拘束圧と体積ひずみの関係は求まりました。次回は、軸差応力と軸ひずみの関係を検討しますが、軸ひずみは弾性的に挙動しているという、驚きの仮説を紹介します。