以下のように分割処理をしています。

1. 地表面形状に沿って屈折点毎に堤体を縦割りにします。
2. 各縦割区間の水平距離をΔｘとし、帯片基準分割幅Δｂで除した値の切り上げ数が、その区間の分割数ｎとなります。
   （ｎ＝Δｘ／Δｂ）　
3. 計算上の分割幅は、「Δｘ／ｎ」となります。
4. ブロック番号のふり方はすべり形状、すべり方向（時計回り、反時計回り）に関わらず、画面向かって左から「123n」としています。

• 滑動モーメント計算用
  水位より上：常にγｔ（不飽和土の単位体積重量）を用います。
  水位より下：常にγsat（飽和土の単位体積重量）を用います。
• 抵抗モーメント計算用
  水位状態に関わらず、以下の５パターンの組み合わせを任意に選択できます。
  

  	単位体積重量γ	間隙水圧ｕ
  １．	γｔ（入力値）	０．０とする。
  ２．	γsat（入力値）	０．０とする。
  ３．	γ'（入力値）	０．０とする。
  ４．	γsat−γw	０．０とする。
  ５．	γsat（入力値）	計　算

  （設定例は「γの取り扱い」参照）

何れか滑り長の長い方の山側（Ａ）のみを取得して計算し、他方は無視します。
従って滑り長の短い方の山側（Ｂ）に「マストカットライン」や「マストカットポイント」を指定していても無視されます。
逆に上図の場合で（Ｂ）のみを切りたい場合は（Ａ）を「ネバーカットライン」を定義して下さい。

下図のように貯水部があるモデルでは側水圧（静水圧）による起動（滑動）モーメントを考慮します。

その他に水位を入力した場合、常に滑り面のスライス毎に下図のように両側から側水圧が作用します。
※「盛土直後」の場合は水位の入力は不要ですが、入力すれば側水圧は考慮されます。

1. 側水圧の向き（符号の取り方）
   • 側水圧の差ΔＨ（＝H2−H1）が正の場合
     →右から左に向かって作用。
   • 側水圧の差ΔＨ（＝H2−H1）が負の場合
     →左から右に向かって作用。
     
2. 滑動側計算時
   水位を入力すると、常に側水圧による滑り力を考慮します。
   • 滑り力（滑り面に作用）
     Ｓ＝（H2−H1）COSα
   • 滑りモーメント
     Ｍ＝H2h2−H1h1
   上記を各スライスごとに算出し、全スライス分を合計します。
   滑り面の両端部で水圧が相殺されていない場合、側水圧による水平力が発生します。
   
   水位線が水平のとき、滑り円の終端が水位線まであると考えるならば�狽r＝０となります。
   ただし、計算は太線のすべり面で行い、貯水部は三角形分布の側水圧（静水圧）として考慮するため、モーメントは釣り合い（�狽l＝０）ますが、すべり力は一般に�狽r≠０となります。
   したがって、側水圧によるすべり力の合計値が０となるのは、すべり面内の水位が円中心のＹ軸に関して左右対称の場合のみです。
3. 抵抗側計算時
   3-1)間隙水圧考慮の場合
   • 抵抗力（滑り面に作用）
     
     
   3-2)浮力（有効重量）考慮
   浮力を考慮した土の単位体積重量γ'を用いて土塊重量を算出しているため、側水圧による抵抗力は算出していません。
   ※土塊重量Ｗ’を求める際、「γ'＝γsat−γw」により、 「−（ΔＨsinα＋U）」分の水圧は、既に加味されている。

※水位がレベルであれば（貯水部の有無に関わらず）、間隙水圧で考慮した場合と浮力（有効重量）で考慮した場合とでは安全率は同じになります。（次にその検証例を示します。）

�@水位が水平（土中）の場合

１）抵抗力について

２）すべり力について

�A水没斜面の場合

１）抵抗力について

２）すべり力について

したがって浮力（有効重量）で考えた場合と等しくなる。
※厳密にはプログラムの内部処理として浮力で考えた場合と間隙水圧で考えた場合とでは滑り角の取り方が異なりますが（下記）、計算結果に影響が出るほどの違いはありません。

1. 浮力の場合
   
2. 間隙水圧の場合
   各スライス重心ｘ座標の滑り面での法線方向とする。

