復(fù)合材料力學(xué)

1、目錄TOC o 1-5 h z HYPERLINK l bookmark2 o Current Document 復(fù)合材料細(xì)觀力學(xué)1 HYPERLINK l bookmark22 o Current Document 簡支層合板的自由振動9 HYPERLINK l bookmark50 o Current Document 不同條件下對稱層合板的彎曲分析14 復(fù)合材料細(xì)觀力學(xué)混凝土細(xì)觀力學(xué)一、研究背景復(fù)合材料細(xì)觀力學(xué)復(fù)合材料細(xì)觀力學(xué)是20世紀(jì)力學(xué)領(lǐng)域重要的科學(xué)研究成果之一，是連續(xù)介質(zhì)力學(xué)和材料科學(xué)相互衍生形成的新興學(xué)科。近20年來，我國科技工作者應(yīng)用材料細(xì)觀力學(xué)的理論和方法，成功研究了許多復(fù)合

2、材料的增強(qiáng)，斷裂和破壞問題，給出了一些特色和有價值的研究成果?；炷良?xì)觀力學(xué)混凝土作為一種重要的建筑材料已有百余年的歷史，它廣泛應(yīng)用于房屋、橋梁、道路、礦井、及軍工等諸多方面。在水工建筑方面，混凝土也被大量使用，特別是大體積混凝土，它是重力壩和拱壩的主要組成部分，對混凝土各項力學(xué)性能的準(zhǔn)確把握及應(yīng)用，在一定程度上決定了水工建筑物的質(zhì)量和安全性能。二、研究目的長期以來，在混凝土應(yīng)用的各個領(lǐng)域里，人們對混凝土的力學(xué)特性進(jìn)行了大量的研究。如何充分的利用混凝土的力學(xué)性能，建造出更經(jīng)濟(jì)、更安全和更合理的建筑物或工程結(jié)構(gòu)，一直都是結(jié)構(gòu)工程設(shè)計領(lǐng)域研究的重要課題。三、研究現(xiàn)狀混凝土是由粗骨料和水泥砂漿組成的

3、非均質(zhì)材料，它的力學(xué)性能受到材料的品質(zhì)、組分、施工工藝和使用條件等因素的影響。過去，人們對混凝土力學(xué)性能的研究很大程度上是依靠實驗來確定的。隨著實驗技術(shù)的發(fā)展，混凝土各種力學(xué)性能被揭示出來。但由于實驗需要花費大量的人力、物力和財力，而且所得到的實驗成果往往由于實驗條件的限制也是很有限的?，F(xiàn)代科學(xué)的一個重要的思維方式與研究方法就是層次方法，在對客觀世界的研究中，當(dāng)停留在某一層次，許多問題無法解決時，深入到下一個層次，問題就會迎刃而解。對混凝土斷裂問題的研究歸納為如下四個研究層次：宏觀層次：混凝土這種非均質(zhì)材料存在著一個特征體積，經(jīng)驗的特征體積相應(yīng)于34倍的最大骨料體積。當(dāng)混凝土體積大于這種特征體

4、積時，材料被假定為均質(zhì)的，當(dāng)小于這種特征體積時，材料的非均質(zhì)性將會十分明顯。有限元計算結(jié)果反映了一定體積內(nèi)的平均效應(yīng)，這個特征體積的平均應(yīng)力和平均應(yīng)變稱之謂宏觀應(yīng)力和宏觀應(yīng)變。細(xì)觀層次：在這個層次中，混凝土被認(rèn)為是一種由骨料、砂漿和它們之間的粘結(jié)帶組成的三相非均質(zhì)復(fù)合材料，細(xì)觀內(nèi)部裂隙的發(fā)展將直接影響混凝土的宏觀力學(xué)性。細(xì)觀層次的模型一般是毫米或厘米量級。微觀層次：在這個層次上，認(rèn)為砂漿的非均質(zhì)性是由漿體中的孔隙所產(chǎn)生的。由于砂漿中孔隙很小而且量多，隨機(jī)分布，水泥砂漿可近似看作細(xì)觀均質(zhì)損傷體。同配合比、同條件的砂漿試件，通常其力學(xué)性能也比較穩(wěn)定，可由試驗直接測定。納觀層次：混凝土的納觀層次認(rèn)為

