このプロンプトは、机の端から突き出した定規を弾いた際に生じる「振動」と、それが空気抵抗や内部摩擦によって停止するまでの「減衰」を、物理エンジンや数値計算シミュレーターで再現するための定義式および計算パラメータのセットである。 ### 1. 物理モデルの定義 定規を「片持ち梁（Cantilever Beam）」としてモデル化する。以下のパラメータを順次定義し、計算エンジンに読み込ませること。 1. **幾何学的パラメータ** * L (長さ): 定規の全長 [m] * l (突き出し長さ): 机の端から出ている部分 [m] * w (幅): [m] * t (厚さ): [m] * I (断面二次モーメント): 矩形断面の場合 $I = \frac{w \cdot t^3}{12}$ [m^4] 2. **材料特性パラメータ** * E (ヤング率): 材質に応じた定数 [Pa]（プラスチックなら約2-4GPa、金属なら約70-200GPa） * ρ (密度): [kg/m^3] * m (質量): $m = \rho \cdot l \cdot w \cdot t$ [kg] 3. **動力学パラメータ** * k (バネ定数): 片持ち梁の理論値 $k = \frac{3EI}{l^3}$ [N/m] * ωn (固有角振動数): $\omega_n = \sqrt{\frac{k}{m_{eff}}}$ （ここで $m_{eff} \approx 0.24m$ を使用） ### 2. 減衰モデルの設定（Rayleigh減衰の適用） 定規の揺れが止まる現象は、主に「空気抵抗（粘性減衰）」と「内部摩擦（ヒステリシス減衰）」の合成である。計算には以下の減衰定数 ζ（ゼータ）を導入する。 * **減衰比 (ζ) の設定目安** * ζ = 0.01 ～ 0.02: 非常に硬い金属定規（減衰が少ない） * ζ = 0.05 ～ 0.08: 一般的なプラスチック定規 * ζ = 0.10 以上: 柔らかい樹脂やゴム製定規 * **減衰振動の微分方程式** $m \cdot \frac{d^2x}{dt^2} + c \cdot \frac{dx}{dt} + k \cdot x = 0$ ここで、減衰係数 $c$ は $c = 2\zeta\sqrt{mk}$ と定義する。 ### 3. 数値計算アルゴリズム（Python/疑似コード） 以下のコードスニペットは、Runge-Kutta法を用いて定規の先端の変位 $x(t)$ を算出する。 ```python import numpy as np def simulate_ruler_vibration(t_max, dt, zeta, omega_n, x0): """ t_max: シミュレーション時間 zeta: 減衰比 omega_n: 固有角振動数 x0: 初期の弾き幅 """ t = np.arange(0, t_max, dt) # 減衰振動の公式: x(t) = x0 * exp(-zeta * omega_n * t) * cos(omega_d * t) omega_d = omega_n * np.sqrt(1 - zeta**2) x = x0 * np.exp(-zeta * omega_n * t) * np.cos(omega_d * t) return t, x # 例: プラスチック定規 params = { 'zeta': 0.06, 'omega_n': 2 * np.pi * 5.0, # 5Hzで振動と仮定 'x0': 0.05 # 5cm弾いた } ``` ### 4. 実行時の入力用テンプレート シミュレーターに対して以下のプロンプトを入力し、出力を調整せよ。 --- **[入力プロンプト]** 以下の条件で片持ち梁の自由減衰振動をプロットせよ。 - 突き出し長さ: [入力値: 例 0.15m] - 材料ヤング率: [入力値: 例 3.0e9 Pa] - 減衰比 (ζ): [入力値: 例 0.05] - 初期変位: [入力値: 例 0.03m] - 計算時間: 2.0秒 - 出力要求: 時間(t)に対する変位(x)のグラフを描画し、振幅が初期値の1/eになる時間（時定数）を算出すること。 --- ### 5. 物理的整合性を保つためのチェックリスト 計算結果が現実の感覚と乖離している場合は、以下の変数を確認せよ。 1. **質量の分布**: 梁の質量を先端に集中させすぎると周期が短くなる。分布質量モデル（オイラー・ベルヌーイ梁）として計算できているか。 2. **非線形性の無視**: 定規を過度に大きく（定規の長さの20%以上）弾いた場合、線形バネ定数 $k$ の仮定が崩れる。その場合は、変位に応じた復元力の非線形項を追加せよ。 3. **空気抵抗の非線形**: 高速で振動する場合、減衰力は速度の2乗に比例する項（ドラッグ）が支配的になる。その際は $c \cdot \frac{dx}{dt}$ ではなく $c_d \cdot (\frac{dx}{dt})^2$ を導入すること。 このプロンプトは、力学的な構造解析の基礎となる。まずはこのモデルで基本挙動を再現し、必要に応じて材料の異方性や形状因子を加えて精緻化すること。物理の計算は、変数の定義さえ正確であれば、必ず現象を正確に写し出す。