朱 波,周 叮,劉偉慶
(南京工業(yè)大學(xué) 土木工程學(xué)院,南京 210009)
復(fù)合材料具有比強度高、比模量高、抗疲勞性好、抗震能力強、結(jié)構(gòu)可設(shè)計性等優(yōu)點,已在航天航空、汽車、建筑等各個領(lǐng)域得到了廣泛的應(yīng)用。復(fù)合材料夾層結(jié)構(gòu)由面板和芯材通過樹脂固化而成,在芯材上開鑿貫穿厚度的圓孔,采用真空導(dǎo)入工藝,芯材中含有的圓柱狀樹脂柱能夠更好地連接面板和芯材,提高面板與芯材之間的抗剝離和協(xié)調(diào)工作能力,起到類似于螺栓的作用,同時亦能提高芯材的受壓受剪性能,稱之為樹脂柱點陣增強型復(fù)合材料。
過去對復(fù)合材料的研究大多集中于靜力學(xué)分析和試驗,近年來對復(fù)合材料的屈曲及動力學(xué)性能的研究已取得不少成果,如:Chen和Yang[1]用一階剪切變形理論分析了層合梁的屈曲和自由振動,并用有限元法計算了層合梁的自振頻率;文獻[2]用解析法分析了考慮轉(zhuǎn)動和剪切效應(yīng)時復(fù)合材料夾層結(jié)構(gòu)的自振頻率。先進復(fù)合材料結(jié)構(gòu)精細(xì)、工藝復(fù)雜周期長,以傳統(tǒng)的方法研制費時費力,且難以得到準(zhǔn)確而系統(tǒng)的科學(xué)結(jié)論,有限元數(shù)值模擬技術(shù)具有試驗方法和理論解析不可比擬的優(yōu)勢,是研究結(jié)構(gòu)力學(xué)性能的重要分析工具,ANSYS是目前有限元算法中較為流行的商業(yè)軟件之一,已廣泛應(yīng)用于各種工程領(lǐng)域。
針對復(fù)合材料,ANSYS提供了一些特殊的殼單元及實體單元,目前對復(fù)合材料夾層結(jié)構(gòu)主要采用這類單元進行模擬和分析[3-9]。但是,對于點陣增強型復(fù)合材料夾層結(jié)構(gòu),需考慮芯層樹脂柱對整個結(jié)構(gòu)的影響,這在已有單元中未能得到很好的考慮。文[10]將芯層樹脂柱和泡沫分開劃分單元,采用單胞概念,對含樹脂柱的夾層板進行靜力學(xué)分析,研究了樹脂柱的分布和材性對層間界面應(yīng)力分布及板豎向位移的影響。
本文對含有樹脂柱的點陣增強型復(fù)合材料夾層板建立單胞分析模型,將芯層樹脂柱與芯材分開,采用8節(jié)點SOLID45實體單元,應(yīng)用ANSYS有限元軟件研究樹脂柱對夾層板的屈曲以及動力學(xué)性能的影響。
圖1所示為點陣增強型復(fù)合材料夾層板模型示意圖,上下兩層為面板,厚度為h1,中間芯層厚h2。對于點陣增強型輕質(zhì)芯材,可認(rèn)為芯材是由一系列單胞在空間排列組成,單胞由芯材基體及樹脂柱共同構(gòu)成,樹脂柱呈圓柱狀,在宏觀結(jié)構(gòu)中呈規(guī)則正交分布。每個單胞的邊長為b,單胞中所含樹脂柱直徑為d。假設(shè)面板、芯材、樹脂柱三者之間無粘結(jié)滑移,且未發(fā)生剝離現(xiàn)象,因此面板單元、芯材單元、樹脂柱單元之間可采用位移協(xié)調(diào)式,各單元共用節(jié)點。
考慮一矩形復(fù)合材料夾層板,其尺寸為:長L1=800 mm,寬L2=320 mm,上下面板厚h1=5 mm,芯材厚h2=50 mm,單胞邊長b=40 mm,樹脂柱直徑d=8 mm。
圖1 點陣增強型復(fù)合材料夾層板示意圖Fig.1 The schematic diagram of a sandwich plate with grooved perforation
纖維增強樹脂基面板以平面內(nèi)受力為主,采用四軸向[0/45/90/-45]的準(zhǔn)正交玻璃纖維布,將面板定義為各向同性材料,并與一般鋪設(shè)情況進行比較,對比結(jié)果見表1,誤差在0.8%以內(nèi),可以忽略,故本文在后續(xù)分析中均將面板定義為各向同性材料以方便建立模型。芯材選用我國量大面廣的泡桐木,面板、芯材及樹脂柱的材料性質(zhì)列于表2。本文研究不同邊界條件下的夾層板:① 四邊簡支;② 兩短邊簡支,另兩邊自由;③ 一短邊固支,其余三邊自由。采用8節(jié)點SOLID45三維實體單元,劃分網(wǎng)格后的復(fù)合材料夾層板模型以及單胞局部網(wǎng)格劃分示意圖分別如圖2、圖3所示。
表1 各向同性與正交鋪設(shè)時夾層板計算結(jié)果的比較Tab.1 Comparisons of isotropic and orthotropic results for a sandwich plate
表2 復(fù)合夾層板各組分材料性質(zhì)Tab.2 Material properties of each component of the composite sandwich plate
圖2 復(fù)合材料矩形板ANSYS模型Fig.2 The meshing model of a sandwich plate
圖3 單胞局部網(wǎng)格劃分示意圖Fig.3 The meshing model of one cell
考慮四邊簡支夾層板,承受平面內(nèi)軸向載荷。保持單胞和樹脂柱尺寸不變,改變樹脂柱的彈性模量,表3給出了不同樹脂柱彈模下夾層板的屈曲臨界載荷。
表3 不同樹脂柱彈模下夾層板的屈曲臨界載荷Tab.3 Critical buckling loads of sandwich plates for different elastic modulus of resinic columns
從表中可以看出,樹脂柱彈模越高,板的屈曲載荷值越大,然而,樹脂柱彈模對復(fù)合材料板屈曲性能的影響并不顯著。由表可知,樹脂柱的彈模提高3倍(1500 MPa到4500 MPa),屈曲臨界載荷只提高了約2.6%,影響有限。
保持芯材樹脂柱直徑與單胞邊長之比d/b=1/5不變,即樹脂柱在芯層中的體積含量一定。改變芯層中樹脂柱分布的疏密程度,即改變單胞底面邊長b,表4給出了不同樹脂柱間距對屈曲載荷的影響,其中b=∞表示芯層無樹脂柱。從表中可以看出,當(dāng)樹脂柱間距變密(即b由160 mm變?yōu)?0 mm時),板的臨界載荷提高了11.8%,比芯層無樹脂柱時增大了13.7%,說明樹脂柱的添加可提高復(fù)合材料夾層板的軸向承載能力,且樹脂柱越密集,軸向承載力越大。
表4 不同樹脂柱間距下夾層板的屈曲臨界載荷Tab.4 Critical buckling loads of sandwich plates for different cell sizes
材料性質(zhì)及模型尺寸如前所述,研究兩對邊簡支另兩對邊自由夾層板的無阻尼自由振動。改變樹脂柱的彈性模量,表5給出了樹脂柱彈模對板前五階固有頻率的影響??梢钥闯?,樹脂柱彈模越高,夾層板的固有頻率越大,彈性模量的變化對各階頻率的影響大致相等,如:樹脂柱的彈模提高3倍(1500 MPa到4500 MPa),一至五階頻率均提高了約1.4%,可見樹脂柱彈性模量的變化對夾層板自由振動特性的影響較小,且自振頻率隨樹脂柱彈模的增加而呈緩慢線性增長趨勢。
表5 不同樹脂柱彈模下夾層板的前五階自振頻率Tab.5 First five natural frequencies of the sandwich plate for different elastic modulus of resinic columns
保持樹脂柱在單胞中所占比例d/b=1/5不變,即樹脂柱在芯層中的體積含量一定。改變芯層中樹脂柱的疏密程度,即改變相鄰樹脂柱的間距b。表6給出了不同樹脂柱疏密度下夾層板的前五階固有頻率,其中b=∞表示芯層無樹脂柱。由表可知,隨著樹脂柱間距的變密,板的各階固有頻率均隨之增大,且低階頻率的變化更為顯著,如:當(dāng)b由160 mm變?yōu)?0 mm時,板的一階頻率提高了28.1%,比芯層無樹脂柱時增大了29.2%;第五階頻率提高了24.5%,比芯層無樹脂柱時增大了25.9%。可見,樹脂柱分布疏密程度對夾層板自振頻率有明顯影響。圖4-圖6給出了b=40 mm時夾層板的前三階自由振動模態(tài)。
表6 不同樹脂柱間距下夾層板的前五階自振頻率Tab.6 First five frequencies of sandwich plate for different cell sizes
考慮一短邊固支,其余各邊自由的夾層板,自由短邊的中點承受在有限上升時間內(nèi)的平面外瞬態(tài)沖擊力,從t=0時刻開始,載荷初始值為0,當(dāng)t=0.075 s時達到最大值20 kN,加載結(jié)束時間為t=1 s,加載過程如圖7所示。
圖4 夾層板自由振動的第一階模態(tài)Fig.4 Fundamental mode shape of the sandwich plate
圖5 夾層板自由振動的第二階模態(tài)Fig.5 Second mode shape of the sandwich plate
圖6 夾層板自由振動的第三階模態(tài)Fig.6 Third mode shape of the sandwich plate
圖7 瞬態(tài)荷載1秒內(nèi)的加載過程Fig.7 Transient loading process within 1 second
圖8 夾層板懸臂端中點處豎向位移隨時間的變化曲線Fig.8 Vertical displacement at cantilever end of the plate
圖9 不同樹脂柱間距下板中線下沿的最大振幅曲線Fig.9 Maximum amplitude of the plate for different cell sizes
考慮樹脂的粘彈性性質(zhì)及周圍空氣阻力的影響,引入阻尼,取等效粘性阻尼系數(shù)c=8 kN·s/m。圖8給出了b=40 mm及b=80 mm時,懸臂端中點處豎向位移隨加載時間的變化曲線,可以看出,當(dāng) t=0.088 s時位移達到最大值。圖9給出了不同樹脂柱間距下,t=0.088 s時夾層板中線下沿的豎向振幅曲線。從圖中可以看出,懸臂端振幅最大,且隨樹脂柱間距變密而減小,當(dāng)b由80 mm變?yōu)?0 mm時,懸臂端振幅減小了約6.9%,比無樹脂柱時降低了約8.7%??梢?,在瞬態(tài)荷載作用下,樹脂柱對于降低夾層板的豎向動力位移有一定的作用。
本文應(yīng)用ANSYS有限元軟件分析了樹脂柱點陣增強輕木芯材復(fù)合面板屈曲載荷、自由振動以及瞬態(tài)響應(yīng)的影響。分析結(jié)果表明,樹脂柱對板的力學(xué)性能有一定的貢獻,其中材料特性對夾層板屈曲及動力學(xué)特性的影響有限,可以忽略;而芯層樹脂柱分布疏密程度的變化對點陣增強型夾層板的屈曲及動力學(xué)特性有較大影響;此外,樹脂柱的添加對于減小瞬態(tài)荷載作用下夾層板的豎向動力位移有一定的積極作用??梢?,為提高夾層板的力學(xué)性能,合理安排樹脂柱的分布要比單純增大其彈性模量效果更好,研究結(jié)果可以為點陣增強型復(fù)合材料夾層板的參數(shù)設(shè)計提供依據(jù)。
[1]Chen A T,Yang T Y.Static and dynamic formulation of symmetrically laminated beam finite element for a microcomputer[J].J Comp Mat,1995,19:459 -475.
[2]Chandrashekhara K,Krishnamurthy K,Roy S.Free vibration of composite beamsincluding rotary inertia and shear deformation [J]. Composite Structures, 1990, 14:269-279.
[3]況 祺,夏凌輝,常春偉.ANSYS在復(fù)合材料夾層結(jié)構(gòu)屈曲分析中的應(yīng)用[J].科學(xué)技術(shù)與工程,2005,5(32):1836-1838.
[4]石建軍,吳東輝,遲 波,等.Shell91單元在復(fù)合材料蜂窩夾層結(jié)構(gòu)分析中的應(yīng)用[J].纖維復(fù)合材料,2006,3:40-42.
[5]李 濤,樊慶文.ANSYS在復(fù)合材料應(yīng)力分析中的應(yīng)用[J].機械,2006,12(33):47-48.
[6]廖英強,蘇建河,柯善良,等.ANSYS在復(fù)合材料仿真分析中的應(yīng)用[J].玻璃纖維,2006,4:21-26.
[7]廖 原.ANSYS在復(fù)合材料梁設(shè)計中的應(yīng)用[J].廣東建材,2007,12:96-98.
[8]Zhang S H,Chen H L.A study on the damping characteristics of laminated composites with integral viscoelastic layers[J].Composite Structures,2006,1(74):63 -69.
[9]Barbero E J, Lonetti P, Sikkil K K. Finite element continuum damage modeling ofplain weave reinforced composites[J].Composites Part B:Engineering,2006,2-3(37):131-147.
[10]姚秀冬,周 叮,劉偉慶.復(fù)合材料夾層板樹脂柱的影響[J].材料科學(xué)與工程學(xué)報,2009,27(6):1-6.