ひずみエネルギー計算

'///////////////////////////////////////////////

'     被積分関数(ひずみエネルギー計算)

'///////////////////////////////////////////////

Sub Energy(cE As CComplex)

    

    Set Gogh = Femtet.Gogh

    

    Dim Sgm() As CComplex

    Dim Eps() As CComplex

    

    ' 応力テンソルとひずみテンソル取得

    Gogh.Galileo.Tensor = GALILEO_STRESS_C

    Gogh.Galileo.Integral.GetTensor Sgm

    Gogh.Galileo.Tensor = GALILEO_STRAIN_C

    Gogh.Galileo.Integral.GetTensor Eps

    

    ' E = 0.5 * σ * ε （ひずみのせん断成分は2倍する必要がある）

    cE.Real = 0

    cE.Real = cE.Real + Sgm(0).Real * Eps(0).Real

    cE.Real = cE.Real + Sgm(1).Real * Eps(1).Real * 2

    cE.Real = cE.Real + Sgm(2).Real * Eps(2).Real * 2

    cE.Real = cE.Real + Sgm(4).Real * Eps(4).Real

    cE.Real = cE.Real + Sgm(5).Real * Eps(5).Real * 2

    cE.Real = cE.Real + Sgm(8).Real * Eps(8).Real

    cE.Real = 0.5 * cE.Real

    

End Sub