基于并行策略的熱固耦合多材料柔順機(jī)構(gòu)設(shè)計(jì)*

張憲民 胡凱 王念峰 郭浩

(華南理工大學(xué) 機(jī)械與汽車工程學(xué)院∥廣東省精密裝備與制造技術(shù)重點(diǎn)實(shí)驗(yàn)室, 廣東 廣州 510640)

基于并行策略的熱固耦合多材料柔順機(jī)構(gòu)設(shè)計(jì)*

張憲民 胡凱 王念峰?郭浩

(華南理工大學(xué) 機(jī)械與汽車工程學(xué)院∥廣東省精密裝備與制造技術(shù)重點(diǎn)實(shí)驗(yàn)室, 廣東 廣州 510640)

熱效應(yīng)在結(jié)構(gòu)設(shè)計(jì)領(lǐng)域受到廣泛的重視，文中利用拓?fù)鋬?yōu)化的方法，對(duì)熱驅(qū)動(dòng)柔順機(jī)構(gòu)設(shè)計(jì)進(jìn)行了研究，提出了相應(yīng)的基于并行策略的求解模型.該模型的核心是將一個(gè)復(fù)雜的多材料熱固耦合問(wèn)題離散成為單材料熱固耦合子問(wèn)題，然后對(duì)這些子問(wèn)題并行求解，并根據(jù)整體目標(biāo)的需要，對(duì)所有子問(wèn)題的解進(jìn)行調(diào)整以得到原始問(wèn)題的解.對(duì)所有的單材料子問(wèn)題，基于SIMP(固體各向同性材料懲罰模型)材料插值模型，建立了熱固耦合拓?fù)鋬?yōu)化求解模型，然后利用順序耦合方法進(jìn)行熱-結(jié)構(gòu)多物理場(chǎng)耦合分析，獨(dú)立并行求解.該模型在理論和實(shí)踐上都比較簡(jiǎn)單，可以處理任意多種材料，并且可以避免零碎的拓?fù)浣Y(jié)構(gòu)，因而有利于制造.文中還通過(guò)數(shù)值算例說(shuō)明了此方法的有效性.

拓?fù)鋬?yōu)化；柔順機(jī)構(gòu)；熱固耦合；敏度分析

拓?fù)鋬?yōu)化一直是結(jié)構(gòu)優(yōu)化領(lǐng)域的熱點(diǎn)問(wèn)題，采用拓?fù)鋬?yōu)化的方法設(shè)計(jì)柔順機(jī)構(gòu)具有較大的設(shè)計(jì)自由度，因而受到人們的重視[1- 5].常用的拓?fù)鋬?yōu)化設(shè)計(jì)方法主要有變密度法[6- 8]、基礎(chǔ)結(jié)構(gòu)法[9- 10]、水平集法[11- 12]等.

多場(chǎng)耦合是指在一個(gè)系統(tǒng)中，熱、電、磁和結(jié)構(gòu)場(chǎng)等相互作用和影響的現(xiàn)象，多場(chǎng)耦合下的結(jié)構(gòu)拓?fù)潢P(guān)系比較復(fù)雜，與直觀感覺和經(jīng)驗(yàn)相差較大，需要一種能夠綜合考慮多場(chǎng)關(guān)系來(lái)設(shè)計(jì)結(jié)構(gòu)布局形式的方法，因而涉及耦合場(chǎng)結(jié)構(gòu)分析的拓?fù)鋬?yōu)化應(yīng)用研究成為一個(gè)重要的發(fā)展方向[13].Sigmund[14- 15]研究了電熱場(chǎng)耦合下的雙材料插值模型問(wèn)題和驅(qū)動(dòng)器的拓?fù)鋬?yōu)化建模、設(shè)計(jì)和制造問(wèn)題；李冬梅等[16]將可靠性分析用于熱固耦合拓?fù)鋬?yōu)化中，提出了一種基于概率約束的可靠性拓?fù)鋬?yōu)化設(shè)計(jì)方法；Takezawa等[17- 18]研究了拓?fù)鋬?yōu)化中結(jié)構(gòu)強(qiáng)度和熱傳導(dǎo)約束問(wèn)題，在優(yōu)化的后續(xù)階段采用PFM法更新密度函數(shù)，從而避免中間密度單元，得到清晰的結(jié)構(gòu)邊界；Koga等[19]采用拓?fù)鋬?yōu)化的方法，在穩(wěn)態(tài)熱模型中綜合考慮傳導(dǎo)和強(qiáng)制熱對(duì)流的影響，設(shè)計(jì)了完整的循環(huán)對(duì)流性能最佳的散熱器.然而迄今為止，基于多材料的多場(chǎng)耦合問(wèn)題的研究相對(duì)較少，尚有待深入研究.

