非線形梁要素のモデル化（材端集中・ファイバー・MS）

非線形解析（プッシュオーバー・時刻歴応答）では、線形の弾性ティモシェンコ梁（ティモシェンコ梁要素（12×12））に代えて、部材の弾塑性挙動を表す 3 種類の梁要素モデルを用います。材端集中ばねモデルは弾性ティモシェンコ梁を土台とし、ファイバー・MS モデルは Timoshenko 適合内挿（φ 依存形状関数、ファイバーモデル）によりせん断変形を曲げと直列に合成します（剛域なしでは弾性ティモシェンコ梁と厳密に一致し、剛域ありでは軸方向を除き一致します）。

材端集中ばねモデル

弾性ティモシェンコ梁の両端に**回転ばね（M–θ 履歴則、部材履歴則の設定と集中ばね履歴則）**を直列接続し、塑性化を材端の集中回転に代表させます（構造動力学で標準的な集中ばね（one-component / two-component）モデル）。

算定式

• 回転ばねは強軸曲げ面（要素局所系の \( r_z \) まわり。要素局所 y 軸＝断面のせい方向とする規約による）に入れます。一軸曲げ（ForceRegime::UniaxialBendingShear）のモデルであり、骨格の降伏モーメント（強軸 \( Z_p \sigma_y \)）・初期剛性 \( 6EI/l \)（\( I \) は強軸）も強軸で与えます。
• 要素剛性は、可撓長（部材長から剛域を除いた長さ）で組んだ弾性ティモシェンコ梁に、端部条件（ピン・半剛）と材端回転ばねを直列に静縮約し、最後に剛域変換で節点自由度へ移して得ます（弾性梁 ティモシェンコ梁要素（12×12） と同じ土台）。ばねは直列のため接合部ばねと塑性ばねの順序は結果に影響しません。
• ばねの変形量は相対回転 \( \gamma_k = \theta_{n,k} - \theta_{b,k} \)（節点回転 − 可撓端回転）とし、履歴則 \( M_s(\gamma) \) を適用します。可撓端回転 \( \theta_b \) は内部自由度として保持し、各トライアルで両端の内部釣合い

\[ R_k = \bigl[K_{\text{flex}}\,\hat{\boldsymbol u}\bigr]_{d_k} - M_{s,k}(\gamma_k) = 0 \]

（\( d_k \) は材端回転自由度の番号、\( \hat{\boldsymbol u} \) は回転スロットを \( \theta_b \) で置き換えた可撓端変位）を要素内 Newton 反復（2 自由度連成）で解きます。履歴則は区分線形のため通常数回で収束します。固定端（節点回転 0）でも可撓部端モーメントが降伏モーメントへ達すればばねが塑性回転し、柱脚ヒンジが形成されます。

• 復元力は弾性可撓部の \( K_{\text{flex}}\,\hat{\boldsymbol u} \) とばねの履歴モーメント \( M_s(\gamma) \) から経路整合に評価します（節点の回転自由度にはばねを介してモーメントが伝わるため、回転スロットはばね側の履歴力で置き換えます。内部釣合いの解では両者は一致します）。接線剛性はばねの接線 \( k_t(\gamma) \) を直列に静縮約した \( K^{*} \) を用います。
• 履歴則は部材個別指定（履歴則の設定（増分用・時刻歴用）の 2 系統）→ 構造種別ごとの既定表で解決します。RC/SRC/CFT 梁は武田型トリリニア、S 梁は**標準型（kinematic バイリニア）**を材端ばねに用います（増分・時刻歴共通の既定）。
• 軸-曲げ相互作用（N–M 相関）は、現在軸力 \( N \) に応じて回転ばねの降伏モーメントを

\[ M_{\text{lim}} = M_{y0} \cdot \left(1 - \frac{|N|}{N_{\text{allow}}}\right) \]

で更新します（下限 \( 0.02 \cdot M_{y0} \)）。この相関は降伏モーメントを更新できる材料（標準型・辻・山田型・座屈考慮型）に適用し、武田型等の骨格固定の履歴材料は対象外とします。

