位相最適化?

Size: px

Start display at page:

Download "位相最適化?"

へいぞうのあき
5 years ago
Views:

1 均質化設計法藤井大地 ( 東京大学 )

2 位相最適化?

3 従来の考え方境界形状を変化させて最適な形状位相を求める Γ t Ω b Γ D

4 境界形状を変化させる問題点解析が進むにつれて, 有限要素メッシュが異形になり, 再メッシュが必要になる位相が変化する問題への適応が難しい Γ Γ t t Ω b Ω b Γ D Γ D

5 領域の拡張と特性関数の導入 χ Ω ( x) = f 0 f x Ω x Ω 985 Murt χ Ω = χ Ω = 0 t Γ Ω Γ D b D

6 材料の ON/OFF 問題整数計画問題微分不可能問題

7 位相最適化の手法連続緩和法 ( 整数条件の緩和 ) 均質化設計法 ( ミクロ的な材料の ON/OFF) 密度法 ( 材料定数が密度のべき乗に比例 ) 離散最適化問題として解く方法遺伝的アルゴリズムを利用した方法セルオートマトンを利用した方法 ESO 法 (Evolutonry Structurl Optmzton)

8 均質化設計法 (HDM) マクロ構造ミクロ構造 x 2 b x y 2 y θ Unt cell 質量制約条件下で外力の仕事量の最小化

9 均質化設計法の目的関数従来の最適化問題応力の制約条件の下での重量の最小化均質化設計法重量制約条件の下でのコンプライアンス ( 外力の仕事量 ) の最小化さらに一般的には, 最小ポテンシャルエネルギーの最大化

10 なぜ応力の制約を用いない? 連続体ではいたるところで応力は発散! 境界角点集中荷重点支持点

11 なぜコンプライアンスの最小化? コンプライアンスを最小化することで構造物の剛性が最大化される応力は局所的に発散してもコンプライアンスは有限値となるコンプライアンスと応力との間には関係があり, コンプライアンスが最小化されれば, 最大応力も抑えられる感度解析が非常に容易になる

12 最小ポテンシャルエネルギーの最大化の方がより一般的外力が与えられる問題では, コンプライアンスの最小化変位が与えられる問題では, コンプライアンスの最大化外力が与えられる問題も, 変位が与えられる問題も, 最小ポテンシャルエネルギーの最大化問題となる

13 均質化設計法の最適化問題 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

14 均質化された弾性マトリックス (,, ) ( ) 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

15 有限要素法による離散化 { α 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

16 均質化弾性マトリックスの離散化 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 ユニットセル

17 計算の効率化均質化弾性マトリックス D H のデータベース化穴の大きさ,b に関してデータベースを作成するデータベースの値を利用して, その間の値は補間関数を用いて補間するユニットセルの角度 θ は, その要素の応力の主軸方向に一致させる

18 均質化設計法の解グレースケールチェッカーボード

19 フィルタリング法境界線長さの制約法 (Permeter control) チェッカーボードを防ぎ, シンプルな位相形状を求めることができるグレースケールを防ぐ処理を付加する必要がある重力制御法 (rvty control) チェッカーボードもグレースケールも同時に防ぎ, シンプルな位相形状を求めることができる

20 フィルタリングの効果

21 重力制御関数 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.75 = = 0.5 = 0 = = 2 = 3 = 4 = , 0.5, rvty control = 2.5 =.75 =

22 フィルタリングの導入法 = ( α ) = 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

23 最適化問題の解法逐次線形計画法 (SLP) 制約条件の扱いが易しくロバスト性が高い対称問題を解いても対称な位相が求まらない最適性規準法 (OC) 制約条件の数は限られるが収束が速い対称問題では対称な最適位相が得られる凸線形化法 (CONLIN) 多数の制約条件を扱え, 収束も速い対称問題では対称な最適位相が得られる

24 最適性規準法による解法ラグラジアンの定義 N L C b W G G = ( α ) = ( α ) Λ ( ) Λ ( α ) 歪みエネルギーの 2 倍 ( 平均コンフライアンス ) 質量制約条件 ( ) 重力制御関数制約条件 α= {,,,, b, b,, b } 2 N 2 N b : 各要素のミクロ構造ユニットセルの穴の大きさ穴の角度 θ は要素中心の応力の主軸方向 θ

25 最適性規準から得られる更新式ラグラジアン最小化の条件 ( α ) ( α ) ( α ) ( α ) 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)

26 更新式の導出法 ( ) ( ) 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 β Λ Λ + = α α β: 更新幅を制御するべき乗係数

27 変数制約条件の考慮 ( ) 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

28 設計変数のムーブリミット ( 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, ζ : 設計変数の変動幅を制約するムーブリミット

29 各ステップ内での更新 s s ( k+ ) ( k+ ) b b = = ( k) ( k + ) ( α ) ( k ) G Λ + Λ + C ( α ) G ( α ) ( k) ( k+ ) ( k ) Λ + Λ + b C ( α ) b 制約条件をアクティブにするただし, 感度係数は更新しない β β b

30 最適性規準法のパラメータ更新式のべき乗係数 β 設計変数のムーブリミット ζ 設計変数, ラグランジェ乗数の更新回数感度係数を更新する外側ループの繰返し数感度係数を更新しない内側ループの繰返し数 β = 0.25, ζ = 0. 外側ループの繰返し数 40 回内側ループの繰返し数 Λ に関しては00 回 Λ に関しては5 回

31 解析例 (MBB はり ) L /6 L /50 L

32 フィルタリングの効果 () 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

33 () 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

34 無償公開ソフト ( 連続緩和法 3 兄弟 ) 骨組構造の位相最適化ソフト (Otto) グランドストラクチャー法 2 次元連続体の位相最適化ソフト (Isler) 均質化設計法 3 次元連続体の位相最適化ソフト (Gud) 密度法

35 Otto の位相最適化手法 - グランドストラクチャー法 - 節点配置グランドストラクチャー最適位相

36 Isler の位相最適化手法 - 均質化設計法 - マクロ構造ミクロ構造 x 2 b x y 2 y θ Unt cell

37 Gud の位相最適化手法 - 密度法ー Element ρ =.0 ρ = 0 ρ = 0.5

38 Isler の解析例 (Snle lod)

39 Isler の解析例 (Snle Lod)

40 Isler の解析例 (Mult Lod)

41 Isler の解析例 (Mult Lod)

42 Isler の解析例 (Dstrbuted lod)

43 Isler の解析例 (Dstrbuted lod)

44 Isler の解析例 (Dstrbuted lod)

45 Isler の解析例 (Snle lod) 周辺単純支持中央集中荷重

46 Isler の解析例 (Dstrbuted lod) 周辺単純支持鉛直等分布荷重

47 まとめ公開ソフトは, 下記のホームページからダウンロード可能実務での利用を考えておられる方は下記ソフトをご購入下さい MSC Nstrn Opt-shpe Alte Opt-struct

構造物のトポロジー最適化について東京大学藤井大地

構造物のトポロジー最適化について東京大学藤井大地構造物のトポロジー最適化について東京大学藤井大地講演内容力学教育, 設計教育に用いることを目的として開発した骨組, 連続体の位相最適化ツールについて位相最適化における問題点とその解決法について位相最適化に関するその他の研究について本日の発表資料は下記に掲載 http://www.nasl.t.u-tokyo.ac.jp/~dfuj/homepage.htm 書き物の中にあります力学系科目の衰退