質点系解析

本節では、各階を質点に縮約したせん断型モデルの固有値と時刻歴を扱います。 立体 FEM の固有値・時刻歴とは独立した解析です。

立体のモードは部材の曲げや局部変形を含むため、質点系の次数と並べて読むことはできません。 周期だけを比較してください。 曲げをばねとして分離する変種は、自由度の置き方が立体とさらに食い違うため使いません。

実行は右バー「質点系」です。 表とグラフは結果タブ「質点系」、球とばねは 3D ビューアの「質点モード」「質点時刻歴」です。

前提

3 次元のねじりばね（KR）は、線形・非線形のいずれでも線形ばねです。 非線形の並進だけが増分解析のトリリニアを使い、ねじりは偏心率精算の KR のままにするためです。

既定値

波形はサンプル波か波形ライブラリから選びます。 立体時刻歴にある「波形 CSV を開いて実行」は、質点系パネルにはありません。 プロジェクトファイル（.ovika）には設定・固有値・時刻歴結果を同梱し、波形の実体はライブラリ参照です。

サンプル波は立体時刻歴と同じ式です。

\[ \ddot{x}_g(t) = A \sin(\omega t)\, e^{-0.3 t},\quad \omega = 2\pi / T \]

\( A \) は振幅、\( T \) は周期、\( dt \) は刻みです。

質量

各層の質量は上端階の物理質量相当重量 \( W_{\text{mass}} \) を重力加速度で割った値です（\( m = W_{\text{mass}}/g \)）。 W_mass は準備計算の階生成が節点の物理質量相当重量から階ごとに集計した Story::dynamic_mass.mass_equiv_weight_n です。 質量方式（固有値解析の質量モデル）に依らず全量を用います（串団子は部材分布質量を持たないため）。

物理質量相当重量は物理質量（鋼材は物理質量密度 7.85 t/m³）に基づき、設計地震用重量 seismic_weight（鋼材は 78.5 kN/m³）とは一致しません（地震用重量・層データの生成の「動的解析用の質点質量の設定」）。 手入力の地震用重量（weight_override）は設計地震力のみを変え、質点系の質量は変えません。

未算定なら準備計算で階を生成する必要があり、質点系は明示エラーで止まります（設計重量での代用はしません）。

3 次元の回転慣性 \( J \) [t·mm²] は、同じ物理質量分布から質量重心まわりに直接算定した値（Story::dynamic_mass.inertia_t_mm2）です。 質点の平面位置はその質量重心 center_xy_mm です。 従来の剛床マスター RZ 質量の比例換算は使いません。

\[ J = \sum \frac{w_{\text{mass},i}}{g} , \lvert r_i - G \rvert^2 \]

\( r_i \) は節点 \( i \) の平面位置、\( G \) は質量重心です。 2 次元の 3D 表示では、剛床マスター（なければ原点）を代表点にします。

記号・単位

層並進剛性

層 Q/δ（既定、線形時）

地震静的 EX の層せん断 \( Q \) と層間変位 \( \delta \) から \( K_x = |Q|/|\delta| \)、EY から同様に \( K_y \) を算定します。 層せん断は当該層の柱せん断の合計です。 柱せん断はすでに当該層より上の水平力を含むので、層間で足し直しません。

2 次元の層間変位は、階の代表節点（剛床マスター、なければ質量最大の節点）の水平変位差です。 3 次元の層間変位は、その層の剛心の平面位置における剛床変位の上下差です。 3×3 の並進ばねは剛心位置のばねなので、\( \delta \) も剛心で取ります。 上下とも同じ平面位置で取るため、セットバックでマスターがずれても層間になります。

剛心座標とねじり剛性 KR [N·mm/rad] は、偏心率の精算と同じ柱剛性 \( k_i = Q_i/\delta_i \)（雑壁の寄与を含む）から算定します。

柱 ki（線形時）

\( K_x \)・\( K_y \)・剛心・KR をすべて柱の \( k_i \) 合計（偏心率精算と同じ）から定めます。

