皮浩，石艷

（四川輕化工大學(xué) 機(jī)械工程學(xué)院，四川 自貢 643000）

關(guān)鍵字：聚丙乙烯；高分子拉伸失效；瞬態(tài)分析；有限元分析

在現(xiàn)代工程技術(shù)領(lǐng)域中，聚丙烯因其優(yōu)良的性能、廉價(jià)的成本，常制成容器、管材和各種機(jī)械零部件，使用量逐年增加。因此，有必要對它的性能進(jìn)行嚴(yán)格的測試和確認(rèn)，以表征材料的性質(zhì)和狀態(tài)［1］。但是，聚丙烯這類高分子材料具有微觀、細(xì)觀和宏觀層次的特征，導(dǎo)致與之相對應(yīng)的強(qiáng)度及破壞行為比較復(fù)雜，不僅涉及材料學(xué)、固體力學(xué)、化學(xué)和物理學(xué)等學(xué)科領(lǐng)域，還具有多種學(xué)科互相滲透的特點(diǎn)［2-3］。總的說來，分析難度較大，常不能很好的預(yù)估聚丙烯這類高分子材料的使用年限和剩余壽命。

國內(nèi)對聚丙烯這類高分子材料失效過程的研究重點(diǎn)，多集中在失效機(jī)理的分析上。蘭州交通大學(xué)的王軍璽等［4］人做了 “大體積混凝土溫度場仿真分析在ANSYS上的實(shí)現(xiàn)”，實(shí)現(xiàn)了大型工程問題分析過程參數(shù)化，具有借鑒意義。南京理工大學(xué)的梁亮偉等人對低溫注塑進(jìn)行了相關(guān)的研究，提出了仿真與實(shí)驗(yàn)存在誤差主要是邊界條件難以精確確定這一觀點(diǎn)［6］。國外復(fù)合泡沫塑料領(lǐng)域力學(xué)行為是其研究的熱點(diǎn)［5］，對于高分子材料失效的仿真研究較少，常以直接實(shí)驗(yàn)的方式測試材料性能。研究方向大多集中在如何合成優(yōu)質(zhì)材料上［6-8］?？偟恼f來，學(xué)界對聚丙烯這類高分子的失效過程的理論研究比較豐富，對失效過程仿真研究相對偏少，研究成果也僅限于高質(zhì)量材料的發(fā)現(xiàn)與制備方面。

本文參照拉伸實(shí)驗(yàn)得到的位移與時(shí)間等參數(shù)指標(biāo)，通過有限元瞬態(tài)分析法，調(diào)試載荷步仿真拉伸失效過程。使用拉伸試驗(yàn)機(jī)，在不同拉伸速率下，分別測試聚丙烯試件拉伸失效過程。建立試件物理模型和體單元數(shù)值模型，研究拉伸速率、單元類型、網(wǎng)格密度對拉伸失效過程仿真效果的影響。并分析了60mm／min時(shí)，實(shí)驗(yàn)組與仿真結(jié)果組出現(xiàn)差異的原因。

1 有限元仿真研究

1.1 建立試件物理模型

以命令流的方式建立了物理模型模型，模型尺寸參考 《塑料拉伸標(biāo)準(zhǔn)2006》。以代碼方式建立的模型兼容性更好。本次實(shí)驗(yàn)試件選用典型的高分子材料，聚丙烯。試件全新無磨損。其性能參數(shù)如表1所示。

表1 聚丙烯試件性能參數(shù)

1.2 網(wǎng)格劃分

定義單元格類型為結(jié)構(gòu)單元SOLID226。SOLID226是三維單元格，滿足應(yīng)力應(yīng)變以及溫度場的耦合分析。采用智能網(wǎng)格劃分，以sweep的方式完成，生成了六面體網(wǎng)格。對1mm、1.5mm、2.5mm、3mm尺寸的仿真對比后，得出2.5mm在滿足精度要求下，計(jì)算時(shí)間最優(yōu)，且計(jì)算機(jī)不會(huì)受到過高的負(fù)載，適合一般計(jì)算機(jī)。

1.3 施加邊界條件

