非均質(zhì)材料熱力耦合及彈塑性損傷分析的多尺度方法研究.pdf

1、非均質(zhì)材料廣泛應(yīng)用于航空航天、化工、微電子及日常生活中，在復(fù)雜熱力學(xué)環(huán)境中，這些材料呈現(xiàn)出與傳統(tǒng)均質(zhì)材料不同的熱力學(xué)特性。在這些材料的生產(chǎn)制造中，可通過(guò)調(diào)整微結(jié)構(gòu)的設(shè)計(jì)參數(shù)，進(jìn)一步優(yōu)化材料性能。然而，在非均質(zhì)材料結(jié)構(gòu)中，非均勻性特征長(zhǎng)度l和結(jié)構(gòu)整體尺寸L存在差異，由于基于連續(xù)介質(zhì)理論的傳統(tǒng)方法（如有限元法）在描述材料非均質(zhì)性上有明顯的不足，需要足夠細(xì)的劃分網(wǎng)格才能給出期望的精度，這意味著會(huì)產(chǎn)生較大規(guī)模自由度。因此，如何對(duì)非均質(zhì)材料在多物

2、理場(chǎng)耦合作用下的行為進(jìn)行高效模擬，成為眾多學(xué)者研究的熱點(diǎn)。2014年，美國(guó)負(fù)責(zé)科研與工程（R&E）的助理國(guó)防部長(zhǎng)(ASD) Bob Baker在國(guó)防科技重點(diǎn)報(bào)告中指出，工程超材料是六大顛覆性基礎(chǔ)研究領(lǐng)域之一，其中超大規(guī)模復(fù)雜結(jié)構(gòu)的建模與設(shè)計(jì)方法又是該領(lǐng)域的主要挑戰(zhàn)，需要在復(fù)雜人工結(jié)構(gòu)計(jì)算工具與快速算法上取得重大突破。從該報(bào)告中可以看出非均質(zhì)材料的數(shù)值方法研究是對(duì)未來(lái)國(guó)防科技發(fā)展具有重大影響和深遠(yuǎn)意義的基礎(chǔ)研究領(lǐng)域之一。
　　本文首

3、先針對(duì)非均質(zhì)材料熱彈性問(wèn)題，提出了一種耦合擴(kuò)展多尺度有限元方法(CEMsFEM)?？紤]傳熱和力學(xué)分析的耦合效應(yīng)，基于求解標(biāo)量場(chǎng)問(wèn)題的多尺度有限元法和材料彈性力學(xué)性能的擴(kuò)展多尺度有限元法，對(duì)溫度場(chǎng)和位移場(chǎng)分別構(gòu)造反映材料非均質(zhì)特性的數(shù)值基函數(shù)，并在此基礎(chǔ)上對(duì)熱力耦合問(wèn)題進(jìn)行多尺度分析。在熱力耦合問(wèn)題中，數(shù)值基函數(shù)的構(gòu)造忽略了粗單元域內(nèi)微觀熱流或熱膨脹載荷的影響，需通過(guò)在粗單元域施加邊界約束，求解微觀載荷對(duì)局部域的影響。此外還進(jìn)一步對(duì)多尺度

4、方法的存儲(chǔ)需求和復(fù)雜度進(jìn)行了理論分析。最后討論了溫度場(chǎng)和位移場(chǎng)數(shù)值基函數(shù)構(gòu)造的線性、周期性邊界條件及超樣本技術(shù)對(duì)處理不同類(lèi)型夾雜的適用性，同時(shí)為驗(yàn)證本文算法的精度和對(duì)復(fù)雜實(shí)際結(jié)構(gòu)模擬的有效性，對(duì)含有不同夾雜的非均質(zhì)材料結(jié)構(gòu)進(jìn)行了多尺度熱彈性分析，其中對(duì)于復(fù)雜邊界模型，采用了多尺度單元與常規(guī)單元混合建模方法，提高了建模的靈活性。
　　其次，考慮非均質(zhì)材料非線性熱彈性問(wèn)題分析中材料性質(zhì)的溫度相關(guān)性，提出了一種預(yù)修正多尺度方法。對(duì)微觀溫

5、度和位移進(jìn)行了加法分解，即數(shù)值基函數(shù)表征部分和微觀擾動(dòng)部分，基于此構(gòu)造了材料性質(zhì)溫度相關(guān)性問(wèn)題的多尺度分析方法，通過(guò)把微觀擾動(dòng)帶來(lái)的內(nèi)部載荷項(xiàng)預(yù)先考慮進(jìn)入宏觀尺度計(jì)算，提高了分析精度。此外還研究了微觀擾動(dòng)構(gòu)造的局部放松邊界條件及超樣本技術(shù)，有效地消除邊界效應(yīng)，提高算法精度。最后通過(guò)數(shù)值算例測(cè)試了本文方法對(duì)多種非均質(zhì)材料模型（如功能梯度材料、冷卻結(jié)構(gòu)）多尺度非線性熱彈性分析的計(jì)算精度和效率，以及對(duì)邊界效應(yīng)處理的有效性。
　　此外針對(duì)

