Excelで定积分微分方程式-有限要素法やら.pptVIP

８．４ 偏微分方程式の数値解法 ［偏微分の復習（１）］ そこで お互いに相手の変数を定数とみなして微分する。 偏微分の例 を定数とみなして ［偏微分の復習（２）］ ［ベクトル解析の復習①］ イメージをつかむために ［ベクトル解析の復習③］ ［ベクトル解析の復習④］ ［ベクトル解析の復習⑤］ ［ベクトル解析の復習⑥］ ［ベクトル解析の復習⑦］ ８．４．１　数値解法における差分の方法（１）差分の種類 常微分の場合と同様，以下のような方法がある。 　　 　　①前進差分（forward difference） 　　②中心差分（centered difference） 　　③後退差分（backward difference） （２）前進差分 微小格子を考え，以下のように近似する。 （３）後退差分 微小格子を考え，以下のように近似する。 （４）中心差分 微小格子を考え，以下のように近似する。 （５）差分式使用上の注意点 どの方法を選択するかは， 　　①与えられた条件， 　　②問題の性格， 　　③モデル化のしやすさ 等に応じて決めること。 ８．４．２　物理現象と偏微分方程式（１）考え方 ①偏微分方程式は，物理現象を支配する法則等により体積無限小の極限で導出される。 ②有限な体積要素における法則により，基礎式が直接導出されることも多い。 　（例：コントロールボリューム法，直接差分法） （２）２階偏微分方程式の分類 ２階偏微分方程式の一般式 （３）楕円型のプログラム例 ［シートの定義］ VBAによる定義　①データ宣言 Private ポアソン As Boolean Private KX As Integer Private KY As Integer Private NN As Integer Private EPS As Double Private U() As Double Private ID1 As Integer Private ID2 As Integer Private F() As Double Private UMax As Double Private UMin As Double Private UDX As Double ②初期条件の設定と境界条件設定 Private Sub 初期条件(KX, KY, MX, MY, DX, DY) For j = 0 To KY For i = 0 To KX U(0, i, j) = 0: U(1, i, j) = 0: F(i, j) = 0 If ポアソン Then ‘　ポアソンの方程式のとき右辺設定 　　　　F(i, j) = 80 * (DX * (MX - i) + DY * (MY - j)) 　　End if Next Next End Sub Private Sub 境界条件(KX, KY, DX, DY) For k = 0 To 1 For i = 0 To KX: U(k, i, 0) = 0: U(k, i, KY) = 16: Next For j = 0 To KY: U(k, 0, j) = 0: U(k, KX, j) = 8: Next Next End Sub ③計算本体（その１） Private Sub 計算(KX, KY, NN, EPS) Dim MX As Integer: MX = KX + 1 Dim MY As Integer: MY = KY + 1 Dim DX As Double: DX = 1# / KX Dim DY As Double: DY = 1# / KY ReDim U(1, KX, KY), F(KX, KY) Dim F1 As Double: F1 = 1 / (DX * DX) Dim F2 As Double: F2 = 1 / (DY * DY) Dim F3 As Double: F3 = 0.5 / (F1 + F2) 初期条件 KX, KY, MX, MY, DX, DY 境界条件 KX, KY, DX, DY ④計算本体（その２） 収斂計算 ID1 = 1: ID2 = 0 For N = 0 To NN - 1 ID = ID1: ID1 = ID2: ID2 = ID: ER = 0 For j = 1 To KY - 1 For i = 1 To KX - 1 U(ID2, i, j) = F3 * (F1 * (U(ID1, i + 1, j) + U(ID1, i, j + 1)) + _ 　 F2