模擬試件拉伸時(shí)一端固定、另一端在拉伸實(shí)驗(yàn)機(jī)作用下拉伸變形的實(shí)際情況。對模型端面施加了固定位移約束。對模型表面施加空氣自然對流散熱。對流系數(shù)根據(jù)一般情況設(shè)置為6.7，環(huán)境溫度設(shè)置為20℃。前處理結(jié)果如圖1所示。

圖1 模型前處理

1.4 施加載荷

試樣在拉伸的過程中一般要經(jīng)歷彈性、上屈服、下屈服、強(qiáng)化、頸縮5個(gè)階段。這里簡化模型，只考慮彈性、上屈服、下屈服、強(qiáng)化這四個(gè)特征變化明顯的階段。由本次設(shè)計(jì)實(shí)驗(yàn)組得到仿真數(shù)據(jù)，每個(gè)階段有著各自對應(yīng)的平均位移與平均時(shí)間，這里為了貼近試驗(yàn)，采用試驗(yàn)實(shí)測數(shù)據(jù)。拉伸試驗(yàn)機(jī)以60mm／min拉伸時(shí)，記錄的時(shí)間位移參數(shù)整理如表2所示。

表2 拉伸階段時(shí)間位移對應(yīng)表

對應(yīng)地模擬拉伸過程，將此次仿真設(shè)為4個(gè)載荷步，分別代表上表所示四個(gè)階段?？紤]到計(jì)算時(shí)間的影響，每個(gè)載荷步驟分為20個(gè)子步。受載模型如圖2所示。

圖2 施加載荷

2 后處理階段與分析對比

高分子聚合物內(nèi)部由高分子鏈段組成，錯(cuò)綜復(fù)雜交聯(lián)成網(wǎng)，不完全符合連續(xù)性假設(shè)。聚丙烯試件材料內(nèi)會(huì)形成隨機(jī)的塊狀結(jié)晶團(tuán)或類似物，這些應(yīng)力集中物隨機(jī)排列放置，破壞了均勻性假設(shè)。并且，由于分子鏈段隨機(jī)排列或糾纏，不可避免的導(dǎo)致各向異性。高分子聚合物受載時(shí)會(huì)調(diào)節(jié)分子內(nèi)部構(gòu)象，拉伸速率快的組，高分子重構(gòu)和調(diào)整的時(shí)間短；拉伸速率慢的組，能夠有時(shí)間調(diào)整分子內(nèi)部構(gòu)象，實(shí)現(xiàn)一定的強(qiáng)度增強(qiáng)。以上因素，會(huì)影響仿真數(shù)值結(jié)果。以60mm／min拉伸時(shí)的情況進(jìn)行說明。

2.1 應(yīng)變分析

試件經(jīng)拉伸試驗(yàn)機(jī)拉伸失效時(shí)的應(yīng)變尺度與仿真效果如圖3、圖4所示。

圖3 拉伸斷裂的試件

圖4 仿真變形與實(shí)際變形

對比試驗(yàn)數(shù)據(jù)與仿真結(jié)果可知，仿真計(jì)算結(jié)果最大位移為317.307mm，小于實(shí)驗(yàn)組平均值330.05mm。拉伸過程中，試件固定端的應(yīng)變較小，被拉端拉伸端應(yīng)變較大，且沿軸向成遞增態(tài)勢，符合實(shí)際情況試件的變形。聚丙烯試件本身是彈塑性材料，拉伸失效實(shí)驗(yàn)過程中帶有一定蠕變變形，而蠕變變形程度受溫度、時(shí)間和載荷的影響，這一分析與盧子興等人的研究符合［9］。而仿真默認(rèn)模型為線彈性體，未能考慮到此影響，產(chǎn)生了偏差；試件拉伸的失效的過程中，高分子鏈段受到作用，不斷調(diào)整內(nèi)部鏈段構(gòu)象以應(yīng)對變形，這種性質(zhì)具有不確定性，對試件材料性能有影響［10-11］。因此本次拉伸試驗(yàn)中，出現(xiàn)了強(qiáng)度隨拉伸而升高的現(xiàn)象。

2.2 應(yīng)力分析

拉伸試驗(yàn)機(jī)可以獲取到拉伸過程的數(shù)據(jù)。整個(gè)拉伸過程大致劃分為彈性、上屈服、下屈服、變形等階段，直至試件突然斷裂而試驗(yàn)終止。由拉伸試驗(yàn)機(jī)得到的應(yīng)力-時(shí)間變化趨勢圖如圖5。