在給定約束條件和負(fù)載下，以一定的設(shè)計(jì)目標(biāo)，尋求多種材料在設(shè)計(jì)域內(nèi)的最優(yōu)分布，即為多材料的拓?fù)鋬?yōu)化.相比于單材料機(jī)構(gòu)，多材料機(jī)構(gòu)能夠充分發(fā)揮各種材料的優(yōu)良屬性，擁有更大的設(shè)計(jì)自由度，受到了許多的關(guān)注與研究.Guo等[20]提出基于水平集的變分一致策略，研究了應(yīng)力約束下多材料拓?fù)鋬?yōu)化問(wèn)題；Wang等[21]提出了基于水平集的多材料拓?fù)鋬?yōu)化模型，每一相材料由不同的水平集函數(shù)組合而成，優(yōu)化得到光滑的邊界和清晰的表面結(jié)構(gòu)；Tavakoli[22]基于多相體積約束的金茲堡-朗道能量泛函，提出一種解決多材料拓?fù)鋬?yōu)化問(wèn)題的模型，采用分步投影的最速下降法對(duì)設(shè)計(jì)變量進(jìn)行更新.上述情形中沒有涉及多材料與多場(chǎng)耦合組合拓?fù)鋬?yōu)化的問(wèn)題，同時(shí)對(duì)機(jī)構(gòu)及耦合驅(qū)動(dòng)也尚未有深入的研究.

文中重點(diǎn)研究了多材料柔順機(jī)構(gòu)的熱固耦合拓?fù)鋬?yōu)化問(wèn)題.首先，將原始的多材料多場(chǎng)耦合優(yōu)化問(wèn)題按照材料的種類離散成單材料多場(chǎng)耦合優(yōu)化子問(wèn)題；其次，對(duì)所有的單材料子問(wèn)題，基于SIMP(固體各向同性材料懲罰模型)建立了熱固耦合拓?fù)鋬?yōu)化模型，利用順序耦合方法進(jìn)行熱-結(jié)構(gòu)多物理場(chǎng)耦合分析，獨(dú)立并行求解；接著調(diào)整所有子問(wèn)題的解以得到符合原始的多場(chǎng)耦合問(wèn)題的解；最后通過(guò)算例對(duì)所提方法的有效性進(jìn)行了驗(yàn)證.

1 材料插值模型及并行優(yōu)化策略

1.1 基于SIMP模型的變密度法

各向同性材料的SIMP法又稱變密度法，其基本思想為：設(shè)計(jì)域被離散成為微結(jié)構(gòu)單元，設(shè)計(jì)變量定義為離散單元的相對(duì)密度，可取0到1之間任意值.單元的材料屬性，比如楊氏模量與單元相對(duì)密度呈指數(shù)變化關(guān)系.設(shè)E0為材料原始的彈性模量，E是單元相對(duì)密度值等于xe時(shí)的彈性模量值，則存在關(guān)系式：

E=(xe)pE0

(1)

其中，p為懲罰因子，其目的在于對(duì)中間密度單元項(xiàng)進(jìn)行懲罰，以盡量減少結(jié)構(gòu)中間密度單元的數(shù)目，使結(jié)構(gòu)單元密度盡可能靠近1或0兩極.通常p取值大于3[14].

1.2 并行優(yōu)化策略

