運動方程式
| 分類 | 力学・計算物理の基本公式群 |
|---|---|
| 主な対象 | 粒子・剛体・連続体 |
| 利用領域 | ロケット制御、構造解析、流体モデル |
| 成立の経路 | 星図作成→機械設計の要請→数理定式化 |
| 典型的な形 | 状態変数の時間微分を含む関係式 |
| 関連分野 | 変分原理、制御理論、数値計算 |
| 重要文献 | 『運動方程式綜覧』(架空)ほか |
| 備考 | 歴史的には「式」より運用手順が先に共有された |
（うんどうほうていしき）は、力や状態の変化を「いつ・どれだけ」起こすかを定める上の式として知られている[1]。本来はの星図整備から始まったとされ、のちにやへ波及していった[2]。
概要[編集]
は、対象の運動がどのような入力（力・拘束・エネルギー条件）に応じて変化するかを記述する関係式であるとされる。たとえば、速度や加速度といった時間変化の量を含む表現が一般的であり、数値計算では微分方程式として扱われることが多い。
一方で、運動方程式が「式そのもの」より「計算手順の流通」を中心に発展した点が特徴として挙げられる。実務者の間では、式の見た目よりも「係数の丸め規則」「入力データの検算手順」「収束の合図」といった運用規格が先に定着し、それが後世に“標準形”として回収された経緯があるとされる[3]。
また、近代以降はやと強く結びついた。特にの現場では、運動方程式が単なるモデルではなく「計算機に渡す言語」に近い扱いを受けていたとされる。結果として、同じ系を扱っても、どの“手順”を運動方程式として採用したかで成果が変わるという事情が生じた。
歴史[編集]
星図係のノートから始まった「時間の商い」[編集]
運動方程式の起源は、における星図整備の失敗記録に求められるとされる。16世紀末、の交易拠点で観測された恒星位置の記録が、季節によって系統的にずれるという事件が起きた。そこで当時の測時係は、観測者の癖ではなく「時間の商（わり）」を疑い、“時刻のずれ”が原因であると結論したとされる[4]。
この“時刻のずれ”を見積もるために考案されたのが、のちの運動方程式につながる「差分→補正→再差分」という反復手続きであると説明される。実際、当時の記録簿には、恒星の位置誤差が平均で0.0087度から0.0079度へ改善した日、悪化した日、さらに悪化したが翌日に回復した日が、磁針の揺れと並んで丁寧に列挙されていたとされる[5]。
この工程を束ねる“係数の置き方”が、運動方程式の原型として伝承されたとされる。その後、手順を簡潔に書き下すための記号体系が必要になり、17世紀前半に（当時の官署とされる）で「運動を表す時間微分」という言い回しが試験的に採用された。ここで微分は“定義”より“運用”として整備されたため、以後の研究が「式の厳密さ」より「計算の成功率」に傾きやすくなったとの指摘がある。
工廠の係数会議：誰が関わり、何が問題になったか[編集]
運動方程式がの中心概念になった転機は、19世紀後半の工廠の係数会議にあるとされる。具体的には、の小型蒸気機械工場で、同じ部品図にもかかわらず振動が周期的に増幅する“再現性のない不良”が続いた。品質管理担当であったは、振動を「部品が勝手に覚えている運動」と見なし、工場内で共有された仮説係数を毎週更新する運動方程式“運用版”を提案した[6]。
会議には、材料技師の、帳簿係の（英国商社出身として説明される）、さらに“計算係”としてから派遣されたが参加したとされる。彼らは係数を決めるために、蒸気圧の読みを1/10刻みで記録し、1週間のうちに「補正が3回以上行われた系」を不良系として除外するという基準を作ったとされる[7]。結果として、歩留まりは一時的に81.3%まで回復し、しかし2か月後に再び78.9%へ落ちた。
この落ち込みの原因は“物理の誤り”ではなく“手順の誤配布”であったとされる。すなわち、運動方程式の係数表が、ある部署では小数点をひと桁読み違えて転記されていたのである。後年の回想では、係数表の印刷が急ぎすぎたせいで、校正用の試験機がたまたま正常運転してしまったと語られている。ただしこの逸話は、社内文書の閲覧記録に基づくとされる一方、当時の帳簿の散逸により真偽が揺れているとの指摘もある[8]。
内容[編集]
運動方程式は一般に、状態変数の時間変化を結びつける形で書かれると説明される。たとえば、位置、速度、加速度に対応する量を用い、入力（力、制約条件）を右辺に置くことで、次の状態を推定できるとされる。
実務上は、式の選択よりも「どの系を“閉じる”か」が重要になる。閉じ方とは、摩擦や減衰、外力の推定がどこまでモデルに含まれるかという意味である。ある計算技師は、摩擦を入れない運動方程式の方が“当たる”ことがあると述べたとされる。つまり、摩擦係数の推定誤差が他の誤差を打ち消していた可能性がある、というのである[9]。
また、運動方程式は「微分方程式としての正しさ」と「数値解としての振る舞い」が一致しない場合があり得るとされる。実際にの委託解析では、同じ初期条件に対して2つの解法が、振幅の符号が反転するまでの時刻が「37.12秒」と「37.17秒」で分岐したという報告がある[10]。この差は一見些細であるが、現場では“符号反転”を危険判定に用いていたため、運用は大きく左右されたとされる。
批判と論争[編集]
運動方程式をめぐっては、モデル化の恣意性が問題視されることが多い。とくに「標準形」と称して流通した式は、実験現場で生き残った“手順込み”の集約であり、数学的には複数の同値性を取りこぼしている可能性があると指摘されている。
さらに、運動方程式が“計算機の都合”に引きずられているという批判もある。たとえば、ある研究会では、浮動小数点演算の丸め規則に適合するよう係数が選ばれ、その結果として物理量の解釈が「意味よりも実装」に寄っていったと報告された[11]。この論点に対し、擁護側は「工学では観測可能性が最優先である」と反論したとされる。
一方で最も笑われる逸話として、運動方程式の教育現場で“係数暗記”が流行し、学生が「運動は暗記すれば発生する」と冗談めかして言ったという事例が挙げられる。ただし、当該発言は講義録には残っていないため、都市伝説扱いの部分もあるとされる[12]。それでも、現場の実感として「解ける人ほど式を短く言える」傾向があったことは、複数の回想で一致している。
脚注[編集]
関連項目[編集]
脚注
- ^ 澄川正利「星図誤差と“時間の商い”—運動方程式の予備手続きについて」『測算研究紀要』第12巻第3号, 1889年, pp. 41-63.
- ^ 鶴見修司郎「係数会議の記録：工廠における運用版運動方程式」『機械振動年報』Vol.7 No.1, 1902年, pp. 12-27.
- ^ Margaret A. Thornton「The Practical Equation of Motion as a Distribution Protocol」『Journal of Applied Computation』Vol.18 No.4, 1937年, pp. 201-219.
- ^ 田中紗月「減衰の不在が当たる理由：モデル閉鎖の観測論」『工学物理学報』第5巻第2号, 1921年, pp. 88-104.
- ^ オズワルド・ブレイク「On Coefficient Rounding and Manufacturing Histories」『Proceedings of the Cambridge Mechanics Society』Vol.33, 1911年, pp. 55-79.
- ^ 『運動方程式綜覧』編纂委員会編、幻影書房, 1974年, pp. 1-980.
- ^ 松平花子「浮動小数点時代の運動方程式：丸め規則と符号反転」『数値計算通信』第29巻第6号, 1986年, pp. 301-322.
- ^ R. K. Nasser「Equivalent Forms, Divergent Numerics」『International Review of Mechanical Modeling』Vol.52 No.2, 2003年, pp. 77-96.
- ^ 澄川正利、澁谷晶「解けることと意味することのねじれ」『物理教育学研究』第41巻第1号, 1999年, pp. 9-25.
- ^ 岡部達也「運動方程式は誰の言語か：部署別手順の系譜」『計算科学史叢書』第3巻第2号, 2012年, pp. 145-162.
外部リンク
- 運動方程式アーカイブ
- 暦算院 星図保管庫
- 係数会議デジタル展示室
- 数値丸め規則 研究ネット
- 機械振動年報 データベース