戻る

5 代数方程式のマトリックス表示


式(2-20)の連立方程式は一般に以下のようなマトリックス表示で表すことができる。

         (2-21)

ここで、は時刻における節点での未知関数ベクトルであり、各マトリックスを

式(2-22)と式(2-24)のように求めると、以下に説明するように式(2-20)と式(2-21)とは

等価になる。 まずは次式で表される。

                             (2-22)

ここで 各のマトリックスは要素マトリックスとよばれ、またVe:要素()の領域(の範囲),

,:要素(e)での形状関数 である。 またM:総要素数  であり、狽ヘ各要素マトリックスを

足し加えて全体マトリックスを作成することを意味している。

式(2-22)の要素マトリックスでは、たとえば一行一列目の成分は図2-3と式(2-13)より次式で表される。


結局、 を考慮して上式を積分すると次式が得られる。

                                     (2-23)

同様に、 マトリックスの他の成分や他のマトリックスについての定式化の結果は次式で表される。


 ,


    ,


                          (2-24)
 

また各要素マトリックスの成分を次のように足し加えて全体マトリックスを作成すると

式(2-21)と式(2-20)とは等価になる。たとえば式(2-24)の中の要素(図2-3、

要素(図2-3、 に対する各マトリックスの成分は図2-4のように表される。
 

                                         

            

  (a) Element             (b) Element

       Fig. 2-4 Element matrices and corresponding nodes


ここで、要素( の一行目の成分は当該要素の一番目の節点()に関する方程式、

すなわち  の()行目に寄与する。また二行目の成分は当該要素の二番目の節点()に

関する方程式、すなわち行目に寄与する。 一方、二行一列目の成分(h/6)は の()列目に、

二行二列目の成分(h/3)は()列目に加える。  次に要素()の要素マトリックスに

おける一行一列目の成分(h/3)を行()列目に、一行二列目の成分(h/6)を行()列目に

加える。二行目の成分は の()行目に寄与し、それ以外の要素からの第()行目への寄与は

ゼロである。その結果、の第行目は となる。


またマトリックス についても同様に作成し、両辺に  をかけて無次元係数を

導入すると行目の代数方程式として式(2-20)を得る。以上の節点に関する議論をすべての内部節点

( 節点に対して節点 ()と()が存在 )に適用し境界条件を導入すると未知関数値の個数と

同数の連立方程式が得られるため、これを解いて各節点値を求めることができる。


上述のガレルキン法は有限要素法の最も基本的な一種であるが、流体解析においては

一般に数値振動が発生するため、前述の各種の安定化解法が適用されている(1),(2)。  

また移流項に関するマトリックス は非対称となり、連立方程式の解法では対称マトリックス

のみの場合に比べてより多くの計算時間と計算容量を必要とする。これらの問題を解決するため

著者らは次章で述べる変形ガレルキン法を提案した。