原点を中心とした復元力がはたらくとき，運動方程式は次のように表せる。

ただし， は物体の質量， は復元力の定数（ >0 ）とする。

この微分方程式を解くと，変位 を時間 の関数として表すことができる。ただし， とする。

本研究では，(1)式をもとにしたシミュレーションを行うことで，微分方程式を解かずに，(2)式のような変位と時間の関係が得られることを示したい。ただし，本内容は「表計算ソフトを利用したシミュレーション１〜モンキーハンティング〜」を学習していることを前提とする。

差分法を用いて，復元力を受けた物体の運動をシミュレーションする。

時刻 での加速度を ，変位を とすると，(1)式は以下のように変形できる。

時刻 での速度 が，微小時間 の間に一定の加速度 で変化すると近似すると， 後の速度 は

と表せる。

同様に，時刻 での変位 が， の間に一定の速度 で変位すると近似すると， 後の変位 は

と表せる。

ここで，経過時間を微小時間 ごとに区切り， 番目の時刻 に対応する物理量を下の表１のように表す。

	時刻	加速度	速度	変位
第項

表１　時刻 に対応する物理量

このとき，(3)〜(5)式より， 番目の各項は次のような漸化式で示せる。

(6)〜(9)式をもとに，表計算ソフトのセルに数式を入力する。表２は入力の一例である。ただし，２行目と５行目の淡黄色のセルは初期設定値の入力欄で，セルの数値を変更することで運動の条件を変更できる。

	A	B	C	D
1	時間間隔 （s）	正の定数	質量 （kg）
2	0.0005	1.0000	0.1000
3
4	時間 （s）	加速度 （m/s2）	速度 （m/s）	変位 （m）
5	0	= - $B$2 / $C$2 * D5	0	1.0000
6	= A5 + $A$2	= - $B$2 / $C$2 * D6	= C5 + B5 * $A$2	= D5 + C5 * $A$2
7	= A6 + $A$2	= - $B$2 / $C$2 * D7	= C6 + B6 * $A$2	= D6 + C6 * $A$2
8	= A7 + $A$2	= - $B$2 / $C$2 * D8	= C7 + B7 * $A$2	= D7 + C7 * $A$2

表２　表計算ソフトへの入力例

表２のように入力して，６行目を下方のセルにコピーしていくと，下の表３のような結果を得ることができる。

	A	B	C	D
1	時間間隔 （s）	正の定数	質量 （kg）
2	0.0005	1.0000	0.1000
3
4	時間 （s）	加速度 （m/s2）	速度 （m/s）	位置 （m）
5	0.0000	-10.00000000	0	1.0000
6	0.0005	-10.00000000	-0.005000000	1.000000000
7	0.0010	-9.999975000	-0.010000000	0.999997500
8	0.0015	-9.999925000	-0.014999988	0.999992500
・ ・ ・	・ ・ ・	・ ・ ・	・ ・ ・	・ ・ ・
3362	0.9345	9.849675589	-0.587557105	-0.984967559
3363	0.9350	9.852613375	-0.582632267	-0.985261337
3364	0.9355	9.855526536	-0.577705960	-0.985552654
・ ・ ・	・ ・ ・	・ ・ ・	・ ・ ・	・ ・ ・

表３　表計算結果

表３の結果を用いて，時間と変位，速度，加速度の関係をグラフに示すと次のようになる。

図１　復元力（変位と時間）

図２　復元力（速度と時間）　　　　　　　　　　　　　　　　

図３　復元力（加速度と時間）

表３や図１から復元力がはたらくとき物体は単振動することが分かる。また，各図から変位，速度，加速度の変化の周期が等しく，そのとき，変位と速度はπ/2，変位と加速度はπだけ位相がずれることも分かる。理論式(2)から得られる上記入力例での振動の周期は1.99秒，振幅1.0m，最大速度3.16m/s，最大加速度10m/s²であるので，表計算の結果は現象をよく再現している。

物体に復元力がはたらくとき，振動の中心付近では速度が最大になる一方で加速度の変化が極微小である。差分法で計算する場合，中心付近で各物理量に若干の誤差が発生する。計算時間間隔を小さくすることで精度を高めることができる。

本教材は、高等学校学習指導要領の「第２章　各学科に共通する各教科」における

・第４節　数学　　第６　数学活用　　(2) 社会生活における数理的な考察　　ウ　データの分析

・第５節　理科　　第３　物理　　(1) 様々な運動　　ウ　円運動と単振動

・第５節　理科　　第10　理科課題研究

・第10節　情報　　第２　情報の科学　　(2) 問題解決とコンピュータの活用　　ウ　モデル化とシミュレーション

という項目で利用できる。

本研究で紹介した表計算のサンプルファイルと，本内容を応用して作成した弾性力のシミュレーションソフト（下図）も用意した。下記リンク「simulation04.zip(88.0KB)」より取得できる。なお，回転運動を含んだ減衰振動までシミュレーションできるソフトになっているので，単振動をさせるには初期設定を以下の例のように変更する必要がある。
例　初期座標 x成分：2.5→0 y成分：0→-1　初速度 y成分：-5→0　バネ定数：50→9.8　自然長：1.5→1　抵抗係数：0.1→0

simulation04.zip(88.0KB)