5、水泥漿體本身也是一種隨機(jī)的復(fù)合材料，其組分是未水化的熟料顆粒、水化硅酸鈣、氫氧化鈣晶體、毛細(xì)孔和其他的化學(xué)成分，該層次毛細(xì)孔的尺寸是微米或亞微米量級。在過去長期的研究中，人們?yōu)榱搜芯可系姆奖?，往往忽略混凝土?fù)雜的內(nèi)部結(jié)構(gòu)，在對材料宏觀本構(gòu)關(guān)系的描述中，這種材料是“黑匣子”，采用唯象學(xué)的觀點，把它們平均化和均勻化為宏觀均勻連續(xù)體，并以實驗結(jié)果為基礎(chǔ)發(fā)展了彈性、彈塑性以及黏彈塑性的混凝土本構(gòu)模型。這種方法能夠使人們對工程結(jié)構(gòu)的均勻化狀態(tài)有一個總體上的認(rèn)識，可以作為工程設(shè)計的依據(jù)和參考。但是，這些模型無法深入了解混凝土在外力作用下內(nèi)部微裂紋萌生、擴(kuò)展及其貫通，直至宏觀裂紋形成，促使試樣失穩(wěn)破裂的整個

6、過程，更無法反映混凝土斷裂過程中表現(xiàn)出來的變形局部化和應(yīng)力重分布等基本特征。四、研究方法事實上，材料的宏觀斷裂過程必然與其細(xì)觀的非均勻結(jié)構(gòu)是密切相關(guān)的，因而，進(jìn)行材料微、細(xì)觀結(jié)構(gòu)的模擬對于了解混凝土宏觀斷裂機(jī)理是非常有意義的。進(jìn)行混凝土材料的斷裂過程以及破壞機(jī)理的研究，必須從混凝土的細(xì)觀結(jié)構(gòu)入手，利用細(xì)觀力學(xué)的研究方法，抓住混凝土材料及其力學(xué)性質(zhì)的非均勻性，借助統(tǒng)計學(xué)和數(shù)值計算方法，以實驗結(jié)果為基礎(chǔ)建立數(shù)值模型，進(jìn)行混凝土斷裂過程的數(shù)值模擬。從混凝土試樣的拋光斷面上（如下圖1），很容易鑒別的兩個相是不同尺寸、形狀的骨料顆粒以及分布紊亂的硬化水泥漿體物質(zhì)所組成的膠凝性介質(zhì)，因此從宏觀水平而論，混

7、凝土可視為由骨料顆粒分散在水泥砂漿基質(zhì)中所組成的材料。圖1混凝土試樣的拋光斷面上從細(xì)觀水平而言，結(jié)構(gòu)的兩個相既不是彼此均勻分布，其本身也不是均勻的，顯示出一定的非均質(zhì)復(fù)雜性。事實上，在混凝土實驗中觀察到的很多現(xiàn)象以及混凝土在應(yīng)力作用下行為的許多方面只能將水泥砂漿與粗骨料的界面視為混凝土結(jié)構(gòu)的第三相才能作出解釋。因此在細(xì)觀層次上，混凝土被認(rèn)為是由粗骨料、砂漿及二者間的粘結(jié)帶組成的三相非均質(zhì)復(fù)合材料，其中粘結(jié)帶為圍繞大骨料周圍存在的一層薄殼。在細(xì)觀層次上，把混凝土看作是由粗骨料、砂漿及二者問的粘結(jié)帶組成的三相非均質(zhì)復(fù)合材料，研究混凝土內(nèi)裂縫形成、擴(kuò)展，追溯混凝土細(xì)觀單元從損傷至斷裂的全過程，預(yù)測材

8、料的破壞抗力等力學(xué)指標(biāo)，建立混凝土的細(xì)觀結(jié)構(gòu)與宏觀性能之間的關(guān)系。其工作主要分為以下幾個部分：1)分析混凝土細(xì)觀結(jié)構(gòu)組成，采用蒙特卡羅隨機(jī)取樣原理，編程產(chǎn)生一個在形狀、尺寸及粗骨料的空間分布上都非常類似于真實混凝土的混凝土隨機(jī)骨料結(jié)構(gòu)。2)結(jié)合VisualFortran、AutoCAD及ANSYS等計算機(jī)軟件工具，對產(chǎn)生的隨機(jī)骨料結(jié)構(gòu)劃分有限元網(wǎng)格，生成可控制粘結(jié)帶厚度的隨機(jī)骨料結(jié)構(gòu)有限元網(wǎng)格。3)用無量綱化指標(biāo)建立細(xì)觀層次均質(zhì)水泥砂漿損傷體的單軸拉壓損傷本構(gòu)關(guān)系，將損傷耦合到混凝土受拉和受壓的本構(gòu)方程中，把損傷力學(xué)和斷裂力學(xué)結(jié)合起來進(jìn)行分析和計算。4)考慮到細(xì)觀層次分析的特殊性，提出結(jié)合抗拉