• 材端ばねの骨格（トリリニア）の既定値は次のとおりです。終局点は \( M_u = 1.1 \cdot M_y \)、終局回転角は降伏回転角の 4 倍（塑性率 4）とします。ひび割れ点は RC 系が \( M_c = 0.56 \cdot \sqrt{F_c} \cdot Z_e \)（梁（曲げ・せん断・軸）。\( M_y \) の 10〜90% にクランプし、折れ点の逆転を防ぎます）、それ以外は \( M_y/3 \) とします。降伏時剛性低下率 \( \alpha_y \) は RC 矩形梁（水平材）を菅野式、それ以外を 0.3 とします。武田型の除荷剛性低下指数は 0.4 とします。
• 辻・山田型を指定した材端ばねは、二次勾配 \( K_2 = 0.01 \cdot k_{rot} \)・混合硬化比 \( \beta = 0.5 \) のバイリニア＋混合硬化とし、N–M 相関を適用します。座屈考慮型を指定した材端ばねは終局点を \( M_u = 1.1 \cdot M_y \) とし、N–M 相関を適用します。標準型・降伏モーメント不定の材端ばねは kinematic バイリニア（降伏後剛性比 0.01）とし、N–M 相関を適用します。

記号・単位

ファイバーモデル

積分点（ガウス点）の断面を多数のファイバーに分割し、平面保持仮定のもとで各ファイバーが独立した履歴状態を持つ一軸材料として応答することで（構造力学のファイバーモデル）、軸力と曲げの連成（N–M 相関）を断面積分から自然に得ます。断面応答（M–φ 積分）の詳細は ファイバー断面（M–φ 積分）を参照してください。

算定式

