Microsoft Word - 付録A,Bとその図

Transcription

1 付録 A 1 自由度系 ( 自由振動 ) の解法 はじめに振動現象を解明するのに基本となる 1 自由度不減衰系 ( 自由振動 ) の運動方程式の作成方法とその微分 ( あるいは偏微分 ) 方程式の解法を説明する. 1 自由度系モデルには, 単振動のばね 質量モデルと数学振子を用いる. A.1 運動方程式 ( 微分方程式 ) を立てる A.1.1 ばね 質量の場合 ( 1) 単振動の運動から運動方程式を求める k x 図 A.1 単振動 ( ばね 質量 )( x 軸 ) m 単振動を想定して, 次の式 (A.1) とおく. sin ( A.1) したがって cos ( A.) sin ( A.3) いま, ニュートンの運動方程式により次式が求められる. ( A.4) 式 (A.4) に式 (A.3) を代入すると ( A.5) よって - 1 -

2 ( A.6) を求められる. 解 ( 式 (A.1)) を想定して始めたので, 運動方程式 ( 微分方程式 ) の作成 ( 式 ( A.4)) と同時に解 ( 式 (A.6)) が得られる. ( ) レイリーの方法 ( 最大の運動エネルギー = 最大の位置エネルギー ) 運動のエネルギーと位置のエネルギーをそれぞれ T, U とおくと ( A.7) sin ( A.8) ( A.9) とおけば ( A.10) ( A.11) 式 (A.10)= 式 (A.11) より ( A.1) の不減衰固有角振動数が式 (A.6) と同様に, 式 (A.1) が求められる. この場 合には 力学的エネルギー保存の法則 を用いている. 運動方程式は, 次に述べ るラグランジュの方程式から求める. ( 3) ラグランジュの運動方程式から求める 後述の式 (A.3) に詳述するラグランジュの運動方程式を次式に示す. 0 (A.13) 前述の式 (A.7) と式 (A.8) から, 次の式 (A.14) と式 (A.15) を得る., ( A.14) - -

3 ( A.15) 式 (A.14) と式 (A.15) を式 (A.13) へ代入すると 0 ( A.16) このように, 式 (A.4) と同じく微分方程式, 式 (A.16) が求められる. ラグランジュの運動方程式の座標は直線座標でも角座標 ( 一般座標 ) でもよい. A.1. 単振子 ( 数学振子 ) の場合 ( 1) 一般座標の円弧の接線方向のつり合いから求める 0 糸の長さ l 図 A. 単振子 ( 一般座標 ) 円弧 s の接線方向において, ニュートンの運動方程式から, 次式が求められる. sin ならば しかるに sin 0 ( 0 ( A.17) A.18) ( A.19) 結局, から式 (A.4) と同じく微分方程式, 式 (A.19) を得る. ( ) ラグランジュの方程式から求める 運動のエネルギー (A.0) 位置のエネルギー cos (A.1) ここでラグランジュの運動方程式を用いる. 0 (A.) - 3 -

4 いま 0 ( A.3) であるから, ( A.4) sin ( A.5) 式 (A.4) と式 (A.5) を式 (A.3) へ代入すると sin0 (A.6) sin の場合 0 ( A.7) を得る. しかるに, 式 (A.4), 式 ( A.19) と同じく微分方程式, 式 (A.7) が求められる. なお, ラグランジュの式は非保存的力学系にも次の式で示される. (A.8) D は散逸関数. ( 3) 単振子の運動をデカルト座標 (, の直交座標 ) から求める 0 糸の張力 F t, 図 A.3 単振子 ( デカルト座標 ), 座標からニュートンの運動方程式より, 式 (A.9), 式 ( A.30) が成立する. sin cos (A.9) ( A.30) - 4 -

