基于神經(jīng)網(wǎng)絡(luò)的橡膠混凝土抗凍性預(yù)測

徐金花， 陳四利， 董冠男， 岳子建

(沈陽工業(yè)大學(xué) 建筑與土木工程學(xué)院， 沈陽 110870)

橡膠混凝土作為一種新型的綠色建材，具有自重輕、變形能力強(qiáng)等優(yōu)勢，有較好的抗裂性和抗沖擊性，尤其在抗?jié)B和抗凍等耐久性方面具有良好的性能[1-3].為了將這種混凝土更好地在我國季節(jié)性冰凍地區(qū)推廣使用，很有必要對其凍融破壞性能和抗凍性進(jìn)行預(yù)測研究.有關(guān)普通混凝土抗凍性的預(yù)測，很多學(xué)者建立了以經(jīng)典理論為基礎(chǔ)的理論模型和以試驗數(shù)據(jù)做支撐的經(jīng)驗?zāi)Ｐ?，但目前學(xué)術(shù)界對以上模型均未完全達(dá)成共識[4].一般認(rèn)為，混凝土抗凍性是一個受多因素影響的復(fù)雜性和爭議性問題，涉及到許多不確定的影響因素，建立理論模型要用到材料學(xué)、力學(xué)、化學(xué)等多個學(xué)科的理論知識，還要結(jié)合亞微觀和宏觀兩個層面來研究，就目前來說，建立理論模型的難度相當(dāng)大.而目前的經(jīng)驗?zāi)Ｐ投嘁詢鋈跈C(jī)理為宏觀指導(dǎo)，考慮影響混凝土抗凍性的主要因素，以試驗數(shù)據(jù)為基礎(chǔ)進(jìn)行回歸分析，建立經(jīng)驗?zāi)Ｐ?，但模型的實用性有待商榷[5].相對于普通水泥混凝土而言，橡膠混凝土中由于摻入了橡膠粉，混凝土的各種性質(zhì)都發(fā)生了不同程度的改變，使其抗凍耐久性問題也變得更為復(fù)雜，建立理論模型的難度就更大[6].人工神經(jīng)網(wǎng)絡(luò)模型由于強(qiáng)大的非線性映射能力和自適應(yīng)學(xué)習(xí)、記憶的特點，使其非常適合解決混凝土抗碳化性能、抗化學(xué)侵蝕性能和抗凍性預(yù)測等一系列耐久性問題[7-8]，因此，本文將運用人工神經(jīng)網(wǎng)絡(luò)對橡膠混凝土的抗凍性進(jìn)行分析.為獲得足夠數(shù)量的訓(xùn)練樣本和測試樣本進(jìn)行了橡膠混凝土的凍融試驗，以相對動彈性模量作為衡量抗凍性的重要指標(biāo).選用MATLAB神經(jīng)網(wǎng)絡(luò)工具箱中的BP網(wǎng)絡(luò)模型，以試驗所得數(shù)據(jù)為基礎(chǔ)，建立橡膠混凝土的相對動彈性模量BP網(wǎng)絡(luò)預(yù)測模型，對其抗凍性進(jìn)行預(yù)測，為工程應(yīng)用提供參考.

1 橡膠混凝土凍融試驗

1.1 試驗原材料

本次凍融試驗采用的橡膠集料為沈陽天禹橡膠粉廠生產(chǎn)的50目橡膠粉，表觀密度為1 020 kg/m3.所用水泥為沈陽產(chǎn)PO42.5級普通硅酸鹽水泥.細(xì)骨料采用普通河砂，細(xì)度模數(shù)為2.55，表觀密度為2 600 kg/m3.粗骨料采用最大粒徑為25 mm的碎石，進(jìn)行5～25 mm連續(xù)級配，表觀密度為2 650 kg/m3.

1.2 試驗配合比

在橡膠混凝土抗凍性試驗中，配合比設(shè)計采用正交試驗設(shè)計方法.以C40級普通混凝土的配合比為基礎(chǔ)，運用常規(guī)配合比設(shè)計方法確定每立方米混凝土中各組成材料的質(zhì)量，然后分別用不同摻量粒徑為50目的橡膠粉等體積替換細(xì)骨料砂.采用橡膠摻量WR、水灰比W/B和砂率Sp3個重要因素的4個水平進(jìn)行設(shè)計，共得到16組不同配合比試件.橡膠摻量WR取5%、10%、15%和20%，分別用A1、A2、A3和A4來表示；水灰比W/B取0.35、0.40、0.45和0.50，分別用B1、B2、B3和B4來表示；砂率Sp取0.33、0.34、0.35和0.36，分別用C1、C2、C3和C4來表示.具體的配合比信息詳見文獻(xiàn)[6].

1.3 試驗方法

