Q&A　（製品評価や導入の際に役立つQ&Aです） 　１．機能・概要 Ｑ１−１． トンネル掘削による周辺地盤の変形解析において「掘削解放率」の設定はできるか？また、併設トンネルの掘削で、先行トンネルの「掘削解放率」が100%未満の状態で後行トンネルの掘削を入力することはできるか？ Ａ１−１． 掘削解放率の設定は可能です。例えば、 Stage1で自重計算 Stage2で先行トンネルを掘削する。 この先行トンネルの解放率を40％、次のStage3で残りの60％に設定。 Stage3で後行トンネルを掘削する。 この後行トンネルの解放率を100％に設定。 というような段階施工検討ができます。 本プログラムでは、[ステージ設定]ウィンドウの[要素定義]タブの右側に並んでいるボタン群の一番下から2番目に[応力解放率の設定]ボタンがありますので、そちらより設定していただけます。 　 Ｑ１−２． 「弾塑性地盤解析(GeoFEAS)2D」の斜面の安定解析と「斜面の安定計算」の違いについて教えてほしい。 Ａ１−２． 「弾塑性地盤解析(GeoFEAS)2D」の斜面の安定解析は、せん断強度低減法のオプションを選択して弾塑性解析を行い、プログラム内の計算上でせん断強度を少しずつ低減していく過程で、斜面内に大きなせん断ひずみが発生するので、その面（画面上では線として見ます）をすべり面とし、その時に低減した割合を全体安全率とする解析です。したがって、有限要素解析のモデルですべり面と安全率を求めることができます。 一方、「斜面の安定計算」では、極限平衡理論をもとに力の釣り合いだけで、円弧により仮想的に設けたすべり面に対して安定計算をし、仮定したすべり円について安全率の計算を繰り返し、最終的に最小安全率となるすべり円を求める計算となります。 既往の研究では両者の結果がほぼ一致するとの報告があります。 　 Ｑ１−３． 粘性土で圧密沈下が生じ側方変位が予想される。橋台に対する影響を解析したいが、どのような土のモデル（構成則）を用いたらよいか。 Ａ１−３． 圧密計算を有限要素解析で理論的に行う場合にはCam Clayモデルや関口太田モデルが組み込まれたプログラムを必要としますが、本プログラムにはそうしたモデルはありません。 しかしながら、単純に弾性論上の範囲で周辺地盤の側方変形を考慮するならば線形弾性もしくはモールクーロンのような土のモデル（構成則）を用いることで解析することが可能です。 　 Ｑ１−４． 斜面の計算で、地下水位を考慮した場合としない場合の２通りを計算したいが、別々に実行したほうのがよいか、ステージを分けて実行した方がよいか？ Ａ１−４． 地下水位を考慮した場合としない場合について、一つのデータの中でステージを分けた場合、前のステージの応力等を引き継ぐので注意が必要 です。むしろ別々のデータを作成したほうが、目的に沿った解析ができます。 　 Ｑ１−５． 堤防の地震時を解析したい場合に、要素モデルに「液状化モデルや低減モデル」を設定して地震による水平慣性力をかけて解析するようなことは可能か？ Ａ１−５． プログラム上のメニューとしては設定が可能です。しかしながら、液状化時の計算は「河川構造物の耐震性能照査規準」に則り特殊な取り扱いをしております。同規準では水平慣性力としてかけることはないので、液状化モデルでは避けたほうがよろしいでしょう。 本来、液状化現象は地震によって間隙水圧が発生する事で、多少のタイムラグがあって地盤の剛性が小さくなり自重により沈下します。実際の現象から見ても、地震がかかっているときに瞬時に（同時に）地盤の剛性が低減するわけではないので、液状化モデルによる剛性の低減と水平慣性力とは同時にかけないのがよろしいでしょう。 　 Ｑ１−６． メッシュをはきだすことは可能か？ Ａ１−６． メッシュをデータとして見ることは可能です。拡張子がMhsのファイルをご覧ください。 　 Ｑ１−７． 要素シミュレーションは可能か？（土質試験にあるような応力ひずみ曲線の結果をみてパラメーターを検討できるか） Ａ１−７． 製品の中に、要素試験シミュレーションプログラムが添付されておりますのでご覧ください。 　 Ｑ１−８． ALID相当の液状化時の変形解析および体積圧縮による沈下量の算出方法について教えてほしい。 Ａ１−８． たとえば、 　　Stage1　：　液状化前 　　Stage2　：　液状化時 　　Satge3　：　体積圧縮時 とします。 変形による沈下量は、Stage2-Stage1で、すなわち、Stage2と1の相対変位で得られます。 同様に、体積圧縮による沈下量は、Stage3-Stage2で、すなわち、Stage3と2の相対変位で得られます。 また、合計値は、Stage3-Stage1で、すなわち、Stage3と1の相対変位で得られます。 Page Top 　２．使い方 Ｑ２−１． ダイレタンシー角の入力方法について、教えてほしい。 Ａ２−１． ダイレタンシー角の入力方法としては、変形解析のときは内部摩擦角と同じとし、安定解析のときはゼロとする方法が推奨されています。 Page Top 　３．要素・構成則 Ｑ３−１． 「地盤と構造物の接触面の指定（ジョイント要素）」は可能か？ Ａ３−１． 可能です。弾性ジョイント要素、MCジョイント要素を用意しています。ただし、いわゆる剥離現象を再現できるものではありません。 　 Ｑ３−２． ジョイント要素を複数設定したいが、可能か？ Ａ３−２． ジョイント要素は地盤と地盤（ソリッド要素）との間に設けます。 設定する場所と隣合うソリッド要素を規定するため一箇所ごと（一直線上の線分の固まり）に設定します。 Page Top 　４．解析理論 Ｑ４−１． 「斜面安定解析」とはどのような解析か？通常、FEM計算で斜面安定解析を行わないと思うが、どのような使い方をするか教えて欲しい。 Ａ４−１． せん断強度低減法の機能を用いると、各要素の最大ひずみ増分が求まります。これをコンタ図で描くと、すべり面の発生状況をみることができます。 同時に全体安全率も求めることができます。つまり、最小安全率のすべり面を自動検索してくれるというイメージです。 　 Ｑ４−２． 「水圧の変動が地盤に及ぼす影響の検討」とは、実際に、可能なのか。どのような使い方をするか教えてほしい。 Ａ４−２． GeoFEAS2Dは、全応力解析（地盤の透水現象を考慮しない解析）を行うプログラムですが、水位線を与えてプログラム内部で水圧を節点荷重として設定することにより、水圧の上昇による地盤の強度低下を考慮した解析（全体安全率など計算）することができます。 水圧が上昇するとなぜ地盤強度が低下するかについては、水平の地下水面を持つ地盤を考えてください。地下水面が上昇すれば、浮力が大きくなり、（有効）応力が小さくなります。従って、地盤の強度が低下します。 　 Ｑ４−３． 「応答震度法」とは何をイメージしているのか教えてほしい。 Ａ４−３． 静的解析手法で、動的解析に近似の結果を得ようとする解析法と考えられます。 具体的には、 水平方向加速度（鉛直方向に一様でない水平加速度も対応できる）を地盤にかけて、地盤には水平方向の慣性力を作用します。そして、地中構造物などに地震の影響を検討することができます。 　 Ｑ４−４． 軸対称解析を行う場合の適用モデル並びにその有効性、優位性についてほしい。 Ａ４−４． LNG地下タンクのために、円形掘削による地盤の変形解析 タンク基礎の沈下 三軸試験結果のシミュレーション 等。 つまり、地盤や構造物は軸対称であり、荷重も軸対称である場合は、軸対称解析を行うことが必要です。 軸対称解析を行うと、入力は2次元断面ですが、くるりと一周した積分解析となりますので、3次元的な影響を考慮できることになります。 同じ2次元断面を入力した場合、平面ひずみ解析は、奥行き方向に半無限に同じ断面があることになるので、奥行きの影響は考慮されませんが、軸対称は全周にくるりと同じ断面があることになるので、奥行きの影響を考慮したことになり、結果的に、3次元的な解析になるとご理解下さい。 Page Top 　５．モデル化 Ｑ５−１． メッシュ分割の際、メッシュの大きさをどのくらいに設定したらよいか。 Ａ５−１． 何メートルメッシュでなければならないというような規定はありません。そのため、これぐらいの大きさでメッシュ分割しなさいという助言も致しかねる状況です。 メッシュ分割を非常に細かくすれば、FEMの計算結果は厳密解に近づくと言われていますが、その確認は簡単ではないということのようです(はじめて学ぶ有限要素法：地盤工学会より)。 群馬大学の蔡先生のセミナー資料によると、要素数を３倍にしても、安全率のオーダーとしては、小数第三位に違いがあるだけのようです。よって、あまり、細かく分割する必要はないのではないかと思われます。 メッシュ分割において、中間節点の指定を有(2次要素)にすることも計算精度を向上させる手法だと思います。メッシュ分割の定義タブの右側最下段のボタンで設定してみてください。 　 Ｑ５−２． 弾完全塑性の物性値のサンプルはないか？ Ａ５−２． 弾完全塑性の物性値としては、強度特性とダイレタンシー角が主なパラメータとなります。 強度特性については、「山岳トンネルにおける 模型実験と数値解析の実務」トンネルライブラリー第１６号 土木学会などを参照すると、さまざまな機関における一般値が掲載されています。三軸試験など室内試験結果があれば、そのＣ，φを用い て結構です。 ダイレタンシー角についてはゼロとする場合や内部摩擦角と同じにする場合もあります。三軸試験結果を参考にすることも考えられます。　「弾塑性有限要素法がわかる」、地盤工学会および「地盤の変形解析: 基礎理論から応用まで （39）」、地盤 工学会， 2002.7(地盤工学・基礎理論シリーズ）などが参考となります。 サンプルデータとしては、GeoFEAS2Dに添付されているデータ（slope， embankmentなど）をご参考にしてください。 　 Ｑ５−３． ソリッドのプロパティの設定での構成則や方式の使い分けを教えてほしい。 Ａ５−３． 概ね、以下のようにお考え下さい。 積層弾性モデル 直交異方性を表すために用います。積層ゴムのように縦方向は剛で水平方向には柔な構造を対象としています。地盤においても堆積過程や地中応力状態において異方性が大なり小なりがあると言われています。 非線形弾性モデル Ｄｕｎｃａｎ方式：応力ひずみ関係を双曲線で表します。ひずみのレベルや拘束圧（対象とする地層の深さ）によって変形係数が変わることをモデル化します。Ｄｕｎｃａｎ方式１と２の違いはポアソン比を定数（方式１）とするか、変化する（方式２）かの違いにあります。 非線形弾性モデル 破壊接近度法：主にトンネルの掘削にともなう地盤の緩みを解析するときに用いられます。応力状態が破壊基準に近づくに従って変形係数が減少（地盤が緩んで変形しやすくなる状態）を表現します。 非線形モデルＲＯモデル せん断弾性係数とせん断ひずみを双曲線でモデル化して地盤のせん断変形（地震時の変形など）に対してひずみ依存性を表現するときに用います。 非線形モデル　ＨＤモデル ＲＯモデルと同様にせん断変形に対してひずみ依存性を考慮するときに用います。ＲＯモデルとは関数の形が異なります。 非線形モデル　ＵＷclayモデル せん断応力とせん断ひずみの骨格曲線を双曲線で表しますが、強度をパラメータと取り入れていることで変形と強度の関係を簡易に表すことができます。斜面の残留変形解析に用いた事例があります。 弾完全塑性モデル　Mohr-Coulomb方式（ＭＣ） 破壊基準（Mohr-Coulomb)を規定し、応力状態を表す応力円（Ｍｏｈｒ円）との関係を判別して破壊基準に達していなければ弾性領域、破壊基準に達すれば塑性領域に分けて計算します。塑性領域に達した場合は大きな変形が生じます。 弾完全塑性モデル　Drucker Prager方式（ＤＰ） 破壊基準の規定のしかたがMＣと異なりますが、２次元（平面ひずみ）の場合は違いはありません。 弾完全塑性モデル　MC-DP 破壊基準としてＭＣ，塑性ポテンシャルにはＤＰにより規定します。破壊に達した状態で、せん断にともなう体積変化であるダイレタンシーが収縮するのか（負のダイレタンシー）膨張する（正のダイレタンシー）のか地盤によって規定したいときに用います。 弾塑性モデル　粘土のPastor-Zienkiewiczモデル 弾塑性モデルの一般化したモデルで、パラメータによってひずみ軟化やひずみ硬化、または硬化係数（破壊に達したあとの変形係数）などを細かく規定できるモデルです。多くのパラメータを規定する必要があります。 No-tension材料 線形モデルあるいは積層線形モデルと基本は同じですが、引っ張り応力に対して抵抗しない材料を表現するモデルです。 　 Ｑ５−４． 覆工と吹付けの間にジョイント要素を設定しても問題ないか？ Ａ５−４． 実現象を再現しようとして細かく設定し過ぎると意図しない挙動を示すモデル化になる場合もありますので、ジョイント要素を用いる場合は注意が必要です。 　 Ｑ５−５． 薬液注入工がある場合、変形係数を原地盤より大きくすることは可能か？ Ａ５−５． 薬液注入する場合は、最初のステージ１では土、次のステージ２で薬液注入後の剛性を入れ直して、プロパテイ設定をすることができます。薬液注入後の剛性で、その後の変形解析をします。ただし、ステージ２でステージ１と比べ荷重の変化がなければ、変形等は発生しません。 　 Ｑ５−６． トンネルの解析後、沈下防止のため地盤改良を実施したい。どの時点でのステージを用いたらよいか。 Ａ５−６． 沈下防止のために地盤改良ということであれば、沈下の原因（解放応力もしくは荷重がかかる）のあるステージと同じステージでプロパティ設定を変更します。すなわち、トンネル掘削に伴い沈下が生じる場合は、同じステージにおいてプロパティ設定を変更すれば、そのステージのプロパテイを基に掘削時の変形を解析します。 　 Ｑ５−７． 風化に伴う物性値の変更は可能か？ Ａ５−７． 別のデータをプレ部から作成して別の物性を与えた場合に、材料の変化（風化）を考慮して比較できます。 　 Ｑ５−８． 推進工法で破壊接近度法を用いて検討したい。緩み係数ｋの妥当値について教えてほしい。 Ａ５−８． 破壊接近度は本来、岩石を対象として提唱された理論であり、パラメータについては岩石試験などの試験結果から決めます。お問い合わせ内容では、土かぶり厚が４ｍ程度で砂質土とのことなので、一般に提唱されている最も低い（非線形の強い材料）をもとにパラメータを決めればよいでしょう。 　 Ｑ５−９． 地中部に存在する埋設配管上を重機(クレーン、トラック等)が走行する場合に埋設配管に作用する応力(地中応力)の算定は可能か？ Ａ５−９． 解析可能と思われます。本製品は、全応力静的FEM解析ソフトになります。 埋設配管に作用する荷重(地表面上の重機荷重)は、たとえば、分布荷重として与条件として、ご自身でご用意して頂く必要がございます。地盤をソリッド要素、埋設管をはり要素でモデル化して頂ければ、埋設配管の断面力、変位などを直接得ることができます。 　 Ｑ５−１０． 橋台と地盤との間にジョイント要素を用いたい。使い方を教えて欲しい。 Ａ５−１０． ジョイント要素のバネは周辺地盤よりも大きな剛性を設定して、荷重レベルが小さい段階では２重節点を維持し、ある設定強度（モールクーロン方式ではＣ，φ）に達すると不連続を再現します。本来、岩盤の亀裂などを再現するために作られたので、通常の土を対象にする場合は不用意に用いないほうが安定した解が求まります。むしろ橋台と地盤との間にダミーの境界面を意味する薄い層を設けたほうがよい場合もございます。 　 Ｑ５−１１． 断層はどのようにモデル化するのか。 Ａ５−１１． ジョイント要素になります。 Page Top 　６．連動 Ｑ６−１． 土留との連動内容について、教えてほしい。 Ａ６−１． 「土留め工の設計」にて「有限要素法(FEM)」 による地盤のみモデル化し、別途弾塑性法により計算した、あるいは掘削時に計測した壁体変位を入力し、地盤変形を計算する方法について「弾塑性地盤解析(GeoFEAS2D)」用の入力データの作成に対応しております。  >> サポートページ　弾塑性地盤解析（GeoFEAS） 2D　Ｑ＆Ａ集