PowerPointプレゼンテーション -金学会99春.ppt

計算材料科学ー量子力学と熱統計力学の基礎ー 工学研究科?材料工学専攻 助教授 西谷滋人 金属材料の特性 高強度，高剛性 高信頼性 耐熱性 気密性(宇宙機では特に重要) 加工性(成型，溶接，切削) 軽量(Al系，Ti系) 役立つTi 航空機部品 ロケット部品(H-IIA) 自転車，眼鏡のフレーム，時計，ラップトップパソコン H-IIロケットエンジンLE-7 H-IIロケット8号機エンジンLE-7の破損状況 http://www.nasda.go.jp/Home/Press/Press-p/199912/h28_991207_f05_j.jpg FTPインデューサ翼の破損状況 インデューサ欠損部の1つの破面(a2破面：写真2.1参照)が、疲労破壊であると判断された。疲労破面に特有に見られるストライエーション等の破面観察結果から、金属材料技術研究所では変動応力範囲と破断繰り返し数を以下のように算定した。 変動応力範囲：数百MPa 破断繰り返し数：105回程度http://yyy.tksc.nasda.go.jp/Home/Press/Press-p/200003/ h28_000317_f01_j.jpg 材料設計の指針 実験データ 知識データ 周期律表 状態図 周期律表 密度, 融点 Mg(1.74) 649K Al (2.70) 660K Ti (4.5) 1670K Fe (7.86) 1536K Cu (8.93) 1083K Au (18.86) 1063K 状態図(水) 状態図(炭素) -同素体- Tiに現れる(a) bcc, (b) hcp, (c) ω格子 状態図 (a) 水，(b) Ti 計算材料科学 量子力学と熱統計力学 量子力学 前期量子論 シュレディンガー方程式 簡単な例(ポテンシャルとエネルギー準位) 化学結合の起源 ばねモデル Tiの第一原理計算結果 前期量子論 黒体輻射, 光電効果, コンプトン散乱 波動ー粒子の二重性(duality)(de Broglie,1923) 電子線回折(1927 Davisson Germer) 物質波 Schr?dinger方程式(1926) Schr?dinger方程式 導出(各自トライ) 簡単な例(一次元調和振動子) ポテンシャル エネルギー準位 化学結合の起源(ポテンシャルと準位) Linear Combination of Atomic Orbitals原子軌道の線形結合 固有値，永年方程式の解 結合，反結合準位 結合準位の距離依存性 ばねモデル Tiの電子状態 Tiの結合エネルギーの体積依存 p=-dE/dV 熱統計力学 熱伝導，熱容量 分子動力学 モンテカルロシミュレーション アインシュタインモデル 比熱 熱伝導，熱容量 実験 ビタークラフト 分子動力学 Newtonの運動方程式 Lagrange力学Hamilton力学 Nose-Hoover力学 モンテカルロ(MC)シミュレーション 初期状態から少しのエネルギー差を仮定して，出現する可能性の高い状態を次々と生成させていく 位相空間 応用（焼きなまし法） 巡回セールスマン 問題 CPU配置，配線 位相空間 MCのアルゴリズム (1) 配置a を仮定しE (a )を求める． (2) a からすこし違った配置a +d?a を作る． (3) DE = E (a ) - E (a )を求める． (4) DE 0なら新たな配置を採用する． (5) DE 0なら新たな配置をexp(- DE/T)の確率で受け入れる． (6) 手順2以下を適当な回数繰り返す． アインシュタインモデル ポテンシャルの焼き直し Energy計算 C=dE/dT Einsteinモデルの計算結果 残された課題 bcc-Tiの安定性 相変態 組織形成 機械的性質 おまけ(数式処理ソフトについて) Dulong-Petitの法則 Diamond(θE=1320K) 物理工学総論，22/5/02　京都大学 H-IIロケットエンジンLE-7