5 ただし, 糸の張力 F t とおく. 式 (A.9) cosθ, 式 (A.30) sinθ により cos sin cos sin sin sin cos ( A.31) (A.3) 式 (A.31)- 式 (A.3) より cos sin sin ( A.33) ~~~~~~~~~~~~~~~ 参考 : 公式, なるとき, また, ~~~~~~~~~~~~~~~~ いま,sin であるから, 時間 t で 回微分すると ( 上記公式参照 ) cos ( A.34) cos cos sin ( A.35) また,cos であるから, 同様に sin ( A.36) sin 式 (A.35) と式 (A.37) を書き直すと sin cos ( A.37) sin cos cos sin ( A.38) ( A.39) 式 (A.38) と式 (A.39) を式 (A.33) へ代入 cos sin cos sin cos sin sin cos sin sin ( A.40) ( A.41) - 5 -

6 ここで,sin ならば sin sin 0 (A.43) 0 ( A.4) ( A.44) 結局, 式 (A.4) と同様に微分方程式, 式 (A.44) を得ることができる. [ 参考 1] 糸の張力を求めるここで参考のため, 図 A.3 の糸の張力を求める. 糸の張力は F t である. 1 法線方向 ( 遠心力とのつり合い ) cos (A.45) いま, より cos ( A.46) 接線方向 sin (A.47) より, 左辺に, 右辺に を乗じる sin ( A.48) よって両辺に dt を乗じて両辺それぞれ積分すると sin (A.49) ここで 0 のとき とする cos ( C は積分定数 ) ( A.50) (A.51) したがって 1 cos ( A.5) 3 糸の張力 - 6 -

7 座標, において 1 cos ( A.53) したがって ( A.54) これにより, 式 (A.46) は cos cos 0 3cos 3 ( また, 公式 1 cos sin より ) 16sin ( A.55) として得られる. ただし,θ> π/ のときは糸がたるむので, 式 (A.55) の成立は θ π/ の場合である. ちなみに, 振子を水平 (θ = π/) の位置で静止させ, 静かに放した場合の最下点 ( y 軸上 ) に来たときの糸の張力は θ=0 を式 (A.55) に代入し, であるから 3 となる. [ 参考 ] 二重振子 図 A.4 二重振子 長さ, 質量 の単振子を 個直列に接続した系の運動方 - 7 -

8 程式と固有角振動数, 固有モードを求める. まず, 運動エネルギー と位置エネルギー は ( A.56) これより, 式 (A.56) を本書の式 (.119) ラグランジュの運動方程式に当ては めて, ラグランジュの運動方程式は次のとおり 0 1,,, 次の運動方程式が得られる. 0 0 ( A.57) よって = (A.58) 1 次モード 次モード (A.59) ( A.60) A. 運動方程式 ( 微分方程式 ) を解く 運動方程式が立てられたので, これらの方程式を解く. ( 1) 微分方程式を, 演算子を用いて解く単振動の微分方程式を解く. はじめに, 階の微分方程式である次の演算子の説明をする ( A.61)

9 第 階微分方程式を解くと次式の公式に当てはまる. ( A.6) さて, ここで とおいて, これを演算子と呼ぶ. 式 (A.61) について演算子を用いて書くと 0 ( A.63) である. ( A.64) 演算子を用いて 0 ( A.65) よって ( A.66), したがって式 (A.6) に当てはめると ( A.67) ただし, 1 の虚数単位とする. 式 (A.67) は次のオイラーの公式を利用して正弦波または余弦波に変換する. ~~~~~~~~~~~~~ オイラーの公式 cossin ~~~~~~~~~~~~~ 式 (A.67) の実時間 t に対する実現象の変位 x は実数値なので, 虚数部が零であ - 9 -