試驗采用《普通混凝土長期性能和耐久性能試驗方法標(biāo)準(zhǔn)》(GB/T50082-2009)中的快凍法[9]，試件尺寸為100 mm×100 mm×400 mm.通常來說，動彈性模量對混凝土結(jié)構(gòu)在凍融循環(huán)作用下的內(nèi)部損傷尤為敏感，而且可以運用動彈儀進(jìn)行無損檢測來測定，可減少大量的試驗工作量，節(jié)約試驗試件和成本[10].因此，該試驗中以相對動彈性模量作為分析和評價橡膠混凝土抗凍性能的重要參數(shù).圖1為運用動彈儀對凍融循環(huán)每隔25次的試件進(jìn)行動彈性模量的測定.

圖1 凍融循環(huán)試驗后測定試件的動彈性模量Fig.1 Dynamic elasticity modulus of samples after freeze-thaw cycle tests

1.4 試驗結(jié)果

經(jīng)過不同凍融循環(huán)次數(shù)后，按照上述試驗方法測定并計算得到各組橡膠混凝土試件的相對動彈性模量，如表1所示.通過對上述試驗結(jié)果進(jìn)行極差分析和方差分析可知，在選取的橡膠摻量WR、水灰比W/B和砂率Sp這三個試驗因素中，橡膠摻量和水灰比對動彈性模量影響顯著，砂率的影響很小.

表1 凍融循環(huán)作用下橡膠粉混凝土相對動彈性模量Tab.1 Relative dynamic elasticity modulus of rubber powder concrete subjected to freeze-thaw cycles %

2 橡膠混凝土樣本

2.1 樣本的收集和分類

以凍融試驗獲得的摻入50目橡膠粉的混凝土相對動彈性模量試驗數(shù)據(jù)為研究依據(jù)，考慮不同橡膠摻量、水灰比和凍融循環(huán)次數(shù)等因素變化帶來的影響，對其相對動彈性模量進(jìn)行分析和預(yù)測，共收集到139組試驗數(shù)據(jù).將這些數(shù)據(jù)組進(jìn)行隨機(jī)排列，得到整個數(shù)據(jù)樣本集.確定前112組作為網(wǎng)絡(luò)模型的訓(xùn)練樣本(約占整個數(shù)據(jù)樣本的80%)，后27組作為網(wǎng)絡(luò)模型的測試樣本(約占整個數(shù)據(jù)樣本的20%).考慮到各組試件采用的砂率均在合理范圍內(nèi)，且其變化對各組試件的動彈性模量影響不顯著，所以忽略砂率這一次要因素，確定BP網(wǎng)絡(luò)的輸入?yún)?shù)分別為橡膠摻量x1、水灰比x2和凍融循環(huán)次數(shù)x3，輸出參數(shù)為相對動彈性模量y.各參數(shù)取值為：橡膠摻量x1為5%、10%、15%和20%；水灰比x2為0.35、0.40、0.45和0.50；凍融循環(huán)次數(shù)x3為0、25、50、…、200；相對動彈性模量y為50%～100%.

2.2 樣本數(shù)據(jù)的歸一化

由于凍融循環(huán)次數(shù)和相對動彈性模量分布范圍較大，為了避免因某個凈輸入絕對值過大而使BP神經(jīng)元輸出飽和，使權(quán)值調(diào)整進(jìn)入平坦區(qū)，也為了避免某些數(shù)值很小但卻對輸出影響較大的向量被忽略，導(dǎo)致因輸入?yún)?shù)不均勻而引起的輸出失真和網(wǎng)絡(luò)麻痹現(xiàn)象，需要對原始數(shù)據(jù)進(jìn)行尺度變換處理，即數(shù)據(jù)的歸一化處理.本文采用式(1)對輸入變量x3和輸出變量y進(jìn)行歸一化處理，將輸入及輸出參數(shù)變換在[0.2，0.8]的區(qū)間內(nèi)，即

(1)

式中：m′為歸一化處理后的值；m為原始數(shù)據(jù)真實值；mmin、mmax為原始數(shù)據(jù)真實值中的最小值和最大值.

3 相對動彈性模量BP網(wǎng)絡(luò)模型

3.1 模型建立

運用MATLAB神經(jīng)網(wǎng)絡(luò)工具箱中的BP網(wǎng)絡(luò)模型，設(shè)置網(wǎng)絡(luò)模型參數(shù)如下：

1) 輸入層.輸入節(jié)點數(shù)為3，即橡膠摻量x1、水灰比x2和凍融循環(huán)次數(shù)x3.

2) 輸出層.輸出節(jié)點數(shù)為1，即相對動彈性模量y.

3) 訓(xùn)練允許誤差.均方差MSE=0.000 1.

4) 隱含層.隱含層為1層，隱含層節(jié)點數(shù)為18，是通過反復(fù)試算后確定的.

5) 傳遞函數(shù).確定輸入層與隱含層之間、隱含層與輸出層之間的傳遞函數(shù)均為Tansig.

