ため 2), 供試体の締固め度を 85% としている 表 1 土試料の物理特性 土試料の物理特性 土粒子密度 [g/cm 3 ] 最大間隙比 最小間隙比.761 平均粒径.354 [mm] 均等係数 76.7 細粒分含有率 26. 最適含水比 9.2 最大乾燥密度 1.898

Transcription

1 第 59 回地盤工学シンポジウム平成 26 年度論文集 19-5 不飽和砂質土の三軸圧縮および等方圧縮時の挙動と空気 - 水 - 土連成解析法によるシミュレーション Mechanical behavior of unsaturated sand under triaxial and isotropic compression, and simulation using Air-Water-Soil coupled method 赤木俊文 *, 森本恭弘 **, 石川椋 ***, 木元小百合 **** Toshifumi AKAKI, Yasuhiro MORIMOTO, Ryo ISHIKAWA, Sayuri KIMOTO 本研究では, 不飽和砂質土の強度 - 変形特性を把握するため, 木津川堤防で採取した砂質土を用いて非排気 - 非排水条件下で三軸圧縮試験を行った また, サクション一定条件で等方圧縮試験を行い, 正規圧密線および限界状態線について骨格応力を用いて整理した その結果, 実験条件によらず限界状態における応力比はほぼ同一の値を示した また, 限界状態における間隙比 - 骨格応力関係において初期サクションの増加に伴い曲線の傾きが大きくなる傾向が見られた さらに, 弾粘塑性構成式および間隙比依存性を考慮した水分特性モデルを用いた空気 - 水 - 土連成解析手法により, 非排気 - 非排水条件の三軸試験の再現解析を行い, 応力 -ひずみ挙動や供試体内のひずみ等の分布について検討した キーワード : 不飽和土骨格応力限界状態サクション Unsaturated soil, skeleton stress, critical state, suction 1. はじめに土 水 空気の三相からなる不飽和土の力学挙動は土の水分量に大きく影響され, 降雨や地下水位の変動により容易に変動するなど, 飽和土の力学挙動に比較して複雑である 不飽和土は自然界における一般的な土の存在形態であり, 地盤工学的問題が発生する実地盤は多くの場合不飽和土である 加えて近年, 我が国において集中豪雨や巨大台風が多発し, それに伴い斜面の崩壊, 道路盛土の変形, 河川堤防の盛土における浸透変形問題など, 不飽和土に関わる重大な地盤工学的問題が発生している こうしたことから, 不飽和土の力学特性の研究が関心を集めるとともに, その重要性を増している状況にある そこで, 本研究では非排気 - 非排水条件下での三軸圧縮試験および不飽和等方圧縮試験を行い, 不飽和土の力学特性について調べた また, 非排気 - 非排水三軸圧縮試験の結果に対して, 空気 - 水 - 土連成解析手法による再現解析を試み, 試験中における応力やひずみの挙動および供試体内 のひずみ等の分布について検討した 特に水分特性曲線に間隙比依存性を表現可能な Gallipoli らによるモデル 1) を導入し, 非排気 - 非排水条件下でのせん断時における有効性を検討している 2. 試験方法 2.1 試験試料本実験で使用した試料は, 京都府城陽市の木津川右岸堤防で採取した土を 2 mm 以下の粒形にふるい分けしたものである 表 1 に土試料の物理特性の一覧を示す 木津川堤防砂は日本統一分類では SFG に分類される, 細粒分を 26 % 含有した砂質土である 2.2 供試体作成方法試料に最適含水比 1.5% となるように蒸留水を加え, 静的締固め装置を用いて成型した 供試体は高さ 1cm, 直径 5cm の円柱形である 河川堤防の場合, 締固め度を平均で 9% とするという規定があるが, 実際には, 既存堤防の土の締固め度はそれを下回っている事例が少なくない * 京都大学大学院工学研究科 学生 Student., Graduate school of Engineering, Kyoto University ** 鉄道 運輸機構 Japan Railway Construction, Transport and Technology Agency *** 京都大学大学院工学研究科 学生 Student., Graduate school of Engineering, Kyoto University **** 京都大学大学院工学研究科 准教授 Associate professor., Graduate school of Engineering, Kyoto University

