Feel＆Think

今回からは、地層科学研究所が開発中のG-SUPRAを使って、凍結膨張のシミュレーションを行っていきます。なお、以下の解析では圧力の単位としてMPa、長さの単位としてm、温度の単位として℃、熱量の単位としてJ、時間の単位としてhを用いています。
解析に用いた物性値を表-1に示します。

表-1　解析に用いた物性値一覧

用いたモデルは、図-1に示したような各辺1mの立方体モデルです。最初に、自重状態を作ります。

将来行う凍結膨張解析では、x、y、zの3方向への変形量を調べたいことから、変形の境界条件として原点を含む側面と底面の法線方向の変位を0とするとともに、モデルの上面と側面は表面力境界とし変形が生ずるようにしました。
ただし、自重解析による変形が不均一とならないように、側面では自重による変形が0となるような表面力を加えました。これは次のように求めました。
解析モデルの間隙水圧分布は、自重解析の結果次式となるはずです。φ0はモデル上面の間隙水圧です。

また、水平面内の全応力分布は次のようになるはずです。

水平変位を0とするためには、水平方向全応力と間隙水圧の差（有効応力）にv/(1－v)を乗じ、間隙水圧を加えた大きさの表面力を加えることになります（図-2）。

水に関する境界条件としては、側面からの水の出入りを禁止し、上面と下面で間隙水圧が一定に保たれるものとしました（図-3）。また、熱に関する境界条件としては、上面と側面で熱の出入りを禁止し、底面は温度変化が可能なように温度固定境界としました（図-4）。

解析結果を、解析モデルを単一要素とした場合と要素分割した場合について以下に示します。
図-5と図-6が、解析により得られた自重による沈下です。狙いどおり、水平方向の変位が0となっていることがわかります。要素分割した場合でも同様です。要素分割した場合では、図-7と図-8に示すように、間隙水圧や鉛直応力が深さ方向に線形に分布します。次回からは、凍結膨張の解析に取り掛かります。