由于原始的問(wèn)題為多材料多場(chǎng)耦合拓?fù)鋬?yōu)化問(wèn)題，當(dāng)按照材料的種類離散為單材料子問(wèn)題后，所有的子問(wèn)題都成為單材料多場(chǎng)耦合問(wèn)題.對(duì)于單材料多場(chǎng)耦合拓?fù)鋬?yōu)化問(wèn)題，國(guó)內(nèi)外學(xué)者已經(jīng)做了許多研究[23- 24].文中提出一種并行優(yōu)化策略，作為一種用于同時(shí)解決多材料的多場(chǎng)耦合問(wèn)題的新方法.

所謂并行優(yōu)化，就是對(duì)所有的單材料子問(wèn)題獨(dú)立并行求解.不難看出，這種優(yōu)化策略的核心與難點(diǎn)在于后續(xù)如何統(tǒng)籌所有子問(wèn)題的解，以得到原始優(yōu)化問(wèn)題的解.

(2)

2 熱固耦合理論分析與建模

2.1 耦合場(chǎng)有限元分析

耦合場(chǎng)分析是指在有限元分析的過(guò)程中綜合考慮了多種物理場(chǎng)的交叉作用和相互影響.常見的耦合分析方法有直接耦合與間接耦合方法[13].文中采用間接耦合中的順序耦合分析方法，首先進(jìn)行熱場(chǎng)有限元分析，計(jì)算得到熱場(chǎng)作用下的溫度分布矩陣，然后將熱分析得到的節(jié)點(diǎn)溫度用于計(jì)算熱載荷，作為體力施加給彈性場(chǎng)，最后求解熱載荷驅(qū)動(dòng)下的結(jié)構(gòu)響應(yīng).

根據(jù)SIMP模型，單元材料的屬性值依賴單元相對(duì)密度值，熱和結(jié)構(gòu)耦合系統(tǒng)可以用下面兩個(gè)有限元方程描述：

(3)

其中：x=(x1,x2,...,xM),表示多種材料相對(duì)密度向量；K1表示熱剛度矩陣；K2表示結(jié)構(gòu)剛度矩陣；U1為溫度向量；U2為位移矢量；P1、P2分別為熱載荷向量和結(jié)構(gòu)外載荷矢量.

2.2 耦合場(chǎng)的敏度分析

基于SIMP法的優(yōu)化問(wèn)題采用基于梯度算法更新設(shè)計(jì)變量，因而需要求解結(jié)構(gòu)響應(yīng)對(duì)設(shè)計(jì)變量的敏度信息.敏度分析一般有3種方法：直接法、半解析法和伴隨矩陣法.此處采用伴隨矩陣法求解目標(biāo)函數(shù)對(duì)各種材料設(shè)計(jì)變量的敏度.

對(duì)于一個(gè)熱固耦合柔順機(jī)構(gòu)的優(yōu)化問(wèn)題可描述為

(4)

其中Ve和V0分別為單元體積和最大許用體積,xmin是最小允許相對(duì)密度值向量.結(jié)構(gòu)響應(yīng)規(guī)定為給定點(diǎn)的輸出位移uout，可寫為

uout=LTU2

(5)

LT是一個(gè)矢量，除了輸出自由度方向的值為1，其他分量全部為0.假定有限元方程(3)已經(jīng)求得，則式(5)可重寫為

(6)

(7)

(8)

這可以通過(guò)設(shè)置如下伴隨向量的值實(shí)現(xiàn)：

(9)

(10)

目前適用于拓?fù)鋬?yōu)化的算法主要有優(yōu)化準(zhǔn)則法和數(shù)學(xué)規(guī)劃算法兩大類[25]，文中采用優(yōu)化準(zhǔn)則算法對(duì)設(shè)計(jì)變量進(jìn)行求解更新.它是一種簡(jiǎn)單方便、收斂迅速的算法，非常適合于文中并行求解單材料子問(wèn)題.采用SIMP模型求解的過(guò)程中，不可避免地會(huì)出現(xiàn)棋盤格現(xiàn)象，文中采用Sigmund[26]提出的敏度過(guò)濾技術(shù)抑制棋盤格式問(wèn)題.

