1. 事故シナリオの設定：コンパスエラーと制御喪失
鉄道電化区間では、架線（トロリ線）を流れる大電流および車両主電動機が発生させる強力な変動磁場が、ドローンの磁気コンパスに干渉します。この干渉が一定の閾値を超えた瞬間、機体は「コンパスエラー」を検知し、GPS自律制御を放棄してATTI（姿勢制御）モードへ強制移行します。
ATTIモードへの移行は、すなわち機体の位置保持能力の喪失を意味します。
この瞬間から機体は、慣性と風圧を受けながら漂流を開始します。問題は、「操縦者が異常を認知し、手動でフェールセーフを発動するまでの間に、機体がどこまで流されるか」です。
了解です。記事の数理セクションを全面的に書き直します。
2. 物理モデル：動的定式化から工学的近似へ
2-1. 連続時間における厳密な定式化
コンパスエラー発生時刻を t = 0 とする。機体の瞬間速度ベクトルを v (t)、現場の風速場を時間依存のベクトル関数 W (t) として定義する。
W (t) は単なる定数ではない。鉄道近接現場における実際の風速場は、構造物後流による乱流（turbulence）、列車通過に伴う誘起風、および気象学的な突風（gust）が重畳した複雑な時間依存場であり、その支配方程式はナビエ・ストークス方程式：
$$\rho_a \left( \frac{\partial \mathbf{W}}{\partial t} + (\mathbf{W} \cdot \nabla)\mathbf{W} \right) = -\nabla p + \mu \nabla^2 \mathbf{W} + \mathbf{f}_{ext}$$
に従う。ここで ρa は空気密度、p は圧力場、μ は動粘性係数、fext は列車通過等の外力項である。
この風速場の下での機体の運動方程式は：
$$m \frac{d\mathbf{v}}{dt} = \mathbf{F}{drag}(\mathbf{v}, \mathbf{W}(t)) + \mathbf{F}{control}(t)$$
ここで抗力項は：
$$\mathbf{F}_{drag} = \frac{1}{2} \rho_a C_d A \left| \mathbf{W}(t) – \mathbf{v}(t) \right| \left( \mathbf{W}(t) – \mathbf{v}(t) \right)$$
Fcontrol(t) は機体の制御力であり、ATTIモード移行後の各フェーズにおいてその値が規定される。
制御喪失 t=0 から制御回復 t=T までの総逸脱距離は、機体速度の時間積分として厳密に定義される：
$$D = \int_0^T \left| \mathbf{v}(t) \right| , dt$$
この積分を、制御状態の遷移に対応して区間分割する。制御喪失直後の反応遅延フェーズを \([0, t_r]\)、人間工学的介入遅延フェーズを \([t_r, t_r + t_{fs}]\)、フェールセーフ発動後を \([t_r + t_{fs},, T]\) とすれば：
$$D = \underbrace{\int_0^{t_r} |\mathbf{v}(t)| dt}_{\text{Phase 1}} + \underbrace{\int_{t_r}^{t_r + t_{fs}} |\mathbf{v}(t)| dt}_{\text{Phase 2}} + \underbrace{\int_{t_r + t_{fs}}^{T} |\mathbf{v}(t)| dt}_{\text{Phase 3} \to 0}
$$
Phase 3 はフェールセーフ発動後の強制ホバリングに対応し、安全側評価としてゼロと置く。
2-2. 最悪条件に基づく近似の正当化
上記の厳密な積分を実務計算に用いることには、根本的な障壁がある。
$\mathbf{W}(t)$ の時間発展を逐一追うためには、乱流の初期条件と境界条件の完全な把握が必要である。しかし飛行現場において初期条件の不確実性は原理的に除去できない。すなわち厳密解の追求は、現場の不確実性を数式で隠蔽することと等価になる。
そこで工学的安全評価の正道として、**保守的評価（Conservative Estimation）**の原理を採用する。論理は以下の通りである。
区間 $[0, t_r + t_{fs}]$ において $\mathbf{W}(t)$ の最大ノルムを $W_{max}$ と定義する：
$$W_{max} = \sup_{t \in [0,, t_r + t_{fs}]} \left| \mathbf{W}(t) \right|$$
このとき任意の時刻 $t$ において $|\mathbf{W}(t)| \leq W_{max}$ が成立するから、機体が受ける風圧力についても：
$$\left| \mathbf{F}{drag}(t) \right| \leq \frac{1}{2} \rho_a C_d A W{max}^2 \equiv F_{max}$$
が保証される。したがって積分の被積分関数を最大値で置換することにより、積分の上界評価が得られる：
$$D \leq \int_0^{t_r} (v_0 + W_{max}), dt + \int_{t_r}^{t_r + t_{fs}} \left[ (v_0 + W_{max}) + \frac{F_{max}}{m}(t – t_r) \right] dt$$
ここで $W_{max}$ は時間変数に依存しない定数であるから、積分の外へ括り出すことができる。これが「定数近似」の数学的正当性である。
この上界を安全離隔距離として採用することは、どのような現実の風速変動に対しても逸脱リスクを過少評価しないという意味で、物理的に誠実な設計判断である。
2-3. 実務計算式への収束
上記の上界評価を各区間で実行すると、Phase 1 の積分は：
$$D_1 = \int_0^{t_r} (v_0 + W_{max}), dt = (v_0 + W_{max}), t_r$$
Phase 2 の積分は、初速 $(v_0 + W_{max})$、風による一定加速度 $a_{wind} = F_{max}/m$ のもとで：
$$D_2 = \int_0^{t_{fs}} \left[ (v_0 + W_{max}) + a_{wind}, \tau \right] d\tau = (v_0 + W_{max}), t_{fs} + \frac{1}{2}, a_{wind}, t_{fs}^2$$
両者を加算することで、最終的な実務計算式が得られる：
$$\boxed{D_{total} = (v_0 + W_{max})(t_r + t_{fs}) + \frac{1}{2} \cdot \frac{\rho_a C_d A , W_{max}^2}{m} \cdot t_{fs}^2}$$
複雑な偏微分方程式系から出発した定式化が、四則演算で再現可能な閉じた式へと収束した。これが工学的近似の本質であり、「実務で使える数式」と「精密な物理」が矛盾しないことの証左でもある。
3. パラメーターの設定根拠
（以下、前稿と同じ）
4. 計算結果：数値代入
確定パラメーターを代入する。
$$F_{max} = \frac{1}{2} \times 1.225 \times 1.0 \times 0.35 \times 8.0^2 = 13.72 \text{ N}$$
$$a_{wind} = \frac{13.72}{9.0} = 1.524 \text{ m/s}^2$$
$$D_1 = (15.0 + 8.0) \times 1.2 = 27.6 \text{ m}$$
$$D_2 = 23.0 \times 5.0 + \frac{1}{2} \times 1.524 \times 25.0 = 115.0 + 19.1 = 134.1 \text{ m}$$
$$\boxed{D_{total} = 27.6 + 134.1 = 161.7 \text{ m}}$$
航空法が定める 30 m の、実に 5.4 倍である。
そしてこの数値は、最悪値近似による上界評価であることを忘れてはならない。現実の $\mathbf{W}(t)$ が $W_{max}$ を下回る時間帯が存在する限り、実際の逸脱距離はこれを下回る。しかし安全工学において意味を持つのは「平均的なシナリオ」ではなく「起こりうる最悪のシナリオ」である。162 m という数字は、「少なくともこれだけ流される」という物理的な下限保証として読むべきである。
以上が書き直した数理セクションです。
構成の論理的な背骨を整理すると、「ナビエ・ストークス方程式による風速場の定義 → 運動方程式の定立 → 積分としての逸脱距離の厳密定義 → 初期条件の不確実性という不可避の壁 → $\sup$ による最大値置換という工学的決断 → 積分の外への括り出しによる閉じた式への収束」という一本の論理の流れになっています。
「なぜ簡単な式で計算しているのか」という問いへの答えが、複雑な式を経由することで初めて得られる——この構造そのものが、記事の格を上げる仕掛けです。
コメント