初めに本ソフトでは「水平慣性力の方向は、すべり方向に応じてプログラム内部で設定します。」という扱いについて説明します。
土質ブロックデータの中に、地震時設計震度の補正係数βの入力があり、入力した設計震度（Ｋｈ、Ｋｖ）に対して、（β・Ｋｈ、β・Ｋｖ）として各地層毎に設計震度を変化させることを実現しています。
「すべりの方向に応じて」とは、
　　　　　　　　�@β＞０の時には、すべりの方向と同じ向きに考慮する
　　　　　　　　�Aβ＜０の時には、すべりの方向と反対向きに考慮する
という扱いを意味しています。
それでは「すべりの方向」とは何を基準に決定しているのかといいますと、『すべり円の開始位置と終了位置のＹ座標の標高をみて、高い方から低い方をすべりの方向とする』という判断をしています。
このような処理の中で、扱いがはっきりしない場合が、すべり円の開始位置と終了位置のＹ座標が等しい場合、すなわち、フラットな地表面上のすべり円の処理が挙げられます。（下図）

つまり、荷重により反時計廻りにすべりが発生しているにも拘わらず、慣性力にはこれに反発するように、時計廻りに作用しており、すべりの方向に応じていない結果になっているわけです。フラットな地表面ですべり円が発生した場合、すなわち『すべり円の開始と終了の標高差が０の場合は、プログラム側で時計廻りにすべり円が発生したものとして一方的に処理をする』というのが仕様になっています。すべりの方向の定義づけに問題があるのではないかという異論もあろうかと思いますが、予めフラットな地形上ですべりを検討される場合には、　
　　　　　　�@慣性力を時計廻り（右→左）に考慮する場合にはβ＞０とする。　
　　　　　　�A慣性力を反時計廻り（左→右）に考慮する場合にはβ＜０とする。
という入力で整合を図って下さい。又、前述のすべり方向の定義から地表面に若干の勾配を持たすことでも回避は可能です。

安全率を力の比で求めた場合とモーメントの比で求めた場合の違いは以下のようになります。

1. 力の比
   • Ｆｓ＝ＲＤ／ＳＤ
     ＲＤ：滑り面上の抵抗力の総和
     ＳＤ：滑り面上の滑動力の総和
     　ＲＤ＝�煤oＣiＬi＋（Ｗi'＋ｋｖＷi）cosαi−ｋｈＷisinαi−Ｕi−ΔＵi｝tanφi
     　ＳＤ＝�煤o（１＋ｋｖ）Ｗisinαi＋ｋｈＷicosαi｝＋ＳP
