領域シンポ発表

Size: px

Start display at page:

Download "領域シンポ発表"

ふじきみこうじょう
5 years ago
Views:

1 1 次元の減衰運動の中の強制振動 ) ( f d d d d d e f e ce ) ( si ) ( 1 ) ( cos ω =ω -γ とおくと一般解は外力 f()=f siω の場合 f d d d d si f ce f ce si ) cos( cos si ) cos( この一般解は 1

2 φ は外力と変位との間の位相差で a 時間が経つと第 1 項は無視できるこの場合の振幅を Dとすると D これに対して一定の外力 (Ω=) の場合の変位をD とすると D 1 D 4 1 h f ただし h,

3 外力の振動数に対する振幅外力の振動数に対する位相のずれ 3 h=,ξ =1( 粘性がなく振動数が同じ ) 場合共鳴が起こり時間と共に発達する位相に関して外力の振動数が小さい時は変位はそれについてゆくしかし Ω が大きくなると外力の変位とは逆 (φ =π ) になるこの境界はΩ =ω である時間がたった時の変位は Dsi 速度は Dcos で変位に対してπ /だけ位相が進む

4 パラメータ励振振り子の振動数 ω の倍で支点を上下させる 4 m d d si ブランコをこぐときにはブランコの周期の半分の周期で腰を上下させるこうして支点から重心までの距離を周期的に変えるとブランコの振幅を大きくすることができる単振り子の糸を周期的に長くしたり短くしたりすると単振り子の振幅が大きくなる

5 パラメータ励振するのは d m si d 5 時間 d の間に加えられる仕事 ( 角 θ を変位をみたときの仕事 ) は dw d si d si d d と常に正であり振り子はたえずエネルギーを増しふれは次第に大きくなる例 : ブランコ身近な振動と共鳴の数々ブランコ

6 今までは時間変動だけの振動であった 6 1 次元進行波 - 時間と距離 - の記述の仕方を眺める 1 次元波動の記述の仕方

7 一つの単色波の伝搬 η=re( A ep {i ( -ω)} ) 7 = 1 A = ( > 1 ) π/ A cos( ) where A A r A i, a Ai A r

8 二つの単色波の重ね合わせ ) )cos( si ( 1 1 si 1 1 si 搬送波 (carrier waves) 包絡線振幅を表しエネルギーに関係する (evelopes) 8

9 仮定 9 : 位相速度二つの速度が定義される : : 群速度 = 1, = 複数の波動の重ね合わせ :ω=f(): 分散関係 c g d : d に依存しない非分散に依存する分散性がある c 群速度 ( エネルギーの速度 ) : 位相速度

10 空間的に孤立した波にかたまりについても波数空間では局在したところに集中して存在する 1

11 次元波動 z z z S z λ λ λ z 物理空間波数空間三角形 ABC の面積 S A C B,, K K z z 方向に平行波数空間で考える方がベクトルとしてわかりやすい! K 11

12 地球の自転が小さい場合の方程式系 1 大気の状態を記述する方程式系 ( ただし地球回転と大気粘性は無視 ) dv d d d 1 C T p g v dq d p p p p R 完全系 p g : 重力加速度, T p R : 乾燥空気の気体定数, p : 気圧, (A1) ( 運動方程式 ) (A) ( 連続の式 ) (A3) ( 熱力学第一法則の式 ) (A4) ( 温位の定義式 ) 二次元の場合 u, w, p, ρθ, の 5変数に対してつの式がある 5 v ( u, w) : 風の二成分, R / C p, : 密度, : 鉛直方向の単位ベクトル, : 温位, : 定圧比熱, : 代表的な気圧例えば ( 1 hpa), p / p : エクスナー関数 ( 無次元の圧力 ) p C p

13 浮力 13 同高度で二つの異なる空気塊がある場合どちらが重いかどうかは空気塊の重さつまりその密度 ρ で判断する軽い密度 ( 高温 ) の空気塊の方が重い ( 低温 ) 空気塊より上向きの力を得る ( その力を浮力という ) 気体の状態方程式より p ρ T ρ 1/T (p が同圧の場合 ) 異なる高度にある空気塊間の浮力は? 空気塊を断熱的に別の空気塊と同じ高度まで持ってきて同圧にして密度あるいはその時の温度で比較する. すべての空気塊を規準気圧まで持ってきてその高度における密度 ( あるいは温度の逆数 ) で比較する. 温位が適している空気塊が異なる高度であっても密度とも簡単に結びつき浮力を相互比較できる便利な量

14 鉛直方向の運動方程式静的安定度 (1) dw d = - g - 1 静力学平衡の ( 空気塊が鉛直方向に運動をしない ) 状態 = - g - 1 p z p 気圧が周囲のものと変わらないものと仮定し空気塊が平衡状態から僅かに断熱的に上方に移動したときの鉛直方向の運動方程式は w = d z d = - g - 1 p z = - g - 1 p z (- g) = g - = g - 14

15 ここで = T p = RT 静的安定度 () p p - R Cp ( 温位と気温の関係式 ) ( 状態方程式 ) 15 (z) + ( d dz )z よって ( 温位減率で成層状態を決める ) d z d = - g (d dz )z - N z, N : Bru-Vaisala 振動数 N g d d z ブラント-バイサラ振動数時間の逆数の二乗

16 Bru-Vaisala 振動数の意味 (1) 16 浮力により自分で上下方向に振動する振動数 --- 標準的な大気で N = 1 1 ( /s) N < : 絶対不安定大気 N = : 中立 ( 等温位 ) 大気 N > : 安定大気

17 Bru-Vaisala 振動数の意味 () 17 断熱的に持ち上げ w = 浮力 w = -w ma w = 浮力 w = -w ma 浮力 w = 1/N ime 等温位大気 (N = ): 気塊を持ち上げても浮力を得ないので振動しない安定度 ( / z) が高い程浮力が大きくなり振動数 (N) が大きくなる

18 内部重力波 1 方向に復元力がある場合 18

19 1 z 19 α 波面の進む方向 N 1 N どうしてこうなる? N cos

20 u u ) w w (z (z) (z) u w z b u w z b wn N g z ( C ) p (1) () (3) (4) ( b g ) 大気の安定度 : 地球大気は正値を持つ安定成層をする

21 w w z w N δ = 1 w Aep[ i( z )] Aep[ i( z )] N i b Aep[ i( z )] δ =1 u Aep[ i( z )] N 図非静力学系 ( 実曲線 ) と静力学系 ( 実直線 ) における無次元の波数 (/) と振動数 (ω/n) の関係の比較 N が与えられた時の (,,ω) の関係分散関係

22 3/ 3/, N N. N 3 / 3 / N N c g ベクトルのと c g は直交する!

23 cos ) / ( ) / ( / z z z N 波数を波長に直すと K 群速度の方向波数空間 Λ Λ z 物理空間波面の進む方向 α z 3

24 物理的解釈 4 内部重力波 : 流体粒子の振動ある波面斜面における質点の運動 : 斜面方向に g cosα の力が働く g ω ω N

剛体過去問解答例 2 1.1) 長さの棒の慣性モーメントは公式より l I G = Ml /12 A 点のまわりは平行軸の定理より 2 2 I A = Ml /12 + M ( l / 2) = Ml 2 / 3 B y 2) 壁からの垂直抗力を R, 床からの垂直抗力と摩擦力を N,f とすると

剛体過去問解答例 2 1.1) 長さの棒の慣性モーメントは公式より l I G = Ml /12 A 点のまわりは平行軸の定理より 2 2 I A = Ml /12 + M ( l / 2) = Ml 2 / 3 B y 2) 壁からの垂直抗力を R, 床からの垂直抗力と摩擦力を N,f とすると剛体過去問解答例. 長さの棒の慣性モーメントは公式より l G l A 点のまわりは平行軸の定理より A l l l B y 壁からの垂直抗力を R, 床からの垂直抗力と摩擦力を N,f とすると運動方程式は方向 : R f, y 方向 : y N l 回転 : G { N f R cos } A 静止しているとき方向の力と力のモーメントがつり合うので y ~ より R ' また摩擦力が最大静止摩擦力より大きいとはしごは動き出すので