ryujimiyaの日記

C#を使って数値解析したい

IvyFEM開発日誌(4).NET有限要素法CAEライブラリIvyFEM Navier-Stokes方程式など各種問題に対応(バージョン0.0.0.5)

IvyFEM.dll バージョン0.0.0.5をリリースしました。(まだデバッグ版のみです)

 GitHub - ryujimiya/IvyFEM: .NET向け有限要素法(FEM)を用いたCAEライブラリ

 

いまできる解析は次の通りです。太字が今回新しく対応したものです。

    弾性体力学:線形弾性体
          超弾性体
          Saint Venant Kirchhoff
          Mooney-Rivlin (非圧縮、微圧縮)
          Ogden (非圧縮、微圧縮)
          多点拘束(Multipoint Constraint, MPC)(直線)
          剛体との接触(直線、円)
          弾性体二体接触
    流体力学 :Navier-Stokesの方程式
    電磁気学 :H面導波管の伝達問題
    各種方程式:Poisson方程式
          熱拡散方程式
          移流拡散方程式
          Helmholtz方程式

 

(1)Navier-Stokesの方程式

(1-1)空洞問題

上の境界に平行な速度の流体を付加したときの対流を計算しました。

定常状態、時間領域の両方に対応しています。

f:id:ryujimiya:20190428094803j:plain



 

(1-2)バックステップ

左から流体を流入させたときの分布を計算しました。

右の境界の手前だけ粘性係数を大きくしています。

f:id:ryujimiya:20190428094837j:plain



 

(2)Poissonの方程式

中央の円領域を励振したときのポテンシャルの分布を計算しました。

左下側にポテンシャル0の壁を設けています。

f:id:ryujimiya:20190408221503j:plain

(3)Helmholtzの方程式

音などの周波数領域の解析です。中央に波源を置いたときの分布を計算しました。

開放系の解析なので周囲の境界条件には二次の吸収境界条件を用いました。

f:id:ryujimiya:20190408221935j:plain

(4)熱拡散方程式

上の境界を低温度、下の境界を高温度とし、右上の熱源から熱を放射したときの時間変化を計算しました。

f:id:ryujimiya:20190408222500j:plain

 

(5)移流拡散方程式

反時計回りの速度で対流しているところに熱源を置いたときの対流を計算しました。

定常状態、時間領域の両方に対応しました。

f:id:ryujimiya:20190408222735j:plain