2 ため 2), 供試体の締固め度を 85% としている 表 1 土試料の物理特性 土試料の物理特性 土粒子密度 [g/cm 3 ] 最大間隙比 最小間隙比.761 平均粒径.354 [mm] 均等係数 76.7 細粒分含有率 26. 最適含水比 9.2 最大乾燥密度 [g/cm 3 ] 2.3 試験装置の概要 図 1 に試験装置の概略を示す 間隙空気圧を供試体上部 から, 間隙水圧を供試体下部からそれぞれ加えている 供 試体上部にはポリフロンフィルターを, 供試体上部には A.E.V が 1kPa あるいは 3kPa のセラミックディスクを 用いている = + = (1) ここで は全応力, は骨格応力, は平均間隙圧, は飽和度, は間隙水圧, は間隙空気圧, はクロネッカーのデルタである 本研究の結果の整理には骨格応力を用いる 2.6 非排気 - 非排水三軸圧縮試験の三軸圧縮試験の試験条件試験条件と試験手順非排気 - 非排水三軸圧縮試験の試験条件を表 2 に示す 本実験ではセル圧は 3 kpa, 間隙空気圧を 2 kpa とし, 間隙水圧を変化させることで異なるサクションを与え, 試験を行っている 試験手順の概要は次のようである 供試体を設置後セル圧と間隙空気圧の差を 2kPa に保ちながらセル圧 間隙空気圧を上昇させ, 間隙水圧の上昇が収束するのを待つ 間隙水圧を所定の値として排水経路を開け, 1 時間後にセル圧を 3kPa まで上げて排水量が落ち着くまで放置する ( 圧密排水過程 ) その後, コックを閉じ非排気 - 非排水状態にしてせん断を開始する 間隙空気圧 (8) (4) 荷重 (2) (6) (1) (3) (5) (1) セル (2) キャップ (3) 間隙空気圧計 (4) 空気作動式バルブ (5) ポリフロンフィルター (6) 供試体 (7) セラミックディスク (A.E.V=1 or 3 kpa) (8) ギャップセンサー (9) アルミホイル Case 表 2 非排気 - 非排水三軸圧縮試験の試験条件 Initial conditions ( = 3 kpa, = 2 kpa, =.1 %/min) Us Us (9) Us (7) Us 間隙水圧 コック 図 1 試験装置の概要 Us : cell pressure, : pore air pressure, : pore water pressure; : suction, : void ratio, : moisture content, : saturation, : axial strain rate. 2.4 体積ひずみ測定方法本研究では軸変位計と側方に設置した非接触変位計 ( ギャップセンサー ) を用いた方法及び, 試験中に撮影し, 画像処理を施した供試体の写真から計測する方法の 2 通りの方法で体積ひずみを測定している 後者の方法では, 軸ひずみ 2% ごとに供試体の写真を撮影している 2.5 骨格応力 2.7 等方圧縮試験の試験条件と試験手順試験条件と試験手順サクション一定条件つまり, 排気 排水条件下で等方圧縮試験を行った 試験条件を表 3 に示す 試験はサクション 1 kpa, 5 kpa の2ケースについて実施した 初期セル圧を 12kPa とし, セル圧と間隙空気圧の差が 2 kpa 4 kpa 8 kpa 16 kpa 32 kpa 64 kpa となるようにセル圧のみを上昇させた 本研究では応力変数として骨格応力 3) を用いる 骨格応 力は全応力と平均間隙圧と以下の関係にある