3 設(shè)計(jì)變量并行更新

當(dāng)單材料子問(wèn)題全部求解之后，得到包含M種單材料的解向量xi.按照前述的SIMP模型，在一個(gè)單元中，密度值是一個(gè)0到1之間的連續(xù)值.則對(duì)于一個(gè)多材料問(wèn)題，任意種材料均可在任意單元存在，則單元密度之和可能超過(guò)1，這種單元稱為過(guò)飽和單元.顯然，這種單元的存在違背了SIMP模型的原理.出于對(duì)模型有效性的考慮，需要對(duì)過(guò)飽和單元進(jìn)行處理，從而保證所有的單元密度之和均不超過(guò)1.此處為每種材料預(yù)先設(shè)定一個(gè)優(yōu)先級(jí)，優(yōu)先級(jí)最低的材料最先從過(guò)飽和單元中去除，假設(shè)材料的優(yōu)先級(jí)順序?yàn)镾1,S2,…,SM(從高到低)，調(diào)整后的過(guò)飽和單元密度值為

(11)

需要注意的是，這里處理過(guò)飽和單元的過(guò)程是一個(gè)人工選擇的過(guò)程.它根據(jù)人為設(shè)定的方法確定出優(yōu)先度順序?qū)^(guò)飽和單元的材料進(jìn)行取舍.這種更新方法的流程如圖1所示.

圖1 設(shè)計(jì)變量更新流程

該方案是基于單材料拓?fù)鋬?yōu)化的理論研究提出的，其本質(zhì)是通過(guò)特殊策略將多材料問(wèn)題分解為多種材料對(duì)應(yīng)的單材料子問(wèn)題進(jìn)行更新，并保證設(shè)計(jì)結(jié)果的可行性.因此，該方案下除了可以使用擅長(zhǎng)處理復(fù)雜約束問(wèn)題的移動(dòng)漸近線法等進(jìn)行更新外，還可以使用對(duì)簡(jiǎn)單約束問(wèn)題具有較優(yōu)收斂結(jié)果以及速度的優(yōu)化準(zhǔn)則法進(jìn)行更新，從而使得多材料拓?fù)鋬?yōu)化問(wèn)題的解決效率得到提高.此外，這種方案的執(zhí)行過(guò)程具有很高的分立性與并行性，因此可以很容易地?cái)U(kuò)展到更多種材料的拓?fù)鋬?yōu)化問(wèn)題.

文中的收斂準(zhǔn)則定義為任意種材料的設(shè)計(jì)變量最后相鄰兩步的改變量不大于某一極小值.假定迭代進(jìn)行了t步，修正后的設(shè)計(jì)變量的t-1與t步之間最大差值為

ε=max|(xi)t-(xi)t-1|

(12)

其中,ε是收斂算子，通常取ε≤0.01.

4 數(shù)值算例及討論

本節(jié)選取幾個(gè)基于多材料多場(chǎng)耦合拓?fù)鋬?yōu)化的的例子，說(shuō)明前述理論模型的有效性.為不失一般性，在給出雙材料拓?fù)浣Y(jié)果之后，三材料的例子也已給出.設(shè)計(jì)域內(nèi)溫度向量為恒量，所有算例均不施加外在負(fù)載，只需要考慮純熱驅(qū)動(dòng)下柔順機(jī)構(gòu)的設(shè)計(jì)問(wèn)題.

4.1 雙材料熱驅(qū)動(dòng)器

本例(算例1)采用一個(gè)比較常見的熱驅(qū)動(dòng)器作為優(yōu)化對(duì)象，兩種材料的設(shè)計(jì)參數(shù)初始值如表 1所示.設(shè)計(jì)域與邊界條件如圖2所示，上、下、左3個(gè)壁均絕熱，在右邊中間區(qū)域獲得位移輸出.由于對(duì)稱性，只需要對(duì)上半部分進(jìn)行優(yōu)化就可得到整體機(jī)構(gòu)結(jié)果.圖3(a)-3(d)分別為兩種材料(紅：綠)體積比為5∶1的熱驅(qū)動(dòng)器設(shè)計(jì)中間結(jié)果，圖3(e)、圖4分別為熱驅(qū)動(dòng)器最終結(jié)果以及優(yōu)化迭代曲線.