6、非均質(zhì)材料的熱彈塑性問(wèn)題，提出了一種增強(qiáng)數(shù)值基函數(shù)多尺度分析方法。該方法基于單位分解的基本思想，在擴(kuò)展多尺度有限元法框架下構(gòu)造了新的增強(qiáng)數(shù)值基函數(shù)，該基函數(shù)由兩部分構(gòu)成，即基本部分和擴(kuò)充部分。其中基本部分采用擴(kuò)展多尺度有限元法中的基函數(shù)，在多尺度分析之前預(yù)先構(gòu)造，在分析過(guò)程中保持不變。擴(kuò)充部分采用微觀位移擾動(dòng)和單位分解函數(shù)的乘積構(gòu)成，通過(guò)在每一增量載荷步首次迭代過(guò)程中重新計(jì)算微觀擾動(dòng)位移，更新擴(kuò)充基函數(shù)。為了使本文算法更高效，進(jìn)一步提出

7、了擴(kuò)充基函數(shù)的自適應(yīng)更新準(zhǔn)則，在計(jì)算過(guò)程中只對(duì)滿足準(zhǔn)則的數(shù)值基函數(shù)進(jìn)行更新。由于本文方法通過(guò)擴(kuò)充基函數(shù)把基本基函數(shù)難以表征的微觀擾動(dòng)帶入了有限元試探函數(shù)空間，可通過(guò)降尺度計(jì)算直接準(zhǔn)確的獲得微觀位移和應(yīng)力，不需要對(duì)結(jié)果進(jìn)行修正。數(shù)值算例結(jié)果驗(yàn)證了增強(qiáng)數(shù)值基函數(shù)的精度及熱彈塑性問(wèn)題多尺度分析方法的計(jì)算效率。
　　最后，采用增強(qiáng)數(shù)值基函數(shù)的多尺度分析方法，進(jìn)一步研究了非均質(zhì)材料彈塑性損傷問(wèn)題。與上述考慮應(yīng)變硬化的熱彈塑性分析不同的是，損

8、傷會(huì)導(dǎo)致材料軟化，從而導(dǎo)致應(yīng)變局部化問(wèn)題。因此所采用的迭代策略是在每一個(gè)載荷步收斂時(shí)，更新擴(kuò)充數(shù)值基函數(shù)，并對(duì)這個(gè)載荷步重新計(jì)算，直到所有的數(shù)值基函數(shù)均滿足精度需求。通過(guò)對(duì)含有圓形和橢圓夾雜的非均質(zhì)結(jié)構(gòu)模型數(shù)值分析驗(yàn)證了非線性多尺度算法的精度和有效性。
　　本論文中所研發(fā)的程序是在自主產(chǎn)權(quán)的開(kāi)放式計(jì)算力學(xué)軟件SiPESC上展開(kāi)的，在其工程數(shù)據(jù)庫(kù)管理系統(tǒng)和結(jié)構(gòu)有限元子系統(tǒng)的基礎(chǔ)上，研發(fā)了傳熱及熱力耦合問(wèn)題的面向?qū)ο笥邢拊治鲎酉到y(tǒng)，

9、并實(shí)現(xiàn)了本文提出的多尺度算法。程序系統(tǒng)利用插件及擴(kuò)展的管理機(jī)制，采用面向?qū)ο蟮腃++程序設(shè)計(jì)語(yǔ)言和Factory及Builder等軟件設(shè)計(jì)模式，具有良好的開(kāi)放性、擴(kuò)展性和模塊化的特點(diǎn)。多尺度單元的數(shù)值基函數(shù)通過(guò)多點(diǎn)約束的方法引入到了平臺(tái)的節(jié)點(diǎn)控制陣中，總體矩陣和載荷矢量的集成以及節(jié)點(diǎn)和單元數(shù)據(jù)的計(jì)算（即升尺度和降尺度）均與常規(guī)有限元方法相同，不需要重新編寫(xiě)代碼。多尺度算法的實(shí)現(xiàn)是基于平臺(tái)“算法+模型”的設(shè)計(jì)模式，可實(shí)現(xiàn)動(dòng)態(tài)組裝，并能夠完