10 るためには,, は複素数でなければならない. したがって,, を複素 数とおいて, オイラーの公式を用いると sin ( A.68) ただし ( A.69) は初期条件により決まる位相, は初期条件により決まる振幅として求められる. 不減衰角固有振動数は式 (A.6), 式 ( A.1), と同じく, 式 (A.69) が得られる. ******************** 公式 (A.6) の考え方の説明をする. 階の微分を行う式 (A.61) が成立するには, その解 x は正弦波 (sin) または余弦波 (cos) あるいは指数導関数以外に存在しない. 式 (A.61) には / の項が無いので周期振動が確実であり x は正弦波または余弦波と置くことができるが,/ の項がある場合にも対応できるように, 解 x を指数関数と置く. さらに, 解は 次方程式から求めるため, 根であるからは と の重ね合わせとする. したがって, 解は式 (A.6) と置くことができる. ********************* ( ) 微分方程式をラプラス変換により解く 0 (A.70) の場合について, 両辺をラプラス変換を行う ʆ : ( ラプラス変換記号 ) ʆ 0 ( A.71) すなわち ʆ ʆ 0 ( A.7) ここで

11 ʆ ʆ 0 0 (A.73) 初期条件 0, 0 を代入し ʆ ʆ ( A.74) よって ʆ ʆ 0 ( A.75) しかるに ʆ ( A.76) ここでラプラス逆変換を行う. ʆ - 1 : ( ラプラス逆変換記号 ) ʆ ʆ ʆ ʆ ʆ ( A.77) ここでラプラスの原関数 と像関数 の関係 ( ラプラス変換 逆変換 ) より, ~~~~~~~~~~~~~~~~ 参考 : ラプラス変換と逆変換 ʆ, ʆ cos, sin, ~~~~~~~~~~~~~~~ cos sin sin ( A.78) ただし

12 (A.79) とする. tan ( A.80) このように式 (A.78) において変位 x の時間による変化が式 (A.68) と同様に 求められる. 不減衰角固有振動数は式 (A.6), 式 ( A.1), 式 ( A.69) と同じく, 式 (A.79) が得られる

13 付録 B いろいろな連続体の振動モードと固有振動数 はじめに 連続体の固有角振動数の算出方法に両端単純支持の梁の横曲げを取り上げ, 円 板や長方形板などの連続体の固有モードと固有振動数をまとめた. B.1 梁の横振動 細い梁の横振動を考える. x 軸上の梁のたわみの変位 y は, 距離 x と時間 t の 関数として表されるので, これを x t y, とおく. 梁の長さ l, 断面積 A, 材料の 縦弾性係数 E, 断面二次モーメント I とおく.dx 部分の左側に作用するせん断力 F, 曲げモーメント M とすると ( 図 B.1 参照 ) 0 x dx x y y F F df y Adx t 図 B.1 棒の横振動 dx 部分の右側に作用するせん断力 F M y F EI ( B.1) x x x F df は y df EI x x x y EI dx x ここで材料の密度を ρ とおくと, dx 部分の質量は Adx ( B.), これにせん断力の作用

14 する方向で運動方程式を立てると より したがって 質量 加速度 = 外力 y y y y ( B.3) Adx EI EI dx EI t x x x x x x 均一な棒の場合は,EI= 一定. ここで, 曲げ振動の基本方程式が求まった. 運動方程式の解 x t X x, 変数 t だけを含む関数 Rt 関数 で表す. 基準関数を y y A 0 EI ( B.4) t x x 4 y y A EI 0 ( B.5) 4 t x y, を変数分離法によって求めるために, 変数 x だけを含む y として x, t X x Rt X x; 振動モードを決定する. t A 基準座標を R ;, B 0 0 を任意定数, 固有角振動数 ω とすると R t A0 cost B0 sint したがって y x t X x A cost B sint, 0 0 (B.6) 式 (B.6) を式 (B.5) の基礎方程式に代入して固有角振動数等を求めることにな る. 代入すると 0 より 4 d X EI dx x 4 AX ( x) 0 ここで A 4, EI EI A とおくと