3 表 3 排気-排水条件での等方圧縮試験の試験条件 速度.1%/min の非排気-非排水三軸試験の応力径路では初 期サクション kpa のケースでは飽和土の非排水三軸試験 Test conditions ( = 1 kpa) ௪ s1 9 1 s5 5 5 Case と同様な傾向を示しているのに対し 初期サクション ݓ ݏ 1kPa, 5kPa, 8kPa, 2kPa のケースではせん断に伴う平 均骨格応力の減少が小さいかほぼ見られず 初期サクショ ンが大きくなるにつれて非排気-非排水三軸圧縮試験の結 : pore air pressure, ௪ : pore water pressure, : suction, : 果に近づいている void ratio, ݓ : moisture content, ݏ : saturation. 3. 試験結果 3.1 非排気-非排水 非排気 非排水三軸圧縮試験結果 非排水三軸圧縮試験結果 図 2 に応力-ひずみ関係 図 3 に軸ひずみ-体積ひずみ 関係 図 4 に軸ひずみサクション関係を示す なお体積 Deviator stress (kpa) 4 Us Us1 Us5 Us8 Us ひずみは写真撮影によるものを示している 応力-ひずみ 関係では初期サクションが大きいケースほど 強度が増加 ずみ軟化挙動が見られるが その他のケースではひずみ軟 最初期には体積圧縮が終了し その後体積ひずみが一定値 を示している 他のケースでははじめ せん断に伴い体積 4 Volumetric strain (%) ずみ関係では 初期サクション kpa のケースではせん断 Us Us1 Us5 Us8 Us 圧縮する その後初期サクション 1kPa, 8kPa のケースで は体積が一定となる傾向が見られ 初期サクション 5kPa, 25 い 軸ひずみ-サクション関係を見ると初期サクション 2 Suction (kpa) と初期サクションの大小の間には明確な関係は見られな 加し その後一定となる 初期サクション 5kPa, 8kPa の 2 25 図 3 軸ひずみ-体積ひずみ関係 2kPa のケースでは体積膨張を示している 体積変化量 kpa, 1kPa のケースではせん断の初期にサクションが増 25 図 2 軸ひずみ-応力関係 する傾向が見られる 初期サクション kpa では明確なひ 化挙動を示さないかその程度は小さい 軸ひずみ-体積ひ 2 Us Us1 Us5 Us8 Us ケースではせん断初期にサクションが低下した後増加し 5 ていき その後一定値をとっている 初期サクション 5 2kPa のケースではせん断に伴いてサクションが減少し 図 4 軸ひずみ-サクション関係 ている 7 実験とは別に実施した 同じ試料を用いた飽和土の非排水 6 および排水せん断試験結果 CU1,CD1 初期サクシ ョン 1kPa, 5kPa でセル圧 4kPa の非排気-非排水三軸試 験結果 Us1-4 Us5-4 初期サクション 1kPa, 5kPa でひずみ速度.1%/min の非排気-非排水三軸試験結果 Us5-.1 排水排気三軸試験結果 Ds1-.1 Ds5.1 も合わせて示している これらの試験結果より 試 験条件に依らない限界応力比の存在が示唆される ひずみ Deviator stress (kpa) 図 5 に骨格応力を用いて整理した応力径路を示す 本 5 M= CD1 CU1 Us Us1 Us5 Us8 Us2 Us1-4 Us5-4 Us5-.1 Ds1-.1 Ds Mean skeleton stress (kpa) 図 5 骨格応力径路

4 3.2 不飽和等方圧縮試験結果図 6 に等方圧縮試験より得られた間隙比 - 平均骨格応力関係を示す また図中に初期サクション 1kPa, 5kPa での非排気 - 非排水三軸圧縮試験における限界状態のプロットおよび飽和非排水三軸試験での限界状態線 (CSL) を合わせて示している 平均骨格応力が増加するに従い, 圧縮指数が大きくなる 間隙比 - 平均骨格応力の最終的な傾きは非排気 - 非排水三軸圧縮試験から求めた限界状態線の傾きにほぼ等しい また, 初期サクションが大きいほどその傾きが大きくなる傾向も確認され, サクション 1 kpa で λ=.8, サクション 5 kpa で λ=.92 となっている ただし, 現時点では試験回数が限られているため, 今後データの蓄積が必要である.65 かけの相対流速である 同様に気相の質量保存則と固相の 質量保存則から気相の連続式が導かれる , = (3) ここで は空気の密度, は空気の土粒子相に対する 見かけの相対流速である 水相と空気相の見かけの相対流 速, はダルシー則により与えられるとする 4.3 混合体の運動量保存則 強い非線形性を示す問題を扱うために増分型境界値問 題として定式化するため, 混合体の運動量保存則を速度型 に変形して用いる 速度型の運動量保存則は以下で与えら れる = (4) ここで は現在配置における公称応力速度である Void ratio s1 s5 Critical states (s1) Critical states (s5) CSL for saturated soil Mean skeleton stress (kpa) 図 6 間隙比 - 平均骨格応力関係 4.4 不飽和弾粘塑性構成式土骨格に対する構成式として岡らによる弾粘塑性構成式を用いる 4) 本構成式は Perzyna の超過応力型流れ則を用いた弾粘塑性構成式 5) を基本とし, 構成式中の応力変数として骨格応力テンソルを用い, 過圧密降伏面および静的降伏関数にサクションによる硬化 - 軟化を導入することで, 不飽和土の構成式に拡張したものである 4.5 水分特性曲線不飽和土では飽和土とサクションの間に密接な関係があり, その関係性を表す式として van Genuchten 式 6) が広 4. 空気 - 水 - 土連成解析手法 4.1 多相混合体理論地盤を土粒子, 水, 空気の 3 相から構成されるとし, 各相の仮想連続体の重ね合わせとして連続体力学の枠組み く用いられているが, 本研究では間隙比変化の影響を考慮するため Gallipoli ら 1) による方法を導入する Gallipoli らは van Genuchten 式に飽和度の間隙比依存性を導入し, 水分特性曲線を次式で与えている で地盤をモデル化する 岡らは多相混合体理論に基づき不 飽和土の解析手法として空気 - 水 - 土連成解析手法を提案 = 1 + (5) している 4) 以下では本手法について概説する 4.2 水相および空気相の連続式水が非圧縮流体であると仮定して, 土粒子相の質量保存則と水相の質量保存則から水相の連続式が以下のように導かれる + +, = (2) ここで, は土粒子相のストレッチングテンソル, は飽和度, は間隙率, は水相の土粒子相に対する見 = (6) ここで は有効飽和度,,, は材料パラメータ, はそれぞれ最大飽和度と最小飽和度, は = 1 で与えられる 4.6 透水 透気係数の間隙比, 飽和度依存性 透水係数, 透気係数は間隙比依存性および飽和度依存性 を有するとして次式で与えられる