6) 訓(xùn)練函數(shù).采用Trainlm.BP網(wǎng)絡(luò)的訓(xùn)練函數(shù)包括：梯度下降BP算法訓(xùn)練函數(shù)Traingd、動量反傳遞的梯度下降BP算法訓(xùn)練函數(shù)Traingdm、動態(tài)自適應(yīng)學(xué)習(xí)率的梯度下降BP算法訓(xùn)練函數(shù)Traingda、動量反傳遞和動態(tài)自適應(yīng)學(xué)習(xí)率的梯度下降BP算法訓(xùn)練函數(shù)Traingdx與Levenberg-marquardt的BP算法訓(xùn)練函數(shù)Trainlm.本文采用Trainlm是因為該函數(shù)通過自適應(yīng)調(diào)整可以獲取最快的調(diào)整量，收斂速度較快.

7) 學(xué)習(xí)函數(shù).采用帶動量項的Learngdm.

橡膠混凝土相對動彈性模量的BP網(wǎng)絡(luò)結(jié)構(gòu)如圖2所示.

圖2 橡膠混凝土相對動彈性模量的BP網(wǎng)絡(luò)結(jié)構(gòu)Fig.2 BP network structure of relative dynamic elasticity modulus of rubber concrete

3.2 模型訓(xùn)練

運用圖2中的BP網(wǎng)絡(luò)模型，訓(xùn)練結(jié)束后得到的訓(xùn)練誤差為7.12e-05，達(dá)到了預(yù)定的MSE=0.000 1的精度要求.將該模型訓(xùn)練結(jié)果與訓(xùn)練樣本實測結(jié)果進(jìn)行對比，結(jié)果如圖3所示.由圖3可知，該網(wǎng)絡(luò)模型的訓(xùn)練結(jié)果與實測結(jié)果比較吻合.

圖3 相對動彈性模量訓(xùn)練結(jié)果和實測結(jié)果對比Fig.3 Comparison between trained and measured results for relative dynamic elasticity modulus

該橡膠混凝土相對動彈性模量BP網(wǎng)絡(luò)的權(quán)值和閾值分別為

b2=[-1.034 6]

3.3 模型測試

運用該模型對27組橡膠混凝土相對動彈性模量預(yù)測樣本進(jìn)行測試，其預(yù)測結(jié)果如圖4所示.預(yù)測結(jié)果和實測值的變化趨勢基本一致，但兩者之間存在一定的誤差，為了更清晰地了解該網(wǎng)絡(luò)模型的性能，將預(yù)測結(jié)果反歸一化處理，其表達(dá)式為

圖4 相對動彈性模量預(yù)測結(jié)果和實測結(jié)果對比Fig.4 Comparison between predicted and measured results for relative dynamic elasticity modulus

(2)

由式(2)可得各組橡膠混凝土試件的相對動彈性模量預(yù)測值.將預(yù)測值與實測值進(jìn)行對比，其絕對誤差和相對誤差如表2所示.

表2 相對動彈性模量的BP網(wǎng)絡(luò)預(yù)測值和誤差Tab.2 Predicted values and errors of relative dynamic elasticity modulus %

由表2中橡膠混凝土相對動彈性模量的預(yù)測結(jié)果可以看出，大部分測試樣本的誤差都較小，有88.9%的測試樣本相對誤差均低于3%，只有11.1%的測試樣本相對誤差超過了3%，但均未超過5%，能夠滿足實際工程的精度要求.因此，運用該BP網(wǎng)絡(luò)模型可以較好地對凍融循環(huán)作用下橡膠粉混凝土的相對動彈性模量進(jìn)行預(yù)測.

4 結(jié) 論

本文通過分析得出以下結(jié)論：

1) 橡膠混凝土在凍融循環(huán)作用下的抗凍性受多種因素影響，本文所建立的橡膠混凝土抗凍性BP網(wǎng)絡(luò)模型，以橡膠摻量、水灰比和凍融循環(huán)次數(shù)為輸入變量，以相對動彈性模量為輸出變量，避免了建立精確的數(shù)學(xué)理論模型所面臨的困境.

2) 該BP網(wǎng)絡(luò)模型輸入層與隱含層之間、隱含層與輸出層之間的傳遞函數(shù)均為Tansig，訓(xùn)練函數(shù)為Trainlm，學(xué)習(xí)函數(shù)為Learngdm.模型收斂速度非?？欤?xùn)練結(jié)束后得到的訓(xùn)練誤差為7.12e-05.

3) 經(jīng)測試樣本驗證，該BP網(wǎng)絡(luò)模型具有較強(qiáng)的預(yù)測能力，其預(yù)測結(jié)果與試驗結(jié)果吻合較好，精度較高，88.9%的測試樣本相對誤差均低于3%，所有測試樣本的相對誤差均未超過5%.隨著訓(xùn)練樣本數(shù)量的增加，通過超強(qiáng)的自我學(xué)習(xí)能力，預(yù)測能力和精度將進(jìn)一步提高.能夠滿足實際工程的精度要求，可用于橡膠混凝土工程的抗凍性評估和預(yù)測.