45 研究ノート 空気抵抗を受ける鈍い物体の落下運動に関する解析と実験 * 防衛大学校名誉教授五十嵐保 神奈川工科大学自動車システム開発工学科 石綿良三 空気抵抗を受ける物体の落下運動は, 簡単に実験ができ中学生の理科教育のデモンストレーション用の教材として適している. また, 高校生や大学生の物理 数学の演習問題および基礎実験の教材として使える. そこで, 名刺サイズのカード, コースターやストロー, 円柱, 球などの鈍い物体が姿勢を変えずにほぼ鉛直に落下する場合に対し, 物体指数を導入し落下運動の解析を行った. 次に, 各物体を水平な状態で m の高さを落下させ, それぞれの落下時間を測定し比較した. 得られた各実測値は, それぞれの解析値と % 以内で一致した. Analysis and Experiments on the Falling Motion of Bluff Bodies Affected by Air Resistance Tamotsu IGARASHI, Professor Emeritus, National Defense Academy Ryozo ISHIWATA, Kanagawa Institute of Technology (Received, September, 5; in revised form 8, December, 5) A falling motion of bluff bodies affected by air resistance is performed a simple experiment, which is one of a demonstration teaching material for science room of junior high school students. This phenomenon is fitted to the practical exercise of physics and mathematics for high school and foundation course of university students. The theoretical analysis on the falling motion of various cards, a coaster, a straw, a cylinder, a sphere which are plumb down in a same position were carried out introducing the body index. Next, these bluff bodies are dropped horizontally at the height of m above the floor. The falling times of various bluff bodies were measured and compared with theses data. The experimental values obtained were agreed within % with the theoretical results. (KEY WORDS): falling motion, bluff bodies, air resistance, drag coefficient, gravity force, acceleration force, body index, second order linear differential equation まえがき 高等学校や大学の教養の物理では質点の落下運動が取り上げられているが, 多くは重力以外の力は無視する自由落下である. 空気抵抗を考慮した解析はわずかで, 速度に比例した空気抵抗, あるいは速度の 乗に比例した空気抵抗が働く場合を扱っている ) ~ 5). 実際の現象としては, 小さな雨滴の落下で, 終端速度を求めるに留まっている. 一方, 理論解析と共に空気中を落下する球状物体の落下速度を測定する物理学実験 6) や自由落下法による重力加速度の測定の実験報告 7), 8) もある. * 36-46 横浜市金沢区釜利谷西 3-56-5 E-mail: tamoigarashi@yahoo.co.jp 五十嵐 9)~4) は高校生 大学生の理科離れ対策の一つとして, 身近な物理現象を取り上げ, 高等学校の数学 物理を活用した解析を行い, 簡単な実験を行い報告してきた. 石綿らの 流れのふしぎ 5) は身のまわりにあるもので, 遊びを通して流れに親しみ, そのふしぎさを楽しく学ぶ一般書である. 本報では空気抵抗を受ける鈍い物体の落下運動を取り上げた. 物体の落下の運動方程式から, 微分方程式を解いて, 落下速度および落下距離の一般解を求めた. 次に, 鈍い物体の一例として身近にある名刺サイズのカード, コースター, ストロー, 円柱, ピンポン, 球等を m の高さを落下させて, その着地時間を測定した. この時間差から空気抵抗の大きさを簡単に体験できる. 得られた各物体の実測値はそれぞれの理論解と % 以内で一致した. 6 日本流体力学会
46 記号 a, b : 矩形の縦と横の長さ B : 物体指数 = (/)C D (S/m) C D : 物体の抗力係数 D : 物体に働く抗力 d : 円板, 円柱および球の直径, 代表長さ g : 重力加速度 H : 床からの高さ L : 円柱の長さ m : 物体の質量 S : 物体の前方投影面積 Re : レイノルズ数 = d/ T* : 終端速度に達する時間 T : m 落下するのに要する時間 : 物体の落下速度 : 物体の終端速度 x : 鉛直方向下向きの座標 X : 落下距離 : 空気の密度 : 空気の動粘性係数 理論解析 名刺サイズのカードを水平にして, 床からの高さ H の位置から落下させる. カードはほぼ真下に落下する. そのほか, コ - スター, 円柱や球などの任意形状の鈍い物体が空気中を落下するときの落下距離 X と落下速度 を理論解析する. 物体に働く力は図 に示すように, 重力 mg, 浮力 F およびニュートンの粘性法則による流体抵抗 D である. 流体抵抗 D は, 物体の抗力係数 C D に空気の動圧.5 と物体の前方投影面積 S の積である. ただし, 浮力 F は無視する. D C.5 S () D は空気の密度で. kg/m 3 ( ) である. 抗力係数 C D は物体形状により異なるが, Re 数, d/ が十分大きい場合は, ほぼ一定である. D = C D.5 S = dx/dt D mg 図 空気中を落下する物体に働く力 x S. 落下の運動方程式 鉛直下向きに x 軸をとり, 運動方程式を立てる. d x m mg D dt () 上式の D に式 () を代入し, 両辺を m で割る. d x CD g S (3) dt m C D が一定の場合, 微分方程式 (3) は解析的に解くことができる. 式を簡単化するため B = (/) C D (S/m) (4) g/b =, g / B (5) とする. ここで, B は物体の質量 m, 物体の前方投影面積 S および形状による抗力係数 C D を含んだ関数で, 物体指数と定義する. なお, 式 (3) で d x/dt = のとき, 右辺は g B =, すなわち, = となる. 式 (5) の は空気抵抗を受ける物体の落下における終端速度である. 落下速度 は dx/dt で, 式 (3) は次の 階の微分方程式となる. d dt B g B( ) (6) 初期条件 :t = で =, x= (7) 微分方程式 (6) は変数分離して, 部分分数分解法を使って解析的に解くことができる. d B dt (8) ( ) t d Bdt これを積分して, 次式が得られる. tanh Bt C 初期条件 (7) より, C = である. 次式が得られる. tanh ( ) Bt 結局, 落下速度 は次式で与えられる. e e Bt Bt tanh( Bt) () () (9) ここで, 式 (5) を用い を消去する. 落下速度 は次式で与えられる.
47 / g / B = tanh Bg t () 次に, 時間 t における落下距離 X を求める t X dt X = ( g / B t tanh Bgt dt g / B / Bg ) ln [ cosh Bg t ] (3) 従って, X は次の無次元式で与えられる. BX = ln [ cosh ( Bg t) ] (4) なお, cosh x, tanh x は双曲線関数であり, 次式で定義される. x x x x sinh x ( e e ) /, cosh x ( e e ) /, x x sinh x e e tanh x cosh x e x e x (5) 空気抵抗を受ける物体の落下運動は, 物体指数 B と重力加速度 g, および時間 t で記述される.. 解析結果のまとめ (m/s) X (cm) B cm.3..6.3.. t (s) 図 物体指数による落下速度の変化 5 4 3 B cm.3..6.3... 終端速度に達する時間 T * と落下距離 X * 落下速度 が終端速度 の.995 に達する時間 T* は式 () より求まり, 次式 (6) で与えられる. tanh Bg T* =.995, Bg T*= 3,.5.5 t (s) 図 3 鈍い物体の落下距離 T* = 3/ Bg (6) また, この間の落下距離 X* は式 (4) から求まる. 結局, 落下距離 X* は次式 (7) で与えられる. X* = ln (cosh 3) /B =.39/B (7).. 物体指数による落下速度 と落下距離 X 次に, 式 () を用いて, 物体指数 B に対する落下速度 (t) の変化を図 に示す.B の値が小さいほど, 落下速度 は大きく, 終端速度 に達する時間は長くなる. 式 (4) を用いて, 物体指数 B に対し任意時間 tの落下距離 X(t) 示したのが図 3 である. B =.~.3 cm の範囲で, X = cm の落下に要する時間 T は.65~.43 s である...3 落下時間 t と落下速度 の求め方落下時間 T は次の手順で求める. まず, 式 (4) より X = cm のときの B T の関係を求め, グラフにしておく. 次に, 物体指数 B の値を求め, グラフから落下時間 T を予測する. これが, 式 (4) を満たすまで繰り返し計算して T が決定する. さらに, 式 () に代入して, 落下速度 (T ) が求まる. 3 実験 3. 供試物体と落下実験実験に用いた供試物体は身近にあるものを利用した. その形状は矩形板, コースター, 円柱, 球で 6) ある. 以上の鈍い物体の抗力係数 C D の値を表 に示す. 本実験レイノルズ数範囲は Re = d/ = 3 ~ 4 で, 抗力係数は一定と見なす. 表 鈍い物体の抗力係数 形状 矩形 円板 円柱 球 CD.5...4 表 に供試物体の形状, サイズ, 前方投影面積 S および質量 m を記した. 矩形 ~5 は Suica, キャッシュカードである. 円板と正方形はコースターである. 円柱 はストロー, 円柱, 3 は硬質発泡スチロール製とトイレットペーパーの芯である. 球 はピンポン玉, 球, 3, 4 は発泡スチロール球である. 実験は物体を床から高さ c m の位置から落下させた. 物体の m の落下運動はカメラを長時間露光にして, ストロボスコープで 3 Hz で発光して
48 表 供試物体 形状 サイズ cm a b, d l, d 投影面積 S cm 質量 m g 矩形 5.4 8.57 46.3 5.3 5.38 8.56 46. 5. 3 5.4 8.54 46. 4.8 4 5.4 8.55 46. 5. 5 5.4 8.6 46.4. 円板 9. 63.6 3.7 9. 63.6.85 3 9. 63.6 8. 正方形 9. 8. 3. 9.3 86.5 4.9 円柱.6 6.5 9.9.7 3. 6. 6. 3 3.8. 38.8 5.4 球 d = 4..6.7 d = 4.87 8.6.5 3 d = 7.36 4.8 4.6 4 d =. 79. 8.49 連続重ね撮りした画像から落下時間を読み取った. この画数 N は整数ではなく, 小数点以下 桁まで読み取り, 落下時間は T = N/3 (s) より求めた. 3. 落下の様子鈍い物体の落下の様子を図 4 (a), (b) に示す. なお, 物体は揺動の少ない重いものを選んだ. カード 3 円板 円柱 図 4 (a) カード, 円板, 円柱の落下の様子 d = 4. cm, d = 7.36 cm, d =. cm 図 4 (b) 球,, 3 の落下の様子 3.3 実験結果と解析との比較物体がほぼ真下に落下した場合の ~3 回のデータの平均値を着地時間とした. 測定誤差は % 以下で, 実測値 T exp を解析で得られた第 章の式 (4) の所要時間 T th および物体指数 B の値を表 3 に示す. 表 3 鈍い物体の着地時間の実測値 形状 B cm T th s T exp s 矩形.6.77.767.635.779.783 3.663.785.777 4.63.774.767 5.6.983.983 円板.6.89.96.3..5 3.534.756.74 正方形.75..3 正方形..95.9 円柱..863.85.75.699.69 3.56.739.733 球..663.66.389.74.76 3.45.69.69 4..687.68
49 図 5 に実測値と解析解との比較を物体指数 B に対して示す. 各物体の実測値は式 (4) より求めた解析解 T th との誤差は % 以下である. T (s).5.5 図 5 鈍い物体の着地時間 : 解析解と実測値の比較 物体指数が B.67 cm では着地時間 T は.8 s 以下で, B. cm では.85~.5 s に増加している. 4 まとめ 空気抵抗を受ける鈍い物体の落下運動に関し, 理論解析を行った. 落下速度 および落下距離 X は, 物体指数 B = (/) CD S/m (cm - ), 重力加速度 g, 時間 t によって, 次式で記述される. = Tth 矩形円形正方形円柱球... B (cm - ) g / B ( tanh Bg t ) X = (/B) ln [ cosh ( Bg t )] ここで, CD, S, m は物体の抗力係数, 前方投影面積, 質量, は空気の密度である. m の高さからの物体の落下所要時間の実測値は解析解と誤差 ±% 以内で一致する. 物体指数 B は, 落下時間 t と落下距離 X を一義的に表している. また, 落下速度が終端速度 の.99 に達する時間 T* および落下距離 X* は, それぞれ式 (5), 式 (5) で与えられる. [ 付録 ] パチンコ玉の落下パチンコ玉 (d =. cm, m = 5.4 g) の物体指数は B = 4. 5 cm である. これを高さ X = m から落下させと時の着地時間は T =.64 s である. 一方, 空気抵抗のない自由落下の場合, h =.5gt より着地時間は.639 s で, 両者の差異は.4 % である. この時間測定から 得られる重力加速度は g = 97 cm/s である. したがって, パチンコ玉の落下実験は, おおよその重力加速度 g を求める実験に使える. 謝辞本研究の実験に協力いただいた神奈川工科大学学部生の清水直道, 藤本航, 八田翼 3 君 ( 石綿研究室 ) に謝意を表します. 引用文献 ) 小口高, [ 解明 ] 新物理, 文英堂, (987), 47-49. ) 伊東敏雄, なーるほど! の力学 ( 第 版, 第 5 刷 ) 学術図書出版, (999), 5-. 3) 原島鮮, 質点の力 [ 復刊 ] 基礎物理学選書 (), 裳華房, (), 84-86. 4) 岡真, 質点系の力学, ニュートンの法則から固体の回転まで, 共立出版, (3), 9-5. 5) 小西克亨, 終端速度の求め方, 埼玉工業大学, 機械工学学習支援セミナー用教材 J3 (3). 6) 東海大学物理学実験連絡協議会編, 中村誠太郎監修, 物理学実験, 東海大学出版会, 改訂 3 版 7 刷, (6),83-9. 7) 石川俊明, 日本物理教育学会誌, 9-4 (98), 97-3. 8) 天野卓治 内田直, 日本物理教育学会誌, 3-4 (983), 8-. 9) 五十嵐保, 静水に浮く角材の姿勢, 日本流体力会誌, ながれ 9-4 (), 53-6. ) 五十嵐保 中村元, 飲料水の入った缶が斜めに立つ安定解析, 日本流体力学会誌, ながれ 4-(5), 95-3. ) 五十嵐保 中村元, 静水に浮く矩形断面柱の姿勢解析, 日本流体力学会誌, ながれ 6-6 (7), 393-4. ) 五十嵐保, 静水に浮く有限長円柱の姿勢解析, 日本流体力学会誌, ながれ 8-5 (9), 49-4. 3) 五十嵐保, 水に浮く比重 以上の円盤, 日本流体力学会誌, ながれ 8- (9), 49-57. 4) 五十嵐保, 容器の縁から盛り上がる水面の高さに関する実験と解析, 日本流体力学会誌, ながれ 33- (4), 6-65. 5) 日本機械学会編, 石綿良三, 根本光正著, BLE BACKS 流れのふしぎ, 講談社 (4), 54-55. 6) 日本機械学会編, 機械工学便覧基礎編 α4, 流体工学 (6), 83-85.