5 = exp 1 1 (7) 表 4 初期条件 Initial conditions for simulation ( = 1, =.1 [%/min]) = exp 1 1 (8) Case σ ここで,, は材料パラメータである 4.7 空気の構成式本研究では空気を理想気体として扱う 理想気体の状態方程式は次式で与えられる = (9) は空気の平均分子量,は気体定数,は温度である 5. 計算条件図 7 に空気 - 水 - 土連成有限要素解析で用いた有限要素メッシュおよび境界条件を示す 円柱供試体の 1/4 をモデル化している 六面体 2 次要素を用いている 変位境界条件は非排気 - 非排水三軸圧縮試験の試験条件に基づいて設定している 底面はすべて固定, 上面は水平方向固定で鉛直方向は変位速度一定, 供試体中心側面は鉛直方向自由で水平方向固定である 試験条件である軸ひずみ速度.1 (%/min) となるようにモデル上面に強制変位を与える 表 4 に初期条件を示す 本研究では, 非排水 - 非排気三軸圧縮試験のうち初期サクション 5kPa のケースについて再現解析を試みた 加えてそれぞれの計算条件に対して水分特性曲線を van Genuchten モデルとした計算と間隙比依存性を考慮した Gallipoli らによるモデル 1) とした計算の両方を行っている 初期応力については等方応力状態としている Us Us is initial pore air pressure; is the axial strain rate; is initial pore water pressure; is initial suction; is initial void ratio; is the coefficient of earth pressure at rest. 表 5 に解析で用いた材料パラメータを示す 圧縮指数 λ は不飽和等方圧縮試験から得られた間隙比 - 平均骨格応力 関係 ( log 曲線 ) より求めた 初期せん断剛性 は 排水排気三軸圧縮試験で得られた応力 - ひずみ曲線の微小 ひずみ領域における勾配から, 変相応力比 は限界状態 線からそれぞれ求めた その他の材料パラメータは, 排水 - 排気三軸圧縮試験結果の要素シミュレーションにより別 途決定した 表 5 材料パラメータ Elasto-viscoplastic model Compression and swelling index,.8,.9 Initial shear modulus Quasi-overconsolidation ratio 1.2 Viscoplastic parameters 39. Viscoplastic parameters, , [1/s] Structural parameter /, β.241, 1. Stress ratio at maximum compression 1.17 Suction parameters,,.3,.25, 5 Parameters for water content characteristic curve Shape parameters, 2., 1.84 Parameters for void ratio dependency, , cm Maximum and minimum saturations 77.6,..2 cm, Parameters for permeability Water and air permeability for full , 図 7 有限要素メッシュと境界条件 saturation, [m/s] Shape parameters,, 3., 2.3,

