数値解析

Python

Python♪FEM:配列の次元数別の複製速度(list、NumPy配列)

list、NumPy配列の1~3次元配列について、参照渡し、浅いコピー、深いコピーなどの複製時間を計測しました。速度の差を感覚的に知っておくことは重要だと思います。私自身、予想外の計測結果になったものもあり、楽しむことができました。
Python

Python♪FEM:部材座標系の要素剛性マトリクス

バネの式は中学生に習いますが、これを行列で表現すると聞くと急に難しく感じてしまいます。でも、順を追って理解すれば難しいものではありません。F=k・xの式から要素剛性マトリクスを求めるまでの過程を説明したいと思います。
Python

Python♪FEM:2015年に高校数学から行列計算が消えた!

ご存じの通り(?)2015年から行列が高校数学から消えました。何考えているの!行列計算って情報処理の基本でしょ。中学生とPythonを勉強する記事を書いている私が言うのもなんですが、小中学生へのプログラムの早期教育よりも大切なのでは。
Python

Python♪FEM:行列の積や転置行列のコード記述のコツ

行列の積や転置行列についてNumpyの関数を使えば簡単ですが、関数を使わない方法も紹介します。PythonのFEMの書籍はまだ少ないため、CやFORTRANのコードを読むときに役に立ちます。熟語を覚えるような感覚でコードを覚えれば簡単です。
Python

Python♪FEM:平面トラスの有限要素法サンプルコード

平面トラスの有限要素法のプログラムをPythonで組んでみました。以前に自分がつくったコードを引っ張り出してきて、Pythonに移植しました。このくらいのボリュームのプログラムにじっくり取り組むと面白いですね。
Python

Python♪FEM:微小変形の仮定と人の感覚の違い

有限要素法のプログラムの自作にチャレンジするなら、「微小変形」は必ず理解しましょう。人間の感覚とずれているので知らなければ苦労します。解析では常識に近い初歩的な内容なので書籍やサイトなどでも説明が省略されていることが多いです。