15 4 d X dx x x 式 (B.7) の一般解は したがって, 上式を改めて X x X x 4 Ce x 0 4 X ( B.7) として式 (B.7) に代入すると 4 4, j ( j は虚数単位 ) より, 一般解は X x Ae jx Be jx Ce x De C1 cosh x C sinh x C3 cosx C4 sin x (B.8) となる. 任意定数 C 1 ~ C 4 は境界条件から求める. いま, 梁の境界条件を両端単純支持 ( 支持 - 支持 ) の場合を想定して解を求める. 単純支持の場合は, たわみと曲げモーメントは零 (0) であるから x y y 0, 0 x EI であるから, x 0 X, d X dx x 0 (B.9) である. 同様に y l においても式 (B.9) が成り立つ. ~~~~~~~~~~~~~~ ここで次の関係式を利用する. e sinh x x e cosh x x e e x x d sinh x dx d cosh x dx cosh x sinh x ~~~~~~~~~~~~~~ これより x0 のとき X x 0, x l のとき X x 0, d X dx d X dx x x 0 0 (B.10)

16 したがって, 式 (B.8) と式 (B.10) より C 1C3 C 1C3 0 0 C1 cosh l C sinh l C3 cos l C4 sin l 0 C1 cosh l C sinh l C3 cos l C4 sin l 0 式 (B.11) から (B.11) C 1 C C3 sin l 0 式 (B.1) の振動数方程式が得られた. 0 (B.1) よって 固有角振動数は l n ( n= 1,, 3, ) n n l EI A 基準関数は任意定数 C(=C 4 ) として ( n= 1,, 3, ) ( B.13) X x x Csinx Csinn l が求められた. 後述の表 B. 1 を参照されたい. ( B.14) B. いろいろな連続体形状の固有モードと固有振動数 連続体の代表例として前述 B.1 に両端単純支持の梁を取り上げた. 次にいろい ろな形状の連続体の固有振動数を列挙する. 図 B.4, 表 B.~ 表 B.4 の数値は 日 本機械学会編 : 機械工学便覧, 基礎編 a 機械力学, 丸善,007 を参照した. なお, f, f : 固有振動数 ( 単位 :Hz あるいは s - 1 ), : 固有角振動数 ( 単位 :rad/s). ( 1) 棒の縦振動

17 f n n l E (B.15) n : 振動系の次数 E : 縦弾性係数 : 密度 ( 単位体積当たりの質量 ) c E : 波の伝播速度 l : 棒の長さ n : 自由 - 自由 ; 1,, 3 3, 固定 - 自由 ;, 3 /, 5 ( ) 丸棒のねじり振動 一端固定, 他端自由 1 / 3 /, f n n l G (B.16) n : 振動系の次数, G : 横弾性係数 n 0,1,,3, : 密度 ( 単位体積当たりの質量 ) l : 棒の長さ n : 自由 - 自由 ; 1,, 3 3, 固定 - 自由 ;, 3 /, 5 ( 3) 梁の横振動 1 / 3 /, f n n EI (B.17) l A n : 振動系のモード次数 ( 図 B. 参照 ) A : 梁の断面積 I : 断面二次モーメント E : 縦弾性係数 n : 表 B. 1 参照 : 密度 ( 単位体積当たりの質量 ) l : 棒の長さ

18 節 節 節 節の数 : 節の数 :3 1 次モード 次モード 節の数 :4 3 次モード ( a) 両端自由 ( 自由 - 自由 ) 腹 1 次モード 次モード 腹 節 腹 3 次モード 腹 (b) 片持ち梁 ( 固定 - 自由 ) 図 B. 梁の曲げ振動モード