9、強(qiáng)度和斷裂韌性的綜合破壞準(zhǔn)則，便于采用連續(xù)介質(zhì)有限元分析細(xì)觀層次的斷裂問題，并提出了適用于細(xì)觀混凝土力學(xué)分析的單元開裂控制法和損傷收斂準(zhǔn)則。開發(fā)軟件的后處理功能，可動態(tài)地模擬試件破壞的全過程，輸出任一時刻的變形、裂縫圖像，使得細(xì)觀斷裂分析全過程可以借助圖像顯示直觀的表現(xiàn)出來。應(yīng)用隨機(jī)骨料結(jié)構(gòu)，采用非線性有限單元法進(jìn)行混凝土單軸受拉和單軸受壓實驗全過程的數(shù)值模擬，得到混凝土數(shù)值試件的單軸抗拉和抗壓強(qiáng)度。對細(xì)觀混凝土數(shù)值試件進(jìn)行雙向受力分析，得到混凝土在雙軸載荷作用下的強(qiáng)度包絡(luò)線，并分析混凝土在雙向受力情況下的破壞模式及應(yīng)力場。對細(xì)觀混凝土骨料結(jié)構(gòu)隨機(jī)性進(jìn)行分析，證明骨料結(jié)構(gòu)的隨機(jī)性與混凝土材料本

10、身的離散性具有一致性，進(jìn)一步說明了隨機(jī)骨料結(jié)構(gòu)用于模擬混凝土實驗的合理性。五、結(jié)論與展望建立細(xì)觀力學(xué)模型，通過小量的實驗研究就可以較好地預(yù)測混凝土的各種宏觀力學(xué)性能，進(jìn)而取代大部分的混凝土材料實驗，特別是實驗中無法得到的大骨料混凝土的力學(xué)性能，從而大大地減少實驗研究的工作量，具有潛在的顯著的經(jīng)濟(jì)和社會效益?；炷两Y(jié)構(gòu)的研究從宏觀到細(xì)觀，雖然在科學(xué)探索的道路上邁出可貴的步伐，但混凝土細(xì)觀結(jié)構(gòu)的一些復(fù)雜現(xiàn)象，如水灰比對界面性質(zhì)的影響、水泥砂漿中氣泡的形成和水泥結(jié)石的形成原理等一些現(xiàn)象，在細(xì)觀層次上還解釋不清楚，要揭示更深層次的機(jī)理，也只有迸入到下一層次的研究中去，人類對混凝土物質(zhì)結(jié)構(gòu)及其規(guī)律的深入

11、了解，必將對改變混凝土上一層次的性能，特別是工程宏觀性能，向著按預(yù)定性能設(shè)計新材料，起到深遠(yuǎn)的影響和巨大的作用。參考文獻(xiàn)邱志章,混凝土細(xì)觀力學(xué)分析與實驗?zāi)M張研，張子明,材料細(xì)觀力學(xué)王宗敏，不均質(zhì)材料（混凝土）裂隙擴(kuò)展及宏觀計算強(qiáng)度與變形于驍中，居襄混凝土的強(qiáng)度和破壞簡支層合板的自由振動這部分主要研究了四邊簡支條件下層合板的自由振動。求出了其理論解和ANSYS數(shù)值解，并做了誤差分析。一、問題四邊簡支的正交鋪設(shè)對稱矩形層合板，鋪設(shè)角度為90/0/90/0/90/0/90單層厚度為0.2e-3m,a=0.6m,b=0.3m。材料參數(shù)E=208GPa，E=18.9GPa,G=5.7GPa,V=0.2