表 1 設(shè)計(jì)參數(shù)初始值(算例1)

圖2 熱驅(qū)動(dòng)器設(shè)計(jì)域及邊界條件(算例1)

Fig.2 Design domain and boundary conditions of heat actuator (Example 1)

圖3 熱驅(qū)動(dòng)器設(shè)計(jì)結(jié)果(算例1)

圖4 對(duì)應(yīng)圖3中熱驅(qū)動(dòng)器的優(yōu)化收斂曲線

Fig.4 Conversion history of heat actuator corresponding to Fig.3

觀察圖3最終的拓?fù)浣Y(jié)果，可以發(fā)現(xiàn)使用文中提出的多材料設(shè)計(jì)方法得到的機(jī)構(gòu)分塊明顯而獨(dú)立，從制造的角度看非常適用，可以先分塊制造，再使用搭接等連接技術(shù)連接成一體以生產(chǎn).

4.2 有約束區(qū)域的雙材料熱驅(qū)動(dòng)器

本例(算例2)設(shè)計(jì)一個(gè)材料分布區(qū)域受到限制的雙材料熱驅(qū)動(dòng)器，設(shè)計(jì)域與邊界條件如圖5所示，設(shè)計(jì)域中間區(qū)域作為主動(dòng)受熱驅(qū)動(dòng)材料分布區(qū)，禁止被驅(qū)動(dòng)材料存在，而周圍區(qū)域作為受動(dòng)材料分布區(qū)，禁止主動(dòng)材料存在，同時(shí)結(jié)構(gòu)上、下兩側(cè)允許水平方向移動(dòng).兩種材料的設(shè)計(jì)參數(shù)初始值與算例1相同，其中中間限制區(qū)域的大小為34 mm×34 mm.圖6(a)、圖7(a)分別為兩種材料的體積比均為5∶1的熱驅(qū)動(dòng)器設(shè)計(jì)結(jié)果以及目標(biāo)函數(shù)收斂曲線；圖6(b)、圖7(b)分別為兩種材料(紅∶綠)體積比分別為2∶1的熱驅(qū)動(dòng)器設(shè)計(jì)結(jié)果以及收斂曲線.

圖5 熱驅(qū)動(dòng)器設(shè)計(jì)域及邊界條件(算例2)

Fig.5 Design domain and boundary conditions of heat actuator (Example 2)

同樣，從圖6中的設(shè)計(jì)結(jié)果可以看出，兩種材料邊界清晰，分塊獨(dú)立.通過(guò)迭代曲線不難發(fā)現(xiàn)，優(yōu)化迭代過(guò)程收斂穩(wěn)定，說(shuō)明對(duì)于復(fù)雜的多材料和多場(chǎng)耦合組合優(yōu)化問(wèn)題，文中所提出的方法具有較好的魯棒性.

圖6 熱驅(qū)動(dòng)器不同約束下的設(shè)計(jì)結(jié)果(算例2)

Fig.6 Design result of heat actuators under different constraints(Example 2)

Fig.7 Conversion history of heat actuators corresponding to Fig.6

4.3 三材料熱驅(qū)動(dòng)器

為不失一般性，此處(算例3)用一個(gè)三材料的熱驅(qū)動(dòng)器說(shuō)明方法的有效性.其設(shè)計(jì)域和邊界條件與算例1相同，3種材料的設(shè)計(jì)參數(shù)初始值如表 2所示，其他參數(shù)與之前算例相同.圖8(a)、圖9(a)為3種材料體積比為4∶1∶1的熱驅(qū)動(dòng)器設(shè)計(jì)結(jié)果及收斂曲線；圖8(b)、 圖9(b)為3種材料體積分別為10∶1∶1的熱驅(qū)動(dòng)器設(shè)計(jì)結(jié)果及收斂曲線.

