均質化設計法 藤井大地 ( 東京大学 )
位相最適化?
従来の考え方 境界形状を変化させて最適な形状 位相を求める Γ t Ω b Γ D
境界形状を変化させる問題点 解析が進むにつれて, 有限要素メッシュが異形になり, 再メッシュが必要になる 位相が変化する問題への適応が難しい Γ Γ t t Ω b Ω b Γ D Γ D
領域の拡張と特性関数の導入 χ Ω ( x) = f 0 f x Ω x Ω 985 Murt χ Ω = χ Ω = 0 t Γ Ω Γ D b D
材料の ON/OFF 問題 整数計画問題 微分不可能問題
位相最適化の手法 連続緩和法 ( 整数条件の緩和 ) 均質化設計法 ( ミクロ的な材料の ON/OFF) 密度法 ( 材料定数が密度のべき乗に比例 ) 離散最適化問題として解く方法 遺伝的アルゴリズムを利用した方法セル オートマトンを利用した方法 ESO 法 (Evolutonry Structurl Optmzton)
均質化設計法 (HDM) マクロ構造 ミクロ構造 x 2 b x y 2 y θ Unt cell 質量制約条件下で外力の仕事量の最小化
均質化設計法の目的関数 従来の最適化問題 応力の制約条件の下での重量の最小化 均質化設計法 重量制約条件の下でのコンプライアンス ( 外力の仕事量 ) の最小化さらに一般的には, 最小ポテンシャルエネルギーの最大化
なぜ応力の制約を用いない? 連続体ではいたるところで応力は発散! 境界角点 集中荷重点 支持点
なぜコンプライアンスの最小化? コンプライアンスを最小化することで構造物の剛性が最大化される 応力は局所的に発散してもコンプライアンスは有限値となる コンプライアンスと応力との間には関係があり, コンプライアンスが最小化されれば, 最大応力も抑えられる 感度解析が非常に容易になる
最小ポテンシャルエネルギーの 最大化の方がより一般的 外力が与えられる問題では, コンプライアンスの最小化 変位が与えられる問題では, コンプライアンスの最大化 外力が与えられる問題も, 変位が与えられる問題も, 最小ポテンシャルエネルギーの最大化問題となる
均質化設計法の最適化問題 H mx mn ( 0) ( 0) ( )( 0) dd 0 dγ 0 Π v = 2 xv D α xv L v K D v t α Γ { α 0, 0 b, ρdd m } S, { 0 0 onγ D} L= K = v v = D x 2 D Γ Γ D D x Globl mcro coordntes y 2 b Unt cell y θ Locl mcro coordntes
均質化された弾性マトリックス (,, ) ( ) H bθ = θ ( b, ) ( θ) H D R D R E ( b, ) = mn ( χ ) ( χ ) H D I y D I y χ Y perodc dy x 2 D Γ D y 2 b y Unt cell θ R ( ) 2 2 cos θ sn θ cosθsnθ = sn2θ sn2θ cos2θ 2 2 θ sn θ cos θ cosθsnθ χ Γ D x Globl mcro coordntes Locl mcro coordntes
有限要素法による離散化 { α b ρdd m } S mn dkd, L= 0, 0, α L 0 D 要素剛性マトリックス x 2 K H 0e e B0 D B0dDee D = e x b y 2 θ y
均質化弾性マトリックスの離散化 H D = DdY X KX Y K X = q e e e K q e = B DB Y e = Y e e B DdY dy χ y χ = BX e = NX e e e Χ j(4) Χ j (3) Χ j() Χ j (2) =,2 j=,2,3 Fnte element ユニットセル
計算の効率化 均質化弾性マトリックス D H のデータベース化 穴の大きさ,b に関してデータベースを作成する データベースの値を利用して, その間の値は補間関数を用いて補間する ユニットセルの角度 θ は, その要素の応力の主軸方向に一致させる
均質化設計法の解 グレースケール チェッカーボード
フィルタリング法 境界線長さの制約法 (Permeter control) チェッカーボードを防ぎ, シンプルな位相形状を求めることができる グレースケールを防ぐ処理を付加する必要がある 重力制御法 (rvty control) チェッカーボードもグレースケールも同時に防ぎ, シンプルな位相形状を求めることができる
フィルタリングの効果
重力制御関数 N N N m N ( ρ ρ ρ ρ ) G = m = + m j j = = = j= = ρ = b, ρ = ρ, 0< G m = 4 = 4 = 3 = 2 = = 0 = 0.25 = 0.5 = 0.75 = = = = = = = 0.875 = 0.75 = 0.625 = 0.5 = 0 = = 2 = 3 = 4 = 3.25 0, 0.5, rvty control = 2.5 =.75 =
フィルタリングの導入法 = ( α ) = UKU α = { } mn C α, 2,, N, b, b2,, b subject to : N ( ) W = b W, 0 α, =,,2N N G G 制約条件として加える ただし, G m N ( α) = ρρ + ( ρ)( ρ ) j j = j= = N m
最適化問題の解法 逐次線形計画法 (SLP) 制約条件の扱いが易しくロバスト性が高い対称問題を解いても対称な位相が求まらない 最適性規準法 (OC) 制約条件の数は限られるが収束が速い対称問題では対称な最適位相が得られる 凸線形化法 (CONLIN) 多数の制約条件を扱え, 収束も速い対称問題では対称な最適位相が得られる
最適性規準法による解法 ラグラジアンの定義 N L C b W G G = ( α ) = ( α ) Λ ( ) Λ ( α ) 歪みエネルギーの 2 倍 ( 平均コンフ ライアンス ) 質量制約条件 ( ) 重力制御関数制約条件 α= {,,,, b, b,, b } 2 N 2 N b : 各要素のミクロ構造ユニットセルの穴の大きさ 穴の角度 θ は要素中心の応力の主軸方向 θ
最適性規準から得られる更新式 ラグラジアン最小化の条件 ( α ) ( α ) ( α ) ( α ) N N N C G C G δl( α ) = + bλ+ Λδ + + Λ+ Λδb ( b ) WδΛ ( G G( α )) δλ = 0 = = b b = G ( α ) ( k) ( k) b Λ + Λ ( k+ ) ( k) = C b ( α ) ( k+ ) β G ( α ) ( k) ( k) Λ + Λ b ( k) = b C ( α ) b Λ ( k+ ) ( k) ( k) ( k) ( ) β N = W = b Λ β ( k+ ) Λ G = G( α ) β Λ ( k)
更新式の導出法 ( ) ( ) 0 C G b Λ Λ + + = α α α α α α α α ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) k k k k k k k k k G b C β Λ Λ + + = α α ( ) ( ) ( ) ( ) k k G b C Λ Λ + = α α ( ) ( ) ( ) ( ) k k G b C β Λ Λ + = α α β: 更新幅を制御するべき乗係数
変数制約条件の考慮 ( ) 0, 0 b =,, N, Λ 0, Λ 0 b Λ Λ { } ( k + ) = { s } mn mx 0,, { } ( k + ) = { sb } mn mx 0,, N = mn 0, ms = ( k+ ) ( k) ( k) ( k) ( ) ( k+ ) β G ( k) = mn 0, Λ G( α ) b β Λ ただし, s s b b = = ( k) Λ G + C ( α ) ( k) ( α ) G ( α ) ( α ) Λ ( k) ( k) Λ + Λ b C b β β b
設計変数のムーブリミット ( k+ ) ( k) s b {( ζ ) } mx,0 = mn +, ( k+ ) ( k) sb {( ζ ) } {( ζ ) b } mx,0 = mn +, {( ζ ) b } k {( ζ ) } {( ζ ) } ( + ζ ) ( k) ( k) {( + ζ ) } s ( k) ( ) f s mx,0 { } ( k) ( k) ( k) f mx,0 s mn, f mn, k {( ζ ) b } {( ζ ) } ( + ζ ) ( k) ( k) {( + ζ ) b } sb ( k) ( ) b f s mx,0 { } ( k) ( k) ( k) b f mx b,0 s mn b, f mn, ζ : 設計変数の変動幅を制約するムーブリミット
各ステップ内での更新 s s ( k+ ) ( k+ ) b b = = ( k) ( k + ) ( α ) ( k ) G Λ + Λ + C ( α ) G ( α ) ( k) ( k+ ) ( k ) Λ + Λ + b C ( α ) b 制約条件をアクティブにする ただし, 感度係数は更新しない β β b
最適性規準法のパラメータ 更新式のべき乗係数 β 設計変数のムーブリミット ζ 設計変数, ラグランジェ乗数の更新回数 感度係数を更新する外側ループの繰返し数感度係数を更新しない内側ループの繰返し数 β = 0.25, ζ = 0. 外側ループの繰返し数 40 回 内側ループの繰返し数 Λ に関しては00 回 Λ に関しては5 回
解析例 (MBB はり ) L /6 L /50 L
フィルタリングの効果 () G = 0, G = 0.73 (c) G = 0.85, G = 0.89, C/ C =.08 0 (b) G = 0.8, G = 0.87, C/ C =.02 0 (d) G = 0.90, G = 0.90, C/ C =.8 0
() m = 0.3, G = 0.85, G = 0.85, C/ C =.05 S 0 (b) m = 0.4, G = 0.85, G = 0.87, C/ C =.08 S 0 (c) m = 0.5, G = 0.85, G = 0.89, C/ C =.08 S 0 (d) m = 0.6, G = 0.8, G = 0.89, C/ C =.03 S 0
無償公開ソフト ( 連続緩和法 3 兄弟 ) 骨組構造の位相最適化ソフト (Otto) グランドストラクチャー法 2 次元連続体の位相最適化ソフト (Isler) 均質化設計法 3 次元連続体の位相最適化ソフト (Gud) 密度法
Otto の位相最適化手法 - グランドストラクチャー法 - 節点配置 グランドストラクチャー 最適位相
Isler の位相最適化手法 - 均質化設計法 - マクロ構造 ミクロ構造 x 2 b x y 2 y θ Unt cell
Gud の位相最適化手法 - 密度法ー Element ρ =.0 ρ = 0 ρ = 0.5
Isler の解析例 (Snle lod)
Isler の解析例 (Snle Lod)
Isler の解析例 (Mult Lod)
Isler の解析例 (Mult Lod)
Isler の解析例 (Dstrbuted lod)
Isler の解析例 (Dstrbuted lod)
Isler の解析例 (Dstrbuted lod)
Isler の解析例 (Snle lod) 周辺単純支持 中央集中荷重
Isler の解析例 (Dstrbuted lod) 周辺単純支持 鉛直等分布荷重
まとめ 公開ソフトは, 下記のホームページからダウンロード可能 http://www.nsl.t.u-tokyo.c.jp/~dfuj/homepe.htm 実務での利用を考えておられる方は下記ソフトをご購入下さい MSC Nstrn Opt-shpe Alte Opt-struct