12、3,密度為1212122700kg/m3,考慮層合板在慣性力作用下的自由振動。二、MATLAB求理論解對中間鋪設(shè)角度為0時的層合板的解進(jìn)行驗證：因為Bj0，*二A26二D16二D26二0，振動頻率和振型由下列振動方程描述：（1）邊界條件為：x二0,a:6w二0,6M二D6wD6w二0 x11,xx12,yy（2）y二0,b:6w二0,6M=D6wD6w=0y12,xx22,yy選取：6（x,y,t）=（Acoset+Bsinet）6（x,y）（3）將此問題分為時間和空間兩部分。為使式（3）滿足方程（1）和邊界條件（2），進(jìn)一步選取sin;m兀x.sinaab(x,y,t)二(Acost+Bsi

13、nt)sin將上式帶入方程（1）得2D11+2(Di2+2%:Ib)2+D22(b)式中，各頻率對應(yīng)于不同振型MATLAB編程求解：4MATLABR2014a繪圖通亞口L碑:件&田詈腳本導(dǎo)入工作區(qū)文件LJ國回口.號清除工作區(qū)血清歸令Simulink庫SIMUUNK二舷徑j(luò)iilParallel環(huán)境布局*宙0律C:UsersLitongDocumentsMATLAB當(dāng)前文件夾命令伍口名稱4程|EllipsoidFitting.m：inputl.txttniheshuju.datoutput.txt12S矩陣.txt亡5果.txt選擇文件以查看詳細(xì)信息t=0.2e-3;a=0.6;b=0.3;p=

14、2700;El=208e9:E2=18.9e9:G12=5.7e9;vl2=0.23;v21=0.021;QU=El/(l-vl2v21):Q12=v21*E2/(1-vl2v21):Q22=E2/仃-vl2*v21):Q66=G12:DI1=(t*3)/12(99*Q11+244Q22):D12=343*(t3)/12*Q12;D22=(t3)/12*(99*Q22+244*Q11);D66=343*(t3)/12*Q66:form=1:10w=sqrt(pi4/p*(Dll*(m/a)4+2*(D12+2*D66)(/a)2*(2/b)*2+D22*(2/b)*4):disp(w):end

15、50.407951.803155.760564.031177.785797.3074122.3714152.6496187.8727fx三、ANSYS求數(shù)值解1、選取單元2、材料屬性ChoosePoissonJsRatio1TemperattiresT1FEXEY1.89E+010EZ1.89E+010PRXY0.23PRYZ0.021pm0.23GXY5.7E+009GYZ5.7E+009GXZ5.7E+0094、鋪層設(shè)置5、劃分網(wǎng)格并選擇分析類型7、列出10階固有頻率ASET,LISTCommandFile*INDEXOFDATASETSONRESULTSFILE*SETTIME/FREQ

16、LOADSTEPSUBSTEPCUMULATIUE155.392111269.1041223104.971334164.011445214.591556221.481667239.721778244.081889276.411998、第8階振型結(jié)果10四、誤差分析對于中間層是0度角的鋪層計算結(jié)果對比w-wn二0X100%二8.9%RiEIci田orcJlECTNGCreateRectanglebyDimensiXI,X2X-coordinatesPolygon洞AreaFilletVoluaesBodesEleaentsContactPairFipingIodelsCircuitTransdu

17、cerserateFe/lodiypy9A旦a開始MicrosoftPower.復(fù)合材料力學(xué)作業(yè)加2ANSYSProgramH；3.bmp畫圖CH?|PickamenuitemorenteranAHSYSCommandCPREP7)對模型進(jìn)行網(wǎng)格劃分。施加邊界條件，四邊固定邊界條件。施加載荷，設(shè)置為200N的均布載荷。FJANSYSlultiphysicsUtilityIenu(esh)FileSelectListPlotPlotCtrlsWorkFleltlsPar：ini&terEMacroMenuCtrlsHelp|Q|q|冒1tiprF|壟阿ANSYSToolbar|SAVE.DBRE

18、S1JM_DB|QUITFOWRGRPHANSYSMainMenu1mces:essorm.rsisTypeleLoadsttingsplyStructxiral田Displaceaent莊1Force/loaentELEMENTSANSYSHApplyPRESonareaspr刀essureOnLinesOnModesOn.ModeCoaponenOnEleaentsOnEleaentCoapoFroaFluid.Anal;On.BeaasTeperature田InertiaHFretnsnSectn囲GenPlaneStrain.OtherField.VoluoeIntrInitialCo

19、nditJnLoadVectorFuuctioueLeteIJjOKIApplyCancelIHelp型型血型昌型！2冒U創(chuàng)011J0JJ0J1AQMAAJH_1I2N開始MicrosoftPower.Ifii2WindowsExpl.c：、ANSYS10.0Outp.QANSYSMultiphys.1ELEMENTS*1ir-_iZL-=nrL-=nir-.iANSYSJAN1Z00700:ZO:07ii廠L-=ni計算求解，并得出結(jié)果云圖。2、短邊固定，長邊簡支這里我們直接使用上一步操作里保存的劃分好網(wǎng)格的模型，對其施加邊界條件，短邊固定，長邊簡支。施加載荷并求解得出結(jié)果圖。1ELEMENTS111LJfit*1L-l