圖8 三材料熱驅(qū)動(dòng)器不同約束下的設(shè)計(jì)結(jié)果

Fig.8 Design result of three-materials heat actuators under different constraints

表2 設(shè)計(jì)參數(shù)初始值(算例3)

圖9 三材料熱驅(qū)動(dòng)器對(duì)應(yīng)圖8中設(shè)計(jì)結(jié)果的收斂曲線

Fig.9 Conversion history of three-materials heat actuators in Fig.8

從圖8的設(shè)計(jì)結(jié)果可以看到，優(yōu)化三材料熱驅(qū)動(dòng)柔順機(jī)構(gòu)時(shí)，依然能夠得到比較理想的拓?fù)浣Y(jié)構(gòu)，各材料之間邊界清晰，分塊獨(dú)立.這說(shuō)明基于并行優(yōu)化策略的模型能夠突破傳統(tǒng)的雙材料限制，進(jìn)行真正意義上的任意多種材料的拓?fù)鋬?yōu)化設(shè)計(jì).

5 結(jié)論

(1)提出了一種新的基于多材料設(shè)計(jì)多場(chǎng)耦合柔順機(jī)構(gòu)的求解方法.將復(fù)雜的多材料問(wèn)題離散為簡(jiǎn)單的單材料子問(wèn)題，然后并行獨(dú)立求解.

(2)基于剛度可加性理論，得到了全局剛度矩陣與各材料之間呈線性關(guān)系，并將之推廣至其他耦合物理場(chǎng).

(3)該方法在理論和應(yīng)用上都比較簡(jiǎn)單，并且能夠處理任意多種材料.數(shù)值算例表明，新方法設(shè)計(jì)的機(jī)構(gòu)合理，各材料之間分塊明顯，利于加工，因而方法是有效的.

(4)重點(diǎn)研究了多場(chǎng)耦合中的熱固耦合效應(yīng)，為后續(xù)的熱、電、結(jié)構(gòu)場(chǎng)耦合理論分析和研究提供了極為重要的參考價(jià)值.

(5)算例部分只給出了依靠純熱驅(qū)動(dòng)作用下的柔順機(jī)構(gòu)設(shè)計(jì)與分析，表明熱效應(yīng)對(duì)結(jié)構(gòu)變形的巨大影響.需要說(shuō)明的是，對(duì)于已有較多的外力和熱場(chǎng)同時(shí)作用下的結(jié)構(gòu)分析，文中所提方法也是有效的，熱載荷與外力同時(shí)作用下的情形類似于多輸入柔順機(jī)構(gòu)設(shè)計(jì)；其次，并行優(yōu)化策略尚未在數(shù)學(xué)上得到嚴(yán)格證明，也無(wú)法說(shuō)明優(yōu)化的結(jié)果是全局最優(yōu)的，但是大量的算例證明了該方法的有效性.

[1] 張憲民.柔順機(jī)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì) [J].機(jī)械工程學(xué)報(bào),2003,39(11)：47- 51. ZHANG Xian-min.Topology optimozation of compliant mechanisms [J].Chinese Journal of Mechanical Engineering,2003,39(11)：47- 51.

[2] HAO G.Towards the design of monolithic decoupled XYZ compliant parallel mechanisms for multi-function applications [J].Mechanical Sciences,2013,4：291- 302.

[3] 于靖軍,宗光華,畢樹生.全柔性機(jī)構(gòu)與MEMS [J].光學(xué)精密工程,2001,9(1):1- 5. YU,Jing-jun,ZONG Guang-hua,BI Shu-sheng.Fully compliant mechanisms and MEMS [J].Opt Precis Eng,2001,9(1)：1- 5.

[4] 陳貴敏,李端玲.平面柔順機(jī)構(gòu)的自由度 [J].機(jī)械工程學(xué)報(bào),2010,46(13)：48- 53. CHEN Gui-min,LI Duan-ling.Degree of freedom of planar compliant mechanisms [J].Journal of Mechanical Engineering,2010,46(13)：48- 53.

