部分改良地盤の液状化強度の簡易評価法、および部分改良地盤の変形量の簡易評価法
【課題】部分改良地盤の液状化強度と変形量を簡易に評価する。
【解決手段】(i)改良前の地盤に地震時に生じるせん断応力τdを算出し、(ii)前記せん断応力τdに基づいて改良前の地盤のせん断ひずみγiを算出し、(iii)改良後の部分改良地盤の等価せん断剛性Ｇeqを算出し、(iv)改良後の地盤における改良体間の未改良地盤に生じるせん断ひずみγiを求め、(v)該せん断ひずみγiに応じた等価せん断剛性Ｇeqを再決定し、(vi)再決定した等価せん断剛性Ｇeqを用いてせん断ひずみγiが一定の値に収束するまで計算を行ってせん断ひずみγiを決定し、(vii)決定したせん断ひずみγiに基づいて未改良地盤に生じる過剰間隙水圧比Δｕ/σν’を求め、該過剰間隙水圧比Δｕ/σν’に基づいて部分改良地盤の液状化強度を評価する。さらに、(viii)改良体に作用する外力を算定し、(ix)前記等価せん断剛性Ｇeqおよび外力に基づいて部分改良地盤の変形量を評価する。
【解決手段】(i)改良前の地盤に地震時に生じるせん断応力τdを算出し、(ii)前記せん断応力τdに基づいて改良前の地盤のせん断ひずみγiを算出し、(iii)改良後の部分改良地盤の等価せん断剛性Ｇeqを算出し、(iv)改良後の地盤における改良体間の未改良地盤に生じるせん断ひずみγiを求め、(v)該せん断ひずみγiに応じた等価せん断剛性Ｇeqを再決定し、(vi)再決定した等価せん断剛性Ｇeqを用いてせん断ひずみγiが一定の値に収束するまで計算を行ってせん断ひずみγiを決定し、(vii)決定したせん断ひずみγiに基づいて未改良地盤に生じる過剰間隙水圧比Δｕ/σν’を求め、該過剰間隙水圧比Δｕ/σν’に基づいて部分改良地盤の液状化強度を評価する。さらに、(viii)改良体に作用する外力を算定し、(ix)前記等価せん断剛性Ｇeqおよび外力に基づいて部分改良地盤の変形量を評価する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、バットレス状改良体や格子状改良体もしくは柱状改良体が造成された部分改良地盤を対象として、その部分改良地盤の液状化強度および変形量を簡易に評価する方法に関する。
【背景技術】
【０００２】
液状化対策や液状化に伴う地盤の側方流動対策として、バットレス状改良や格子状改良あるいは柱状改良を行う場合がある。そのような（全面改良ではない）部分改良地盤については、改良体により拘束された自然地盤（未改良地盤）の液状化の可能性や過剰間隙水圧の上昇度の評価が重要であり、たとえば格子状の部分改良地盤についての液状化強度の評価法は特許文献１で提案されている。
【特許文献１】特開２００１−３５５２２９号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
特許文献１に示されている方法は、評価法自体の論理的根拠が乏しく、格子状改良以外の条件も含めて、対象地盤や改良土自体に対する適用範囲を明確にできないという問題がある。
【０００４】
上記事情に鑑み、本発明は、部分改良地盤の液状化強度と地震時の変形量を簡易に評価し得る有効適切な評価法を提供することを目的とする。
【課題を解決するための手段】
【０００５】
本発明の部分改良地盤の液状化強度の簡易評価法は、バットレス状改良体や格子状改良体もしくは柱状改良体が造成された部分改良地盤の液状化強度を簡易に評価する方法であって、改良前の地盤に地震時に生じるせん断応力τdを算出し、前記せん断応力τdに基づいて改良前の地盤のせん断ひずみγiを算出し、改良後の部分改良地盤の等価せん断剛性Ｇeqを算出し、改良後の地盤における改良体間の未改良地盤に生じるせん断ひずみγiを求め、該せん断ひずみγiに応じた等価せん断剛性Ｇeqを再決定し、再決定した等価せん断剛性Ｇeqを用いてせん断ひずみγiが一定の値に収束するまで計算を行ってせん断ひずみγiを決定し、決定したせん断ひずみγiに基づいて未改良地盤に生じる過剰間隙水圧比Δｕ/σν’を求め、該過剰間隙水圧比Δｕ/σν’に基づいて部分改良地盤の液状化強度を評価することを特徴とする。
【０００６】
本発明の部分改良地盤の変形量の簡易評価法は、バットレス状改良体や格子状改良体もしくは柱状改良体が造成された部分改良地盤の変形量を簡易に評価する方法であって、改良前の地盤に地震時に生じるせん断応力τdを算出し、前記せん断応力τdに基づいて改良前の地盤のせん断ひずみγiを算出し、改良後の部分改良地盤の等価せん断剛性Ｇeqを算出し、改良後の地盤における改良体間の未改良地盤に生じるせん断ひずみγiを求め、該せん断ひずみγiに応じた等価せん断剛性Ｇeqを再決定し、再決定した等価せん断剛性Ｇeqを用いてせん断ひずみγiが一定の値に収束するまで計算を行ってせん断ひずみγiを決定し、決定したせん断ひずみγiに基づいて未改良地盤に生じる過剰間隙水圧比Δｕ/σν’を求め、改良体に作用する外力を算定し、前記等価せん断剛性Ｇeqおよび前記外力に基づいて部分改良地盤の変形量を評価することを特徴とする。
【発明の効果】
【０００７】
本発明によれば、液状化に関する既往の要素試験から得られた知見と、部分改良地盤に対する等価せん断剛性を求める手法から、任意の拘束条件下での等価地盤せん断モデルによって地盤の過剰間隙水圧比を求めることにより、様々な形態の部分改良地盤の液状化強度および変形量を簡易にかつ統一的に評価することができる。
【発明を実施するための最良の形態】
【０００８】
以下、本発明の部分改良地盤の変形量の評価法の実施形態について説明するが、これは本発明者らが先に特願２００８−０５６２７２において提供した「部分改良地盤の変形量の簡易算定法」を基礎とするものであるので、まずその算定法について図１〜図８を参照して説明する。
【０００９】
本算定法は、部分改良地盤を単位周期構造体の集合体として築造し、その単位周期構造体の剛性を数学的均質化理論に基づいて求めてこれを部分改良地盤全体の等価剛性とし、この等価剛性と築造後の部分改良地盤に作用する外力とから部分改良地盤の変形量を事前に予測算定するものである。
そして、そのために、本算定法においては、部分改良地盤の等価剛性と改良率との関係を、単位周期構造体の縦横比と改良体の剛性と原地盤の剛性とをパラメータとして予め求めて簡易チャート化しておき、その簡易チャートを用いて部分改良地盤の等価剛性を簡易に求めることを主眼とする。
【００１０】
図１は本算定方の概要を示すものであり、（ａ）は柱状の改良体による部分改良地盤の場合、（ｂ）は格子状の改良体による部分改良地盤の場合である。
（ａ）は直径lRの柱状の改良体をｘ方向（幅方向）の間隔lx、ｙ方向（長さ方向）の間隔lyをもって配列して造成した場合であって、この場合の単位周期構造体は２辺がlx,lyの矩形（lx＝lyの場合には正方形）の中心に直径lRの柱状の改良体が配置されたものである。この場合の改良率R、すなわち単位周期構造体の面積に対する改良体の面積の比、つまりは部分改良地盤全体の面積に対する改良体全体の面積の比は
【数１】
となる。
【００１１】
（ｂ）は厚さlRの板状の改良体をｘ方向とｙ方向の双方にそれぞれlx、lyの間隔で格子状に組み合わせて造成した場合であって、この場合の単位周期構造体はｘ方向の長さlx、ｙ方向の長さlyの矩形の内側に、厚さlRの改良体がｘ方向およびｙ方向に１枚ずつ、あるいは厚さlRの改良体がｘ方向に1枚、厚さlR/2の改良体がｙ方向に２枚、配置されたものである。この場合の改良率Rはいずれも
【数２】
となる。
【００１２】
そして、単位周期構造体は剛性の異なる２つの弾性体、すなわち未改良で低剛性の原地盤（その剛性をESとする）と、改良により高剛性とされた改良体（その剛性をERとする）の複合体と見なすことができ、さらにその複合体は、この複合体全体の剛性と等価とみなせる剛性（以下、これを等価剛性EHという）を有する単一の均質体と見なすことができ、その均質体の等価剛性EHは、数学的均質化理論に基づき単位周期構造体の特性と形状とをパラメータとして次のように求めることができる。
【００１３】
すなわち、数学的均質化理論によれば、２つの弾性体の複合体と等価の１つの均質体の弾性係数CHをマトリックスを表記すると、次式で表される。
【数３】
【００１４】
上式においてCはミクロ周期構造としての単位周期構造体の弾性マトリックスである。また、Ｘはミクロ周期構造に単位マクロ歪みIを与えた場合の応答変位であり、３次元では次式のように６成分からなるものである。
【数４】
【００１５】
また、均質体の弾性係数CHの逆行列（コンプライアンスマトリックス）は次式で表され、この式から各方向の等価剛性を求めることができる。
【数５】
【００１６】
上式におけるExHは均質体のｘ方向の軸剛性、EyHはｙ方向の軸剛性、EzHはｚ方向の軸剛性、GxyHはｘ−ｙ面内のせん断剛性、GyzHはｙ−ｚ面内のせん断剛性、GzxHはｚ−ｘ面（ｘ−ｚ面）内のせん断剛性であり、上式により求められる均質体の各剛性はすなわち単位周期構造体およびその集合体としての部分改良地盤全体の等価剛性を表すものである。
そして、本算定法においては、上式で求められる各方向の等価剛性と改良率Rとの関係を、単位周期構造体のパターンをパラメータとして予め簡易チャート化しておくことにより、その簡易チャートを用いて部分改良地盤の各方向の等価剛性（軸剛性およびせん断剛性）を簡易に求めるものである。
【００１７】
その簡易チャートは、具体例を図２〜図３に示すように、横軸に改良率Rをとり、縦軸に等価剛性ExH、EyH、EzH、GxyH、GyzH、GzxH（改良体の軸剛性ERあるいはせん断剛性GRにより除して正規化してある）をとり、ES/ERあるいはGS/GR（改良体の剛性に対する原地盤の剛性の比）と、単位周期構造体の縦横比lx/lyをパラメータとして作成したものである。
図２（ａ）〜（ｄ）に示す簡易チャートは柱状の改良体の場合において、lx/ly＝１（つまり正方形配列）の場合のものであり、図３（ａ）〜（ｆ）に示す簡易チャートは格子状の改良体の場合において、lx/ly＝1,2,3とした場合のものあり、いずれもパラメータES/ERあるいはGS/GRを0、0.2、0.4、0.6、0.8としたものである。
【００１８】
このような簡易チャートを予め作成しておくことにより、改良率Rと、改良体の剛性ERあるいはGR、原地盤の剛性ESあるいはGSのみから、部分改良地盤の各方向の等価剛性を直ちに求めることができる。たとえば、柱状の改良体を正方形配列して造成する部分改良地盤における地盤厚さ方向（ｚ方向）の軸剛性を求める場合には、図２（ｂ）の簡易チャートを用いてEzHを求めれば良い。また、格子状の改良体による部分改良地盤における水平面（ｘ−ｙ面）におけるせん断剛性を求める場合には、図３（ｆ）のチャートを用いてGxyHを求めれば良い。
【００１９】
本算定法の具体的な適用例として、液状化が想定される護岸に対して側方流動防止対策として地盤改良を行う場合において、液状化が生じた後における護岸の側方流動量を予測算定する場合の一具体例を説明する。
図４はその全体作業手順の概要を示すフローチャートである。
まず、（１）地盤調査を行った結果から液状化層の深度（液状化層厚）を決定し、それに基づき、（２）液状化後に部分改良地盤に作用する外力を決定する。
一方、（３）改良体の形状（柱状あるいは格子状）とその剛性、単位周期構造体の縦横比（lx/ly）を決定する。また、改良率を仮決定し、その改良率に対応する液状化後の等価剛性を簡易チャートにより求める。
そして、（４）等価剛性と外力とにより変形量を算定し、設計条件を満足しなければ改良率を修正して条件を満足するまで以上の手順を繰り返す。なお、その際に必要であれば改良率の修正に併せて、あるいはそれに代えて、他の条件（単位周期構造体のパターンや改良体の剛性等）の見直しを行っても良い。
以上を設計条件を満足するまで繰り返し、条件を満足すれば改良率を確定させることにより、（５）側方流動対策の決定とする。
【００２０】
以上の基本的な手順（１）〜（５）を図５に示す構造の護岸に適用して、その液状化後に生じることが予想される変形量を算定する場合の具体例を以下に示す。
（１）液状化層厚の決定
図５において層厚H4で示す砂層が液状化対策が必要な層であり、その範囲を地盤改良するとする。図示例の場合には格子状の改良体の造成による部分地盤改良とし、単位周期構造体のlx＝15m、ly＝5mとし、したがってlx/ly＝3とする。
【００２１】
（２）外力の決定
改良体の前面および背面の液状化層が全て液状化するとして、液状化後に地盤改良体にかかる外力を算定する。
護岸前面側（水面を原点とする座標Ｚ1として示す）においては、Ｚ1＝0〜H1までは河川あるいは海の水圧がかかり、Ｚ1＝H1〜HI+H2までは水圧＋土圧がかかる。すなわち、地盤の単位体積重量γt、水の単位体積重量γWとすると、
P1＝γwＺ1 at Ｚ1＝0〜H1
P1＝γwＺ1＋γt（Ｚ1−H1） at Ｚ1＝H1〜H1＋H2
護岸背面側（地下水位を原点とする座標Ｚ2として示す）においては、Ｚ2＝0において土被り圧がかかり、Ｚ2＝0〜H4は土被り圧＋地下水圧がかかる。すなわち
P2＝γtH3 at Ｚ2＝0
P2＝γtH3＋γtＺ2 at Ｚ2＝0〜H4
【００２２】
（３）改良体の決定、等価剛性の算定
上記の外力により地盤改良体にはその横断面に沿う鉛直面（ｚ−ｘ面）に沿ってせん断変形が生じるので、ここでは部分改良地盤のその方向の等価せん断剛性GzxHを図３（ｄ）に示した簡易チャートにより求める。この場合、改良体のせん断剛性GRをセメント系地盤改良を想定してGR＝38.5MN/m2とし、原地盤のせん断剛性GSは液状化によりGS＝0.1MN/m2になると想定して、簡易チャートにおけるパラメータGS/GR＝0とする。
そして、改良率Rの仮決定をR＝50％とすれば、単位周期構造体の縦横比lx/ly＝3であるから、図６に示すように簡易チャートの縦軸GzxH／GR＝0.43、ゆえに、等価剛性GzxH＝38.5×0.43＝16.5MN／m2として求められる。
【００２３】
（４）変形量の算定
上記の外力を負荷した際の部分改良地盤の変形量を算定する。その算定は２次元弾性有限要素法によるか、あるいは、より簡易な手法として、部分改良地盤全体をせん断棒にモデル化することにより行うことができ、いずれの場合も以下に示すようにほぼ同様の結果が得られる。
（４−１）２次元弾性有限要素法による場合
図５に示している各諸元、H1＝3m、H2＝7m、H3＝2m、H4＝10m、lx＝15m、ly＝5m、γt＝17kN/m3、γw＝10kN/m3を用いて、２次元弾性有限要素法により変形量を算定する。その結果、図７に示すように最大変形量が3cmと算定された。
（４−２）せん断棒による場合
図８に示すように、部分改良地盤全体をその等価剛性と同等の剛性を有するせん断棒にモデル化し、そのせん断棒に上記の外力を作用させた際に生じる変形量を算定する。この場合
改良体背面側の外力 P2＝(34+204)×10/2＝1190kN/m
改良体前面側の外力 P1＝30×3/2＋(30+149)×7/2＝671.5kN/m
外力の合力 P＝P2-P1＝1190-671.5＝518.5kN/m
等価せん断剛性 GzxH＝16.5MN／m2＝16.5×103kN／m2
等価せん断バネ k＝GzxHW／l＝16.5×103×15/10＝24.8×103kN／m2
水平方向変形量 ｘ＝P/k＝518.5/24.8×103＝0.02m＝2cm
この場合の算定結果は、２次元弾性有限要素法による場合の算定結果に比べて若干の誤差があるものの、この種の解析においては両者の結果は実質的に同等であるといえるし、少なくとも評価結果に影響しない範囲内の誤差であるといえる。
【００２４】
（５）対策の決定
以上で算定された変形量が設計条件を満足すれば、改良率の仮決定（上記の場合はR＝50％）が妥当であったのでその改良率を最終決定として対策決定とする。変形量の算定結果が設計条件を満足しなければ改良率を変更して以上の手順を繰り返す。すなわち、変形量が過大であれば改良不足であるので改良率を大きくするように変更し、変形量が過小であれば改良過剰であるので改良率を小さくするように変更し、満足すべき結果が得られるまで以上の手順を繰り返せば良い。勿論、その際に必要であれば、すなわち改良率の修正のみでは条件を満足できない場合には、単位周期構造体のパターンや改良体の剛性も併せて見直せば良い。
【００２５】
以上で説明したように、本算定法によれば部分改良地盤の変形量の算定に際してまずその等価剛性を簡易チャートにより求めることにより、従来の３次元有限要素法により変形量を直接算定する場合に比べて全体の算定手順を格段に簡略化することができ、特に厳密なモデル化を行うための手間を大きく軽減することができる。したがって、従来のように多大の手間とコスト、時間を要することなく部分改良地盤全体の変形量算定を簡易に実施することができ、その結果、部分改良地盤による改良効果を事前に確認できてその信頼性を充分に高めることができ、特に上記のような護岸の液状化による側方流動防止対策として部分改良地盤を築造するに際して適用して好適である。
【００２６】
以上で本発明の基礎となった算定法について説明したが、次に本発明の部分改良地盤の変形量の評価法についての実施形態を図９〜図１３を参照して説明する。
本実施形態の評価法は、図９にフローチャートとして示した一連のステップ（ｉ）〜
（ix）により構成されるものであり、以下にその詳細を説明する。
【００２７】
（i）地震時に地盤に生じるせん断応力τdを算出する。せん断応力τdは地震応答解析から算定するか、簡易液状化算定法に用いられている（１）式を用いて求める。
【数６】
（１）式において、Ｍ：マグニチュード、αmax：地表面最大加速度、ｚ：深度、σν'：鉛直有効応力、σν：鉛直全応力、ｇ：重力加速度である。
【００２８】
（ii)上記の（１）式で求めたせん断応力τdと、図１１を用いて、改良前の地盤のせん断ひずみγiを算出する。図１１は（２）式をもとに作成したものである。
【数７】
（２式）において、ｈmaxは最大減衰定数で砂の場合は0.2〜0.25の値となる。初期せん断剛性Ｇsoや基準せん断ひずみγrfは砂の種類や拘束圧によって異なるが、典型的な例としてＧso＝50MPa、γrf＝0.1％、ｈmax＝0.2とした。
【００２９】
（iii）上述した算定法で説明した手法により、図３を用いて部分改良地盤の改良率、改良形状、改良パターン、および地盤と改良体のせん断剛性Ｇso、ＧIに応じた等価せん断剛性Ｇeqを算出する。
【００３０】
（iv）上記(i)で求めたせん断応力τdと、上記(iii)で仮算出した等価せん断剛性Ｇeqを用いて、（３）式により改良体間地盤に生じるせん断ひずみγiを求める。
【数８】
【００３１】
（v）上記(iv)で求めたせん断ひずみγiと図１０を用いて、せん断ひずみγiに応じた等価せん断剛性Ｇeqを再決定する。
【００３２】
（vi）再決定した等価せん断剛性Ｇeqと、（１）式で求めたτdを用いて、（３）式により再度せん断ひずみγiを計算し、以上の計算をせん断ひずみγiが一定の値に収束するまで行う。
【００３３】
（vii）決定したせん断ひずみγiを用いて、図１２により過剰間隙水圧比Δｕ/σν’を決定する。
図１２は（２）式と次の（４）〜（６）式を用いて求めたものである。
【数９】
ここで、Ｒlは繰返せん断応力比、Ｒ20は液状化強度で２０回で液状化に至るせん断応力比として規定している。Ｎlは液状化に至った繰返し回数である。ｋは実験定数で-0.25程度の値をとる。
図１３は液状化強度曲線の模式図を示したものである。
（４）式はＲl１回の繰返しで １／Ｎl だけ液状化に近づいたとみなせるので、２０回の繰返しせん断が生じた際の累積損傷度Ｒnは（５）式のように表せる。この状態で生じる過剰間隙水圧比は、De Alba の提案式である（６）式で表すことができる。αrfは実験定数であるが、緩い砂では0.7程度の値をとる。
【００３４】
(viii)液状化時に改良体に作用する外力を、たとえば上述した算定法に示した手法により算定する。
【００３５】
(ix)上記の外力が作用したときの改良体の変形量を、たとえば上述した算定法に示した手法（等価せん断棒や２次元平面ひずみ解析）を用いて算定する。
【００３６】
なお、上記のステップ(vii)において図１２を用いて過剰間隙水圧比Δｕ/σν’を求めることにより、その過剰間隙水圧比Δｕ/σν’が1.0未満であれば地盤は液状化に至らず過剰間隙水圧がある程度上昇するに留まることになり、そのことから部分改良地盤の液状化強度を評価することができる。そこで、仮に過剰間隙水圧比Δｕ/σν’が1.0以上であって地盤が液状化すると評価された場合には、改良率や改良パターンを再設定し、過剰間隙水圧比Δｕ/σν’が1.0未満になるまで以上の手順を繰り返せば良い。
したがって、以上で説明した本発明の部分改良地盤の変形量の評価法は、その過程で過剰間隙水圧比Δｕ/σν’を求めることから、部分改良地盤の液状化強度を評価するための評価法としても適用できるものである。
【図面の簡単な説明】
【００３７】
【図１】本発明の基礎となった部分改良地盤の変形量の簡易算定法の概要を示す説明図である。
【図２】同、簡易チャートの例（柱状の改良体の場合）を示す図である。
【図３】同、簡易チャートの例（格子状の改良体の場合）を示す図である。
【図４】同、護岸を対象とする側方流動防止対策としての部分改良地盤を築造する場合への適用例を示すもので、全体手順を示すフローチャートである。
【図５】同、護岸の構造図である。
【図６】同、簡易チャートにより等価剛性を求める手順を示す図である。
【図７】同、２次元弾性有限要素法による変形量の算定結果を示す図である。
【図８】同、等価せん断棒による変形量の算定手法を示す図である。
【図９】本発明の部分改良地盤の液状化強度の簡易評価法および変形量の簡易評価法の実施形態を示すもので、その全体手順を示すフローチャートである。
【図１０】同、せん断剛性低下率とせん断ひずみの関係を示す図である。
【図１１】同、せん断ひずみを求めるための図である。
【図１２】同、過剰間隙水圧比を求めるための図である。
【図１３】同、液状化強度曲線の模式図である。
【技術分野】
【０００１】
本発明は、バットレス状改良体や格子状改良体もしくは柱状改良体が造成された部分改良地盤を対象として、その部分改良地盤の液状化強度および変形量を簡易に評価する方法に関する。
【背景技術】
【０００２】
液状化対策や液状化に伴う地盤の側方流動対策として、バットレス状改良や格子状改良あるいは柱状改良を行う場合がある。そのような（全面改良ではない）部分改良地盤については、改良体により拘束された自然地盤（未改良地盤）の液状化の可能性や過剰間隙水圧の上昇度の評価が重要であり、たとえば格子状の部分改良地盤についての液状化強度の評価法は特許文献１で提案されている。
【特許文献１】特開２００１−３５５２２９号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
特許文献１に示されている方法は、評価法自体の論理的根拠が乏しく、格子状改良以外の条件も含めて、対象地盤や改良土自体に対する適用範囲を明確にできないという問題がある。
【０００４】
上記事情に鑑み、本発明は、部分改良地盤の液状化強度と地震時の変形量を簡易に評価し得る有効適切な評価法を提供することを目的とする。
【課題を解決するための手段】
【０００５】
本発明の部分改良地盤の液状化強度の簡易評価法は、バットレス状改良体や格子状改良体もしくは柱状改良体が造成された部分改良地盤の液状化強度を簡易に評価する方法であって、改良前の地盤に地震時に生じるせん断応力τdを算出し、前記せん断応力τdに基づいて改良前の地盤のせん断ひずみγiを算出し、改良後の部分改良地盤の等価せん断剛性Ｇeqを算出し、改良後の地盤における改良体間の未改良地盤に生じるせん断ひずみγiを求め、該せん断ひずみγiに応じた等価せん断剛性Ｇeqを再決定し、再決定した等価せん断剛性Ｇeqを用いてせん断ひずみγiが一定の値に収束するまで計算を行ってせん断ひずみγiを決定し、決定したせん断ひずみγiに基づいて未改良地盤に生じる過剰間隙水圧比Δｕ/σν’を求め、該過剰間隙水圧比Δｕ/σν’に基づいて部分改良地盤の液状化強度を評価することを特徴とする。
【０００６】
本発明の部分改良地盤の変形量の簡易評価法は、バットレス状改良体や格子状改良体もしくは柱状改良体が造成された部分改良地盤の変形量を簡易に評価する方法であって、改良前の地盤に地震時に生じるせん断応力τdを算出し、前記せん断応力τdに基づいて改良前の地盤のせん断ひずみγiを算出し、改良後の部分改良地盤の等価せん断剛性Ｇeqを算出し、改良後の地盤における改良体間の未改良地盤に生じるせん断ひずみγiを求め、該せん断ひずみγiに応じた等価せん断剛性Ｇeqを再決定し、再決定した等価せん断剛性Ｇeqを用いてせん断ひずみγiが一定の値に収束するまで計算を行ってせん断ひずみγiを決定し、決定したせん断ひずみγiに基づいて未改良地盤に生じる過剰間隙水圧比Δｕ/σν’を求め、改良体に作用する外力を算定し、前記等価せん断剛性Ｇeqおよび前記外力に基づいて部分改良地盤の変形量を評価することを特徴とする。
【発明の効果】
【０００７】
本発明によれば、液状化に関する既往の要素試験から得られた知見と、部分改良地盤に対する等価せん断剛性を求める手法から、任意の拘束条件下での等価地盤せん断モデルによって地盤の過剰間隙水圧比を求めることにより、様々な形態の部分改良地盤の液状化強度および変形量を簡易にかつ統一的に評価することができる。
【発明を実施するための最良の形態】
【０００８】
以下、本発明の部分改良地盤の変形量の評価法の実施形態について説明するが、これは本発明者らが先に特願２００８−０５６２７２において提供した「部分改良地盤の変形量の簡易算定法」を基礎とするものであるので、まずその算定法について図１〜図８を参照して説明する。
【０００９】
本算定法は、部分改良地盤を単位周期構造体の集合体として築造し、その単位周期構造体の剛性を数学的均質化理論に基づいて求めてこれを部分改良地盤全体の等価剛性とし、この等価剛性と築造後の部分改良地盤に作用する外力とから部分改良地盤の変形量を事前に予測算定するものである。
そして、そのために、本算定法においては、部分改良地盤の等価剛性と改良率との関係を、単位周期構造体の縦横比と改良体の剛性と原地盤の剛性とをパラメータとして予め求めて簡易チャート化しておき、その簡易チャートを用いて部分改良地盤の等価剛性を簡易に求めることを主眼とする。
【００１０】
図１は本算定方の概要を示すものであり、（ａ）は柱状の改良体による部分改良地盤の場合、（ｂ）は格子状の改良体による部分改良地盤の場合である。
（ａ）は直径lRの柱状の改良体をｘ方向（幅方向）の間隔lx、ｙ方向（長さ方向）の間隔lyをもって配列して造成した場合であって、この場合の単位周期構造体は２辺がlx,lyの矩形（lx＝lyの場合には正方形）の中心に直径lRの柱状の改良体が配置されたものである。この場合の改良率R、すなわち単位周期構造体の面積に対する改良体の面積の比、つまりは部分改良地盤全体の面積に対する改良体全体の面積の比は
【数１】
となる。
【００１１】
（ｂ）は厚さlRの板状の改良体をｘ方向とｙ方向の双方にそれぞれlx、lyの間隔で格子状に組み合わせて造成した場合であって、この場合の単位周期構造体はｘ方向の長さlx、ｙ方向の長さlyの矩形の内側に、厚さlRの改良体がｘ方向およびｙ方向に１枚ずつ、あるいは厚さlRの改良体がｘ方向に1枚、厚さlR/2の改良体がｙ方向に２枚、配置されたものである。この場合の改良率Rはいずれも
【数２】
となる。
【００１２】
そして、単位周期構造体は剛性の異なる２つの弾性体、すなわち未改良で低剛性の原地盤（その剛性をESとする）と、改良により高剛性とされた改良体（その剛性をERとする）の複合体と見なすことができ、さらにその複合体は、この複合体全体の剛性と等価とみなせる剛性（以下、これを等価剛性EHという）を有する単一の均質体と見なすことができ、その均質体の等価剛性EHは、数学的均質化理論に基づき単位周期構造体の特性と形状とをパラメータとして次のように求めることができる。
【００１３】
すなわち、数学的均質化理論によれば、２つの弾性体の複合体と等価の１つの均質体の弾性係数CHをマトリックスを表記すると、次式で表される。
【数３】
【００１４】
上式においてCはミクロ周期構造としての単位周期構造体の弾性マトリックスである。また、Ｘはミクロ周期構造に単位マクロ歪みIを与えた場合の応答変位であり、３次元では次式のように６成分からなるものである。
【数４】
【００１５】
また、均質体の弾性係数CHの逆行列（コンプライアンスマトリックス）は次式で表され、この式から各方向の等価剛性を求めることができる。
【数５】
【００１６】
上式におけるExHは均質体のｘ方向の軸剛性、EyHはｙ方向の軸剛性、EzHはｚ方向の軸剛性、GxyHはｘ−ｙ面内のせん断剛性、GyzHはｙ−ｚ面内のせん断剛性、GzxHはｚ−ｘ面（ｘ−ｚ面）内のせん断剛性であり、上式により求められる均質体の各剛性はすなわち単位周期構造体およびその集合体としての部分改良地盤全体の等価剛性を表すものである。
そして、本算定法においては、上式で求められる各方向の等価剛性と改良率Rとの関係を、単位周期構造体のパターンをパラメータとして予め簡易チャート化しておくことにより、その簡易チャートを用いて部分改良地盤の各方向の等価剛性（軸剛性およびせん断剛性）を簡易に求めるものである。
【００１７】
その簡易チャートは、具体例を図２〜図３に示すように、横軸に改良率Rをとり、縦軸に等価剛性ExH、EyH、EzH、GxyH、GyzH、GzxH（改良体の軸剛性ERあるいはせん断剛性GRにより除して正規化してある）をとり、ES/ERあるいはGS/GR（改良体の剛性に対する原地盤の剛性の比）と、単位周期構造体の縦横比lx/lyをパラメータとして作成したものである。
図２（ａ）〜（ｄ）に示す簡易チャートは柱状の改良体の場合において、lx/ly＝１（つまり正方形配列）の場合のものであり、図３（ａ）〜（ｆ）に示す簡易チャートは格子状の改良体の場合において、lx/ly＝1,2,3とした場合のものあり、いずれもパラメータES/ERあるいはGS/GRを0、0.2、0.4、0.6、0.8としたものである。
【００１８】
このような簡易チャートを予め作成しておくことにより、改良率Rと、改良体の剛性ERあるいはGR、原地盤の剛性ESあるいはGSのみから、部分改良地盤の各方向の等価剛性を直ちに求めることができる。たとえば、柱状の改良体を正方形配列して造成する部分改良地盤における地盤厚さ方向（ｚ方向）の軸剛性を求める場合には、図２（ｂ）の簡易チャートを用いてEzHを求めれば良い。また、格子状の改良体による部分改良地盤における水平面（ｘ−ｙ面）におけるせん断剛性を求める場合には、図３（ｆ）のチャートを用いてGxyHを求めれば良い。
【００１９】
本算定法の具体的な適用例として、液状化が想定される護岸に対して側方流動防止対策として地盤改良を行う場合において、液状化が生じた後における護岸の側方流動量を予測算定する場合の一具体例を説明する。
図４はその全体作業手順の概要を示すフローチャートである。
まず、（１）地盤調査を行った結果から液状化層の深度（液状化層厚）を決定し、それに基づき、（２）液状化後に部分改良地盤に作用する外力を決定する。
一方、（３）改良体の形状（柱状あるいは格子状）とその剛性、単位周期構造体の縦横比（lx/ly）を決定する。また、改良率を仮決定し、その改良率に対応する液状化後の等価剛性を簡易チャートにより求める。
そして、（４）等価剛性と外力とにより変形量を算定し、設計条件を満足しなければ改良率を修正して条件を満足するまで以上の手順を繰り返す。なお、その際に必要であれば改良率の修正に併せて、あるいはそれに代えて、他の条件（単位周期構造体のパターンや改良体の剛性等）の見直しを行っても良い。
以上を設計条件を満足するまで繰り返し、条件を満足すれば改良率を確定させることにより、（５）側方流動対策の決定とする。
【００２０】
以上の基本的な手順（１）〜（５）を図５に示す構造の護岸に適用して、その液状化後に生じることが予想される変形量を算定する場合の具体例を以下に示す。
（１）液状化層厚の決定
図５において層厚H4で示す砂層が液状化対策が必要な層であり、その範囲を地盤改良するとする。図示例の場合には格子状の改良体の造成による部分地盤改良とし、単位周期構造体のlx＝15m、ly＝5mとし、したがってlx/ly＝3とする。
【００２１】
（２）外力の決定
改良体の前面および背面の液状化層が全て液状化するとして、液状化後に地盤改良体にかかる外力を算定する。
護岸前面側（水面を原点とする座標Ｚ1として示す）においては、Ｚ1＝0〜H1までは河川あるいは海の水圧がかかり、Ｚ1＝H1〜HI+H2までは水圧＋土圧がかかる。すなわち、地盤の単位体積重量γt、水の単位体積重量γWとすると、
P1＝γwＺ1 at Ｚ1＝0〜H1
P1＝γwＺ1＋γt（Ｚ1−H1） at Ｚ1＝H1〜H1＋H2
護岸背面側（地下水位を原点とする座標Ｚ2として示す）においては、Ｚ2＝0において土被り圧がかかり、Ｚ2＝0〜H4は土被り圧＋地下水圧がかかる。すなわち
P2＝γtH3 at Ｚ2＝0
P2＝γtH3＋γtＺ2 at Ｚ2＝0〜H4
【００２２】
（３）改良体の決定、等価剛性の算定
上記の外力により地盤改良体にはその横断面に沿う鉛直面（ｚ−ｘ面）に沿ってせん断変形が生じるので、ここでは部分改良地盤のその方向の等価せん断剛性GzxHを図３（ｄ）に示した簡易チャートにより求める。この場合、改良体のせん断剛性GRをセメント系地盤改良を想定してGR＝38.5MN/m2とし、原地盤のせん断剛性GSは液状化によりGS＝0.1MN/m2になると想定して、簡易チャートにおけるパラメータGS/GR＝0とする。
そして、改良率Rの仮決定をR＝50％とすれば、単位周期構造体の縦横比lx/ly＝3であるから、図６に示すように簡易チャートの縦軸GzxH／GR＝0.43、ゆえに、等価剛性GzxH＝38.5×0.43＝16.5MN／m2として求められる。
【００２３】
（４）変形量の算定
上記の外力を負荷した際の部分改良地盤の変形量を算定する。その算定は２次元弾性有限要素法によるか、あるいは、より簡易な手法として、部分改良地盤全体をせん断棒にモデル化することにより行うことができ、いずれの場合も以下に示すようにほぼ同様の結果が得られる。
（４−１）２次元弾性有限要素法による場合
図５に示している各諸元、H1＝3m、H2＝7m、H3＝2m、H4＝10m、lx＝15m、ly＝5m、γt＝17kN/m3、γw＝10kN/m3を用いて、２次元弾性有限要素法により変形量を算定する。その結果、図７に示すように最大変形量が3cmと算定された。
（４−２）せん断棒による場合
図８に示すように、部分改良地盤全体をその等価剛性と同等の剛性を有するせん断棒にモデル化し、そのせん断棒に上記の外力を作用させた際に生じる変形量を算定する。この場合
改良体背面側の外力 P2＝(34+204)×10/2＝1190kN/m
改良体前面側の外力 P1＝30×3/2＋(30+149)×7/2＝671.5kN/m
外力の合力 P＝P2-P1＝1190-671.5＝518.5kN/m
等価せん断剛性 GzxH＝16.5MN／m2＝16.5×103kN／m2
等価せん断バネ k＝GzxHW／l＝16.5×103×15/10＝24.8×103kN／m2
水平方向変形量 ｘ＝P/k＝518.5/24.8×103＝0.02m＝2cm
この場合の算定結果は、２次元弾性有限要素法による場合の算定結果に比べて若干の誤差があるものの、この種の解析においては両者の結果は実質的に同等であるといえるし、少なくとも評価結果に影響しない範囲内の誤差であるといえる。
【００２４】
（５）対策の決定
以上で算定された変形量が設計条件を満足すれば、改良率の仮決定（上記の場合はR＝50％）が妥当であったのでその改良率を最終決定として対策決定とする。変形量の算定結果が設計条件を満足しなければ改良率を変更して以上の手順を繰り返す。すなわち、変形量が過大であれば改良不足であるので改良率を大きくするように変更し、変形量が過小であれば改良過剰であるので改良率を小さくするように変更し、満足すべき結果が得られるまで以上の手順を繰り返せば良い。勿論、その際に必要であれば、すなわち改良率の修正のみでは条件を満足できない場合には、単位周期構造体のパターンや改良体の剛性も併せて見直せば良い。
【００２５】
以上で説明したように、本算定法によれば部分改良地盤の変形量の算定に際してまずその等価剛性を簡易チャートにより求めることにより、従来の３次元有限要素法により変形量を直接算定する場合に比べて全体の算定手順を格段に簡略化することができ、特に厳密なモデル化を行うための手間を大きく軽減することができる。したがって、従来のように多大の手間とコスト、時間を要することなく部分改良地盤全体の変形量算定を簡易に実施することができ、その結果、部分改良地盤による改良効果を事前に確認できてその信頼性を充分に高めることができ、特に上記のような護岸の液状化による側方流動防止対策として部分改良地盤を築造するに際して適用して好適である。
【００２６】
以上で本発明の基礎となった算定法について説明したが、次に本発明の部分改良地盤の変形量の評価法についての実施形態を図９〜図１３を参照して説明する。
本実施形態の評価法は、図９にフローチャートとして示した一連のステップ（ｉ）〜
（ix）により構成されるものであり、以下にその詳細を説明する。
【００２７】
（i）地震時に地盤に生じるせん断応力τdを算出する。せん断応力τdは地震応答解析から算定するか、簡易液状化算定法に用いられている（１）式を用いて求める。
【数６】
（１）式において、Ｍ：マグニチュード、αmax：地表面最大加速度、ｚ：深度、σν'：鉛直有効応力、σν：鉛直全応力、ｇ：重力加速度である。
【００２８】
（ii)上記の（１）式で求めたせん断応力τdと、図１１を用いて、改良前の地盤のせん断ひずみγiを算出する。図１１は（２）式をもとに作成したものである。
【数７】
（２式）において、ｈmaxは最大減衰定数で砂の場合は0.2〜0.25の値となる。初期せん断剛性Ｇsoや基準せん断ひずみγrfは砂の種類や拘束圧によって異なるが、典型的な例としてＧso＝50MPa、γrf＝0.1％、ｈmax＝0.2とした。
【００２９】
（iii）上述した算定法で説明した手法により、図３を用いて部分改良地盤の改良率、改良形状、改良パターン、および地盤と改良体のせん断剛性Ｇso、ＧIに応じた等価せん断剛性Ｇeqを算出する。
【００３０】
（iv）上記(i)で求めたせん断応力τdと、上記(iii)で仮算出した等価せん断剛性Ｇeqを用いて、（３）式により改良体間地盤に生じるせん断ひずみγiを求める。
【数８】
【００３１】
（v）上記(iv)で求めたせん断ひずみγiと図１０を用いて、せん断ひずみγiに応じた等価せん断剛性Ｇeqを再決定する。
【００３２】
（vi）再決定した等価せん断剛性Ｇeqと、（１）式で求めたτdを用いて、（３）式により再度せん断ひずみγiを計算し、以上の計算をせん断ひずみγiが一定の値に収束するまで行う。
【００３３】
（vii）決定したせん断ひずみγiを用いて、図１２により過剰間隙水圧比Δｕ/σν’を決定する。
図１２は（２）式と次の（４）〜（６）式を用いて求めたものである。
【数９】
ここで、Ｒlは繰返せん断応力比、Ｒ20は液状化強度で２０回で液状化に至るせん断応力比として規定している。Ｎlは液状化に至った繰返し回数である。ｋは実験定数で-0.25程度の値をとる。
図１３は液状化強度曲線の模式図を示したものである。
（４）式はＲl１回の繰返しで １／Ｎl だけ液状化に近づいたとみなせるので、２０回の繰返しせん断が生じた際の累積損傷度Ｒnは（５）式のように表せる。この状態で生じる過剰間隙水圧比は、De Alba の提案式である（６）式で表すことができる。αrfは実験定数であるが、緩い砂では0.7程度の値をとる。
【００３４】
(viii)液状化時に改良体に作用する外力を、たとえば上述した算定法に示した手法により算定する。
【００３５】
(ix)上記の外力が作用したときの改良体の変形量を、たとえば上述した算定法に示した手法（等価せん断棒や２次元平面ひずみ解析）を用いて算定する。
【００３６】
なお、上記のステップ(vii)において図１２を用いて過剰間隙水圧比Δｕ/σν’を求めることにより、その過剰間隙水圧比Δｕ/σν’が1.0未満であれば地盤は液状化に至らず過剰間隙水圧がある程度上昇するに留まることになり、そのことから部分改良地盤の液状化強度を評価することができる。そこで、仮に過剰間隙水圧比Δｕ/σν’が1.0以上であって地盤が液状化すると評価された場合には、改良率や改良パターンを再設定し、過剰間隙水圧比Δｕ/σν’が1.0未満になるまで以上の手順を繰り返せば良い。
したがって、以上で説明した本発明の部分改良地盤の変形量の評価法は、その過程で過剰間隙水圧比Δｕ/σν’を求めることから、部分改良地盤の液状化強度を評価するための評価法としても適用できるものである。
【図面の簡単な説明】
【００３７】
【図１】本発明の基礎となった部分改良地盤の変形量の簡易算定法の概要を示す説明図である。
【図２】同、簡易チャートの例（柱状の改良体の場合）を示す図である。
【図３】同、簡易チャートの例（格子状の改良体の場合）を示す図である。
【図４】同、護岸を対象とする側方流動防止対策としての部分改良地盤を築造する場合への適用例を示すもので、全体手順を示すフローチャートである。
【図５】同、護岸の構造図である。
【図６】同、簡易チャートにより等価剛性を求める手順を示す図である。
【図７】同、２次元弾性有限要素法による変形量の算定結果を示す図である。
【図８】同、等価せん断棒による変形量の算定手法を示す図である。
【図９】本発明の部分改良地盤の液状化強度の簡易評価法および変形量の簡易評価法の実施形態を示すもので、その全体手順を示すフローチャートである。
【図１０】同、せん断剛性低下率とせん断ひずみの関係を示す図である。
【図１１】同、せん断ひずみを求めるための図である。
【図１２】同、過剰間隙水圧比を求めるための図である。
【図１３】同、液状化強度曲線の模式図である。
【特許請求の範囲】
【請求項１】
バットレス状改良体や格子状改良体もしくは柱状改良体が造成された部分改良地盤の液状化強度を簡易に評価する方法であって、
改良前の地盤に地震時に生じるせん断応力τdを算出し、
前記せん断応力τdに基づいて改良前の地盤のせん断ひずみγiを算出し、
改良後の部分改良地盤の等価せん断剛性Ｇeqを算出し、
改良後の地盤における改良体間の未改良地盤に生じるせん断ひずみγiを求め、
該せん断ひずみγiに応じた等価せん断剛性Ｇeqを再決定し、
再決定した等価せん断剛性Ｇeqを用いてせん断ひずみγiが一定の値に収束するまで計算を行ってせん断ひずみγiを決定し、
決定したせん断ひずみγiに基づいて未改良地盤に生じる過剰間隙水圧比Δｕ/σν’を求め、
該過剰間隙水圧比Δｕ/σν’に基づいて部分改良地盤の液状化強度を評価する
ことを特徴とする部分改良地盤の液状化強度の簡易評価法。
【請求項２】
バットレス状改良体や格子状改良体もしくは柱状改良体が造成された部分改良地盤の変形量を簡易に評価する方法であって、
改良前の地盤に地震時に生じるせん断応力τdを算出し、
前記せん断応力τdに基づいて改良前の地盤のせん断ひずみγiを算出し、
改良後の部分改良地盤の等価せん断剛性Ｇeqを算出し、
改良後の地盤における改良体間の未改良地盤に生じるせん断ひずみγiを求め、
該せん断ひずみγiに応じた等価せん断剛性Ｇeqを再決定し、
再決定した等価せん断剛性Ｇeqを用いてせん断ひずみγiが一定の値に収束するまで計算を行ってせん断ひずみγiを決定し、
決定したせん断ひずみγiに基づいて未改良地盤に生じる過剰間隙水圧比Δｕ/σν’を求め、
改良体に作用する外力を算定し、
前記等価せん断剛性Ｇeqおよび前記外力に基づいて部分改良地盤の変形量を評価する
ことを特徴とする部分改良地盤の変形量の簡易評価法。
【請求項１】
バットレス状改良体や格子状改良体もしくは柱状改良体が造成された部分改良地盤の液状化強度を簡易に評価する方法であって、
改良前の地盤に地震時に生じるせん断応力τdを算出し、
前記せん断応力τdに基づいて改良前の地盤のせん断ひずみγiを算出し、
改良後の部分改良地盤の等価せん断剛性Ｇeqを算出し、
改良後の地盤における改良体間の未改良地盤に生じるせん断ひずみγiを求め、
該せん断ひずみγiに応じた等価せん断剛性Ｇeqを再決定し、
再決定した等価せん断剛性Ｇeqを用いてせん断ひずみγiが一定の値に収束するまで計算を行ってせん断ひずみγiを決定し、
決定したせん断ひずみγiに基づいて未改良地盤に生じる過剰間隙水圧比Δｕ/σν’を求め、
該過剰間隙水圧比Δｕ/σν’に基づいて部分改良地盤の液状化強度を評価する
ことを特徴とする部分改良地盤の液状化強度の簡易評価法。
【請求項２】
バットレス状改良体や格子状改良体もしくは柱状改良体が造成された部分改良地盤の変形量を簡易に評価する方法であって、
改良前の地盤に地震時に生じるせん断応力τdを算出し、
前記せん断応力τdに基づいて改良前の地盤のせん断ひずみγiを算出し、
改良後の部分改良地盤の等価せん断剛性Ｇeqを算出し、
改良後の地盤における改良体間の未改良地盤に生じるせん断ひずみγiを求め、
該せん断ひずみγiに応じた等価せん断剛性Ｇeqを再決定し、
再決定した等価せん断剛性Ｇeqを用いてせん断ひずみγiが一定の値に収束するまで計算を行ってせん断ひずみγiを決定し、
決定したせん断ひずみγiに基づいて未改良地盤に生じる過剰間隙水圧比Δｕ/σν’を求め、
改良体に作用する外力を算定し、
前記等価せん断剛性Ｇeqおよび前記外力に基づいて部分改良地盤の変形量を評価する
ことを特徴とする部分改良地盤の変形量の簡易評価法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【公開番号】特開２０１０−１３３２０４（Ｐ２０１０−１３３２０４Ａ）
【公開日】平成２２年６月１７日（２０１０．６．１７）
【国際特許分類】
【出願番号】特願２００８−３１２７１３（Ｐ２００８−３１２７１３）
【出願日】平成２０年１２月８日（２００８．１２．８）
【新規性喪失の例外の表示】特許法第３０条第１項適用申請有り 発行者名 社団法人地盤工学会 刊行物名 第４３回地盤工学研究発表会 平成２０年度発表講演集（２分冊の２） 発行日 平成２０年６月１２日
【出願人】（０００００２２９９）清水建設株式会社 (2,433)
【Ｆターム（参考）】
【公開日】平成２２年６月１７日（２０１０．６．１７）
【国際特許分類】
【出願日】平成２０年１２月８日（２００８．１２．８）
【新規性喪失の例外の表示】特許法第３０条第１項適用申請有り 発行者名 社団法人地盤工学会 刊行物名 第４３回地盤工学研究発表会 平成２０年度発表講演集（２分冊の２） 発行日 平成２０年６月１２日
【出願人】（０００００２２９９）清水建設株式会社 (2,433)
【Ｆターム（参考）】
[ Back to top ]