6 6. 解析結果図 8 に軸ひずみ- 偏差応力関係, 図 9 に骨格応力径路, 図 1 に軸ひずみ- 間隙圧力関係, 図 11 に軸ひずみ-サクション関係, 図 12 に軸ひずみ- 飽和度関係, 図 13 にサクション- 飽和度関係を示す なお, 解析結果は全要素の平均 Genuchten モデルでは表現できない 今回の計算では Gallipolo らによるモデルを用いた場合でもこの挙動は表現できていない 25 2 値を示している まず, 図 8 の軸ひずみ- 偏差応力関係について, 水分特性曲線として従来の van-genuchten モデルを用いた場合と Gallipoli らによるモデルを用いた場合を比較すると, せん断初期には両者の挙動に明確な差は見られない 軸ひずみが 1% を超えるころから若干の差が見ら (kpa) 15 1 Experiment Us5 (GS) 5 FEM (van-genuchten) れ, せん断応力の発生量は Gallipoli らによるモデルの方が 図 8 軸ひずみ - 偏差応力関係 ( 解析結果 ) 小さい これは, 軸ひずみ - サクション関係から分かるよ うに,Gallipoli らのモデルを用いた場合にサクションが若 干小さく, 土骨格の強度が低くなるためと考えられる 次 Experiment Us5 (GS) FEM (van-genuchten) に実験結果の再現性について, せん断応力の発生量を実験結果よりも小さく計算している 図 9 の骨格応力径路ではせん断初期における, 直線的な応力径路を定性的によく表現できている 計算結果ではせん断終了時点で実験から得られた限界応力比に到達していない 図 1 に示す間隙水圧については, せん断初期における上昇を定性的に表現できているが, 間隙水圧の上昇量が実 (kpa) Mean skeleton stress (kpa) 図 9 骨格応力径路 ( 解析結果 ) 験よりも大きく出ている せん断に伴う間隙空気圧の上昇については, 軸ひずみ 5% 程度までよく再現できており, 5% 以上の軸ひずみ発生時においても定性的に表現できている 図 11 の軸ひずみ-サクション関係について, 水分特性曲線として従来の van-genuchten モデルを用いた場合と Gallipoli らによるモデルを用いた場合を比較すると, 軸ひずみ-サクション関係では大きな差は見られない 図 12 の軸ひずみ- 飽和度関係については,Gallipoli らによるモデルを用いた場合で実験結果をより良く再現できている 図 13 のサクション- 飽和度関係については,Gallipoli らによるモデルを用いた場合に飽和度の増加が大きい 実験結果の再現について見ると, せん断初期のサクションの減少は定性的に表現できている また, 軸ひずみ 1% 付近からの, サクション一定で飽和度が増加していく挙動については Gallipoli らによるモデルを用いることでより良く表現できている 軸ひずみ 1.% 付近からのサクションの増加について, サクション- 飽和度関係から飽和度とサクションがともに増加していることが分かる この挙動は van Pore water pressure Experiment Us5 (GS) FEM (Van-Genuchten) Pore air pressure Experiment Us5 (GS) FEM (Van-Genuchten) 図 1 軸ひずみ - 間隙圧関係 ( 解析結果 ) Suction (kpa) Experiment Us5 (GS) FEM (van-genuchten) 図 11 軸ひずみ - サクション関係 ( 解析結果 )

7 6 Experiment Us5 (GS) FEM (Van-Genuchten) 図 12 軸ひずみ - 飽和度関係 ( 解析結果 ) Saturation (%) 6 Experiment Us5 (GS) FEM (van-genuchten) =4. =12. =2. 図 14 偏差骨格応力の第二不変量の分布 Suction (kpa) 図 13 サクション - 飽和度関係 ( 解析結果 ) 次に Gallipoli らのモデルを用いたケースについて, 試験 中の応力やひずみの供試体中の分布を検討する 骨格応力 =4. =12. =2. 図 15 蓄積粘塑性偏差ひずみの分布 の偏差成分の第二不変量 2 の分布を図 14 に示す は じめは供試体内部で均一に値が増加していき, 最終的には 供試体中心からやや上部および下部に応力が局所的に分 布している 蓄積粘塑性偏差ひずみ の分布図を図 15 に示す 蓄 積粘塑性偏差ひずみ は次式で定義され, 非可逆的なせ ん断変形の程度を表す指標である = (1) ここで は粘塑性ストレッチングテンソルの偏差成分 である 蓄積粘塑性偏差ひずみは供試体上下端部にはあま り蓄積せず, 供試体中央部で増加していることが分かる 最終的には供試体中央部から側面へとせん断帯が発生し ている 粘塑性体積ひずみの分布を図 16 に, 飽和度の分布を図 17 に示す 粘塑性体積ひずみが供試体中心からやや上部 および下部に集中している. 飽和度も同じ部分で高くなっ ていることから, 圧縮している部分で飽和度が高くなって いることが分かる =4. =12. =2. 図 16 粘塑性体積ひずみの分布 =4. =12. =2. 図 17 飽和度の分布

