増分解析（プッシュオーバー解析）
保有水平耐力を求めるため、漸増静的（プッシュオーバー）解析を行います。画面上の表記は「増分解析」です。 外力分布は Ai 分布（荷重）とし、荷重制御から変位制御へと段階的に制御を切替えます。
終了目標
変位増分をどこまで押し込むかは、次の 2 つの終了目標で指定します。それぞれ独立に有効・無効を切替えられ、有効な目標のいずれかに最初に達した時点で解析を打ち切ります。
| 目標 | 内容 | 既定 |
|---|---|---|
| 目標層間変形角 | 全層の最大層間変形角（\( \max_i |\delta_i| / h_i \)、\( h_i \) は階高＝階標高の隣接差分）が指定値（1/n）に達したら終了 | 有効（1/150） |
| 目標変位 | 頂部（最上階の代表節点）の加力方向変位が指定値 [mm] に達したら終了 | 無効 |
両方を無効にした場合は荷重制御（λ=1、設計地震力レベル）までで終了します。 変位制御フェーズの押込み上限は、目標変位はその値を、目標層間変形角は「全層が一様に目標角へ達した場合の頂部変位」＝目標角 × Σ階高を用います（最大層間変形角は平均層間変形角＝頂部変位/Σ階高 以上であるため、上限到達までに必ず判定が成立します）。
終了目標の判定は、水平力を載荷するフェーズ（荷重制御・変位制御）すべてで行います。 あるフェーズで目標に達した場合、以降のフェーズへは移行しません。 長期荷重の初期載荷（λ=0）は水平力に対する応答ではないため、判定の対象外です。
増分の刻み（均等変位刻み制御）
段階制御では、性能曲線の点間隔（頂部変位軸）が全域で概ね均等になるよう増分を刻みます。 目標刻みを du =（押込み上限）/（ステップ数設定）とし、次の規則で増分を決めます。
- 荷重制御: λ 増分を「頂部変位が du 進む見込みの量」に適応させます。増分は直近の確定ステップ間の荷重−変位勾配（初回は初期接線剛性による弾性勾配）から求めます。剛性が変化しない弾性域は粗い λ 刻みで足り、降伏が進み勾配が増すほど刻みは自動的に細かくなります。
- 変位制御: 同じ du を頂部変位の押込み刻みに用います。
弾性勾配を推定できない場合は、荷重制御を等分割 λ 刻み（1/ステップ数）、変位制御を 10 分割とする固定刻みで解析します。 増分方式「荷重増分のみ」は常に等分割の λ 刻みです。
算定式
外力は標準地震荷重生成の節点荷重を自由度へ載荷します（地上部 Ai 分布と地下部の標準式。載荷先は地震力の「Pi の載荷先」に同じ）とし、周期は略算式によります（\( C_0 = 0.2 \), \( Z = 1.0 \), 第2種地盤）:
この増分載荷経路の Z・地盤種別・C0 は上記の固定値で、画面の地域・地盤設定を伝達しません。対象建物への適用はこの制限を確認してください。地下の直上せん断力加算、指定 Ci の副剛床の別荷重、入力診断は標準地震荷重入口と共通です。PH 指定震度は未検証として拒否します。
\[ T = h(0.02 + 0.01\alpha) \]
Newton 反復は次とし（幾何剛性オプション use_kg、1 ステップあたり最大 50 回）、収束を \( \|r\| < 10^{-6} \cdot \max(\|f_{\text{ext}}\|, 1) \) で判定します:
\[ K_t \cdot \delta u = f_{\text{ext}} - f_{\text{int}} \]
解析経過の最大ベースシアは次とします:
\[ V_{base,max} = \text{性能曲線上の最大ベースシア} \]
記号・単位
| 記号 | 意味 | 単位 |
|---|---|---|
| \( \lambda \) | 荷重係数（設計地震力レベルで 1） | - |
| \( \boldsymbol{q} \) | Ai 分布の荷重パターン | N |
| \( K_t \) | 接線剛性行列 | N/mm |
| \( \boldsymbol{f}_{\text{ext}},\ \boldsymbol{f}_{\text{int}} \) | 外力・要素内力ベクトル | N |
| \( \delta u \) | 変位増分ベクトル | mm |
| \( T \) | 一次固有周期（略算） | s |
| \( h \) | 建物高さ | m |
| \( \alpha \) | 鉄骨造部分の高さ比 | - |
| \( V_{base,max} \) | 解析経過の最大ベースシア | N |
| \( Q_u \) | 保有耐力比較用の採用点における層耐力 | N |
| \( Q_{un} \) | 必要保有水平耐力 | N |
| \( D_s \) | 構造特性係数 | - |
| \( F_{es} \) | 形状係数 | - |
| \( Q_{ud} \) | 地震時の設計用層せん断力 | N |
| \( C_0 \) | 標準せん断力係数 | - |
接線剛性が特異・非正定値になったときの扱い
弾塑性の接線剛性 \( K_t \) は、降伏の進行（崩壊機構の形成・耐力劣化）によって正定値性を失うことがあります。増分解析はこれを増分が大きすぎるときの非収束と同じ扱いとし、水平力の荷重制御・変位制御では次の順で対処します。長期荷重の初期載荷は再試行せず、失敗理由を表示して解析を停止します。
- 確定済みの状態から増分を半減して解き直す（最大 5 回。確定 λ・確定頂部変位を基点とした増分のみを縮めるため、単調載荷が保たれます）。
- それでも進めない場合はそのフェーズを打ち切り、次のフェーズ（荷重制御 → 変位制御）へ引き継ぐ。極限点近傍は変位制御の方が安定に追えるためです。
- 1 ステップも確定できなかった場合のみ、原因を切り分けたメッセージでエラー停止する（性能曲線が空の結果を「解析できた」として返すと、保有水平耐力を 0 と誤認させるためです）。
終了理由の記録
解析がどのように終了したか（終了目標へ到達／荷重係数上限まで載荷／押込みスケジュール完了／非収束による打ち切り）を結果へ記録します。現在の解析では、接線剛性の求解失敗で増分半減を尽くした場合も、非収束による打ち切りとして記録します。目標へ到達しないまま非収束・数値的な特異性などで打ち切られた場合、性能曲線は途中で途切れており \( V_{base,max} \)（性能曲線上の最大ベースシア）はその時点までの最大値です。結果画面はこの場合に警告を表示します（過小評価の可能性があるため、増分ステップ数の増加やモデルの見直しを検討してください）。
接線剛性の数値的な特異性・精度限界は、有限精度で信頼できる求解を確認できない状態を表し、物理的には安定した架構でも極端な剛性比などで生じ得ます。この終了理由だけで物理的な機構形成・耐力喪失を断定せず、降伏ヒンジの分布に基づく崩壊機構の分類とは区別します。途中打ち切りの場合は、解析通知にも終了理由と最大ベースシアがその時点までの最大値であることを表示します。
解析開始前の初期接線剛性チェック：初期（未変形・未降伏）の接線剛性を一度分解し、失敗した場合は増分を刻む前にエラー停止します。初めから特異なモデルは増分をいくら刻んでも解けないため、剛性がゼロの自由度を節点 ID と自由度成分で名指しして停止します（例: 梁だけが一直線に並ぶ節点のねじれ回転）。分解できるが対角が負の自由度がある場合は耐力劣化・機構形成として、いずれも該当しない場合は拘束不足・機構の可能性として案内します。
自由度を名指しする診断は、拘束（剛床・剛リンク・MPC）を縮約した後の接線剛性 \( K_r = T^\top K T \)（＝実際に分解する行列）の対角で判定します。名指しの対象は縮約後に独立自由度として残るものだけなので、剛床・剛リンク・MPC のスレーブ自由度は挙がりません。
長期荷重の初期載荷
水平力の増分に先立ち、保存した架構用長期組合せの terms・係数による外力を載荷し、その応力状態（柱軸力・梁端モーメント）を初期条件とします。標準は固定荷重と架構用積載の係数 1 の組合せであり、地震重量用積載や全ての長期系ケースの総和ではありません。組合せの明示選択と完全性の条件は終局検定の重力参照と共有します。 長期荷重の初期載荷は既定で有効です。節点荷重・部材荷重の全量を、荷重係数 1.0 の固定外力に対する1回の静的 Newton 求解で釣り合わせます。収束した場合だけ変位・応力・材料履歴を確定します。収束しない場合は載荷前の状態へ復元して理由付きで停止し、分割載荷や長期荷重なしの解析へ切り替えません。 この設定を無効にした結果は、一般建物の終局合否には使用できません。 長期応力を初期状態とすることで、軸力に依存する部材耐力（柱の曲げ降伏 \( M_y \)・せん断降伏 \( Q_y \) の軸力項、ファイバー断面の N–M 相関）が長期軸力を含んだ状態で評価されます。 長期荷重は全フェーズを通じて外力 \( f_0 + \lambda \boldsymbol{q} \) の \( f_0 \) として保持され、載荷完了状態は荷重係数 λ=0 の 1 ステップとして記録されます（ヒンジ詳細図の N-M 応答経路の始点）。 初期載荷は右パネル「増分解析」で無効化できます。明示的に重力を載荷しない解析は参考の水平載荷実験として扱い、一般建物の終局合否には自動採用しません。初期重力を載荷する場合は、組合せの未指定や必要ケースの不足を重力なしへ読み替えません。
解析には初期重力の terms・係数・入力識別・載荷有無を記録します。確定内力は初期重力を含む総応答です。増分・確定総応答による終局検定では、同じ解析の重力条件と確定評価ステップを採用し、QL を再加算しません。表示ステップの切替だけでは需要採用ステップを変更しません。
材端軸力の符号
ヒンジ詳細図（ヒンジ図の表示）へ記録する材端軸力は、圧縮正・引張負の符号付きで保持します。材端集中ばねの N–M 相関（非線形梁要素のモデル化）は \( |N| \) を、せん断降伏耐力 \( Q_y \) の軸力項は圧縮側の大きさ（引張は 0）を用いるため、解析の挙動はこの符号に依存しません。
部材応答履歴を記録できない場合
ヒンジ詳細図に用いる部材応答履歴は単精度浮動小数点で保持します。軸力 [N]、両端の局所曲げモーメント [N·mm]、弦からの材端回転と材端集中ばねの回転 [rad] のいずれかが非有限値、または記録時の単精度への変換で非有限値になる場合、部材・ステップ・成分・理由を示して解析をエラーで停止します。その実行の途中までの応答を正常な解析結果として返しません。記録不能な値をゼロや最大値へ置き換えたり、記録できない履歴を除いて正常終了したりする扱いは行いません。
有限値として記録できる応答には通常の単精度丸めを適用し、単位・符号は保持します。極小値が丸めによってゼロとなる場合もあります。この数値記録の検査は、解析モデルや荷重の構造上の妥当性を保証するものではありません。
実装参照：sepika_solver::nonlinear::pushover::member_response::record_member_step（crates/sepika-solver/src/nonlinear/pushover/member_response.rs）と step_record::StepRecorder::record（crates/sepika-solver/src/nonlinear/pushover/step_record.rs）が実装します。
制御方式の切替（荷重制御 → 変位制御）
2 つの制御方式には、それぞれ釣合い経路を安定して辿れる領域があります。 増分解析は、経路の進行に合わせて有効な方式へ乗り換えることで、崩壊機構の形成域まで解析を進めます。
- 荷重制御（λ を 0 → 1 へ増分）: 構造が十分に硬い序盤。λ=1 は \( C_0 = 0.2 \) の設計地震力レベル
- 変位制御（頂部変位を増分し λ を解く）: 耐力の頭打ち・低下の前後
荷重制御は、荷重を入力とするため耐力のピーク（極限点）を越えられません。 ピークを越えた大きさの荷重には、釣合う解がそもそも存在しないためです。
保有耐力比較では、明示採用した確定点の層力を必要保有水平耐力と比較します。 必要保有水平耐力 \( Q_{un} = D_s \cdot F_{es} \cdot Q_{ud} \) の \( Q_{ud} \) は標準せん断力係数 \( C_0 = 1.0 \)（大地震レベル）で算定するため、同じ Z・地盤・略算周期で地下や指定 Ci の副剛床を含まない場合、増分解析の λ に換算すると \( 5 D_s F_{es} \)、\( D_s \) が下限の 0.25 でも λ≒1.25 に相当します。 そのため、λ=1 で終わる荷重制御だけでは判定に必要な耐力域へ届きません。
また、建物が保持できない大きさの荷重を与えても、釣合い解が存在せず反復が収束しないだけで、耐力の頭打ちや低下に至る釣合い経路の情報は得られません。 一方、頂部変位は崩壊機構形成の前後を通じて単調に増え続けます。 そのため、変位をパラメータに取る変位制御へ切替えることで、耐力の頭打ちと低下を釣合い経路上の応答として観測できます。解析経過のピークと、目的別に明示採用する Qu・Ds の評価点は区別します。 健全な建物であれば、変位制御の区間でも λ は単調に増加し、荷重の低下はヒンジ形成による剛性低下が生じてはじめて現れます。
なお、頂部変位そのものが経路上で一時的に戻るスナップバックだけは、変位制御でも追えません。 段階制御は、荷重制御と変位制御の 2 段で構成されます。
増分方式の選択（比較検証用）
増分方式は、既定の段階制御のほかに「荷重増分のみ」を選択できます（右パネル「増分解析」）。 荷重増分のみでは変位制御へ移行せず、終了目標が有効な場合は λ=1 を超えて同じ刻み dλ のまま荷重増分を継続します。 増分を半減しても収束しなくなった時点で解析を打ち切ります。 これ以上の荷重に釣合う解がない、すなわち耐力ピーク近傍に達したことを意味するためです。
荷重を入力とするため、この方式では耐力の頭打ち・低下域は追跡できません。 段階制御との結果の違い（どの λ まで釣合いが取れるか、変位制御が結果に与える影響）を確認する用途を想定しています。 保有水平耐力の判定に用いる結果は、崩壊機構の形成まで追跡できる段階制御（既定）で求めてください。 結果画面と CSV には、どちらの方式で得られた結果かを表示します。
変位制御が仮定するもの・しないもの
変位制御と聞くと、「λ=1 時点の変形形状を拡大していくため、ヒンジ発生による変形形態の変化を表現できないのではないか」という疑問が生じるかもしれません。 ですが、変位制御が拘束するのは頂部（最上階の代表節点）の加力方向変位という 1 自由度だけで、変形の形状は仮定しません。 残りの全自由度は、その時点の接線剛性（ヒンジによる剛性低下を反映した剛性）と外力 λ·q に対する釣合いから解かれます。 そのため、たとえば特定の層の柱にヒンジが入って層崩壊形へ向かう場合、同じ頂部変位の増分でも、その内訳は軟化した層の層間変位へ自動的に集中します。 この変形集中は、結果画面の層別の Q-δ 曲線（層間変位−層せん断力）で確認できます。
一方、固定しているのは外力の分布形状 q（Ai 分布）で、全フェーズを通じて λ でスケールするだけです。 ヒンジ形成により振動特性が変わると実際の地震時の慣性力分布も変わり得ますが、これには追随しません。 これは静的漸増載荷（プッシュオーバー）解析に共通の近似で、告示の保有水平耐力計算も Ai 分布による静的載荷を前提としています。
解析前の入力チェック
増分解析は、部材の終局耐力（曲げ降伏 \( M_y \)、せん断終局 \( Q_{su} \)、耐震壁の \( Q_u \)）で応力を頭打ちにすることで崩壊機構を形成します。耐力を算定できない部材は弾性のまま・降伏しないまま扱われ、押し込むほど際限なく応力を負担するため、崩壊機構が形成されず保有水平耐力を過大評価します（危険側）。
そのため、次の入力不足があるモデルは、既定値で補わずに解析をエラーで停止し、該当する部材 ID と是正内容を表示します。
| 対象 | 停止条件 |
|---|---|
| 耐震壁（壁エレメント） | 終局せん断強度 \( Q_u \) を算定できない（下記の一覧） |
| 線材（梁・柱・ファイバー・マルチスプリング） | 断面に材料が割り当てられていない |
| 同上 | 主材料の区分が鉄筋（構造種別の判定。線材の主材料には鋼材かコンクリートを割り当てます） |
| 同上（配筋を持つ断面: RC 矩形・RC 円形） | 主材料の区分が鋼材（鋼材として検定すると耐力を大きく過大評価します） |
| 同上（コンクリート系断面: RC・SRC・CFT） | 主材料にコンクリート強度 \( F_c \) が設定されていない、または \( F_c \le 0 \) |
| 同上（配筋を持つ断面: RC 矩形・RC 円形・SRC 矩形） | 主筋の材料が未割当、またはその材料に正の \( f_y \) がない |
| 同上（鋼材断面形状: 形鋼・鋼管など） | 鋼材の降伏強度 \( f_y \) を解決できない（SRC は内蔵鉄骨の材料、その他は主材料。ファイバー断面は鋼材の降伏進展を追うため必須） |
| 同上（断面形状が未設定） | 主材料に正のコンクリート強度 \( F_c \) も正の降伏強度 \( f_y \) も設定されていない（\( F_c \) があればコンクリート、なければ鋼材のファイバとして組み立てるため、使われる側の強度が必要） |
| 全部材（線材） | 断面の有効せん断断面積 \( A_s \) が 0（せん断有効断面積 As。全解析共通の入力チェック） |
\( F_c \le 0 \)・\( f_y \le 0 \) のような 0 以下の値は、未設定と同じ「耐力を算定できない」入力として扱います。
耐震壁で耐力を算定できないのは次の場合です。材種によらず次を検出します。
- 壁エレメントとして構築できない（4 節点の指定がない、節点座標が退化している）
- 断面が設定されていない、または壁厚が \( \le 0 \)
- 断面に材料が割り当てられていない
RC 壁では、これに加えて次を検出します。
- 材料にコンクリート強度 \( F_c \) が設定されていない、または \( F_c \le 0 \)
- 付帯柱（側柱）があるのに、その断面から主筋量を取得できない（断面形状が RC 矩形でない、主筋の本数・径が 0）
- 付帯柱がなく、かつ明示縦筋比から求める縦筋面積が 0
- 上記のいずれにも該当しないが \( Q_u \) の算定結果が \( \le 0 \) となる（寸法・開口の入力を確認してください）
以下は耐震壁として非線形化するRC壁の必須入力です。スリット等で既存の耐震条件を満たさないRC壁には、既存の線形・極小剛性モデルを明示採用し、必要入力の不足から非線形壁を線形へ切り替える扱いとは区別します。
初期せん断耐力は、同じ解析入力の長期荷重と壁方向の正負Ai荷重による静的一次参照から、解析開始前に確定します。この参照の適用範囲は単独無開口壁と対称同質の既存対応側柱に限られ、材料入力が揃っているだけでは開始可能になりません。全対象壁の正負耐力が未解決なら増分解析を開始せず、採用した耐力は解析中固定します。
RC壁の方向別配筋は小数で指定します。横筋比は実壁厚に対する横筋比、縦筋比は壁板の総縦筋面積を実壁厚と壁長の積で除した値です。直交最小筋比 ps は数量・質量用の既存入力として保持し、方向比の推定には使いません。横筋比を等価壁厚へ換算する係数 t/te は終局せん断式で一度だけ適用します。縦筋の総面積は縦筋比×t×lwで、壁長方向に20等分して配置します。
横筋には断面のせん断補強筋材料、縦筋には断面の主筋材料を割り当てます。両者は独立した鉄筋材料で、対応材種と有限正の明示fyが必要です。材料名だけの強度補完、コンクリートの強度倍率による鉄筋強度の推定は行いません。材料強度基準では、縦筋材料の倍率を縦筋fyへ一度だけ適用します。横筋のfyは主筋の倍率を適用せず使用します（横筋材料の倍率欄もこのせん断式には用いません）。標準配筋の元部位 main・band・stirrup を壁方向へ読み替えません。vertical・horizontal として明示した方向部位は、位置指定付き部位や異なるnative割当を一つの材料へ統合せず、解決不能として扱います。
方向比の未設定は明示0と区別します。0は有限非負の入力として保持します。側柱主筋で引張成分を与えられる横筋式単体は横0でも評価できますが、側柱のない式で縦筋面積0なら適用条件を満たしません。ファイバー断面は明示縦0を鉄筋なしとして生成できます。縦筋比・材料が欠落する場合は軸曲げファイバーを生成できず、非線形解析を開始しません。縦成分だけ不足した側柱付き壁の横筋式単体の評価と、非線形解析の開始可否は別です。
壁版ごとの採用コンクリート強度は展開後も物理壁版に対応します。同断面でも部材指定のFcが異なれば、その壁ごとの材料で評価します。準備計算は不足理由を「非線形解析不可」と表示し、線形解析は独立したモードとして利用できます。必要な成分が不足するRC耐震壁をQu=0や弾性壁柱へ置換して、増分解析・非線形時刻歴を継続しません。入力や材料の編集後は旧成果を再利用せず再計算します。
非線形RC壁の軸・曲げ初期剛性の鉄筋換算にも明示縦筋比を用います。線形モードに残るps換算剛性は別の調査対象であり、方向値不足の妥当性や安全側を保証するものではありません。
縦筋の等価分布と横筋入力は本文に示すモデルを採用します。引用規準原典への適合と現在の履歴則の実験妥当性は未確認です。方向別入力の採用は、終局式の選択、同一評価状態M/Q/N、側柱の役割、参考Qcの接続の妥当性を保証しません。
鋼板耐震壁では、次を検出します。
- 材料に降伏強度 \( F \) が設定されていない、または \( F \le 0 \)
- \( Q_y \) の算定結果が \( \le 0 \) となる（板厚・壁長さの入力を確認してください）
同じチェックは非線形時刻歴応答解析にも適用します（耐力が定まらない部材を弾性のまま解析すると応答を過小評価するためです）。許容応力度計算（線形弾性解析）は部材耐力を用いないため対象外です。
実装参照：sepika_element::factory::{nonlinear_input_issues, ensure_nonlinear_input} が判定し、crates/sepika-solver/src/nonlinear/pushover/driver.rs と crates/sepika-solver/src/dynamic/timehistory/nonlinear.rs が解析の冒頭で呼び出します。
材料強度（割増）
保有水平耐力計算の材料強度として、降伏強度に次の割増係数を乗じます （H12 建告第2464号の運用・各認定条件）。
| 用途 | 割増係数 |
|---|---|
鋼材（材料名から鋼材グレードを解決できる材料。SS400・SN490B 等） | 1.1 |
| 鋼材のうち 590N 級（SA440・TMCP440） | 1.05 |
| RC 主筋（RC壁は必須の明示fyを使用） | 1.1 |
| せん断補強筋 | 割増しない |
| 名称から解決できない直接入力材料 | 割増しない（下記の直接入力で指定可） |
- 鋼材の自動判定は標準材料の対応済み名称と共通です。前後空白は除き、大小文字を区別し、未知の接尾辞は認識しません。
- 材料ごとに割増係数を直接入力でき（材料タブの「割増」欄）、入力がある場合は 自動判定より優先してその値を用います。直接入力材料に割増を適用したい場合に指定します。
- 割増はプッシュオーバー解析（部材の集中ばね・ファイバーの組み立て、 曲げヒンジ・せん断降伏の判定閾値）にのみ適用し、許容応力度計算（一次設計）・ 時刻歴応答解析には適用しません。
実装参照：sepika_solver::nonlinear::pushover（crates/sepika-solver/src/nonlinear/pushover/mod.rs）が算定します。
材料強度の割増は部材組み立て（sepika_element::factory::build_nonlinear_behavior の
StrengthBasis::MaterialStrength）とヒンジ・せん断降伏閾値で適用し、係数は
sepika_core::material_grade::{material_strength_factor_steel, material_strength_factor_rebar}
が算定します。
変位制御は、Ai 分布の比例荷重パターン \( \lambda \boldsymbol{q} \) を保持したまま、荷重係数 \( \lambda \) を「頂部変位 = 目標値」の拘束条件から決定します（Batoz–Dhatt の変位制御法）。各反復で釣合い残差解 \( \delta u_r = K^{-1}(\lambda \boldsymbol{q} - f_{\text{int}}) \) と荷重パターン解 \( \delta u_q = K^{-1} \boldsymbol{q} \) を解き、
\[ \delta\lambda = \frac{u_{\text{target}} - u_{\text{roof}} - \delta u_r[\text{roof}]}{\delta u_q[\text{roof}]}, \qquad \delta u = \delta u_r + \delta\lambda \cdot \delta u_q \]
として頂部変位拘束と力の釣合いを同時に満たす \( \lambda \) を定めます。収束は外力ノルム基準の相対残差（\( 10^{-6} \)）と頂部変位の目標一致で判定し、収束しないステップは押込み増分を半減して再試行します。荷重制御から変位制御へ切替わっても載荷パターンは同一の Ai 分布のままであるため、荷重−変位曲線は連続で、ベースシアの低下は部材の剛性低下（ヒンジ形成・崩壊機構）によってのみ生じます。これにより解析は、設計地震力レベル（\( \lambda = 1 \)）を超えて崩壊機構形成域まで押し切ることができます。 崩壊機構判定は、静的不静定次数 \( r = 3m - 3n + r_{\text{support}} \)（\( r_{\text{support}} \) は加力方向の平面内拘束自由度数）に対して降伏ヒンジ数 \( \ge r+1 \) をゲートとします。ここで \( m \)・\( n \) は、加力方向の平面骨組に関与する部材とその接続節点の数です。具体的には 2 節点の線材（梁・柱・ブレース）のうち部材軸が加力直交の水平方向へ卓越しない（軸単位ベクトルの直交水平成分の絶対値が 0.707 以下の）ものを \( m \) に数え、それらが接続する節点を \( n \)・\( r_{\text{support}} \) の集計対象とします（3 次元モデルを加力方向の平面骨組へ射影して次数式を適用する扱いで、直交方向の梁や耐震壁は算入しません）。 層崩壊の判定に用いるヒンジの層帰属は、層を上端の階（床）で識別し、ヒンジ位置側と相手端の節点のうち高い側の所属階を採ります。柱脚のヒンジは最下層の柱脚を表すため、高い側（柱頭）により最下層へ帰属します。高い側も基部の階にある部材は基礎梁であり、どの層にも属さないものとして層崩壊の判定から除きます。
耐震壁の軸・曲げ
耐震壁（壁エレメントモデル）の壁柱は、柱と同様に既定でファイバー断面の弾塑性評価とします（コンクリート格子＋縦筋の等価分散配置。壁エレメントモデル）。細長壁の曲げ降伏・軸力変動による N–M 相関が増分解析に反映されます。
耐震壁の面内せん断
耐震壁（壁エレメントモデル）は、面内せん断を終局せん断強度 \( Q_u \) で頭打ちとする弾完全塑性として扱う（壁エレメントモデル、耐震壁）ため、壁が伝達する面内水平力は \( Q_u \) を超えません。\( Q_u \) を算定できない壁があるモデルは、上記「解析前の入力チェック」によりエラーで停止します。
側柱（耐震壁の鉛直辺に取り付く柱）は面内曲げ面の端部回転を静的縮約し、面内せん断を負担しません（許容応力度計算と同じ扱い）。壁のせん断断面には側柱断面を算入します。
鋼板耐震壁を含むモデルでは、面内せん断を鋼板のせん断降伏 \( Q_y = t \, l_w F/\sqrt{3} \) で頭打ちとします。せん断座屈を考慮していないため、解析の開始時に対象の枚数・要素 ID・生成元の壁版 ID とあわせて集約した注意事項を 1 回表示します。目標到達前の打切り注意がある場合も鋼板壁の注意は保持します（壁エレメントモデル）。
曲げヒンジ判定（Mc・My・終局）
曲げヒンジの発生レベルを判定します。 RC 曲げひび割れ（\( \kappa = 0.56 \)、技術基準解説書 P.621-623）:
\[ M_c = \kappa \cdot \sqrt{F_c} \cdot Z_e \]
RC 曲げ降伏（\( d = \text{有効せい} \)、同 P.623）:
\[ M_y = 0.9 \cdot a_t \cdot \sigma_y \cdot d \]
有効せい \( d \) は引張縁から引張筋重心までの距離 \( d_t \) を用いて \( d = D - d_t \) とします。
\( d_t \) はかぶり・せん断補強筋径・主筋径・多段配筋を考慮し、RC 断面検定と同じ
sepika_core::rc_rebar_geom の規約で算定します（曲げばねの降伏モーメントも同一）。
鋼 全塑性:
\[ M_p = Z_p \cdot \sigma_y \]
算定式：\( m_{\max} \ge M_y \) かつ \( \mu \ge {} \)ULTIMATE_DUCTILITY(=4.0) のとき終局、\( M_c \le m_{\max} < M_y \) のときひび割れとします。
記号・単位
| 記号 | 意味 | 単位 |
|---|---|---|
| \( M_c \) | RC 曲げひび割れモーメント | N·mm |
| \( M_y \) | 曲げ降伏モーメント | N·mm |
| \( M_p \) | 鋼の全塑性モーメント | N·mm |
| \( m_{\max} \) | 材端モーメントの最大値（危険断面） | N·mm |
| \( \kappa \) | RC 曲げひび割れの係数（0.56） | - |
| \( F_c \) | コンクリートの設計基準強度 | N/mm² |
| \( Z_e,\ Z_p \) | 弾性断面係数・塑性断面係数 | mm³ |
| \( a_t \) | 引張鉄筋断面積 | mm² |
| \( \sigma_y \) | 鉄筋・鋼材の降伏強度 | N/mm² |
| \( d \) | 断面の有効せい | mm |
| \( \mu \) | 部材塑性率 | - |
| \( Q_y \) | せん断降伏耐力 | N |
| \( A_s \) | せん断有効断面積 | mm² |
| \( f_y \) | 鋼材の降伏強度 | N/mm² |
| \( {}_rQ_{su} \) | SRC の RC 部の終局せん断耐力 | N |
| \( {}_sA_w,\ {}_sF \) | 内蔵鉄骨のせん断有効断面積・基準強度 | mm², N/mm² |
式の選択
配筋を持つ断面（RC 矩形・RC 円形）は、材料の区分に依らず上の RC の式を使います。 鋼の \( M_p = Z_p \sigma_y \) は配筋を無視した素の断面係数を使うため、RC 断面へ当てると \( M_y \) が桁で過大になり、曲げヒンジが検出されなくなるためです。 それ以外の断面は構造種別（構造種別の判定）で分け、 S は鋼の全塑性、複合断面（SRC・CFT）と形状を持たない断面は \( M_y = \sigma_y Z_e \)・\( M_c = \kappa \sqrt{F_c} Z_e \) の簡易値とします。 塑性断面係数 \( Z_p \) はすべての鋼断面形状（H 形・箱形・円形鋼管・組立 H 形・T 形・溝形・リップ溝形・山形・平鋼・中実丸鋼）について算定します（断面諸量）。
判定に用いる材端モーメントの位置
比較する材端モーメント \( m \) は、危険断面＝剛域フェイスの値を用います。要素が返す材端力は グローバル系・節点位置のため、(1) 要素局所座標系へ回転し、(2) 剛体アームのモーメント （材端せん断 × アーム長）を差し引いて求めます（剛域）。 節点位置のモーメントは剛体アームぶん大きく、断面耐力 \( M_y \) と直接比較すると、剛域を持つ 部材のヒンジを過早に検出し、崩壊荷重を過小評価します（両端固定柱では \( m_{\text{node}} = Q \cdot L/2 \) となり、判定が剛域長に依存しなくなります）。 局所成分を用いることで、材軸まわりのねじりを曲げと取り違えることもありません。
材端ごとの記録
ヒンジは i 端・j 端それぞれについて判定・記録します。崩壊機構の運動学的ゲート （降伏ヒンジ数 \( \ge r+1 \)）は (部材, 材端) の重複を除いて数えるため、1 部材につき モーメント最大側の 1 端しか記録しないと、両端が降伏した柱も 1 個としか数えられず、 機構が成立していても部分崩壊形と判定されます。
判定対象
曲げヒンジ・せん断降伏の判定対象は 2 節点の線材（梁・柱・ブレース）です。 耐震壁（壁エレメントモデル、4 節点）は線材の材軸（i 端 → j 端）が定義できない ため対象外とし、壁の非線形性は壁要素自身の弾塑性 （壁柱ファイバーの軸・曲げと、終局せん断強度 \( Q_u \) による面内せん断の 頭打ち。上記「耐震壁の軸・曲げ」「耐震壁の面内せん断」）で表現します。
実装参照：sepika_solver::nonlinear::pushover::hinge::{member_moment_thresholds, track_hinges}（crates/sepika-solver/src/nonlinear/pushover/hinge.rs）が算定します。
危険断面への変換は sepika_solver::nonlinear::pushover::geom::member_end_forces_at_face
（crates/sepika-solver/src/nonlinear/pushover/geom.rs）が
sepika_element::frame::rigid_arm::to_flex_force を用いて行います。
RC 主筋の σy は断面の主筋材料の \( f_y \) から解決します（未割当のモデルは解析前チェックで停止します）。鋼材 σy は断面の主材料の \( f_y \)、ない場合は既定 235 とします。
部材塑性率（ductility）の 3 方式
部材塑性率を 3 方式（基点ひずみ・重み付け平均・初降伏）で算定します。
算定式
方式(1) は基点ひずみによります。RC 引張 0.01、圧縮 0.005、鉄骨 0.01 を超えた曲率を基点とし、次とします:
\[ \mu = \kappa_{\max}/\kappa_{\text{ref}} \]
方式(2) は重み付け平均が次を満たす時点を基点とします:
\[ J_m = \frac{\sum \sigma_y \cdot A \cdot |\varepsilon| \cdot \mu_i}{\sum \sigma_y \cdot A \cdot |\varepsilon|} \ge 1.0 \]
方式(3) は初降伏によります。
実装参照：sepika_solver::nonlinear::pushover::ductility::compute_ductility_refs（crates/sepika-solver/src/nonlinear/pushover/ductility.rs）が算定します。
せん断降伏判定（Qy）
せん断降伏を判定します。 鋼系は次によります（von Mises）:
\[ Q_y = A_s \cdot f_y/\sqrt{3} \]
RC 矩形断面は荒川 mean 式系（rc_qsu_simple、軸力を各ステップで反映）によります。せん断補強筋の \( \sigma_{wy} \) は断面のせん断補強筋材料の \( f_y \) から解決し、未割当の場合のみ SD295 相当（295 N/mm²、材料未割当時の代表値）を用います。対応グレードの材料を割り当てているのに \( f_y \) が未設定（0 以下を含む）の場合は解析前チェックのエラーで停止します。
SRC 矩形断面は、RC 部と内蔵鉄骨の累加式によります（SRC 規準の累加強度の考え方）:
\[ Q_y = {}_rQ_{su} + {}_sA_w \cdot {}_sF/\sqrt{3} \]
- \( {}_rQ_{su} \): RC 部の荒川 mean 式系。鉄骨を控除しない全断面 \( b \cdot D \) と配筋で評価し、軸力の各ステップ反映も RC 矩形と同一です。
- \( {}_sA_w \): 内蔵鉄骨のせん断有効断面積。強軸（局所 y）はウェブ内法 \( t_w (H - 2t_f) \)、弱軸（局所 z）は上下フランジ \( 2 B t_f \)（全塑性評価のため応力分布係数による低減は行いません）。
- \( {}_sF \): 内蔵鉄骨の基準強度。鋼種名の前方一致で当該板厚の区分から解決し、解決できない場合は 235 N/mm² とします。材料強度割増（直接入力係数優先、なければ鋼種名から判定: 鋼材 1.1・590N 級 1.05・鋼種名を解決できない場合は割増なし 1.0）を乗じます。
材料に降伏強度 \( f_y \) が設定されていても（主筋 \( \sigma_y \) の解決用など）、断面形状から精算できる RC 矩形・SRC 矩形は常に上記の式で評価します（\( f_y \) による鋼系の式 \( Q_y = A_s f_y/\sqrt{3} \) は、形状から精算できない断面のみに用います）。
SRC 矩形以外のコンクリート系断面（形状情報がない場合等）は \( Q_y = A_s \cdot 0.7\sqrt{F_c} \)（簡易慣用値）へフォールバックします。
材料強度が未入力で \( Q_y \) を算定できない部材、および有効せん断断面積 \( A_s \) が 0 の断面を用いる部材は、解析前チェックにより解析が停止します（上記「解析前の入力チェック」、せん断有効断面積 As）。
実装参照：sepika_solver::nonlinear::pushover::shear_yield::build_dir_threshold（crates/sepika-solver/src/nonlinear/pushover/shear_yield.rs）が算定します。
RC主筋・せん断補強筋の材種別基準値は「要・原典照合」と明記します。RC壁解析は必須の明示fyを使用し、345・295の欠落補完は行いません。
局所せん断を ey/ez へ射影して方向別に判定します。
確定状態の層力履歴と目的別採用
性能曲線の全支持抵抗と骨組負担せん断力
性能曲線のベースシアは、同じ確定状態で固定支持と地盤支持ばねに伝わる水平反力を合計し、 載荷方向に符号をそろえた全支持抵抗です。固定自由度へ直接作用する節点荷重は、部材に変形を 生じない直接支持作用として控除します。部材端点ちょうどの集中荷重も同じ節点荷重として扱います。 固定した代表節点への既存の剛体接続などによって、変位0へ消去される従属自由度も固定と同じように扱います。 端点から離れた集中荷重と分布荷重は、拘束端へ直接伝わる分も含めて部材荷重として回復します。 例えば長さ3 mの両端固定部材へ水平に2 kN/mを載せた場合、変位は0でも両支持の合計抵抗は6 kNです。
性能曲線の層別値は、上端階に属する節点を上層から累積した範囲で評価します。 その節点群の物理材端力と地盤支持作用に、対象層下端SLより上・最上端SL以下に作用する 部材内部の荷重を加えます。部材端点の集中荷重は節点所属で数え、分布荷重は高さで 切り取った区間を積分します。所属階を持たない節点の荷重や最上端より外の荷重がある場合、 層別値と建物全体のベースシアは一致するとは限りません。
二次設計に採用する骨組負担せん断力は、層内水平切断面を横切る物理部材の力です。 柱のせん断力とブレース軸力の水平成分を含み、切断面上側の地盤支持ばねや固定支持の分担は 上側自由体の外部作用として差し引きます。水平外力16 kNを骨組と支持ばねが8 kNずつ担う場合、 性能曲線の全支持抵抗は16 kN、骨組負担せん断力は8 kNです。
元の節点荷重と部材荷重から上側自由体の外力を求め、同じ確定状態の部材力・支持作用と照合します。 剛床代表節点に作用する荷重も元荷重として含みます。 切断位置ちょうどの部材集中荷重は部材i端側から近づいた断面力の極限に対応させるため、 上下端を反転すると切断位置で選ぶ側は変わります。性能曲線の階帯の境界とは別に扱います。
実装参照：sepika_solver::nonlinear::pushover::physical_response::PhysicalStep（crates/sepika-solver/src/nonlinear/pushover/physical_response.rs）、capacity_response（crates/sepika-solver/src/nonlinear/pushover/capacity_response.rs）、story_response（crates/sepika-solver/src/nonlinear/pushover/story_response.rs）が回復と照合を担います。
収束して確定した各 step について、実要素の線材応答、符号付き層切断面力、参照外力、上側の固定支持・支持ばねの外部作用、実外力累積と釣合い許容差を記録します。壁単体の履歴と同じ run・入力世代・方向・step を使います。性能曲線の最大ベースシアは解析経過の指標であり、二次設計の目的別採用点とは区別します。
Ds 判定と保有耐力比較の点は、保有水平耐力の判定に従って明示します。MCP の増分解析では任意の ds_step／capacity_step（0 始まり）を指定でき、未指定は未判定です。解析条件と目的別識別、集計値又は集計不能理由を出力します。通常の終了目標は解析制御の条件であり、法定の Ds 判定停止点を意味しません。
GUI と CSV は性能曲線全体のピークを「解析経過の最大ベースシア」と表示します。MCP の max_base_shear_kN も同じピーク値（kN）です。保有耐力比較に使う Qu は、明示した capacity_step の capacity_story_evaluation に記録する層切断面力 qu_n（N）から採用します。例えば全体ピーク150 kN、比較点120 kNなら、それぞれ150 kNと120 kNを別の意味で出力します。評価点が未指定のときも解析経過ピークは残りますが、比較用の Qu は確定せず、MCP の目的別集計は null です。
結果画面の「解析経過の層別最大せん断力（各層のピーク）」は、各層の性能曲線全体での絶対せん断力の最大値です。各層のピークは同じ step に生じるとは限らず、この表示を層別の採用 Qu として組み合わせません。設計表の Qu は保有耐力比較用に選定した同一確定点の層切断面力です。