|
カテゴリ:CAE
エクセルVBAで6面体8節点アイソパラメトリック要素の有限要素法ソルバーを作っている。
大苦戦の末にようやくできた。 さんざん見直したところ、積分点の座標の符号をまちがっていた。 以下のような2要素のモデルでテストをした。 ![]() ・拘束条件 節点1 xyz拘束 節点4 xz拘束 節点7 xy拘束 節点10 x拘束 ・荷重条件 節点3, 6, 9, 12にx方向に250N。 材料はヤング率210000MPa, ポアソン比0.3の鉄相当の弾性体。 積分点は8点、(ξ,η, ζ)=(±1/√3, ±1/√3, ±1/√3)としてガウス積分をしている。 変位を100倍にして表示した結果が以下。 ![]() 荷重点の変位0.00952381mmとなった。 ![]() にほんブログ村 お気に入りの記事を「いいね!」で応援しよう
[CAE] カテゴリの最新記事
|