[5] ROZVANY G I N.A critical review of established methods of structural topology optimization [J].Structural and Multidisciplinary Optimization,2009,37(3)：217- 237.

[6] SIGMUND O.On the design of compliant mechanisms using topology optimization [J].Journal of Structural Mechanics,1997,25(4)：493- 524.

[7] 羅震,蒙永立,郭文德,等.分布式柔性機(jī)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)的理論和算法 [J].機(jī)械工程學(xué)報(bào),2006,42(10)：27- 36. LUO Zhen,MENG Yong-li,GUO Wen-de,et al.Theoretical and algorithmic on topology optimization design of distributed compliant mechanisms [J].Chinese Journal of Mechanical Engineering,2006,42(10)：27- 36.

[8] BRUNS T E.A revaluation of the SIMP method with filtering and an alternative formulation for solid-void topology optimization [J].Structural and Multidisciplinary Optimization,2005,30(6)：428- 436.

[9] RAMRAKHYANI D S,FRECKER M I,LESIEUTRE G A.Hinged beam elements for the topology design of compliant mechanisms using the ground structure approach [J].Structural and Multidisciplinary Optimization,2009,37(6)：557- 567.

[10] WANG Nian-feng,ZHANG Xian-min.Multi-material topology optimization of complaint mechanism using ground structure approach [C]∥Proceedings of International Conference on Manipulation,Manufacturing and Measurement on the Nanoscale(3M-NANO).[S.l.]：IEEE,2014：249- 254.

[11] ZHU Ben-liang,ZHANG Xian-min,WANG Nian-feng.Topology optimization of hinge-free compliant mechanisms with multiple outputs using level set method [J].Structural and Multidisciplinary Optimization,2013,47(5)：659- 672.

[12] 羅俊召,陳立平,羅震.基于半隱式格式的水平集法連續(xù)體結(jié)構(gòu)形狀和拓?fù)鋬?yōu)化方法 [J].固體力學(xué)學(xué)報(bào),2008,29(2)：175- 180. LUO Jun-zhao,CHEN Li-ping,LUO Zhen.A level set method for topology optimization using AOS scheme [J].Chinese Journal of Solid Mechanics,2008,29(2)：175- 180.

[13] 李冬梅.多場(chǎng)耦合及多相材料的柔順機(jī)構(gòu)拓?fù)鋬?yōu)化研究 [D].廣州：華南理工大學(xué),2011.

[14] SIGMUND O.Design of multi-physics actuators using topology optimization(Part I)：one-material structures [J].Computer Methods in Applied Mechanics and Engineering,2001,190(49)：6577- 6604.

[15] SIGMUND O.Design of multi-physics actuators using topology optimization(Part II)：Two-material structures [J].Computer Methods in Applied Mechanics and Engineering,2001,190(49)：6605- 6627.

[16] 李冬梅,張憲民,王念峰,等.基于可靠性約束的熱固耦合結(jié)構(gòu)拓?fù)鋬?yōu)化 [J].華南理工大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,39(6)：42- 46. LI Dong-mei,ZHANG Xian-min,WANG Nian-feng,et al.Topology optimization of thermo-mechanical coupling structures based on reliability constraint [J].Journal of South China University of Technology(Natural Science Edition),2011,39(6)：42- 46.

[17] TAKEZAWA A,NISHIWAKI S,KITAMURA M.Shape and topology optimization based on the phase field method and sensitivity analysis [J].Journal of Computational Physics,2010,229(7)：2697- 2718.

[18] TAKEZAWA A,YOON G H,JEONG S H,et al.Structural topology optimization with strength and heat conduction constraints [J].Computer Methods in Applied Mechanics and Engineering,2014,276：341- 361.

[19] KOGA A A,LOPES E C C,NOVA H F V,et al.Deve-lopment of heat sink device by using topology optimization [J].International Journal of Heat and Mass Transfer,2013,64：759-772.

[20] GUO X,ZHANG W,ZHONG W.Stress-related topology optimization of continuum structures involving multi-phase materials [J].Computer Methods in Applied Mechanics and Engineering,2014,268：632- 655.

