トポロジー最適化に挑戦する：フリーFEMソフトとExcelマクロで形状最適化（4）（2/4 ページ）

[高橋良一／ＲＴデザインラボ 代表，MONOist]

密度法と有限要素法

　密度法を採用したときの有限要素法プログラムの活用方法について解説します。話を簡単にするために、平面応力問題とし、図7のような2次元三角形要素の剛性マトリクスで説明します。｛f｝は荷重ベクトル、｛u｝は変位ベクトル、［k］は剛性マトリクスです。

図7　2次元三角形要素［クリックで拡大］

　2次元問題の場合、剛性マトリクス［k］は式5で計算できます。なお、式5の導出は参考文献［4］のWebサイトが参考になります。hは要素の板厚です。

式5

　各節点（i，j，k）に変位があって、その変位量が一致していない場合は要素が引き延ばされているので“ひずみ”が発生します。よって、ひずみは変位ベクトル｛u｝で決まることになり、その関係は式6、式7で表されます。｛ε｝はひずみ、｛u｝は変位ベクトルです。［B］は変位−ひずみマトリクスといい、3行6列のマトリクスです。変位−ひずみマトリクス［B］は要素の形状で決まります。

式6

式7

　ひずみと応力は、図8に示すように「x軸を法線とする面に、x軸方向に作用するひずみをε_xx、応力をσ_xxと表記し、x軸を法線とする面に、y軸方向に作用するひずみをε_xy、応力をσ_xyと表記する」ことが正確な表現です。添え字が「x面を法線とする面」と「x軸方向に作用する応力」というように2つあるので、ひずみと応力はテンソル量として2行2列のマトリクスで表現されます。しかし、有限要素法ソフトではこれをベクトル量として表現してプログラミングされています。よって式7において、ひずみをベクトル量として表記しました。

図8　微小四角形の応力とひずみ［クリックで拡大］

　また、式8の関係があって、2種類あったせん断ひずみとせん断応力は1つの量として表されます。

式8

　材料力学で述べられているせん断ひずみγ_xyと図8で述べたせん断ひずみε_xyは少し違っていて、式9にその関係を示しておきます。

式9

　ひずみの表記を式10のように変更すると、ひずみは式7のようになります。

式10

　応力の表記も式11に変更し、応力もベクトル量として表記します。

式11

式12

　式5の［D］は応力−ひずみマトリクスで、3行3列のマトリクスです。応力とひずみの関係は次式で表されます。

式13