8 7. まとめ本研究では不飽和土の力学特性を把握するために木津川堤防砂を用いた非排気 - 非排水三軸圧縮試験および不飽和等方圧縮試験を行った さらに, 間隙比依存性を考慮し 力とひずみが供試体中央部からやや上部および下部の領域に集中していることが示された 今後実験データの蓄積をすすめるとともに, 解析手法の精度向上が必要である た水分特性曲線を導入した空気 - 水 - 土連成解析手法を用 いて実験結果の再現解析を行い, 解析手法の適用性について検討した 加えて, 非排気 - 非排水三軸圧縮試験中の応力やひずみの分布を明らかにした 主な結果は以下の通りである (1) 非排気 - 非排水三軸圧縮試験結果から初期サクションの増加に伴う土の強度増加の傾向が確認された (2) 骨格応力で整理した限界状態における応力比は飽和 不飽和, 非排気 - 非排水, 排気排水, ひずみ速度, 拘束圧等の条件によらず概ね一定の値が得られた (3) 等方圧縮試験から得られた log 曲線では荷重が 増加するに従い圧縮指数が大きくなる log 曲線の最終的な傾きは非排気 - 非排水三軸圧縮試験から求めた限界状態線の傾きにほぼ等しく, 初期サクションが大きいほどその傾きが大きくなる傾向がある (4) 空気 - 水 - 土連成解析手法による再現解析ではせん断中の飽和度増加は Gallipoli らによる水分特性曲線を用いることでより良く表現できている サクションと飽和度がともに増加する挙動については今回の計算では再現できず, 今後の課題である (5) 軸ひずみ- 応力関係, 骨格応力径路は定性的によく再現できている 一方で応力の発生量を過小評価しており, また骨格応力径路が限界状態線に達していないな 参考文献 1) Gallipoli, D., Wheeler, S. J., & Karstunen, M. (23) : Modelling the variation of degree of saturation in a deformable unsaturated soil. Géotechnique., Vol. 53, No.1, pp ) 財団法人国土技術研究センター. (22). 河川堤防の構造検討の手引き. 3) Jommi, C. (2) : Remarks on the constitutive modelling of unsaturated soils. In experimental evidence and theoretical approaches in unsaturated soils, pp , Rotterdam, Balkema. 4) Oka, F., Kodaka, T., Kimoto, S., Kim, Y., & Yamasaki, N. (26) : An elasto-viscoplastic model and multiphase coupled FE analysis for unsaturated soil. In Unsaturated Soils 26, pp ASCE. 5) Adachi, T. and Oka, F. (1982) : Constitutive equations for normally consolidated clay based on elasto-viscoplasticity, Soils and Foundations, Vol. 22, No. 4, pp ) van Genuchten, M. T. (198) : A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil science society of America journal, Vol. 44, No. 5, pp ど, 精度面で課題が見られた (6) 数値解析により非排気 - 非排水三軸圧縮試験中では応 We have conducted a series of triaxial tests for unsaturated sandy soil under fully undrained conditions in order to study the strength and deformation characteristics. In addition, we have conducted isotropic compression tests under constant suction condition and obtained the normally consolidation line by using the skeleton stress as a stress valuable. From the test results, we have found that stress ratio at critical state is almost constant under the various test conditions. The void ratio-skeleton stress relations showed that the slope increases with increasing initial suction. Furthermore, we have simulated the triaxial tests under fully undrained conditions by an air-water-soil coupled method using a soil-water characteristic curve considering the dependency on the void ratio, and discussed the distribution of stress and strain during the test