非線形時の並進骨格は、該当方向の増分解析の層 Q–δ を等包絡面積則でトリリニア縮約したものです（第 1 折点の割線比は上表。長期載荷点を原点とする補正は 質点系（串団子）非線形時刻歴 と同じです）。 固有値と減衰に使う初期剛性も、この骨格の K1 です。 線形時の定義（層 Q/δ または柱 ki）は使いません。

2 次元（せん断串）

加振方向 1 本の質点列です。層間ばねはせん断型三重対角です。

算定式

固有値は

\[ K x = \omega^2 M x \]

（M は対角、K は三重対角）を標準固有値問題へ帰着して解きます。 モード形状は最上階を 1.0 に正規化します。 刺激係数と有効質量は算定しません。 立体固有値の参加係数は M 正規化に紐づく量で、質点系は縮約モデルの周期を見る解析だからです。

時刻歴は Newmark 平均加速度と Newton 反復、減衰は初期剛性比例 \( a_1 = 2h/\omega_1 \) です。 詳細は 質点系（串団子）非線形時刻歴 です。

3 次元（Ux, Uy, θz）

各階の自由度は質量重心における並進 2 成分と鉛直軸まわり回転です。 剛心相対位置を \( r_x = X_s - X_g \)、\( r_y = Y_s - Y_g \) として、層の 3×3 剛性は次です。

算定式

\[ K = \begin{bmatrix} K_x & 0 & -K_x r_y \ 0 & K_y & K_y r_x \ -K_x r_y & K_y r_x & K_x r_y^2 + K_y r_x^2 + K_R \end{bmatrix} \]

\( K_x \)・\( K_y \) は線形なら静解析の層剛性、非線形なら増分骨格の初期剛性です。 階間は相対変位で直列につながります。 地動は加振方向（X または Y）の 1 方向です。

たとえば偏心のない 1 層では、並進の固有周期は \( T = 2\pi\sqrt{m/K} \)、ねじりの固有周期は \( T = 2\pi\sqrt{J/K_R} \) となり、上の 3×3 を対角化した結果と一致します。

非線形時、X・Y の並進ばねは各方向のトリリニア（最大点指向型）、ねじりは KR の線形ばねのままです。 3 次元のモード形状は各階 \( (U_x, U_y, \theta_z) \) で、頂部水平変位のノルムで正規化します。 刺激係数と有効質量は 2 次元と同じく算定しません。

時刻歴の方向別ピーク

時刻歴の層間・層せん断・塑性率は、X・Y・45°・最大の 4 欄で出します。 地動は X か Y の 1 方向でも、ねじりがあれば剛心位置の層間は平面内を動くためです。

層間は剛心位置の並進 \( (\delta_x,\ \delta_y) \) です。 層ばねの変形がここだからです。 層せん断は同じ位置のばね力 \( (Q_x,\ Q_y) \) です。 頂部変位のグラフは質量重心の加振方向成分です。 質点の自由度がそこだからです。

算定式

各時刻の値から、層ごとに時刻歴の最大絶対値を取ります。

45° は 45° と 135° の投影の大きい方です。 斜めの応答はどちらの対角にも出るためです。

塑性率は次です。 \( \delta_{1x} \)・\( \delta_{1y} \) は各方向骨格の第 1 折点変形です。

45° の分母は、軸平行の降伏矩形 \( |\delta_x|=\delta_{1x} \)、\( |\delta_y|=\delta_{1y} \) に沿う斜め距離です。 最大の塑性率は並進ばねが経験した値の包絡です。 合成変位を \( \min(\delta_{1x},\delta_{1y}) \) で割ると、ばねが感じていない量になるため使いません。

2 次元では加振方向以外の並進は 0 です。 45° は加振方向成分の \( 1/\sqrt{2} \) になります。

3D 表示

3D ビューアの「質点モード」「質点時刻歴」で、各階の質点を球、階間ばねを線で描きます。 球の半径は質量に応じます。 「骨組を重ねる」をオフにすると、質点の基準位置は鉛直 1 列です。 質点モードでは変形前の串を破線と高い透過で重ねます。 表示倍率は立体の変形図と同じく、ピーク並進がモデル対角長の 10% になる自動倍率に手動係数を掛けます。 詳細は 質点系の 3D 表示 です。