トポロジー最適化のチェッカーフラグ対策：フリーFEMソフトとExcelマクロで形状最適化（9）（6/6 ページ）

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

トポロジー最適化の計算手順の続き

　前述したD₁とD₂の考察から、次のρ_iは次式で決定します。ρ_iの変化量が大きいと繰り返し計算が発散するので、ρ_iの変化量にリミッターをかけています。連載第6回と同じです。

式24

式25

　もう1つリミッターをかけます。ρ_iは1以上にはなれず、かつマイナス値にもなれません。次式となります。

式26

式27

　この辺の手順は「Excel」のマクロプログラムには組み込んでおらず、シートのセルに式として設定しています。任意に数値やアルゴリズムを変えられると思います。式17の第2項D₂の分母は前述したことと反対に作用しますが、Dが1に近づけばよいので少し目をつむりましょう。

　λは計算で求まりますが、λ_gは簡単に求まりません。なぜなら、ρ_iの変化にリミッターがかかっているので、そう簡単に式9は成立しません。式9の最初の値はプラスの値なので、世代を繰り返すことで式9の値をゼロに近づけます。ここでのλ_gは、式9を成立させながら目標関数を最小化するための調整定数のようなものです。

　式23を再掲します。最適化の最中で、式9がゼロから程遠かったら、慌てて式9がゼロに近づくようにしなければならないため、後述する理由からλ_gを大きくする必要があります。λ_gを大きくさせる度合いがampで、これはフィードバック系のゲインに相当しますね。

式23の再掲

　式23の第2項は、式9の第2項（重力関数）であって、チェッカーフラグの度合いを数値化したものでした。チェッカーフラグになっていれば、式23の第2項は小さくなり、|λ_g|は大きくなります。また、初期状態は全てのρ_iが0.375［-］などになっていて、のっぺりとしたグレーなので、この場合も|λ_g|は大きくなります。

　式17の第1項D₁はプラス値、第2項D₂はマイナス値となり、正負が拮抗しています。|λ_g|が大きければ、|D₂|が大きくなって、次世代のρ_iを決めるとき、D₂の意向が反映されます。つまり、チェッカーフラグや、のっぺりとしたグレー状態を避けるように次世代のρ_iが変化します。

　式23の第2項が目標値G_targetに近づけば、λ_gがゼロに近づき、第2項D₂もゼロに近づきます。この結果、次の世代のρ_iを決めるときに、D₁の意向が反映されます。つまり、各要素の蓄えられる仕事量が小さくなるように、次世代のρ_iが変化します。以上の操作で世代を繰り返すことによって、式7が最小になるρ_iが求まり、同時に式9が成立することになります。以上が図3の*B、*Cの説明です。少し長かったですね。

　最後に、図3の*Dを説明します。「各ρ_iが0ないし1になったか、かつ、h＜0になったか」の判定と計算終了の判断方法を述べます。重力関数は大きければ大きいほどよいので、“目標値G_targetよりも大きくなってもよし”とします。よって、h＝0ではなくh＜0としました。この判断は人間で行いましょう。

---

　以上が、チェッカーフラグ対策付きトポロジー最適化のアルゴリズムです。この部分は文献を参考にせず、筆者が自分で考えましたので改善の余地があると思います。

　次回は、チェッカーフラグ対策付きトポロジー最適化を実施した例を紹介します。チェッカーフラグ対策のあり／なしで最適化形状が変わりますので、それらを比較しましょう。　（次回へ続く）

⇒ 連載バックナンバーはこちら

Profile

高橋 良一（たかはし りょういち）
ＲＴデザインラボ 代表

1961年生まれ。技術士（機械部門）、計算力学技術者 上級アナリスト、米MIT Francis Bitter Magnet Laboratory 元研究員。

構造・熱流体系のCAE専門家と機械設計者の両面を持つエンジニア。約40年間、大手電機メーカーにて医用画像診断装置（MRI装置）の電磁振動・騒音の解析、測定、低減設計、二次電池製造ラインの静音化、液晶パネル製造装置の設計、CTスキャナー用X線発生管の設計、超音波溶接機の振動解析と疲労寿命予測、超電導磁石の電磁振動に対する疲労強度評価、メカトロニクス機器の数値シミュレーションの実用化などに従事。現在ＲＴデザインラボにて、受託CAE解析、設計者解析の導入コンサルティングを手掛けている。⇒ ＲＴデザインラボ