顆粒增強(qiáng)復(fù)合材料隨機(jī)分布模型的數(shù)值模擬

1、 10期單豪良,等:顆粒增強(qiáng)復(fù)合材料隨機(jī)分布模型的數(shù)值模擬2725實(shí)驗(yàn)測(cè)量是否一致,做如下的假設(shè):(a顆粒和基體只發(fā)生線(xiàn)彈性小變形;(b顆粒、基體為各向同性材料,由于顆粒隨機(jī)分布,在宏觀意義上,復(fù)合材料為各向同性材料;(c顆粒為均勻球體,顆粒直徑符合正態(tài)分布;(d顆粒之間無(wú)相互擠壓,在給定的載荷范圍內(nèi),各相材料不發(fā)生破壞;(e顆粒和基體之間粘合完好,不產(chǎn)生界面脫粘?，F(xiàn)以鋁基碳化硅顆粒增強(qiáng)復(fù)合材料為研究對(duì)象,6066鋁和夾雜物碳化硅的材料性能如表1所示。數(shù)值分析采用通用有限元軟件包ANSYS。表1材料性能2邊界條件圖1顆粒隨機(jī)分布模型在單元胞體承受外載作用,胞體的各個(gè)表面都要在各自的法線(xiàn)方向上保

2、持變形一致,也即一個(gè)界面上的所有節(jié)點(diǎn)允許有位移,但每個(gè)節(jié)點(diǎn)的位移應(yīng)該是相同的。這樣,由胞體重復(fù)疊加的復(fù)合材料在變形上才會(huì)連續(xù)。具體量化要求如下:z=0,“3=0;z=口,盯=100MPa,拈,=以;菇=0,菇=口,“l(fā)與彳無(wú)關(guān);Y=0,Y=口,u2與z無(wú)關(guān)。即在上表面施加100MPa的應(yīng)力,%為單胞上表面的位移量,M。、“:、,分別為單元胞體在菇、幾石三個(gè)方向的位移。3應(yīng)力場(chǎng)分布對(duì)胞體進(jìn)行切片以研究其內(nèi)部的應(yīng)力場(chǎng)分布,圖2是邊長(zhǎng)為100la,m顆粒體積分?jǐn)?shù)15%的單元胞體內(nèi)部的VO/1Mises應(yīng)力分布。圖3是顆粒增強(qiáng)復(fù)合材料拉伸實(shí)驗(yàn)遠(yuǎn)離斷口區(qū)域的金相顯微鏡照片。從圖2可以看出,試件受到豎直方

3、向的拉力,在基體中,顆粒的上下邊緣的基體中,應(yīng)力較大,而在顆粒水平方向上應(yīng)力較小;在增強(qiáng)顆粒中,增強(qiáng)項(xiàng)中的應(yīng)力較基體中的應(yīng)力要大得多,顆粒中應(yīng)力最大處往往處于顆粒的邊緣;當(dāng)顆粒密度較小,顆粒間的距離較大,顆粒對(duì)基體的影響較小,應(yīng)力主要集中在顆粒內(nèi)部邊緣,對(duì)于脆性顆粒而言,其裂紋區(qū)域即為高應(yīng)力區(qū)域,顯微照片下的裂紋區(qū)域圖3(a (b與模擬的高應(yīng)力區(qū)域圖2(a(b中的紅色區(qū)域形貌一致。當(dāng)顆粒間距離較小,顆粒間基體的黃色區(qū)域顯示,顆粒與基體之間有應(yīng)力均化的趨勢(shì),基體的應(yīng)力與顆粒內(nèi)部的應(yīng)力幾乎相等,由于基體的強(qiáng)度比顆粒的強(qiáng)度要低,基體首先屈服,微裂紋萌生并擴(kuò)展。特別是在循環(huán)荷載作用下,圖3(c(d中基

4、體中的裂紋沿著高應(yīng)力區(qū)域方向擴(kuò)展,裂紋的形狀與圖2(c(d處基體中黃色高應(yīng)力區(qū)域的形狀一致。以往文獻(xiàn)中多數(shù)認(rèn)為在顆粒增強(qiáng)復(fù)合材料中,顆粒強(qiáng)度遠(yuǎn)遠(yuǎn)高于基體,復(fù)合材料的破壞主要發(fā)生在基體中。但實(shí)驗(yàn)測(cè)試,以及數(shù)值模擬分析顯示,復(fù)合材料微裂紋萌生的位置,由增強(qiáng)顆粒和基體的相對(duì)強(qiáng)度決定。增強(qiáng)顆粒的強(qiáng)度往往高于基體,但在受拉狀態(tài)下,顆粒內(nèi)部的應(yīng)力也高于基體,因此,顆粒與基體都存在破壞的可能。圖圖3 (a(b中的裂紋與圖2中的顆粒內(nèi)部的高應(yīng)力區(qū)域形狀及位置相同。另外,當(dāng)顆粒之間的距離比較接近,微裂紋產(chǎn)生在基體區(qū)域,如圖3(c(d所示。4彈性模量 Ju等人采用整體體積平均方法推導(dǎo)出顆粒增2726科學(xué)技術(shù)與工程

5、9卷圖2100且m胞體等距切片應(yīng)力分布圖強(qiáng)復(fù)合材料剛度系數(shù)f5。樊建平等人也做了單元胞體的數(shù)值模擬【6J。現(xiàn)運(yùn)用顆粒隨機(jī)分布胞體模型對(duì)顆粒增強(qiáng)復(fù)合材料的彈性模量進(jìn)行模擬。并與實(shí)驗(yàn)結(jié)果"1進(jìn)行比較。從圖5和文獻(xiàn)5中町以得到以下結(jié)果:(1利用顆粒隨機(jī)分布模型模擬的彈性模量與實(shí)驗(yàn)測(cè)試的復(fù)合材料彈性模量偏差很小。所有彈性模量的模擬值與實(shí)驗(yàn)測(cè)試得到的彈性模量各個(gè)數(shù)據(jù)點(diǎn)的平均偏差為2.16%。對(duì)于小于5%的誤差,即可認(rèn)為,計(jì)算模擬的彈性模量與實(shí)驗(yàn)測(cè)試彈性模量吻合比較好;(2隨著顆粒體積分?jǐn)?shù)的增加,單顆粒夾雜胞體單元對(duì)復(fù)合材料彈性模量的模擬與實(shí)驗(yàn)測(cè)試結(jié)果的偏差越來(lái)越大.顆粒體積分?jǐn)?shù)較小時(shí),顆粒問(wèn)的

6、間距較大,顆粒間相互影響很小,顆粒間的相互影響可以忽略,所以顆粒體積分?jǐn)?shù)小時(shí),單顆粒模型和多顆粒隨機(jī)模型與實(shí)驗(yàn)測(cè)試得到的彈性模量相差無(wú)幾;顆粒體積分?jǐn)?shù)較大,顆粒間的間距小,顆粒與顆粒的相互作用大,顆粒模型沒(méi)有考慮顆粒問(wèn)的相互干涉機(jī)制,使得在顆粒體積分?jǐn)?shù)較大時(shí),單顆粒模型模擬的彈性模量較實(shí)驗(yàn)測(cè)試值要小得多;顆粒體積分?jǐn)?shù)越大,單顆粒模型模擬的彈性模量與實(shí)驗(yàn)測(cè)試值的偏離越大;(3顆粒數(shù)量的增多能有效提高材料的彈性模量。與單顆粒單元胞體模型相比,既考慮了顆粒隨機(jī)分布對(duì)材料的影響,又考慮了相鄰顆粒之間的相互作用,計(jì)算模型比較符合實(shí)際情況,這是隨機(jī)顆粒模 型模擬結(jié)果優(yōu)于單顆粒模型的原因。10期單豪良,等:

7、顆粒增強(qiáng)復(fù)合材料隨機(jī)分布模型的數(shù)值模擬2727 存在顆粒界面脫粘等情況。顆粒隨機(jī)分布模型能較好地反應(yīng)顆粒分布的隨機(jī)性,也能充分考慮顆粒間的相互作用,反應(yīng)顆粒增強(qiáng)復(fù)合材料的微觀結(jié)構(gòu)圖3遠(yuǎn)離拉伸斷口區(qū)域的微裂紋 圖5復(fù)合材料彈性系數(shù)隨著顆粒體積分?jǐn)?shù)的增加,在復(fù)合材料加工制造階段,相互粘連、擠壓、團(tuán)結(jié)的顆粒增多,而且,還帆布,較準(zhǔn)確地模擬復(fù)合材料的彈性模量。5結(jié)論顆粒增強(qiáng)復(fù)合材料的隨機(jī)分布模型是符合實(shí)際情況的,其高應(yīng)力區(qū)域與實(shí)驗(yàn)觀測(cè)結(jié)果相吻合。234567參考文獻(xiàn)GuildFJ.Bonfield W.Predictive modeling of hydmxyapatitepolyetIlylene

8、composite.Biomatefials。1993;14:985-993BalacI。UskokoviePS,IgnjatovicN,eta1.Stress analysisinhydroxyapatite/polyLlactide compositebiomaterials.ComputationalMaterialsScience,2001;20:275-283U Z14,Sehmauder S,WannerA,n at.Expression tocharacterize theflow behaviorofparticlereinforcedcompositesbasedonaxis

9、ymmetric unitcell models.Scripta MetallurgicaetMaterialia,1995;33:1289-1294Zhang Jiangtao,“u Lisheng,Zhai Pengcheng.eta1.The predictionof the dynamicresponsesof ceramicparticlereinforcedMMCs byUsingmultiparticlecomputationalmieromechanical method.Composites ScienceandTechnology,2007;67:2775-2785Ju J

10、 W。ChenTM.Mieromechanics andeffectivemoduliof elasticcomposites containingrandomlydispersed ellipsoidal inhomeneitiesActaMechaniea。1994;103:103121樊建平,鄧澤賢.兩個(gè)單位胞體模型數(shù)值模擬結(jié)果比較.華中科技大學(xué)學(xué)報(bào),2003;31:8284柏振海,黎文獻(xiàn),羅兵輝.一種復(fù)合材料彈性模量的計(jì)算方法.中南大學(xué)學(xué)報(bào),2006;37:438443Particle Reinforced Composite Random Numerical SimulationMo

11、delSHAN Haoliang,GUO Rongxin,XIA Haiting,YANG Bangcheng,HE Tiancun (The facultyof ArchitecturalEngineering,Kunming Universityof Science andTechnology,Kunming 650024,P.R.ChinaAbstractUsing three.dimensional finite element method method,makeananalog computation by building thepellet random distributio

12、n compatriot body numerical value analysis model,to the particulate composite.Under the effect stretching load unidirectionally。discuss thatyonMises strain distribution law and elasticity modulus numerical value imitate in the model.With experiment bear fruit compare,discover the pellet random distribution compatriot body model elasticity modulus according with reality fa