運動方程式の基本 ニュートンの第一法則 力 = 質量 加速度 大気や海洋に加わる力を, 思いつくだけ挙げてみよう 重力, 圧力傾度力, コリオリ力, 摩擦力 水平方向に働く力に下線をつけよう. したがって水平方向の運動方程式は 質量 水平加速度 = コリオリ力 + 圧力傾度力 + 摩擦力 流体の運動
|
|
- あやか とどろき
- 6 years ago
- Views:
Transcription
1 2. 浅水方程式系の導出 見延庄士郎 ( 海洋気候物理学研究室 ) 第 1 回まとめ 1/2 二つの変数の関係の強さを表す統計量は相関であり, 最小値は -1, 最大値は +1, 無相関は である. 過去数十年間の ( 気象庁は 3 年 ) 月ごとの平均値を, 月平均データの平年値または気候値という. 観測値から平年値を引いたものが, 偏差である. 連続する n 個のデータを平均して, 中央のデータの値に置き換える平滑化が, 移動平均である. 異なる物理量を比較する際には, それぞれの標準偏差で割る, 規格化がよく用いられる. エル ニーニョは,3~6 年に一度熱帯東太平洋の海洋表面水温が異常に上昇する現象で, 北半球の冬季に最盛期. 特に強かったのは,1982/83 年,1997/98 年, 年 地球温暖化で, 極端に強いエルニーニョが増える予想. 1 2 第 1 回まとめ2/2 エルニーニョの最盛期の特徴 海洋表面水温は, 赤道東太平洋で高くなる. 降水は, 赤道東太平洋で増加する. 気圧は, 熱帯東太平洋で低気圧偏差, 海大陸付近で高気圧偏差, となる. このような東西の気圧のシーソーを, 南方振動と呼ぶ. 東西風速は, 赤道中央太平洋で, 西風偏差 ( 東向きの風速偏差 ) となる. 上空 2Paの等圧面高度は, 熱帯東太平洋で高気圧偏差, そのすぐ南北で低気圧偏差, さらにその先の北米北西部と南大洋の6S, 11W 付近に高気圧偏差を示す. 水温躍層は, 赤道東太平洋で深く, 赤道西太平洋で浅くなる. エルニーニョの9ヶ月前の特徴 水温躍層は, 赤道中央太平洋で深くなる. 本日の目的 本授業の後のほうで, エルニーニョのメカニズムを議論したい. たとえば, 概要は ttps:// に解説がある. ただし,3 年後期にふさわしい授業として扱うには, 数式も必要である. そこで, 支配方程式系を導出しよう. 方程式系は, 未知変数と式の数が同じで, 系 (= システム ) の挙動を記述できる連立方程式. 3 つの式を導出する. そのうち二つは, 東西方向と南北方向の運動方程式で, もう一つの式は, これらの式で必要となる圧力を与える式である. たった 3 つの式だが, 大気海洋の気候力学で重要な流れと波の性質を深く学ぶことができる. 3 4
2 運動方程式の基本 ニュートンの第一法則 力 = 質量 加速度 大気や海洋に加わる力を, 思いつくだけ挙げてみよう 重力, 圧力傾度力, コリオリ力, 摩擦力 水平方向に働く力に下線をつけよう. したがって水平方向の運動方程式は 質量 水平加速度 = コリオリ力 + 圧力傾度力 + 摩擦力 流体の運動方程式は単位体積について立てることにして, 質量は密度 ρ で表される. 座標系と変数を導入 (u,v) 独立変数, : 東向き 北向き座標, t: 時間未知の従属変数 u, v: 東向き 北向き速度, p: 圧力, ρ: 密度 = 単位体積 (1 m 3 ) 当たりの質量 5 6 コリオリ力のまとめ ( 復習資料から ) 角速度 Ω> で回転する回転系で, 速度 v で運動する物体は, 進行方向右向きに加速度 2vΩ で加速するように見える. この時, 物質の質量を m とすると, 力 2 v Ωm が働いたように見え, これをコリオリ力という. 地球の角速度ベクトルの, それぞれの地点での鉛直成分を 2 倍したものを, コリオリ パラメータといい,f で表す. 地球の自転角速度を Ω(=2 π/ 自転周期 ) として, 緯度 θ でのコリオリパラメータは, f = 2 Ω sin θ である. コリオリ力のみが働く場合に, 回転系での運動方程式は, 以下のとおりである. du fv dt = dv dt = f u du/dt と / 1/2 du/dt( ディーユー ディーティー ) は全微分またはラグランジェ微分または全微分と呼ばれる, 運動につき従って測る微分である. 一方, / ( ラウンドユー ラウンドティー ) は偏微分またはオイラー微分と呼ばれる. 固定した地点で測る微分である. 両者の関係は, 合成関数の偏微分で得られる. 一般に,u(,,t), =(ξ), =(ξ), t=t(ξ) であれば, 合成関数の偏微分は, 右の依存関係図より du dξ ξ ξ ξ である. ここで, ξ=tとすれば du dt t t u t ξ 7 8
3 du/dt と / 2/2 一つ前のスライドの du dt t で = u, = v に他ならないので, 上式は du u u u = + u + v dt とかける. 同様に, dv v v v = + u + v dt ただしこれら2 式の, 第二 三項は非線形で厄介なのと, また現象の振幅が小さいなら第一項よりも小さくなるので微小振幅を仮定して du/dt / と近似 ( 線形化 ) する. すると二つ前のスライドの下の, 回転系での運動方程式は以下の通りとなる. u v = fv = f u 圧力経度力のみが働く場合 コリオリ力が働かず ( 慣性系 ), 圧力経度力のみが働く場合の, 単位体積当たりの流体 ( 空気または水 ) についての運動方程式は次のとおり. = = ρ ρ コリオリ力と圧力経度力と摩擦力 ( 文字のまま ) が同時に働く場合は, = fv + 方向の摩擦力 (2.1a) ρ = fu + 方向の摩擦力 (2.1b) ρ 符号に注意 9 1 定常で摩擦力が働かない場合 この場合の, 運動方程式は次のとおり. fv = (2.2a) + fu = ρ ρ (2.2b) 2 力のバランス ( 平衡 ) が成り立っており, その 2 力は, コリオリ力と圧力傾度力である. このように吹く風を地衡風, 海の流れを地衡流という. Q. 右の圧力場が与えられている. 黒点における地衡風のベクトルを書き込め 浅水方程式へ 式の数が 2 つで, 未知変数の数が 4 つ. 閉じていない 式の数と未知変数の数が一致する, 閉じた方程式を得たい. そのためには, 簡素化が不可欠. 密度が一定 ρ である流体を仮定する. 密度が定数である流体を一層流体という. また, 地球流体の場合, 高さよりも水平方向のスケールがはるかに長いので, 一層流体のモデルを浅水モデルという. H 最後に残る従属変数, 残らない従属変数, 係数 u, v, p ρ 流体表面 静止時の流体表面 流体底面 (,,t): 表面変位 H(, ): 静止状態の流体の厚さ ( 水深 ) 11 12
4 圧力, ウォーミングアップ Q. ある場所での圧力は, その上に乗っている物体の質量 重力加速度で与えられる. 地表面に 1 気圧の圧力がかかっている場合に,1 m 2 の面積の上に存在する大気の質量を求めよ. ただし簡単のために, 1 気圧を 1 Pa とし, 重力加速度を 1 m/s 2 とする. なお,( ヘクト ) は 1 を意味する SI 接頭辞で, 1, Pa= 1, Pa = 1, N/m 2 である. A. 1, / 1 = 1, kg = 1 ton 1 層流体の圧力と圧力傾度 Q1. 流体表面から深さ d での圧力を求めよ. ただし流体の密度は一定で ρ, 重力加速度は g とする. A1. p = ρ g d Q2. A1 の圧力を図の ( 水面の静止状態からの高さ ) と,Z( 静止状態の水面からの水深 ) とで表せ. A2. p = ρ g (+Z) Q3. A2の圧力の 微分 ( 水平微分 ) を求めよ A3. p = ( ρg( + Z)) = ρg H Z d 流体表面 p ρ 静止時の流体表面 流体底面 層流体の運動方程式 一つ前のスライドの圧力を,(2.1ab) 式に代入すると, その式のρは今 ρ であることに注意して, = fv g + 方向の摩擦力 = fu g + 方向の摩擦力 となる. これが1 層流体の運動方程式である. だけでなく,u, v も,, tのみの関数 ( 深さによらない ) 次に, 摩擦力について考えよう. (2.3a) (2.3b) 15 摩擦力 摩擦力には2 通りある. 一つは, 大気や海洋が運動して, 周囲の空気や水あるいは地表面や海底から抵抗を受けて, 運動を減速させるように働く摩擦力である. これはしばしば簡潔に, 以下のように表され, レイリー ダンピングと呼ばれる. =... ru, =... rv したがって,(2.3ab) の, 摩擦力の項を置き換えて = fv g ru (2.4a) = fu g rv もう一つは, 海洋上を吹く風によって, 海が引きずられることによって生じる摩擦力である. これは海洋を駆動する上で重要である. (2.4b) 16
5 風応力による水柱の加速 単位面積 (1 m 1 m) で深さHの水柱を考える応力 = 面に働く力 水柱の重さはρ H この水柱にかかる風応力は 加速度 = 力 / 質量なので, この効果による加速 度は u =... + ρh したがって,(2.3ab) の摩擦力の項を置き換えて = fv g + ρh = fu g + ρ H u t (2.5a) (2.5b) 次に について考えよう. 上 17 連続の式 1/2 一層流体では密度が一定なので, 流体が圧縮されない ( 非圧縮 ). この場合, 体積が保存する. 下のように3 辺が,, Ηである直方体への流入流出を考えよう. v(, + / 2) w = t u( + / 2, ) Η v(, / 2) u( / 2, ) 側面からの入出の差の分だけ, 表面が上昇または下降する. これを差分式 ( などが有限の長さを持つ式. これが無限小なら微分式 ) 18 連続の式 2/2 = H u +, u, 2 2 H v, + v, 2 2 したがって, 両辺を微小体積 Ηで割って差分を微分で置き換えると ( ノートに計算してください ), = H + (2.6) という 1 層流体における連続の式が得られる. 結局, 閉じた浅水方程式系として,(2.4ab) と (2.6), または (2.5ab) と (2.6) が得られた. 第二回まとめ f をコリオリ パラメータといい, Ω を地球自転の角速度,θ を緯度として,f=2 Ω sin θ で表される. 密度一様な一層流体の, 運動方程式と連続の式からなる, 閉じた方程式系を, 浅水方程式系という. 摩擦が, 運動に伴う減速の場合, 潜水方程式系の運動方程式は, = fv g ru (2.4a) = fu g rv 摩擦が風応力である場合には, 運動方程式は = fv g + (2.5a) = fu g + H H 連続の式は ρ = H + (2.6) ρ (2.4b) (2.5b) 19
運動方程式の基本 座標系と変数を導入 (u,v) ニュートンの第一法則 力 = 質量 加速度 大気や海洋に加わる力を, 思いつくだけ挙げてみよう 重力, 圧力傾度力, コリオリ力, 摩擦力 水平方向に働く力に下線をつけよう. したがって水平方向の運動方程式は 質量 水平加速度 = コリオリ力 + 圧
2. 潜水方程式系の導出 見延庄士郎 ( 海洋気候物理学研究室 ) minobe@mail.sci.okudai.ac.jp 第 1 回まとめ 1/2 二つの変数の関係の強さを表す統計量は相関であり, 最小値は -1, 最大値は +1, 無相関は である. 過去数十年間の ( 気象庁は 3 年 ) 月ごとの平均値を, 月平均データの平年値または気候値という. 観測値から平年値を引いたものが, 偏差である.
More information3. 重力波と沿岸 赤道ケルビン波 2014 年 9 月 30 日 16:35 見延庄士郎 ( 海洋気候物理学研究室 ) 予習課題 : 以下の you tube のビデオを見ておくこと. 個々のビデオは全部は見ずに, 雰囲気がつかめる程度見ればいい.
3. 重力波と沿岸 赤道ケルビン波 2014 年 9 月 30 日 16:35 見延庄士郎 ( 海洋気候物理学研究室 ) minobe@sci.hokudai.ac.jp 予習課題 : 以下の you ube のビデオを見ておくこと. 個々のビデオは全部は見ずに, 雰囲気がつかめる程度見ればいい. 大気の重力波 : hp://www.youube.com/wach?v=yxnkzecu3be 津波シミュレーション
More information4. ロスビー波と 1.5 層モデル見延庄士郎 ( 海洋気候物理学研究室 ) 予習ビデオ : NASA JPL の海面高度偏差でエルニーニョを見るページ. 赤道ケルビン波が見えるか? の動画 ht
4. ロスビー波と 1.5 層モデル見延庄士郎 ( 海洋気候物理学研究室 ) 予習ビデオ : NASA JPL の海面高度偏差でエルニーニョを見るページ. 赤道ケルビン波が見えるか? hps://sealevel.jpl.nasa.gov/elnino015/inde.hml の動画 hps://sealevel.jpl.nasa.gov/elnino015/1997vs015-animaed-800.gif
More information3. 重力波と沿岸 赤道ケルビン波見延庄士郎 ( 海洋気候物理学研究室 ) 予習課題 を印刷解答の上, 授業冒頭で提出してください. 予習のための課題なので
3. 重力波と沿岸 赤道ケルビン波見延庄士郎 ( 海洋気候物理学研究室 ) 予習課題 hp://www.sci.hokudai.ac.jp/~minobe/class/po_clm/03_e_propagaion.pd を印刷解答の上, 授業冒頭で提出してください. 予習のための課題なので, 冒頭のみ提出を受け付けます. また以下の you ube のビデオを見ておくこと. 個々のビデオは全部は見ずに,
More information浅水方程式 順圧であるためには, 静水圧近似が必要 Dw Dt + コリオリ力 = 1 p + 粘性 g ρ z w が u, v に比べて小さい 運動の水平距離に対して水深が浅い 浅水 海は深いが, 水平はさらに広い 最大 1 万 km 浅水方程式 : u, v, の式 水平 2 次元の解 D D
流体地球科学第 11 回 東京大学大気海洋研究所准教授藤尾伸三 ttp://ovd.aori.u-tokyo.ac.jp/ujio/2015ciba/ ujio@aori.u-tokyo.ac.jp 2016/1/8 順圧流の運動方程式 流体の密度が一様ならば, 圧力 静水圧 の水平勾配は鉛直一様 海面の高さによる水平圧力勾配のみ ηx,y px, y, z = ρ g dz = ρgη z p x
More information4. ロスビー波と 1.5 層モデル見延庄士郎 ( 海洋気候物理学研究室 ) 予習ビデオ : NASA JPL の海面高度偏差でエルニーニョを見るページ. 赤道ケルビン波が見えるか? の動画 ht
4. ロスビー波と 1.5 層モデル見延庄士郎 ( 海洋気候物理学研究室 ) 予習ビデオ : NASA JPL の海面高度偏差でエルニーニョを見るページ. 赤道ケルビン波が見えるか? hps://sealevel.jpl.nasa.gov/elnino015/inde.hml の動画 hps://sealevel.jpl.nasa.gov/elnino015/1997vs015-animaed-800.gif
More information2 図微小要素の流体の流入出 方向の断面の流体の流入出の収支断面 Ⅰ から微小要素に流入出する流体の流量 Q 断面 Ⅰ は 以下のように定式化できる Q 断面 Ⅰ 流量 密度 流速 断面 Ⅰ の面積 微小要素の断面 Ⅰ から だけ移動した断面 Ⅱ を流入出する流体の流量 Q 断面 Ⅱ は以下のように
3 章 Web に Link 解説 連続式 微分表示 の誘導.64 *4. 連続式連続式は ある領域の内部にある流体の質量の収支が その表面からの流入出の合計と等しくなることを定式化したものであり 流体における質量保存則を示したものである 2. 連続式 微分表示 の誘導図のような微小要素 コントロールボリューム の領域内の流体の増減と外部からの流体の流入出を考えることで定式化できる 微小要素 流入
More information流体地球科学第 7 回 力のバランス永遠に回れるバランス ( 以下, 北半球 =コリオリ力は進行方向の右向き ) 慣性振動 : 遠心力 =コリオリ力 地衡風 : コリオリ力 = 圧力傾度力 東京大学大気海洋研究所准教授藤尾伸三
流体地球科学第 7 回 力のバランス永遠に回れるバランス ( 以下, 北半球 =コリオリ力は進行方向の右向き ) 慣性振動 : 遠心力 =コリオリ力 地衡風 : コリオリ力 = 圧力傾度力 東京大学大気海洋研究所准教授藤尾伸三 http://ovd.aori.u-tokyo.ac.jp/fujio/205chiba/ fujio@aori.u-tokyo.ac.jp F C F A 旋衡風 : 遠心力
More informationパソコンシミュレータの現状
第 2 章微分 偏微分, 写像 豊橋技術科学大学森謙一郎 2. 連続関数と微分 工学において物理現象を支配する方程式は微分方程式で表されていることが多く, 有限要素法も微分方程式を解く数値解析法であり, 定式化においては微分 積分が一般的に用いられており. 数学の基礎知識が必要になる. 図 2. に示すように, 微分は連続な関数 f() の傾きを求めることであり, 微小な に対して傾きを表し, を無限に
More information図 7-: コリオリ力の原理 以下では 回転台の上で物体が運動したとき 物体にはたらくみかけの力を定量的に求めてみる 回転台は角速度 で回転していて 回転台に乗っている観測者から見た物体の速度ベクトルの動径方向の成分を u 接線方向の成分を v とする 図 7-3: 回転台の上での物体の運動 はじめ
7 大気の力学 () 7. コリオリ力 水平面内に気圧の差があると風が吹く原因となる 気圧の差によって空気塊 高にはたらく力を気圧傾度力 (pessue gaient foce) という 気圧傾度力は等 圧線と直角に 高圧側から低圧側に向かってはたらく しかし 天気図で見ら れる風向と 等圧線とのなす角は直角ではないことが多い これは 地球の自 高転の影響によって 地球上を運動する空気塊にコリオリ力
More information線積分.indd
線積分 線積分 ( n, n, n ) (ξ n, η n, ζ n ) ( n-, n-, n- ) (ξ k, η k, ζ k ) ( k, k, k ) ( k-, k-, k- ) 物体に力 を作用させて位置ベクトル A の点 A から位置ベクトル の点 まで曲線 に沿って物体を移動させたときの仕事 W は 次式で計算された A, A, W : d 6 d+ d+ d@,,, d+ d+
More information領域シンポ発表
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 φ は外力と変位との間の位相差で a 時間が経つと 第 1 項は無視できる この場合の振幅を
More information外から中心に投げたボールの動画 1 中心に向かってまっすぐ投げる 回転盤でボールをキャッチ 円盤の回転速度とボールの速度を合わせれば, 投げたボールを取れる ( 投げた人にはボールが回ってくるように見える ) 投げてからの時間は, 回転の半周期 円盤の外から見る図斜めに飛んでいく 投げた人が見る図コ
流体地球科学第 6 回 外から中心に投げたボールは? 回転盤の外から見た図 ( ) 期待される位置, ( ) 実際の位置 間違った図 1 間違った図 2 正しい図 東京大学大気海洋研究所准教授藤尾伸三 http://ovd.aori.u-tokyo.ac.jp/ujio/215chiba/ ujio@aori.u-tokyo.ac.jp 215/11/2 最終更新日 215/11/24 ボールは左
More information数値計算で学ぶ物理学 4 放物運動と惑星運動 地上のように下向きに重力がはたらいているような場においては 物体を投げると放物運動をする 一方 中心星のまわりの重力場中では 惑星は 円 だ円 放物線または双曲線を描きながら運動する ここでは 放物運動と惑星運動を 運動方程式を導出したうえで 数値シミュ
数値計算で学ぶ物理学 4 放物運動と惑星運動 地上のように下向きに重力がはたらいているような場においては 物体を投げると放物運動をする 一方 中心星のまわりの重力場中では 惑星は 円 だ円 放物線または双曲線を描きながら運動する ここでは 放物運動と惑星運動を 運動方程式を導出したうえで 数値シミュレーションによって計算してみる 4.1 放物運動一様な重力場における放物運動を考える 一般に質量の物体に作用する力をとすると運動方程式は
More information7 渦度方程式 総観規模あるいは全球規模の大気の運動を考える このような大きな空間スケールでの大気の運動においては 鉛直方向の運動よりも水平方向の運動のほうがずっと大きい しかも 水平方向の運動の中でも 収束 発散成分は相対的に小さく 低気圧や高気圧などで見られるような渦 つまり回転成分のほうが卓越
7 渦度方程式 総観規模あるいは全球規模の大気の運動を考える このような大きな空間スケールでの大気の運動においては 鉛直方向の運動よりも水平方向の運動のほうがずっと大きい しかも 水平方向の運動の中でも 収束 発散成分は相対的に小さく 低気圧や高気圧などで見られるような渦 つまり回転成分のほうが卓越している そこで 回転成分に着目して大気の運動を論じる 7.1 渦度 大気の回転成分を定量化する方法を考えてみる
More information粘性 接する流体の間に抵抗が働き, その大きさは速度差に比例 粘性係数流体と壁の間にも抵抗 ( 摩擦 ) は働く 粘性や摩擦の大きさは, 物質の性質で異なる ネバネバ 粘性が大 流体は一緒に動こうとするサラサラ 粘性が小 ある流体粒子の速度を U, 上側 下側の流体粒子の速度を U U, U L 粘
流体地球科学第 9 回 東京大学大気海洋研究所准教授藤尾伸三 http://ovd.aori.-tokyo.ac.jp/fjio/2015chiba/ fjio@aori.-tokyo.ac.jp 2015/12/11 最終更新日 2015/12/9 北大西洋の中緯度の表層における収支 質量保存 = 0 (q 2 q 1 ) + (0 q 3 ) = 0 南から入ってくる流量 q 1, 北に出て行く流量
More informationMicrosoft PowerPoint - 夏の学校(CFD).pptx
/9/5 FD( 計算流体力学 ) の基礎理論 性能 運動分野 夏の学校 神戸大学大学院海事科学研究科勝井辰博 流体の質量保存 流体要素内の質量の増加率 [ 単位時間当たりの増加量 ] 単位時間に流体要素に流入する質量 流体要素 Fl lm (orol olm) v ( ) ガウスの定理 v( ) /9/5 = =( ) b=b =(b b b ) b= b = b + b + b アインシュタイン表記
More informationÿþŸb8bn0irt
折戸の物理 スペシャル補習 http://orito-buturi.com/ NO.3 今日の目的 : 1 微分方程式をもう一度 三角関数の近似について学ぶ 3 微分の意味を考える 5. 起電力 の電池, 抵抗値 の抵抗, 自己インダクタンス のコイルとスイッチを用いて右図のような回路をつくった 始めスイッチは 開かれている 時刻 t = でスイッチを閉じた 以下の問に答えよ ただし, 電流はコイルに
More informationMicrosoft Word - thesis.doc
剛体の基礎理論 -. 剛体の基礎理論初めに本論文で大域的に使用する記号を定義する. 使用する記号トルク撃力力角運動量角速度姿勢対角化された慣性テンソル慣性テンソル運動量速度位置質量時間 J W f F P p .. 質点の並進運動 質点は位置 と速度 P を用いる. ニュートンの運動方程式 という状態を持つ. 但し ここでは速度ではなく運動量 F P F.... より質点の運動は既に明らかであり 質点の状態ベクトル
More information以下 変数の上のドットは時間に関する微分を表わしている (ex. 2 dx d x x, x 2 dt dt ) 付録 E 非線形微分方程式の平衡点の安定性解析 E-1) 非線形方程式の線形近似特に言及してこなかったが これまでは線形微分方程式 ( x や x, x などがすべて 1 次で なおかつ
以下 変数の上のドットは時間に関する微分を表わしている (e. d d, dt dt ) 付録 E 非線形微分方程式の平衡点の安定性解析 E-) 非線形方程式の線形近似特に言及してこなかったが これまでは線形微分方程式 ( や, などがすべて 次で なおかつそれらの係数が定数であるような微分方程式 ) に対して安定性の解析を行ってきた しかしながら 実際には非線形の微分方程式で記述される現象も多く存在する
More information風成循環その 2 北半球の亜熱帯循環系 (w <, curl z τ < ) ) δ w 2π ζ 相対渦度 ζ = を仮定する ( あるいは, 右辺第 2 項が小さい ) 風が与える渦度と惑星渦度の移流とバランス βv = fw ( スベルドラップ平衡 ) w < であれば, v < ( 南向き
風成循環 流体地球科学第 12 回 東京大学大気海洋研究所准教授藤尾伸三 http://ovd.aori.u-tokyo.ac.jp/fujio/215chiba/ fujio@aori.u-tokyo.ac.jp 216/1/22 北半球の亜熱帯循環系 (w
More information領域シンポ発表
総観気象学 のシラバス ( 第 1 回 ) 日本の四季と気圧配置 ( 第 2 回 ) 大気運動の全体像 ( 第 3 回 ) 大気の運動方程式の導出 (1) ( 第 4 回 ) スケール解析 ( 第 5 回 ) コリオリ力 ( 第 6 回 ) 地衡風 温度風 ( 第 7 回 ) ロスビー波 ( 第 8 回 ) 温帯低気圧 ( 傾圧不安定波 ) ( 第 9 回 ) 前線 ( 第 10 回 ) 乾燥対流
More information木村の物理小ネタ 単振動と単振動の力学的エネルギー 1. 弾性力と単振動 弾性力も単振動も力は F = -Kx の形で表されるが, x = 0 の位置は, 弾性力の場合, 弾性体の自然状態の位置 単振動の場合, 振動する物体に働く力のつり合
単振動と単振動の力学的エネルギー. 弾性力と単振動 弾性力も単振動も力は F = -x の形で表されるが, x = の位置は, 弾性力の場合, 弾性体の自然状態の位置 単振動の場合, 振動する物体に働く力のつり合いの位置 である たとえば, おもりをつるしたばねについて, ばねの弾性力を考えるときは, ばねの自然長を x = とし, おもりの単振動で考える場合は, おもりに働く力がつり合った位置を
More information<4D F736F F D20824F B CC92E8979D814696CA90CF95AA82C691CC90CF95AA2E646F63>
1/1 平成 23 年 3 月 24 日午後 6 時 52 分 6 ガウスの定理 : 面積分と体積分 6 ガウスの定理 : 面積分と体積分 Ⅰ. 直交座標系 ガウスの定理は 微分して すぐに積分すると元に戻るというルールを 3 次元積分に適用した定理になります よく知っているのは 簡単化のため 変数が1つの場合は dj ( d ( ににします全微分 = 偏微分 d = d = J ( + C d です
More informationDVIOUT-SS_Ma
第 章 微分方程式 ニュートンはリンゴが落ちるのを見て万有引力を発見した という有名な逸話があります 無重力の宇宙船の中ではリンゴは落ちないで静止していることを考えると 重力が働くと始め静止しているものが動き出して そのスピードはどんどん大きくなる つまり速度の変化が現れることがわかります 速度は一般に時間と共に変化します 速度の瞬間的変化の割合を加速度といい で定義しましょう 速度が変化する, つまり加速度がでなくなるためにはその原因があり
More informationMicrosoft PowerPoint - 1章 [互換モード]
1. 直線運動 キーワード 速さ ( 等速直線運動, 変位 ) 加速度 ( 等加速度直線運動 ) 重力加速度 ( 自由落下 ) 力学 I 内容 1. 直線運動 2. ベクトル 3. 平面運動 4. 運動の法則 5. 摩擦力と抵抗 6. 振動 7. 仕事とエネルギー 8. 運動量と力積, 衝突 9. 角運動量 3 章以降は, 運動の向きを考えなければならない 1. 直線運動 キーワード 速さ ( 等速直線運動,
More information宇宙機工学 演習問題
宇宙システム工学演習 重力傾度トルク関連. 図に示すように地球回りの円軌道上を周回する宇宙機の運動 を考察する 地球中心座標系を 系 { } 軌道面基準回転系を 系 { } 機体固定系を 系 { } とする 特に次の右手直交系 : 地心方向単位ベクトル 軌道面内 : 進行方向単位ベクトル 軌道面内 : 面外方向単位ベクトル 軌道面外 を取る 特に この { } Lol Horiotl frme と呼ぶ
More informationMicrosoft PowerPoint - 流体力学の基礎02(OpenFOAM 勉強会 for geginner).pptx
~ 流体力学の基礎 ~ 第 2 回 流体静力学 2011 年 10 月 22 日 ( 土 ) 講習会のスケジュール概要 ( あくまでも現時点での予定です ) 流体力学の基礎 第 1 回目 2011.09 流体について 第 2 回目 2011.10 流体静力学 第 3 回目 2011.11/12 流体運動の基礎理論 1 第 4 回目 2012.01 流体運動の基礎理論 2 第 5 回目 2012.02
More informationPowerPoint プレゼンテーション
第 2 回 : 海洋数値モデルの概要 海洋大循環モデル海洋プリミティブ方程式圧力項 粘性 拡散項の取り扱い境界条件初期値 境界値問題海洋モデルの種類地球流体力学的現象地衡流平衡と温度風平衡海洋アイソスタシーとリモートセンシングどこまでダウンスケールできるか 海洋大循環モデル 地球流体力学的な現象を数値的に表現する数値モデル地球流体力学的現象 : 地球回転と成層の効果が支配的な現象 海洋大循環モデルの適用範囲
More informationニュートン重力理論.pptx
3 ニュートン重力理論 1. ニュートン重力理論の基本 : 慣性系とガリレイ変換不変性 2. ニュートン重力理論の定式化 3. 等価原理 4. 流体力学方程式とその基礎 3.1 ニュートン重力理論の基本 u ニュートンの第一法則 = 力がかからなければ 等速直線運動を続ける u 等速直線運動に見える系を 慣性系 と呼ぶ ² 直線とはどんな空間の直線か? ニュートン理論では 3 次元ユークリッド空間
More information<4D F736F F F696E74202D20836F CC8A C58B858B4F93B982A882E682D1978E89BA814091B28BC68CA48B E >
バットの角度 打球軌道および落下地点の関係 T999 和田真迪 担当教員 飯田晋司 目次 1. はじめに. ボールとバットの衝突 -1 座標系 -ボールとバットの衝突の前後でのボールの速度 3. ボールの軌道の計算 4. おわりに参考文献 はじめに この研究テーマにした理由は 好きな野球での小さい頃からの疑問であるバッテングについて 角度が変わればどう打球に変化が起こるのかが大学で学んだ物理と数学んだ物理と数学を使って判明できると思ったから
More information微分方程式による現象記述と解きかた
微分方程式による現象記述と解きかた 土木工学 : 公共諸施設 構造物の有用目的にむけた合理的な実現をはかる方法 ( 技術 ) に関する学 橋梁 トンネル ダム 道路 港湾 治水利水施設 安全化 利便化 快適化 合法則的 経済的 自然および人口素材によって作られた 質量保存則 構造物の自然的な性質 作用 ( 外力による応答 ) エネルギー則 の解明 社会的諸現象のうち マスとしての移動 流通 運動量則
More information自由落下と非慣性系における運動方程式 目次無重力... 2 加速度計は重力加速度を測れない... 3 重量は質量と同じ数値で kg が使える... 3 慣性系における運動方程式... 4 非慣性系における運動方程式... 6 見かけの力... 7 慣性系には実在する慣
自由落下と非慣性系における運動方程式 1 1 2 3 4 5 6 7 目次無重力... 2 加速度計は重力加速度を測れない... 3 重量は質量と同じ数値で kg が使える... 3 慣性系における運動方程式... 4 非慣性系における運動方程式... 6 見かけの力... 7 慣性系には実在する慣性力があるか... 7 1 2 無重力 (1) 非慣性系の住人は無重力を体感できる (a) 併進的な加速度運動をしている非慣性系の住人
More informationOCW-iダランベールの原理
講義名連続体力学配布資料 OCW- 第 2 回ダランベールの原理 無機材料工学科准教授安田公一 1 はじめに今回の講義では, まず, 前半でダランベールの原理について説明する これを用いると, 動力学の問題を静力学の問題として解くことができ, さらに, 前回の仮想仕事の原理を適用すると動力学問題も簡単に解くことができるようになる また, 後半では, ダランベールの原理の応用として ラグランジュ方程式の導出を示す
More information第五回まとめ 2/2 東岸域では, 北風がエクマン流を通じて湧昇をもたらす. これを沿岸湧昇と呼ぶ. 沿岸湧昇域での, 局所的な大気海洋結合変動現象であるカリフォルニアニーニョ, ニンガルーニーニョなどが発見されている. エクマン流と地衡流の関係の仮説 前回学んだエクマン流が, どう地衡流と関係する
6. 海洋風成循環と海面高度 見延庄士郎 ( 海洋気候物理学研究室 ) 復習課題 http://www.sci.hokudai.ac.jp/~minobe/class/po_clm/05_review_ekman.pdf を印刷解答の上, 授業冒頭で提出してください. 冒頭のみ提出を受け付けます. 東岸域でのエクマン流 Q. 大洋の東岸域では, 岸に沿う風が吹くと, エクマン流が岸に直交するので,
More informationMicrosoft Word - 中村工大連携教材(最終 ).doc
音速について考えてみよう! 金沢工業大学 中村晃 ねらい 私たちの身の回りにはいろいろな種類の波が存在する. 体感できる波もあれば, できない波もある. その中で音は体感できる最も身近な波である. 遠くで雷が光ってから雷鳴が届くまで数秒間時間がかかることにより, 音の方が光より伝わるのに時間がかかることも経験していると思う. 高校の物理の授業で音の伝わる速さ ( 音速 ) は約 m/s で, 詳しく述べると
More information本日の内容 6. 海洋風成循環と海面高度 見延庄士郎 ( 海洋気候物理学研究室 ) 予習課題 : 下のニュースに目を通しておくこと. 温暖化進めば今世紀半ばまで年 4 兆円支出米政府監査院 ( 朝日新聞 2017/11/25)
本日の内容 6. 海洋風成循環と海面高度 見延庄士郎 ( 海洋気候物理学研究室 ) 予習課題 : 下のニュースに目を通しておくこと. 温暖化進めば今世紀半ばまで年 4 兆円支出米政府監査院 ( 朝日新聞 2017/11/25) http://www.shi.com/rticles/askbt1t4fkbtuhbi002.html クローズアップ科学 黒潮が 12 年ぶりに大蛇行高潮の被害誘発 シラス漁に打撃
More informationMicrosoft PowerPoint - 熱力学Ⅱ2FreeEnergy2012HP.ppt [互換モード]
熱力学 Ⅱ 第 章自由エネルギー システム情報工学研究科 構造エネルギー工学専攻 金子暁子 問題 ( 解答 ). 熱量 Q をある系に与えたところ, 系の体積は膨張し, 温度は上昇した. () 熱量 Q は何に変化したか. () またこのとき系の体積がV よりV に変化した.( 圧力は変化無し.) 内部エネルギーはどのように表されるか. また, このときのp-V 線図を示しなさい.. 不可逆過程の例を
More information<4D F736F F F696E74202D E94D58B9393AE82F AC82B782E982BD82DF82CC8AEE E707074>
地盤数値解析学特論 防災環境地盤工学研究室村上哲 Mrakam, Satoh. 地盤挙動を把握するための基礎. 変位とひずみ. 力と応力. 地盤の変形と応力. 変位とひずみ 変形勾配テンソルひずみテンソル ひずみテンソル : 材料線素の長さの 乗の変化量の尺度 Green-Lagrange のひずみテンソルと Alman のひずみテンソル 微小変形状態でのひずみテンソル ひずみテンソルの物理的な意味
More information実践海洋情報論 57 ここで mは流体の質量 (kg) u は流体の速度 (m sec -1 ) である しかし 緯度 φにおける角速度を直観的に理解することが困難である そこで 地球儀を北極上空から観察しよう 北極点に十字のマークを置くと 地球儀を反時計回りに回転させると 極点上の十字のマークも反
実践海洋情報論 56 3-3 大気 海洋へ及ぶ力と流れ 3-3-1 コリオリ力 コリオリ力は 自転する惑星上の流体に働く見かけの回転力であり 惑星渦度として定義される 惑星渦度は 惑星の回転を示す角速度に依存し 極域で最大となり 赤道上では0となる 北半球において 低気圧が反時計回りに回るのはコリオリ力によるものと説明されている しかし なぜ反時計回りになるか 十分に説明のなされた図書が少ない ここでは
More informationPowerPoint Presentation
Non-linea factue mechanics き裂先端付近の塑性変形 塑性域 R 破壊進行領域応カ特異場 Ω R R Hutchinson, Rice and Rosengen 全ひずみ塑性理論に基づいた解析 現段階のひずみは 除荷がないとすると現段階の応力で一義的に決まる 単純引張り時の応カーひずみ関係 ( 構成方程式 ): ( ) ( ) n () y y y ここで α,n 定数, /
More information物理演習問題
< 物理 > =0 問 ビルの高さを, ある速さ ( 初速 をとおく,において等加速度運動の公式より (- : -= t - t : -=- t - t (-, 式よりを消去すると t - t =- t - t ( + - ( + ( - =0 0 t t t t t t ( t + t - ( t - =0 t=t t=t t - 地面 ( t - t t +t 0 より, = 3 図 問 が最高点では速度が
More information<4D F736F F D2089F082AF82E997CD8A7796E291E A282EB82A282EB82C8895E93AE2E646F63>
いろいろな運動. 自由落下. 投げ上げ 3. 放物運動 4. 標的にボールを当てる 5. 斜面に向かって投げ上げる 6. ブレーキをかけた自動車 7. 摩擦のある斜面上を滑り落ちる物体 8. ばね振り子 ( 単振動 ) 9. 摩擦を受けるばね振り子. 補足 : 微分方程式の解き方 自由落下質量 の質点を高さ h の地点から初速 で落とした. 鉛直上向きを 軸正 の向き, 地表を原点とし, 重力加速度の大きさを
More informationChap2.key
. f( ) V (V V ) V e + V e V V V V ( ) V V ( ) E. - () V (0 ) () V (0 ) () V (0 ) (4) V ( ) E. - () V (0 ) () V (0 ) O r θ ( ) ( ) : (r θ) : { r cos θ r sn θ { r + () V (0 ) (4) V ( ) θ θ arg( ) : π π
More informationMicrosoft Word - cavitation.doc
音響キャビテーション ヤング ラプラスの式 A B dl g C ヤングは 二つの流体の境界に厚さの無視できる架空の膜が存在し しかもこの膜には張力が作用するというモデルを考えた. すなわち 図に示すように 気泡表面と交わる平面 ABC を考えたとき 表面の接線方向の単位長さ dl 当たり g の力が作用していると考え この力を表面張力とよんだ. 次元は [N/m][J/m] となる. 気相と液相が接する界面には表面張力が働く
More informationMicrosoft PowerPoint - H21生物計算化学2.ppt
演算子の行列表現 > L いま 次元ベクトル空間の基底をケットと書くことにする この基底は完全系を成すとすると 空間内の任意のケットベクトルは > > > これより 一度基底を与えてしまえば 任意のベクトルはその基底についての成分で完全に記述することができる これらの成分を列行列の形に書くと M これをベクトル の基底 { >} による行列表現という ところで 行列 A の共役 dont 行列は A
More informationÿþŸb8bn0irt
折戸の物理 スペシャル補習 http://oritobuturi.co/ NO.5(009..16) 今日の目的 : 1 物理と微分 積分について 微分方程式について学ぶ 3 近似を学ぶ 10. 以下の文を読み,[ ア ]~[ ク ] の空欄に適当な式をいれよ 物体物体に一定の大きさの力を加えたときの, 物体の運動について考え よう 右図のように, なめらかな水平面上で質量 の物体に水平に一定の大きさ
More information多次元レーザー分光で探る凝縮分子系の超高速動力学
波動方程式と量子力学 谷村吉隆 京都大学理学研究科化学専攻 http:theochem.kuchem.kyoto-u.ac.jp TA: 岩元佑樹 iwamoto.y@kuchem.kyoto-u.ac.jp ベクトルと行列の作法 A 列ベクトル c = c c 行ベクトル A = [ c c c ] 転置ベクトル T A = [ c c c ] AA 内積 c AA = [ c c c ] c =
More information伝熱学課題
練習問題解答例 < 第 章強制対流熱伝達 >. 式 (.9) を導出せよ (.6) を変換する 最初に の微分値を整理しておく (.A) (.A) これを用いて の微分値を求める (.A) (.A) (.A) (.A6) (.A7) これらの微分値を式 (.6) に代入する (.A8) (.A9) (.A) (.A) (.A) (.9). 薄い平板が温度 で常圧の水の一様な流れの中に平行に置かれている
More information木村の物理小ネタ ケプラーの第 2 法則と角運動量保存則 A. 面積速度面積速度とは平面内に定点 O と動点 P があるとき, 定点 O と動点 P を結ぶ線分 OP( 動径 OP という) が単位時間に描く面積を 動点 P の定点 O に
ケプラーの第 法則と角運動量保存則 A. 面積速度面積速度とは平面内に定点 O と動点 P があるとき, 定点 O と動点 P を結ぶ線分 OP( 動径 OP という が単位時間に描く面積を 動点 P の定点 O に関する面積速度の大きさ という 定点 O まわりを回る面積速度の導き方導き方 A ( x( + D, y( + D v ( q r ( A ( x (, y( 動点 P が xy 座標平面上を時刻
More information風力発電インデックスの算出方法について 1. 風力発電インデックスについて風力発電インデックスは 気象庁 GPV(RSM) 1 局地気象モデル 2 (ANEMOS:LAWEPS-1 次領域モデル ) マスコンモデル 3 により 1km メッシュの地上高 70m における 24 時間の毎時風速を予測し
風力発電インデックスの算出方法について 1. 風力発電インデックスについて風力発電インデックスは 気象庁 GPV(RSM) 1 局地気象モデル 2 (ANEMOS:LAWEPS-1 次領域モデル ) マスコンモデル 3 により 1km メッシュの地上高 70m における 24 時間の毎時風速を予測し 2000kW 定格風車の設備利用率として表示させたものです 数値は風車の定格出力 (2000kW)
More information<4D F736F F D20824F F6490CF95AA82C696CA90CF95AA2E646F63>
1/15 平成 3 年 3 月 4 日午後 6 時 49 分 5 ベクトルの 重積分と面積分 5 重積分と面積分 Ⅰ. 重積分 と で 回積分することを 重積分 といいます この 重積分は何を意味しているのでしょう? 通常の積分 (1 重積分 ) では C d 図 1a 1 f d (5.1) 1 f d f ( ) は 図形的には図 1a のように面積を表しています つまり 1 f ( ) を高さとしてプロットすると図
More information大阪大学物理 8 を解いてみた Ⅱ. 問 ( g cosq a sin q ) m - 台 B 上の観測者から見ると, 小物体は, 斜面からの垂直抗力 N, 小物体の重力 mg, 水平左向きの慣性力 ma を受け, 台 B の斜面と平行な向きに運動する したがって, 小物体は台 B の斜面に垂直な方
大阪大学物理 8 を解いてみた Ⅰ. 問 g 最高点の座標を y max とすると, 力学的エネルギー保存則より \ y m mgy 補足 max g max 小物体の運動方向に対する仕事は重力 ( 保存力 ) の斜面に沿った成分のみであり, 垂直抗力 ( 非保存力 ) の仕事は である よって, 力学的エネルギー保存則が成り立つ これを確かめてみよう 小物体は重力の斜面に沿った外力を受けながらその運動エネルギーを失っていく
More information横浜市環境科学研究所
周期時系列の統計解析 単回帰分析 io 8 年 3 日 周期時系列に季節調整を行わないで単回帰分析を適用すると, 回帰係数には周期成分の影響が加わる. ここでは, 周期時系列をコサイン関数モデルで近似し単回帰分析によりモデルの回帰係数を求め, 周期成分の影響を検討した. また, その結果を気温時系列に当てはめ, 課題等について考察した. 気温時系列とコサイン関数モデル第 報の結果を利用するので, その一部を再掲する.
More informationQ
埼玉工業大学機械工学学習支援セミナー ( 小西克享 ) 剛体の重心と自由運動 -1/8 テーマ 07: 剛体の重心と自由運動 一般的に剛体が自由に運動できる状態 ( 非拘束の状態 ) で運動するとき, 剛体は回転運動を伴った運動をします. たとえば, 棒の端を持って空中に放り投げると, 棒はくるくる回転しながら上昇してやがて地面に落ちてきます. 剛体が拘束されない状態で運動する様子を考察してみましょう.
More information( 慣性抵抗 ) 速度の 2 乗に比例流体中を進む物体は前面にある流体を押しのけて進む. 物 aaa 体の後面には流体が付き従う ( 渦を巻いて ). 前面にある速度 0 の流体が後面に移動して速度 vとなったと考えてよい. この流体の質量は単位時間内に物体が押しのける体積に比例するので,v に比例
空気抵抗があるときの自由落下 抵抗が速度に比例する場合 1. 絵を描く, 座標と情報, 記号を記入する x F0 v
More informationギリシャ文字の読み方を教えてください
埼玉工業大学機械工学学習支援セミナー ( 小西克享 ) 慣性モーメント -1/6 テーマ 01: 慣性モーメント (Momet of ietia) コマ回しをすると, 長い時間回転させるには重くて大きなコマを選ぶことや, ひもを早く引くことが重要であることが経験的にわかります. 遊びを通して, 回転の運動エネルギーを増やせば, 回転の勢いが増すことを学習できるので, 機械系の学生にとってコマ回しも大切な体験学習のひとつと言えます.
More informationChap3.key
区分求積法. 面積 ( )/ f () > n + n, S 長方形の和集合で近似 n f (n ) リーマン和 f (n ) 区分求積法 リーマン和 S S n n / n n f ()d リーマン積分 ( + ) + S (, f ( )) 微分の心 Zoom In して局所的な性質を調べる 積分の心 Zoom Ou して大域的な性質を調べる 曲線の長さ 領域の面積や体積 ある領域に含まれる物質の質量
More informationMicrosoft PowerPoint - elast.ppt [互換モード]
弾性力学入門 年夏学期 中島研吾 科学技術計算 Ⅰ(48-7) コンピュータ科学特別講義 Ⅰ(48-4) elast 弾性力学 弾性力学の対象 応力 弾性力学の支配方程式 elast 3 弾性力学 連続体力学 (Continuum Mechanics) 固体力学 (Solid Mechanics) の一部 弾性体 (lastic Material) を対象 弾性論 (Theor of lasticit)
More information2009 年 11 月 16 日版 ( 久家 ) 遠地 P 波の変位波形の作成 遠地 P 波の変位波形 ( 変位の時間関数 ) は 波線理論をもとに P U () t = S()* t E()* t P() t で近似的に計算できる * は畳み込み積分 (convolution) を表す ( 付録
遠地 波の変位波形の作成 遠地 波の変位波形 ( 変位の時間関数 ) は 波線理論をもとに U () t S() t E() t () t で近似的に計算できる は畳み込み積分 (convolution) を表す ( 付録 参照 ) ここで St () は地震の断層運動によって決まる時間関数 1 E() t は地下構造によって生じる種々の波の到着を与える時間関数 ( ここでは 直達 波とともに 震源そばの地表での反射波や変換波を与える時間関数
More information剛体過去問解答例 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 ' また 摩擦力が最大静止摩擦力より大きいとはしごは動き出すので
More informationえられる球体について考えよ 慣性モーメント C と体積 M が以下の式で与えられることを示せ (5.8) (5.81) 地球のマントルと核の密度の平均値を求めよ C= kg m 2, M= kg, a=6378km, rc=3486km 次に (5.82) で与えら
5.5 慣性モーメント (5-42) 式で与えられたマッカラーの公式は 扁球惑星体の重力加速度とその主な慣性モーメントを関連づけている その公式を使うことで 探査飛行や軌道上を周回する宇宙船によって 例えば慣性モーメントを束縛している惑星の重力場を計測することができる 慣性モーメントは惑星全体の形や内部の密度分布を反映するため 惑星の内部構造を調べるために慣性モーメントの数値を利用することができる
More information<4D F736F F F696E74202D D488A778AEE B4F93B982CC8AEE A2E707074>
宇宙工学基礎 ( 軌道の基礎 松永三郎 機械宇宙学科 機械宇宙システム専攻 ニュートンの法則 第 法則 力が作用作用しないしない限り 質点質点は静止静止ないしはないしは一定速度一定速度で運動するする ( 慣性の法則 慣性空間 慣性座標系慣性座標系の定義第 法則 慣性座標系におけるにおける質点質点の運動 p F ( pɺ t ( F: 全作用力, pmv: 並進運動量 ( 質量と速度速度の積 慣性系を規準規準としてとして時間微分時間微分を行うことにことに注意第
More information今週の内容 後半全体のおさらい ラグランジュの運動方程式の導出 リンク機構のラグランジュの運動方程式 慣性行列 リンク機構のエネルギー保存則 エネルギー パワー 速度 力の関係 外力が作用する場合の運動方程式 粘性 粘性によるエネルギーの消散 慣性 粘性 剛性と微分方程式 拘束条件 ラグランジュの未
力学 III GA 工業力学演習 X5 解析力学 5X 5 週目 立命館大学機械システム系 8 年度後期 今週の内容 後半全体のおさらい ラグランジュの運動方程式の導出 リンク機構のラグランジュの運動方程式 慣性行列 リンク機構のエネルギー保存則 エネルギー パワー 速度 力の関係 外力が作用する場合の運動方程式 粘性 粘性によるエネルギーの消散 慣性 粘性 剛性と微分方程式 拘束条件 ラグランジュの未定乗数法
More information<4D F736F F D20824F E B82CC90FC90CF95AA2E646F63>
1/1 平成 3 年 6 月 11 日午前 1 時 3 分 4 ベクトルの線積分 4 ベクトルの線積分 Ⅰ. 積分の種類 通常の物理で使う積分には 3 種類あります 積分変数の数に応じて 線積分 ( 記号 横(1 重 d, dy, dz d ( ine: 面積分 ( 記号 縦 横 ( 重 線 4 ベクトルの線積分 重積分記号 ddy, dydz, dzdz ds ( Surface: 1 重積分記号
More information計算機シミュレーション
. 運動方程式の数値解法.. ニュートン方程式の近似速度は, 位置座標 の時間微分で, d と定義されます. これを成分で書くと, d d li li とかけます. 本来は が の極限をとらなければいけませんが, 有限の小さな値とすると 秒後の位置座標は速度を用いて, と近似できます. 同様にして, 加速度は, 速度 の時間微分で, d と定義されます. これを成分で書くと, d d li li とかけます.
More information物理学 II( 熱力学 ) 期末試験問題 (2) 問 (2) : 以下のカルノーサイクルの p V 線図に関して以下の問題に答えなさい. (a) "! (a) p V 線図の各過程 ( ) の名称とそのと (& きの仕事 W の面積を図示せよ. # " %&! (' $! #! " $ %'!!!
物理学 II( 熱力学 ) 期末試験問題 & 解答 (1) 問 (1): 以下の文章の空欄に相応しい用語あるいは文字式を記入しなさい. 温度とは物体の熱さ冷たさを表す概念である. 物体は外部の影響を受けなければ, 十分な時間が経過すると全体が一様な温度の定常的な熱平衡状態となる. 物体 と物体 が熱平衡にあり, 物体 と物体 が熱平衡にあるならば, 物体 と物体 も熱平衡にある. これを熱力学第 0
More information第 5 章 構造振動学 棒の振動を縦振動, 捩り振動, 曲げ振動に分けて考える. 5.1 棒の縦振動と捩り振動 まっすぐな棒の縦振動の固有振動数 f[ Hz] f = l 2pL である. ただし, L [ 単位 m] は棒の長さ, [ 2 N / m ] 3 r[ 単位 Kg / m ] E r
第 5 章 構造振動学 棒の振動を縦振動, 捩り振動, 曲げ振動に分けて考える 5 棒の縦振動と捩り振動 まっすぐな棒の縦振動の固有振動数 f[ Hz] f l pl である ただし, L [ 単位 m] は棒の長さ, [ N / m ] [ 単位 Kg / m ] E は (5) E 単位は棒の材料の縦弾性係数 ( ヤング率 ) は棒の材料の単位体積当りの質量である l は境界条件と振動モードによって決まる無
More information3 数値解の特性 3.1 CFL 条件 を 前の章では 波動方程式 f x= x0 = f x= x0 t f c x f =0 [1] c f 0 x= x 0 x 0 f x= x0 x 2 x 2 t [2] のように差分化して数値解を求めた ここでは このようにして得られた数値解の性質を 考
3 数値解の特性 3.1 CFL 条件 を 前の章では 波動方程式 f x= x = f x= x t f c x f = [1] c f x= x f x= x 2 2 t [2] のように差分化して数値解を求めた ここでは このようにして得られた数値解の性質を 考える まず 初期時刻 t=t に f =R f exp [ik x ] [3] のような波動を与えたとき どのように時間変化するか調べる
More informationFEM原理講座 (サンプルテキスト)
サンプルテキスト FEM 原理講座 サイバネットシステム株式会社 8 年 月 9 日作成 サンプルテキストについて 各講師が 講義の内容が伝わりやすいページ を選びました テキストのページは必ずしも連続していません 一部を抜粋しています 幾何光学講座については 実物のテキストではなくガイダンスを掲載いたします 対象とする構造系 物理モデル 連続体 固体 弾性体 / 弾塑性体 / 粘弾性体 / 固体
More informationMicrosoft PowerPoint - シミュレーション工学-2010-第1回.ppt
シミュレーション工学 ( 後半 ) 東京大学人工物工学研究センター 鈴木克幸 CA( Compter Aded geerg ) r. Jaso Lemo (SC, 98) 設計者が解析ツールを使いこなすことにより 設計の評価 設計の質の向上を図る geerg の本質の 計算機による支援 (CA CAM などより広い名前 ) 様々な汎用ソフトの登場 工業製品の設計に不可欠のツール 構造解析 流体解析
More informationMicrosoft Word _資料2(最終版).doc
資料 -2 今回の潮位の特性 ( 追加 ) 2.1 被災時の潮位特性 (1) 概要被災時の潮位特性について 観測潮位データを基に調べた 先ず 28 年 2 月 24 日を含む 日本海沿岸の検潮データを収集して整理した 先ず 小樽 深浦 佐渡 富山 能登 舞鶴 境 浜田 ( いずれも気象庁 ) の8 地点の2 月の観測潮位と潮位偏差を収集して整理した 8 地点の位置を下図に示した 次に 伏木富山 新湊
More informationMicrosoft Word - 微分入門.doc
基本公式 例題 0 定義式 f( ) 数 Ⅲ 微分入門 = の導関数を定義式にもとづいて計算しなさい 基本事項 ( f( ), g( ) が微分可能ならば ) y= f( ) g( ) のとき, y = y= f( ) g( ) h( ) のとき, y = ( f( ), g( ) が微分可能で, g( ) 0 ならば ) f( ) y = のとき, y = g ( ) とくに, y = のとき,
More informationMicrosoft PowerPoint - zairiki_3
材料力学講義 (3) 応力と変形 Ⅲ ( 曲げモーメント, 垂直応力度, 曲率 ) 今回は, 曲げモーメントに関する, 断面力 - 応力度 - 変形 - 変位の関係について学びます 1 曲げモーメント 曲げモーメント M 静定力学で求めた曲げモーメントも, 仮想的に断面を切ることによって現れる内力です 軸方向力は断面に働く力 曲げモーメント M は断面力 曲げモーメントも, 一つのモーメントとして表しますが,
More information1/12 平成 29 年 3 月 24 日午後 1 時 1 分第 3 章測地線 第 3 章測地線 Ⅰ. 変分法と運動方程式最小作用の原理に基づくラグランジュの方法により 重力場中の粒子の運動方程式が求められる これは 力が未知の時に有効な方法であり 今のような 一般相対性理論における力を求めるのに使
/ 平成 9 年 3 月 4 日午後 時 分第 3 章測地線 第 3 章測地線 Ⅰ. 変分法と運動方程式最小作用の原理に基づくラグランジュの方法により 重力場中の粒子の運動方程式が求められる これは 力が未知の時に有効な方法であり 今のような 一般相対性理論における力を求めるのに使う事ができる 最小作用の原理 : 粒子が時刻 から の間に移動したとき 位置 と速度 v = するのが ラグランジュ関数
More information<4D F736F F D2097CD8A7793FC96E582BD82ED82DD8A E6318FCD2E646F63>
- 第 章たわみ角法の基本式 ポイント : たわみ角法の基本式を理解する たわみ角法の基本式を梁の微分方程式より求める 本章では たわみ角法の基本式を導くことにする 基本式の誘導法は各種あるが ここでは 梁の微分方程式を解いて基本式を求める方法を採用する この本で使用する座標系は 右手 右ネジの法則に従った座標を用いる また ひとつの部材では 図 - に示すように部材の左端の 点を原点とし 軸線を
More informationスライド 1
Q & A Q: 猫ひねりができるのって猫だけですか?(2 人 ) A: 動物が専門でない私にとっては難しい質問です おそらく 猫に近い ヒョウ チーター ヤマネコ等はできるのではないかと思います ちなみに猫とともにペットの代表である犬は 猫ほどうまくないようです 犬を抱っこしていて落としてしまい 怪我をする犬もけっこういるようです 猫はおそらく大丈夫です Q: 空気抵抗は気圧に比例したりしますか?
More information応用数学Ⅱ 偏微分方程式(2) 波動方程式(12/13)
偏微分方程式. 偏微分方程式の形 偏微分 偏導関数 つの独立変数 をもつ関数 があるとき 変数 が一定値をとって だけが変化したとす ると は だけの関数となる このとき を について微分して得られる関数を 関数 の に関する 偏微分係数 略して偏微分 あるいは偏導関数 pil deiie といい 次のように表される についても同様な偏微分を定義できる あるいは あるいは - あるいは あるいは -
More information[ 振動の発生 ] 第 1 章 土木振動学序論 [ 振動の発生 ] 外力と内力内力が釣り合って静止釣り合って静止した状態 :[: [ 平衡状態 ] 振動の発生振動の発生 :[ 平衡状態 ] が破られ 復元力復元力が存在すると振動が発生する つまり (1) 平衡 ( 静止 ) 状態が破られる (2)
[ 振動の発生 ] 第 1 章 土木振動学序論 [ 振動の発生 ] 外力と内力内力が釣り合って静止釣り合って静止した状態 :[: [ 平衡状態 ] 振動の発生振動の発生 :[ 平衡状態 ] が破られ 復元力復元力が存在すると振動が発生する つまり (1) 平衡 ( 静止 ) 状態が破られる (2) 運動が発生する (3) 復元力があると 振動状態になる 自由度 (degree of freedom)
More information2011年度 大阪大・理系数学
0 大阪大学 ( 理系 ) 前期日程問題 解答解説のページへ a a を自然数とする O を原点とする座標平面上で行列 A= a の表す 次変換 を f とする cosθ siθ () >0 および0θ
More information物性物理学 I( 平山 ) 補足資料 No.6 ( 量子ポイントコンタクト ) 右図のように 2つ物質が非常に小さな接点を介して接触している状況を考えましょう 物質中の電子の平均自由行程に比べて 接点のサイズが非常に小さな場合 この接点を量子ポイントコンタクトと呼ぶことがあります この系で左右の2つ
物性物理学 I( 平山 ) 補足資料 No.6 ( 量子ポイントコンタクト ) 右図のように つ物質が非常に小さな接点を介して接触している状況を考えましょう 物質中の電子の平均自由行程に比べて 接点のサイズが非常に小さな場合 この接点を量子ポイントコンタクトと呼ぶことがあります この系で左右のつの物質の間に電位差を設けて左から右に向かって電流を流すことを行った場合に接点を通って流れる電流を求めるためには
More information1/10 平成 29 年 3 月 24 日午後 1 時 37 分第 5 章ローレンツ変換と回転 第 5 章ローレンツ変換と回転 Ⅰ. 回転 第 3 章光速度不変の原理とローレンツ変換 では 時間の遅れをローレンツ変換 ct 移動 v相対 v相対 ct - x x - ct = c, x c 2 移動
/ 平成 9 年 3 月 4 日午後 時 37 分第 5 章ローレンツ変換と回転 第 5 章ローレンツ変換と回転 Ⅰ. 回転 第 3 章光速度不変の原理とローレンツ変換 では 時間の遅れをローレンツ変換 t t - x x - t, x 静止静止静止静止 を導いた これを 図の場合に当てはめると t - x x - t t, x t + x x + t t, x (5.) (5.) (5.3) を得る
More information第 4 週コンボリューションその 2, 正弦波による分解 教科書 p. 16~ 目標コンボリューションの演習. 正弦波による信号の分解の考え方の理解. 正弦波の複素表現を学ぶ. 演習問題 問 1. 以下の図にならって,1 と 2 の δ 関数を図示せよ δ (t) 2
第 4 週コンボリューションその, 正弦波による分解 教科書 p. 6~ 目標コンボリューションの演習. 正弦波による信号の分解の考え方の理解. 正弦波の複素表現を学ぶ. 演習問題 問. 以下の図にならって, と の δ 関数を図示せよ. - - - δ () δ ( ) - - - 図 δ 関数の図示の例 δ ( ) δ ( ) δ ( ) δ ( ) δ ( ) - - - - - - - -
More information.( 斜面上の放物運動 ) 目的 : 放物運動の方向の分け方は, 鉛直と水平だけではない 図のように, 水平面から角 だけ傾いた固定した滑らかな斜面 と, 質量 の小球を用意する 原点 から斜面に垂直な向きに, 速さ V で小球を投げ上げた 重力の加速度を g として, 次の問い に答えよ () 小
折戸の物理 演習編 ttp://www.orito-buturi.co/ N..( 等加速度運動目的 : 等加速度運動の公式を使いこなす 問題を整理する能力を養う ) 直線上の道路に,A,B の 本の線が 5. の間隔で道路に 垂直に交差して引かれている この線上を一定の加速度で運 動しているトラックが通過する トラックの先端が A を通過してか ら後端が B を通過するまでの時間は.8s であった
More informationスライド 1
暫定版修正 加筆の可能性あり ( 付録 ) 準備 : 非線形光学効果 (). 絵解き : 第二高調波発生. 基本波の波動方程式 3. 第二高調波の波動方程式 4. 二倍分極振動 : ブランコ 5. 結合波動方程式へ 6. 補足 : 非線形電気感受率 ( 複素数 ) 付録 43 のアプローチ. 分極振動とは振動電場に誘われて伸縮する電気双極子の集団運動. 電気感受率と波動方程式の関係を明らかにする 3.
More information破壊の予測
本日の講義内容 前提 : 微分積分 線形代数が何をしているかはうろ覚え 材料力学は勉強したけど ちょっと 弾性および塑性学は勉強したことが無い ー > ですので 解らないときは質問してください モールの応力円を理解するとともに 応力を 3 次元的に考える FM( 有限要素法 の概略 内部では何を計算しているのか? 3 物が壊れる条件を考える 特に 変形 ( 塑性変形 が発生する条件としてのミーゼス応力とはどのような応力か?
More informationvecrot
1. ベクトル ベクトル : 方向を持つ量 ベクトルには 1 方向 2 大きさ ( 長さ ) という 2 つの属性がある ベクトルの例 : 物体の移動速度 移動量電場 磁場の強さ風速力トルクなど 2. ベクトルの表現 2.1 矢印で表現される 矢印の長さ : ベクトルの大きさ 矢印の向き : ベクトルの方向 2.2 2 個の点を用いて表現する 始点 () と終点 () を結ぶ半直線の向き : ベクトルの方向
More information2017年度 長崎大・医系数学
07 長崎大学 ( 医系 ) 前期日程問題 解答解説のページへ 以下の問いに答えよ () 0 のとき, si + cos の最大値と最小値, およびそのときの の値 をそれぞれ求めよ () e を自然対数の底とする > eの範囲において, 関数 y を考える この両 辺の対数を について微分することにより, y は減少関数であることを示せ また, e< < bのとき, () 数列 { } b の一般項が,
More informationMicrosoft PowerPoint - 測量学.ppt [互換モード]
8/5/ 誤差理論 測定の分類 性格による分類 独立 ( な ) 測定 : 測定値がある条件を満たさなければならないなどの拘束や制約を持たないで独立して行う測定 条件 ( 付き ) 測定 : 三角形の 3 つの内角の和のように, 個々の測定値間に満たすべき条件式が存在する場合の測定 方法による分類 直接測定 : 距離や角度などを機器を用いて直接行う測定 間接測定 : 求めるべき量を直接測定するのではなく,
More information例 e 指数関数的に減衰する信号を h( a < + a a すると, それらのラプラス変換は, H ( ) { e } e インパルス応答が h( a < ( ただし a >, U( ) { } となるシステムにステップ信号 ( y( のラプラス変換 Y () は, Y ( ) H ( ) X (
第 週ラプラス変換 教科書 p.34~ 目標ラプラス変換の定義と意味を理解する フーリエ変換や Z 変換と並ぶ 信号解析やシステム設計における重要なツール ラプラス変換は波動現象や電気回路など様々な分野で 微分方程式を解くために利用されてきた ラプラス変換を用いることで微分方程式は代数方程式に変換される また 工学上使われる主要な関数のラプラス変換は簡単な形の関数で表されるので これを ラプラス変換表
More information(Microsoft PowerPoint - \221\34613\211\361)
計算力学 ~ 第 回弾性問題の有限要素解析 (Ⅱ)~ 修士 年後期 ( 選択科目 ) 担当 : 岩佐貴史 講義の概要 全 5 講義. 計算力学概論, ガイダンス. 自然現象の数理モデル化. 行列 場とその演算. 数値計算法 (Ⅰ) 5. 数値計算法 (Ⅱ) 6. 初期値 境界値問題 (Ⅰ) 7. 初期値 境界値問題 (Ⅱ) 8. マトリックス変位法による構造解析 9. トラス構造の有限要素解析. 重み付き残差法と古典的近似解法.
More information19年度一次基礎科目計算問題略解
9 年度機械科目 ( 計算問題主体 ) 略解 基礎科目の解析の延長としてわかる範囲でトライしてみたものです Coprigh (c) 7 宮田明則技術士事務所 Coprigh (c) 7 宮田明則技術士事務所 Ⅳ- よってから は許容荷重として は直径をロ - プの断面積 Ⅳ- cr E E E I, から Ⅳ- Ⅳ- : q q q q q q q q q で絶対値が最大 で絶対値が最大モーメントはいずれも中央で最大となる
More information1.民営化
参考資料 最小二乗法 数学的性質 経済統計分析 3 年度秋学期 回帰分析と最小二乗法 被説明変数 の動きを説明変数 の動きで説明 = 回帰分析 説明変数がつ 単回帰 説明変数がつ以上 重回帰 被説明変数 従属変数 係数 定数項傾き 説明変数 独立変数 残差... で説明できる部分 説明できない部分 説明できない部分が小さくなるように回帰式の係数 を推定する有力な方法 = 最小二乗法 最小二乗法による回帰の考え方
More information(Microsoft Word - \216\221\227\277\201i\220\333\223\256\201jv2.doc)
宇宙工学基礎講義資料摂動 ( 松永担当分 ) ベクトル行列演算 ) 微分演算の定義 [ ] ) 微分公式 ( ベクトル記法と行列記法 ) E E ここで E は単位行列 チルダ演算は外積演算と等価の反対称行列を生成する演算 : ( ) ) 恒等演算式 : 次元列ベクトル ( ) ( ) ( ) ( ) ( ) E E ) ( ( )( ) ( )( ) ( )( ) ( ) ( ) ( ) ( )
More informationeq2:=m[g]*diff(x[g](t),t$2)=-s*sin(th eq3:=m[g]*diff(z[g](t),t$2)=m[g]*g-s* 負荷の座標は 以下の通りです eq4:=x[g](t)=x[k](t)+r*sin(theta(t)) eq5:=z[g](t)=r*cos(the
7. 制御設計の例 7.1 ローディングブリッジの制御装置 はじめに restart: ローディング ブリッジは 負荷をある地点から別の地点に運びます 台車の加速と減速は好ましくない振動を発生してしまいます そのため負荷はさらに安定し難くなり 時間もかかってしまいます 負荷がある地点から他の地点へ素早く移動し すみやかに安定するような制御装置を設計します 問題の定義 ローディング ブリッジのパラメータは以下の通りです
More information最速降下問題
最速降下問題 西山豊 533-8533 大阪市東淀川区大隅 --8 大阪経済大学経営情報学部 Tel: 06-638-43 E-Mail: nishiyama@osaka-ue.ac.jp. どの経路が速く到達するか図 のように傾斜面がある. 玉がAからBまで転がるとき最短時間であるのはどの曲線であろうか. 今仮に経路を直線, 次関数, サイクロイドとしよう. AとBを結ぶ最短経路は直線であるので直線がもっとも速く到達するかと思えるが意外と遅い.
More information2. エルニーニョ / ラニーニャ現象の日本への影響前記 1. で触れたように エルニーニョ / ラニーニャ現象は周辺の海洋 大気場と密接な関わりを持つ大規模な現象です そのため エルニーニョ / ラニーニャ現象は周辺の海流や大気の流れを通じたテレコネクション ( キーワード ) を経て日本へも影響
トピックス エルニーニョ / ラニーニャ現象 2009 年 7 月 10 日に気象庁から エルニーニョ現象が発生しているとの発表がありました 本 Express では 日本の気候にも大きな影響を与えるエルニーニョ / ラニーニャ現象 ( キーワード ) のメカニズムと日本への影響およびその予測可能性と温暖化について説明します 1. エルニーニョ / ラニーニャ現象とはエルニーニョ現象とは 太平洋赤道域の日付変更線付近から南米のペルー沿岸にかけての広い海域で
More information今後の予定 6/29 パターン形成第 11 回 7/6 データ解析第 12 回 7/13 群れ行動 ( 久保先生 ) 第 13 回 7/17 ( 金 ) 休講 7/20 まとめ第 14 回 7/27 休講?
今後の予定 6/29 パターン形成第 11 回 7/6 データ解析第 12 回 7/13 群れ行動 ( 久保先生 ) 第 13 回 7/17 ( 金 ) 休講 7/20 まとめ第 14 回 7/27 休講? 数理生物学演習 第 11 回パターン形成 本日の目標 2 次元配列 分子の拡散 反応拡散モデル チューリングパタン 拡散方程式 拡散方程式 u t = D 2 u 拡散が生じる分子などの挙動を記述する.
More information<4D F736F F D20959F93878DC48F4390B E F6E82CC895E93AE95FB92F68EAE82C68AB590AB97CD2E646F63>
ニュートンの運動方程式と慣性力 金沢工業大学基礎教育部 福島國雄 ねらい 力学の問題を解く場合には必ず物体に作用しているすべての力とその性質を知る必要がある. たとえば, 太陽の周りを公転している地球の運動を調べるには, 地球に作用している力を知る必要がある. その力は,( 太陽以外の天体の影響を無視し, 太陽に固定された座標系を慣性系と見なして ) 慣性系からみれば太陽が地球に及ぼす万有引力のみである.
More information3回
30 第 3 章ベクトルの微分法 キーワードベクトル ベクトルの演算 ゼロベクトル マイナスのベクトル ベクトルの定数倍 定数ベクトル 関数ベクトル ベクトルの成分表示 ベクトルの微分法 速度ベクトル 加速度ベクトル 極率 極率半径 ベクトルのスカラー積 ベクトル積 3.1 ベクトルの演算 1kgの質量や m 3 の体積などのように量で与えるものをスカラーと呼ぶ これに対し 北東の風 風速 m/sのように方向と大きさで与えるものをベクトルと呼ぶ
More information