19 表 B. 1 式 (B.17) の n 1 1 次 次 3 3 次 単純支持 ( 支持 - 支持 ) π π 3π 片持ち梁 ( 固定 - 自由 ) 片持ち梁 ( 固定 - 支持 ) 自由 - 自由 固定 - 固定 ( 4) 完全輪の半径方向振動 ( 円形輪の中央線が伸縮する半径方向変位のみの場合 ) f n 1 n R E (B.18) n : 全波形の数 ( モード次数 )( 図 B. 3 参照 ) E : 縦弾性係数, : 密度 ( 単位体積当たりの質量 ) R : 輪の半径 ( 輪の中央線 ) n n 3 n 4 図 B.3 完全輪の振動モード ( 5) 円形膜の横振動 f ns T ns (B.19) R a T : 単位長さ当たりの一定の張力 (N/m) a : 単位面積当たりの質量 ns の n : 周方向の次数, s : 半径方向の次数 ( 節円の数 )( 図 B.4 参照 )

20 図 B.4 円形膜の固有モードと λ ただし, 円周の拘束 ( 境界 ) 条件は外周 R に沿って, 半径 R, 周方向 θ, 時間 t に関して 0 とする. ( 6) 円板の横振動 外内周に沿って一様な境界条件をもつ円板 f ns 1 Par D R ( B.0) a 0.33 ポアソン比 n : 節直径の数 s : 節円の数固有モードは図 B.4 と同じ R : 半径 D : 円板の曲げ剛性 1 E : 縦弾性係数, t: 板厚 a : 単位面積当たりの質量 (kg/m ) Par : パラメータ表 B. 参照 - 0 -

21 s n 表 B. 式 (B.0) の Par 周辺自由 周辺固定 ν に無関係 ( 7) 円筒殻の半径方向振動 平均半径 R をもつ薄い円筒殻 1 Par R E f n ( B.1) h R 0.0, 0.3 m, n : 振動数の次数 ( ここでは円筒殻長さ方向の波 m は考慮しない.) n : 周方向波数 ( 図 B.5 参照 ) 拘束条件 ; 両端に沿って単純支持辺 L : 円筒殻の両端 h : 円筒厚 R : 円筒殻の半径 x 0, : 密度 ( 単位体積当たりの質量 ) Par : パラメータ表 B.3 参照 L - 1 -

22 節直径の数 :6 節直径 円筒 ( 両端に沿って単純支持辺 ) 図 B.5 円筒殻の振動モード (n=5 次モード ) 表 B. 3 式 (B.1) の Par L/R モード ( n) ( 8) 長方形板の横振動 1 周辺単純支持 f mn D m n a a b (B.) m, n : 振動次数周辺単純支持 (x, y 方向の半波数 ) 図 B.6 参照 a : 単位面積当たりの質量 kg m a : 横長さ (x 軸 ) b : 縦長さ (y 軸 ) 3 Eh D : 板の曲げ剛性 1 1 E : ヤング率 ( 縦弾性係数 ) h : 板厚 : ポアソン比 =0.3 長方形の周辺境界条件の違いによる固有周波数を求めると, 次の式となる. - -

23 f mn 1 Par D a (B.3) a Par : パラメータ表 B. 4 参照. a 節 節 b 節 節 1.1 次 1. 次.1 次. 次 図 B.6 長方形板の固有モード 表 B. 4 式 (B.3) の Par 境界条件 モード 0.4 / a b 自由 ( モードは節が対角線, 曲線等存在する ) 単純支持 (a/b=1 を境に m, n の値は逆になる ) (m=1,n=1) 14.6(m=1,n=1) 19.74(m=1,n=1) 3.08(m=1,n=1) 71.55(m=1,n=1) 16.19(m=1,n=) 7.41(m=1,n=) 49.35(m=1,n=) 61.68(m=,n=1) 101.(m=,n=1) (m=1,n=3) 43.86(m=,n=1) 49.35(m=,n=1) 98.70(m=1,n=) 150.5(m=3,n=1) (m=1,n=4) 49.35(m=1,n=3) 78.96(m=,n=) 111.0(m=3,n=1) 19.6(m=4,n=1) 固定 (ν に 無関係 )