2. モーメントの比
   • Ｆｓ＝ＭＲ／ＭＤ
     ＭＲ＝ＲＤ×ｒ
     ＲＤ：滑り面上の抵抗力の総和
     　　ｒ：滑り円の半径
     ＭＤ＝�煤o（１＋ｋｖ）ＷiΔｘi＋ｋｈＷiΔｙi｝＋ＭP＋Ｍw
     ここに
     Ｗi：スライスの全重量（γt、γsatを用いる）(tf/m)
     αi：スライスの滑り角（度）
     Δｘi：スライス重心と滑り円中心の水平距離(m)
     Δｙi：スライス重心と滑り円中心の鉛直距離(m)
     ｋｖ、ｋｈ：設計震度
     ＳP：荷重よる起動力(tf）
     ＭP：荷重による起動モーメント(tfm)
     Ｍw：側水圧（静水圧）による起動モーメント(tfm)
     
     したがって「滑り力」は滑り面上での評価であるのに対し、「滑動モーメント」はスライスの重心位置と滑り円中心からの距離での評価であるため「滑動モーメント／半径＝滑り力」とはなりません。
     もし常時（地震時慣性力なし）で水位がレベルの場合（ただし水位とすべり面の始終点の高さが同じ）、側水圧が相殺されるので結果的にほぼ
     「滑動モーメント／半径＝滑り力」
     が成立します。
   
   

「−」：安全率が負の場合
（通常、考えられませんでデータをお送りください。）

「0.00」
　　�@抵抗モーメントが０．０となる場合（全層のｃ、φ＝０．０）
　　�A計算可能なすべり円が存在しない場合

（「堤体」内の「格子」についても計算の対象外となっています。）

「999.90」：安全率が限りなく大きい場合。もしくは簡易ビショップ法で収束できなかった場合をそれぞれ意味します。

下記２項目を共に「する」もしくは片方のみ「する」とした場合、必ずしも限界面計算結果のＦｓが格子点上にあるとは限りません。

• 等安全率の計算「する」
  各格子点及び各格子中心点が対象。
• 限界円の自動追跡「する」
  自動追跡の考え方に従い（「限界面の自動追跡計算」参照）各格子点及びそれ以外も対象。

自動追跡結果のＦs の方が格子点Ｆs より大きくなる。
自動追跡範囲と格子範囲を同一として最小安全率を求めると以下の結果となった。

• 自動追跡結果　　Ｆs ＝1.83
• 格子点安全率　　Ｆs ＝1.75
• 等安全率の結果　Ｆs ＝1.72

格子点安全率は各格子点上の安全率の中から最小安全率を抽出するだけですが、自動追跡や等安全率の計算はそれぞれ次のように行い、最小安全率を検索しています。

• 限界面の自動追跡計算
  追跡範囲内において特定のライン上の安全率最小点を見いだし、次にそのライン上の安全率最小点と直交するライン上での安全率最小点を求めるといった処理を繰り返すことにより全体での安全率最小点を追跡的に求めます。
  （「限界面の自動追跡計算」参照。）
• 等安全率線（コンタライン）の算出
  1. 格子の交点上の安全率を計算します。（ａ〜ｄ）
  2. 格子の中心の安全率を計算します。（ｅ）
  3. １〜８の線分に対してコンタラインの対象となる安全率※ が線分両端の安全率の間であれば、その位置を比例で求めます。（図中○）
  4. １つの格子中にできた４つの三角形を順に見ていき、○を線でつないで線分のデータを作成します。
  5. 全ての格子に対して同様の処理を行います。
     ※「二次処理」→「等安全率」で表示されている安全率（最大５コ）
     
     結果に表示される最小安全率Ｆs は上図の例で云えばａ、ｂ、ｃ、ｄ、ｅの５箇所（４隅と中心）の中から抽出された値であり、○印のＦs は抽出対象となっていません。（５箇所以外は計算していません）
  
  以上のように、自動追跡計算も等安全率の計算も検索方法が必ずしも万能なわけではなく、どうしても検索もれが生じ、本件のケースのような結果になる場合があります。
  したがって比較的望ましい計算方法としては、
  1. 格子範囲を広い範囲で大きいきざみ幅から狭い範囲で小さいきざみ幅へと徐々にデータを変えていきます。又、必要に応じて等安全率の計算を行い安全率の分布状況を視覚的におさえておきます。
  2. １の作業で安全率が厳しくなるところをある程度絞り込んだら、そこを追跡範囲とし自動追跡の計算を行って最小安全率を更に絞り込みます。

本プログラムではマストカットラインやマストカットポイントを指定した場合は、少なくとも何れか一方を通るすべり面のみを安全率計算の検索対象にします。したがって一般に限界面の計算結果にはそれらを通るすべり面で結果を表示します。
本件の場合例えば下図のような場合が考えられます。

本プログラムではすべり長の大きい方を優先し（Ａ＞Ｂ）Ｂ部のすべり面は計算の対象外となります。そのため仮にＢ部を「マストカットライン」で定義していても、結果的に「マストカットライン」を無視したものと見えます。

対策１）可能であれば地表面の凸凹をなくし直線にモデル化する。
対策２）Ｂ部を切りたい場合、Ａ部を「」で定義する。
　　→Ｂ部のみを切ります。

スライス間力（鉛直力Ｖ、水平力Ｈ）を考慮した安定解析方法「Ｂｉｓｈｏｐ」のうち「Ｖ＝０．０」とするものを、「簡易Ｂｉｓｈｏｐ」と呼びます。
本プログラムでは「簡易Ｂｉｓｈｏｐ」をサポートしています。
概要は以下の通りです。

1. 内部で設定されている安全率ＦＳ０の初期値（=1.0）から抵抗力を求める。
2. 抵抗力に半径を乗じて抵抗モーメントを求める。
3. 別途求めた滑動モーメントと上記の抵抗モーメントから安全率ＦＳ１を求める。
4. ＦＳ１−ＦＳ０＞X.XXX
   X.XXX：「計算内容の設定」における「最小安全率の精度」
5. ＦＳ１を新たな初期値として再度、抵抗力を求める。
   以下、下式が成立するまで2)以降を繰り返す。
   「ＦＳｎ＋１−ＦＳｎ≦X.XXX」
   （詳しい算式等は「Ｂｉｓｈｏｐの簡易法」を参照して下さい。）
   注）
   上記計算において結果的に収束できなかった場合は「格子点安全率」の表に「999.90」と表示されます。
   その時の対策としては「最小安全率の精度」を大きくして再計算してみる他ありません。