圖5 實(shí)驗(yàn)組應(yīng)力-時(shí)間變化

仿真計(jì)算完成后，使用時(shí)間歷程后處理器，沿試件軸向選取的四個(gè)節(jié)點(diǎn)的應(yīng)力-時(shí)間變化圖如圖6。

圖6 仿真組中抽取的四個(gè)節(jié)點(diǎn)應(yīng)力-時(shí)間變化

根據(jù)各仿真階段的應(yīng)力云圖和變化趨勢圖，聚丙烯試件拉伸失效過程總體遵循拉伸失效準(zhǔn)則，應(yīng)力與時(shí)間呈現(xiàn)出與實(shí)驗(yàn)組吻合的趨勢，即先增大后減低，然后維持直至斷裂。仿真的結(jié)果相對實(shí)驗(yàn)結(jié)果偏大，原因是，此次仿真默認(rèn)整個(gè)試件拉伸過程為線彈性變形，而實(shí)際情況是只有前期很短位移內(nèi)滿足胡克定律［12］，而之后屈服、頸縮階段并不滿足，而且實(shí)際拉伸時(shí)失效情況還更復(fù)雜。

2.3 溫度場分析

使用熱成像儀實(shí)時(shí)記錄拉伸失效過程溫度場變化情況［13］。選取典型階段的溫度值作為對比依據(jù)。結(jié)果如圖7所示。

圖7 彈性階段末期

彈性階段末期，聚丙烯試件中間段變形最為劇烈，能量最多，最高溫度顯示為23℃，并向兩端呈現(xiàn)擴(kuò)散態(tài)勢。見圖8。

圖8 上屈服階段末期

上屈服階段末期，聚丙烯試件最高溫度達(dá)到50℃，變化劇烈，說明此階段拉伸做功效果顯著。并且可得溫度作用范圍擴(kuò)大，說明變形程度劇烈的段長擴(kuò)大。見圖9。

圖9 下屈服階段末期

下屈服階段末期，聚丙烯試件溫度值變化不再明顯，峰值溫度達(dá)到52℃。溫度作用范圍進(jìn)一步擴(kuò)大，達(dá)到最大。見圖10。

圖10 彈性階段末期

強(qiáng)化階段末期，聚丙烯試件溫度峰值變化不再明顯，維持在51.827℃左右，溫度作用范圍擴(kuò)大。仿真計(jì)算結(jié)束后，截取典型階段時(shí)刻的溫度值進(jìn)行對比，如表3。

表3 實(shí)驗(yàn)組與仿真結(jié)果值對比

選取熱源上的一點(diǎn)，使用Time Hist Postpro后處理器，顯示其時(shí)間與溫度的對應(yīng)關(guān)系，如圖11。

圖12 熱面源中心點(diǎn)時(shí)間-溫度

由以上對比可知，仿真數(shù)值大小與實(shí)驗(yàn)結(jié)果趨勢吻合。在彈性階段，試件通過形變來應(yīng)對拉伸試驗(yàn)機(jī)的做功，因此拉伸實(shí)驗(yàn)機(jī)對試件作的做功主要轉(zhuǎn)化為試件的彈性勢能，溫度變化不大，僅僅增加了3℃左右，但是這個(gè)過程的時(shí)間極短。隨后的屈服階段，聚丙烯試件繼續(xù)受到拉伸試驗(yàn)機(jī)加載，被拉試件無法存儲(chǔ)過多的彈性勢能，伴隨宏觀上試件形變量增大和微觀上高分子鍵的破壞，試件溫度峰值增加明顯，達(dá)到50℃左右。拉伸失效末期，溫度變化不太明顯了，穩(wěn)定在52℃左右，直至斷裂。

3 結(jié)論

以較高速率拉伸聚丙烯試件時(shí)，相對于較低速率，其斷裂失效過程需要的時(shí)間短，試件變形程度小，而溫度峰值更高。以熱源面方式模擬拉伸做功多少，會(huì)影響溫度場隨拉伸的熱源的擴(kuò)散效果，需要尋找更好的替代方案。另外，聚丙烯試件失效瞬間，拉伸試驗(yàn)機(jī)會(huì)較劇烈的抖動(dòng)并伴隨聲響，這種現(xiàn)象還未能被仿真出來。