• 非線形解析で構築するファイバー要素は塑性化域考慮モデルとし、「全可撓長の弾性梁 \( K_{el} \) ＋ 端部の塑性増分ヒンジ」の直列モデルとして解きます。基準長は可撓長 \( L' \) です（剛域があれば節点間長から剛域長を除いた長さ。非線形梁要素の剛域）。塑性化域長 \( L_p \) は部材の指定値、未指定なら断面せいの 0.5 倍（マルチスプリング（MS）モデルと共通の既定。断面せいが取れない場合は 200 mm を仮定）とし、可撓長の 45% までにクランプします。したがって塑性化を担うファイバー断面は可撓部両端の 2 断面（\( \xi = \mp 1 \)、すなわち可撓部の材端＝剛域があればその剛域フェイス。非線形梁要素の剛域）であり、その断面の状態のみを追跡します。
• 弾性梁 \( K_{el} \): \( B(\xi)^{\mathsf T}\,\mathrm{diag}(EA, EI_y, EI_z)\,B(\xi) \) を 2 点ガウス（\( \xi = \mp 1/\sqrt{3} \)）で厳密積分し、一定せん断ひずみ場（後述）と Saint-Venant ねじり剛性を加算して組み立てます（\( EA, EI_y, EI_z \) は公称断面諸元）。弾性状態（ヒンジ回転 \( \gamma = 0 \)）の要素剛性はこの \( K_{el} \) に厳密一致し、弾性ティモシェンコ梁（ティモシェンコ梁要素（12×12））の剛性と一致します。
• 塑性増分ヒンジ: 可撓部両端のファイバー断面（\( \xi = \mp 1 \)、両端×2 軸（\( \kappa_y, \kappa_z \)）で最大 4 自由度）の曲率 \( \kappa \) を内部未知数とし、各トライアルで両端の内部平衡

\[ m_B = D_{\text{nom}} \cdot B(\xi_{\text{end}}) \cdot \hat{\boldsymbol u} = m_{\text{sec}}(\varepsilon_0, \kappa) \]

（\( D_{\text{nom}} = \mathrm{diag}(EA, EI_y, EI_z) \) は公称断面諸元、\( m_{\text{sec}} \) は端部断面のファイバー積分によるモーメント）を、要素内 Newton 反復（整合接線を用いた区分線形問題のため通常数回で収束）で解きます。

• ヒンジ回転（節点回転と可撓端回転の差）は、断面の弾性線を超える塑性超過分のみを持ちます。

\[ \gamma = s \cdot L_p \cdot \left( \kappa - \frac{m_{\text{sec}}(\kappa)}{EI_{\text{sec}}} \right) \]

（\( s \) は端の向き: i 端 +1・j 端 −1）。弾性では \( \gamma = 0 \) となり、\( EI_{\text{sec}} \)（端部ファイバー断面の弾性曲げ剛性。ファイバー離散化後の値）による弾性整合が成立します。

• 可撓端変位 \( \hat{\boldsymbol u} \) は、有効端（剛接端）の回転スロットを可撓端回転 \( \theta_b = \theta_n - \gamma \)（\( \theta_n \) は節点回転側の値）で置き換えた変位です。復元力は \( \boldsymbol f = K_{el}\,\hat{\boldsymbol u} \) から履歴に整合して評価します。

• 整合接線は

\[ K^{*} = K_{el} - \bigl(K_{el}[:,\text{slots}]\cdot G\bigr)\cdot J^{-1}\cdot(D\cdot B\text{行}) \]

（\( G = \partial \theta_b/\partial \kappa \) はヒンジ回転の曲率感度、\( J \) は内部平衡のヤコビアン）です。内部変数（\( \kappa \)）消去による整合接線は一般に非対称になり得るため、対称ソルバ（Cholesky）前提の全体組立に合わせて対称化して用います（復元力は消去前の厳密値のため収束解は変わりません）。

• 軸ひずみ \( \varepsilon_0 = (\hat{\boldsymbol u} \text{ の軸方向差})/L' \) は弾性のまま断面へ渡します（軸方向の塑性変形は考慮しません）。断面はこの \( \varepsilon_0 \) と曲率 \( \kappa \) から各ファイバーのひずみを決めて積分するため、軸力と曲げの連成（N–M 相関）は断面のファイバー積分で自然に考慮されます。

• ヒンジは剛接端（EndCondition::Fixed）のみ有効です。ピン・半剛端は材端解放（非線形梁要素の材端条件（ピン・半剛））の内部自由度で扱い、塑性ヒンジは形成しません。

• 塑性化域を設けない全長ファイバー積分モデル（2 点ガウス積分、\( \xi = \mp 1/\sqrt{3} \)、重み 1）も要素として構築できますが、非線形解析の要素生成では用いません。

• ファイバーは断面形状どおりの領域に配置します（MN 相関曲面と同じ配置規則を共用）。H 形はフランジ＋ウェブ、角形鋼管・鋼管は管壁のみ、CFT は管壁＋充填コンクリート、RC 円形は円形領域、SRC はコンクリート＋主筋＋内蔵鉄骨に配置するため、中空・薄肉断面の断面積・剛性・耐力が正しく表現されます。非対称断面（山形・溝形・T 形等）は断面積重心へ座標補正します。

• RC 断面はコンクリートファイバー格子に加え、主筋を点ファイバー（Menegotto–Pinto 鋼材）として分離配置します。主筋点の縁からの距離・多段間隔は断面検定の有効せいと同じ rc_rebar_geom 規約です（ファイバー断面）。コンクリートは Fc≤60 では材料から解決した E を初期接線へ反映した NewRC、\( f_c > 60 \) では NewRC 固有 Ec の暫定外挿として扱います。鋼材領域（形鋼・鋼管・内蔵鉄骨）・主筋はいずれも Menegotto–Pinto モデル（降伏後剛性比 \( b = 0.01 \)）です。

ファイバー材料の既定値と強度の解決順

ファイバー断面の各領域には、次の既定の構成則・パラメータを割り当てます。 強度は表の解決順で決め、解決できないモデルは非線形解析の入力チェックが解析前に エラーで停止します（弾性で代替しません。ファイバー断面は降伏・ひび割れの進展を 追うことが目的のため、弾性代替は耐力の過大評価＝危険側になります）。

• コンクリートの初期接線は、Fc≤60 では材料から解決した E、Fc>60 では NewRC 固有 Ec の暫定外挿を用います。
• 引張強度は両コンクリート構成則とも既定 \( f_t = 2.0 \) N/mm² とします。
• コンクリートの除荷則は部材ごとに増分用・時刻歴用を選択できます（履歴則の設定（増分用・時刻歴用））。柱等のファイバー要素・MS 要素の自動（未指定）の既定は、増分解析=逆行型、時刻歴応答解析=Karsan–Jirsa 型です。壁柱ファイバーの自動既定は、増分解析=原点指向型、時刻歴応答解析=Karsan–Jirsa 型です。いずれも部材ごとに指定できます。 Fc≤60 は材料から解決した E を初期接線へ反映し、Fc>60 は NewRC 固有 Ec を用いて外挿します。逆行型・原点指向型は ConcreteNewRc、Karsan–Jirsa 型は ConcreteCyclic::newrc に対応します。
• Karsan–Jirsa 型を用いる場合の骨格・引張の既定は次のとおりです: Fc によらず NewRC 骨格（終局ひずみ \( \varepsilon_{cu} = 0.01 \)、 超過は \( \varepsilon_{cu} \) の曲線値を残留として保持）を用います。\( f_c > 60 \) は NewRC の適用範囲外への外挿です。引張は \( f_t = 2.0 \) N/mm²、 軟化勾配 \( E_{ts} = E_0/10 \)（\( E_0 \) は初期接線剛性）です。
• Menegotto–Pinto の既定パラメータは \( b = 0.01,\ R_0 = 20,\ a_1 = 18.5,\ a_2 = 0.15,\ a_3 = 0 \)（等方硬化なし）\( ,\ a_4 = 1 \) です（Menegotto–Pinto 鋼材・鉄筋モデル）。
• 材料強度の割増は保有水平耐力計算（材料強度基準）でのみ適用し、鋼材文脈・RC 主筋文脈それぞれの材料グレード別割増係数を強度に乗じます。時刻歴応答解析等は割増なし（1.0）です。

ファイバー分割の既定

• 断面形状がある場合: 目標ファイバ寸法 = 断面最大寸法 / 16（下限 1 mm）。円形・円環断面は周方向 24 分割、径方向は薄肉 2 分割・中実 8 分割。
• 断面形状がない場合（壁エレメントの壁柱以外の矩形近似）: 幅 12 × せい 20 の中実格子。
• せん断変形は Timoshenko 適合内挿（φ 依存の曲率形状関数＋一定せん断ひずみ場による変位法内挿。interdependent interpolation 系の定式化）で曲げと直列に合成します。せん断変形係数

\[ \varphi = \frac{12EI}{GA_s L^2} \]

は公称断面諸元（Section の \( I \)・\( A_s \) と Material の \( E \)・\( G \)）から算定して凍結するため、降伏後も内挿は弾性時の配分を保ちます（曲げの Hermite 内挿と同型の近似）。公称値を用いるのは線形解析（ティモシェンコ梁要素（12×12）の弾性 Timoshenko 梁）の φ と一致させるためであり、RC 断面ではファイバーの実効初期剛性が公称値と僅かに乖離するため、φ は公称値ベースの近似となります。

• 曲率場は回転項を \( (1 \pm 3\xi) \to (1 \pm 3\xi + \varphi) \)、全体を \( 1/(1+\varphi) \) 倍に補正します。せん断ひずみは曲げ面ごとに要素内一定の

\[ \gamma_y = \frac{\varphi_y}{2(1+\varphi_y)} \left( \theta_{z,i} + \theta_{z,j} - \frac{2\,\Delta u_y}{L} \right), \qquad \gamma_z = \frac{\varphi_z}{2(1+\varphi_z)} \left( \theta_{y,i} + \theta_{y,j} + \frac{2\,\Delta u_z}{L} \right) \]

（\( \Delta u \) = 端部並進差。符号は各曲げ面の曲率規約に従い面ごとに異なる点に注意してください。いずれも剛体回転で恒等的にゼロとなる客観的な測度）とし、弾性せん断剛性 \( GA_s \) を作用させます。曲率剛性とせん断剛性の和は、断面剛性が φ の算定基礎と一致する一様弾性断面で、弾性 Timoshenko 梁（ティモシェンコ梁要素（12×12））の剛性と厳密に一致します（被積分関数は \( \xi \) の 2 次のため 2 点ガウスで厳密です）。

• \( GA_s \le 0 \)（せん断有効断面積が未設定等）の場合は \( \varphi = 0 \)（せん断剛直 = Euler-Bernoulli）へフォールバックするため、幾何剛性・整合質量は Euler 系の係数を用います（ティモシェンコ梁要素（12×12）の弾性梁と同じ扱いです）。
• 整合質量の質量特性は、断面の材料領域から求めた SectionMassProperties を弾性梁と共有します。単位長さ当たり質量と断面回転慣性を用います。Lumped 質量は既存契約どおり、ファイバー面積の総和に主材料密度を掛けて求め、整合質量の入力とは分離します。
• 軸方向と曲げ面の質量は Timoshenko 形状関数から積分し、材軸まわりの回転自由度には質量用極二次モーメント \( \mu I_p = \mu I_y + \mu I_z \) を与えます。ねじり剛性のねじり定数 J は質量用回転慣性へ流用しません。
• 剛域がある場合は、可撓部の質量を剛体アームで節点自由度へ変換し、剛域の分布質量を加えて全長の質量を保存します。集中質量を並進自由度へ分配する扱いは弾性梁と同じです。
• 端部解放を含む整合質量は、塑性状態ではなく初期弾性剛性から固定した変換を用います。剛性側と同じ拡大自由度で ( S=-K_{bb}^{-1}K_{ba} )、( R=[I;S] ) を求め、( M=R^T\bar{M}R ) として 12×12 へ縮約します。回転ばねは無質量です。K_{bb} が特異な場合は解放条件を満たす質量を構成できないため、解放なし質量へフォールバックせず、Beam/Fiber の質量行列入口で明示的にエラー通知します。

マルチスプリング（MS）モデル

部材端の塑性化領域（長さ \( L_p \)、既定は断面せいの 0.5 倍）の断面を、少数の軸方向ばね群で置換し、中央は弾性材で連結します（RC 柱・梁端部の塑性ヒンジを軸ばね群で表す MS モデル）。軸ばね群の合力が \( N \)、図心まわりの偶力が \( M \) となるため、N–M 相関を自然に表現できます。ねじりばねは別系統で付加します。

算定式

• 端部断面のばね配置は幅方向 × せい方向 = 2 × 5 = 10 本の 2 次元配置とし、二軸曲げに対応します。
• 実体は塑性化域考慮ファイバー要素（ファイバーモデルの build_plastic_zone）と同一の定式化であり、端部断面の分割数だけが粗いです（ばね = 粗いファイバー）。したがって trial/commit/rollback・チェックポイントもファイバー要素と同一機構に乗ります。
• 妥当性は、一定軸力下の繰り返し曲げにおける N–M 相関・履歴ループがファイバー要素と整合することで担保します（相互検証）。
• せん断変形は ファイバーモデルの Timoshenko 適合内挿と同一です（実体がファイバー要素のため定式化も共有します）。

モデルの自動選択（ForceRegime）

ElementKind::Beam の要素は、部材の ForceRegime に基づいて非線形要素種別を選択します。

選択規則

• ForceRegime::UniaxialBendingShear → 材端集中ばね
• ForceRegime::AxialBendingInteract → ファイバー
• ForceRegime::Auto → トポロジから判定します（ヒューリスティック）。剛床に所属する水平材（軸力変動が小）→ 材端集中、それ以外（柱などの鉛直材）→ ファイバー。鉛直材判定は部材軸単位ベクトルの \( |e_z| > 0.707 \)（水平から 45° 超。線形静解析の柱判定・偏心率算定と共通の規約）によります。

MS モデルはこの自動選択には乗らず、要素種別 ElementKind::MultiSpring を明示指定した場合にのみ用います。

線形弾性解析は常に弾性 BeamElement を用い、非線形解析でのみ上記の非線形要素を構築します。両者を分けるのは、Auto 判定が剛床に乗らない梁もファイバーへ振り分けるため、共通化すると線形解析の弾性梁まで非線形要素へ置き換わってしまうためです。

非線形梁要素の剛域

ファイバーモデル・MS モデル（いずれも実体は FiberBeam）も、弾性梁（剛域）と同じ剛体アームで剛域を扱います。断面積分・せん断・幾何剛性を、仕口パネル要素の剛域長への反映で定義する可撓長 \( L' \) で組み立て、可撓端自由度を \( T_r \) で節点自由度へ写します。

\[ \boldsymbol{u}_{\text{flex}} = T_r\, \boldsymbol{u}_{\text{node}}, \qquad K_{\text{node}} = T_r^{T}\, K_{\text{flex}}\, T_r, \qquad \boldsymbol{f}_{\text{node}} = T_r^{T}\, \boldsymbol{f}_{\text{flex}} \]

\( T_r \) の成分と剛域・剛体アーム長の解決規則は 剛域 と同一です。

算定式

• ひずみ－変位行列 \( B(\xi) \)・せん断変形係数 \( \varphi \)・一定せん断ひずみ場・弾性梁 \( K_{el} \) の積分は、いずれも \( L' \) を基準長とします。
• 端部のファイバー断面（\( \xi = \mp 1 \)、塑性増分ヒンジを担う断面）は可撓部の材端、すなわち剛域フェイスに位置します。塑性化域 \( L_p \) も剛域フェイスから材内側へ測り、\( L' \) の 45% までにクランプします（ファイバーモデル）。塑性ヒンジは節点位置ではなく部材フェイスに形成されます。
• 断面ひずみは節点変位ではなく可撓端変位 \( \boldsymbol{u}_{\text{flex}} \) から評価します。剛体回転（節点回転＝弦回転）では可撓端に相対たわみが生じないため、内力はゼロになります（客観性）。
• ねじりは断面積分と連成しない独立項であり、剛体アーム変換がねじり自由度に作用しないため、弾性梁と同じく節点間長基準 \( GJ/L \) とします（剛域では増大させません）。

軸方向の扱い（弾性梁との違い）

弾性梁（剛域）は軸断面積を \( A \cdot (L'/L) \) に補正して軸剛性を \( EA/L \)（節点間長基準）へ戻し、剛域を曲げのみに効かせます。ファイバー要素ではこの補正を行わず、軸剛性は \( EA/L' \)（剛域を軸方向にも剛とする剛体オフセットの扱い）になります。

理由は、ファイバー断面の積分が軸力と曲げを連成させることにあります。断面は軸ひずみと曲率から各ファイバーのひずみ \( \varepsilon = \varepsilon_0 - \kappa_z\, y + \kappa_y\, z \) を決め、その応力を積分して \( (N, M_y, M_z) \) を同時に返します。軸成分だけを別の基準長へ補正すると、断面が返す軸力が断面のひずみ状態と整合しなくなり、N–M 相関（ファイバーモデルの要点）が崩れます。したがって軸方向は剛体アームの運動学にそのまま従わせます。

この結果、剛域をもつ部材では、ファイバー要素の弾性状態の剛性は軸自由度を除いて弾性 Timoshenko 梁と厳密に一致し、軸剛性のみ \( L/L' \) 倍だけ大きくなります。

非線形梁要素の材端条件（ピン・半剛）

ファイバーモデル・MS モデルも、弾性梁（端部条件（剛接・ピン・半剛）の反映）と同じ定式化で材端の回転解放を扱います。剛接端は「節点回転 = 要素端回転」を厳密に満たすため内部自由度を作りません。ピン・半剛の端のみ、要素端回転を内部自由度へ分離し、節点回転との間に回転ばね \( k_s \) を挟んで静縮約します。

• ピン: \( k_s = 0 \)（厳密なモーメント解放）
• 半剛: \( k_s = k_\theta \)（接合部の回転剛性 [N·mm/rad]）

算定式

• 縮約は可撓長で組んだ要素剛性に対して行います。組立順は「可撓長で要素剛性 → 材端解放の静縮約 → 剛体アーム変換 → 全体座標変換」で、弾性梁（局所自由度と組み立ての全体像）と同一です。縮約は \( K^{*} = K_{aa} - K_{ab} K_{bb}^{-1} K_{ba} \) により節点自由度 12 へ戻します。
• 内部自由度 \( \boldsymbol{u}_b \)（解放した要素端回転）は状態として保持し、各ステップで内部釣合い

\[ R_k = f_{\text{elem}}[d_k] + k_{s,k}\,\bigl(u_{b,k} - u_{\text{flex}}[d_k]\bigr) = 0 \]

を満たすよう要素内 Newton 反復で解きます（\( d_k \) は解放した回転自由度の番号、\( \boldsymbol{f}_{\text{elem}} \) は要素内力）。弾性域では線形のため 1 反復で厳密に収束します。断面が降伏した後も内部釣合いを解き直すため、ピン端のモーメントゼロは非線形域でも保たれます。

• 解放端の節点回転が受け持つのは回転ばねのモーメントのみであり、ピン端では 0 になります。
• ねじり剛性を持たない部材（\( J \le 0 \)）の rx は、解放しても縮約行列 \( K_{bb} \) の対角がゼロになり特異化するだけでモーメント解放としての意味を持たないため、解放しません。
• 線材（梁・柱）の i 端のねじれは、端条件が剛接でも既定で解放します（ねじり剛性を期待しないモデル化。適用判定・例外は 部材のねじり剛性の既定モデル化（i 端ねじれ解放） と同一で、弾性梁と非線形梁で同じ規則です）。

非線形解析での部材内力

弾塑性要素の断面内力は、部材端の復元力を材軸方向へ釣合いで分配して求めます（線形解析の内力回収 線形静解析 と同じ分配式）。復元力は、接線剛性 × 全変位では降伏後の内力が実際より小さく評価されるため、接線剛性と全変位の積ではなく各積分点の断面応答（ファイバーの履歴状態）から算定した値を用います。

評価断面（危険断面）は線形解析と同じ規則（節点芯・部材中央・柱フェース位置、および部材付帯情報の追加位置。6.2.3）で与え、線形解析と同じ断面で応力を比較できるようにします。

適用範囲：ファイバー・MS・材端集中ばね（sepika_element::frame::concentrated）・弾性材・トラス（軸材）に適用します。