[21] WANG Y,LUO Z,KANG Z,et al.A multi-material level set-based topology and shape optimization method [J].Computer Methods in Applied Mechanics and Engineering,2015,283：1570- 1586.

[22] TAVAKOLI R.Multimaterial topology optimization by volume constrained Allen-Cahn system and regularized projected steepest descent method [J].Computer Methods in Applied Mechanics and Engineering,2014,276：534- 565.

[23] 左孔天,錢勤,趙雨?yáng)|,等.熱固耦合結(jié)構(gòu)的拓?fù)鋬?yōu)化設(shè)計(jì)研究 [J].固體力學(xué)學(xué)報(bào),2006,26(4)：447- 452. ZUO Kong-Tian,QIAN Qin,ZHAO Yu-dong,et al.Research on the topology optimization about thermo-structural coupling field [J].Acta Mechanica Solida Sinica,2005,26(4)：447- 452.

[24] 莊春剛.基于水平集的多材料結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)方法與應(yīng)用 [D].上海：上海交通大學(xué),2007.

[25] 羅震.基于變密度法的連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)技術(shù)研究 [D].武漢：華中科技大學(xué),2005.

[26] SIGMUND O.A 99 line topology optimization code written in Matlab [J].Structural and Multidisciplinary Optimization,2001,21(2)：120- 127.

Design of Thermal-Mechanical Coupling Multiple-Material Compliant Mechanisms Based on Parallel Strategy

ZHANGXian-minHUKaiWANGNian-fengGUOHao

(School of Mechanical and Automotive Engineering∥Guangdong Provincial Key Laboratory of Precision Equipment and Manufacturing Technology,Guangzhou 510640, Guangdong, China)

In the field of structural design, thermal effects have been extensively discussed. In this paper, based on the topology optimization method, the design of the thermally-driven compliant mechanism is investigated, and a corresponding parallel penalization model is constructed. The core of this model is to decompose a complex multiple-material thermal-mechanical coupling problem into a series of single-material thermal-mechanical coupling sub-problems. Then, the sub-problems are solved concurrently, and according to the requirement of the overall system, the solutions to all the sub-problems are coordinated so as to achieve the solution to the original problem. For all the sub-problems, the thermal-mechanical coupling optimization models are constructed based on the solid isotropic material with penalization(SIMP) model, and then the sequential coupling method is adopted to conduct the thermal-structural coupling physics analysis, thus achieving concurrent and independent solutions. The constructed model is simple in concept and practice and can be used to handle an arbitrary number of candidate materials. In addition, it can avoid fragile topologies and is, therefore, helpful in manufacturing. Finally, the effectiveness of the constructed model is verified by several numerical examples.

topology optimization; compliant mechanism; thermal-mechanical coupling; sensitivity analysis

2015- 12- 07

國(guó)家自然科學(xué)基金資助項(xiàng)目(51205134,51575187,91223201)；廣州市產(chǎn)學(xué)研協(xié)同創(chuàng)新重大專項(xiàng)(2014Y2-00217)；廣東省自然科學(xué)基金資助項(xiàng)目(S2013030013355)；黃埔區(qū)重大專項(xiàng)(20150000661)；華南理工大學(xué)中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金資助項(xiàng)目(2015ZZ007) Foundation items: Supported by the National Natural Science Foundation of China(51205134,51575187,91223201) and the Natural Science Foundation of Guangdong Province(S2013030013355)

張憲民(1964-)，男，教授，博士生導(dǎo)師，主要從事機(jī)構(gòu)學(xué)及精密制造裝備研究.E-mail：zhangxm@scut.edu.cn

? 通信作者: 王念峰(1978-)，男，教授，博士生導(dǎo)師，主要從事結(jié)構(gòu)優(yōu)化及機(jī)器人技術(shù)研究.E-mail：menfwang@scut.edu.cn

1000- 565X(2016)10- 0022- 07

TB 112；TB 114

10.3969/j.issn.1000-565X.2016.10.004