以下のサイトを参考にしながら、2次元ソリッド四角形4節点アイソパラメトリック要素のFEMコードをエクセルVBAで組むことに挑戦中。
http://www.geocities.jp/moridesignoffice/FEM-sum-4.html
前に、2次元ソリッド三角形要素は作ったことがあるので、これを使えばすぐにできるかと思いきや、どっかで間違ってるらしく苦戦中。
http://plaza.rakuten.co.jp/takupin/diary/201209300000/
プログラムではなくて、エクセルスプレッドシートで一通り、ヤコビアン、Bマトリクス、Dマトリクス、要素合成マトリクス、全体剛性マトリクス、拘束条件と荷重条件で列と行をソート、逆行列計算、荷重が既知の拘束されていない節点変位計算、拘束節点の荷重計算を確認してみたけれど、やっぱり間違っている。
すっかりドツボにはまってしまった。。。
にほんブログ村