龍武劍，周 波，梁沛堅(jiān)，孫柔嘉

深圳大學(xué)土木工程學(xué)院，廣東省濱海土木工程耐久性重點(diǎn)實(shí)驗(yàn)室，廣東深圳 518060

【土木建筑工程 / Architecture and Civil Engineering】

顆粒堆積模型在混凝土中的應(yīng)用

龍武劍，周 波，梁沛堅(jiān)，孫柔嘉

深圳大學(xué)土木工程學(xué)院，廣東省濱海土木工程耐久性重點(diǎn)實(shí)驗(yàn)室，廣東深圳 518060

評述顆粒堆積模型的基本理論、發(fā)展及其最新應(yīng)用. 對比分析多種顆粒堆積模型間的差異，指出顆粒堆積模型的主要優(yōu)點(diǎn)、局限性及其在復(fù)合顆粒體系材料領(lǐng)域中的應(yīng)用范圍. 分析表明，利用顆粒堆積模型優(yōu)化計(jì)算混凝土的堆積密實(shí)度，可以降低復(fù)合顆粒體系間的空隙體積，使顆粒系統(tǒng)達(dá)到最緊密堆積狀態(tài)，減少混凝土中膠凝材料的用量，實(shí)現(xiàn)低膠凝材料用量綠色混凝土優(yōu)化設(shè)計(jì).

混凝土；顆粒堆積模型；堆積密實(shí)度；顆粒體系；松動效應(yīng)；附壁效應(yīng)；空隙率；配合比

近年來，綠色生態(tài)混凝土的研究及制備在新型建筑材料研發(fā)領(lǐng)域備受關(guān)注. 混凝土作為應(yīng)用最廣泛且用量最大的建筑材料，雖然給人類提供了很大的便利性，但其生產(chǎn)制備過程也給環(huán)境帶來極大的負(fù)面影響. Gao等[1]研究發(fā)現(xiàn)，全球僅水泥產(chǎn)業(yè)每年就產(chǎn)生約1.8 ×1011kg的二氧化碳排放量，占世界二氧化碳總排量的7%左右，導(dǎo)致了生態(tài)環(huán)境被嚴(yán)重破壞.

實(shí)現(xiàn)綠色生態(tài)混凝土的制備及生產(chǎn)途徑有許多，其中之一是通過降低水泥用量，達(dá)到減少能源消耗，降低二氧化碳排放量的目的[2-5]. 從宏觀而言，混凝土混合料可被認(rèn)為是由不同尺寸顆粒堆積而成的顆粒系統(tǒng)，并被膠凝材料水化產(chǎn)物包裹膠結(jié)而成的一個整體. 因此，基于顆粒堆積模型可以優(yōu)化混凝土混合料級配，使混凝土混合物盡可能達(dá)到最密實(shí)狀態(tài)，從而將顆粒間的空隙率達(dá)到最低；同時，在保持原水膠比不變的前提下，通過減少填充空隙之間的用水量，最終降低其水泥用量，達(dá)到環(huán)保的目的.

本研究評述了顆粒堆積模型的基本理論和應(yīng)用研究，旨在為綠色混凝土材料設(shè)計(jì)奠定基礎(chǔ). 目前顆粒堆積模型的應(yīng)用技術(shù)仍處于實(shí)驗(yàn)階段，尚未形成完善的適用于混凝土材料設(shè)計(jì)的應(yīng)用理論. 系統(tǒng)研究顆粒堆積模型在綠色生態(tài)混凝土中的應(yīng)用理論及設(shè)計(jì)方法，有利于低水泥用量綠色混凝土在建材行業(yè)中的健康、可持續(xù)發(fā)展.

1 顆粒堆積理論

顆粒有宏觀和微觀之分. 從宏觀角度來看，它可以指砂、石、土，甚至是大的混凝土試件等；從微觀角度而言，它指的是分子和原子等. 同時，顆粒根據(jù)本身的不同特征可有不同的分類形式，例如，根據(jù)表觀形態(tài)，可以分為規(guī)則形態(tài)(球體和立方體等)和不規(guī)則形態(tài)；根據(jù)堆積方式，可分為有序堆積和隨機(jī)堆積(在重力影響下的一種隨機(jī)堆積方式). 對于宏觀顆粒的堆積，一般指在自身重力作用下的一種堆積方式. 微觀顆粒的堆積影響因素相對較為復(fù)雜，因還需考慮靜電作用和范德華力.

由于顆粒本身形貌及組成的復(fù)雜性，近一個世紀(jì)以來顆粒堆積問題一直是數(shù)學(xué)家、化學(xué)家、材料科學(xué)家及其他相關(guān)行業(yè)領(lǐng)域?qū)＜页掷m(xù)關(guān)注的焦點(diǎn)之一. 顆粒堆積理論主要是在20世紀(jì)30年代提出，可分為連續(xù)分布的顆粒體和不連續(xù)分布的顆粒體. 連續(xù)分布的顆粒體由某一粒級范圍內(nèi)所有尺寸顆粒組成，不連續(xù)分布顆粒體則由該粒級范圍的有限尺寸顆粒所組成.

Westman等[6]提出了的顆粒堆積理論及計(jì)算多元顆粒最大堆積率的方法，表觀體積為

(1)

其中， Va為單位實(shí)際體積顆粒的表觀體積； ε為單位實(shí)際體積空隙所占的空隙率； ω為單位實(shí)際體積顆粒所占的容積率.

Westman等[6]認(rèn)為當(dāng)粗細(xì)顆粒尺寸的比值達(dá)到某一個非常大的值時，可得出以下結(jié)論：① 當(dāng)單一容積內(nèi)粗顆粒所占組成部分接近100%時，粗細(xì)顆?；旌狭系谋碛^體積主要由粗顆粒決定，而細(xì)顆粒填入粗顆粒之間的空隙，可以忽略其所占容積；② 當(dāng)單一容積內(nèi)細(xì)顆粒所占組成部分接近100%時，細(xì)顆粒形成氣孔并環(huán)繞粗顆粒堆積，這時混合料的表觀體積除包括細(xì)顆粒的表觀體積之外，還包含粗顆粒的實(shí)際體積.

Andreason[7]認(rèn)為，用同一粒徑顆粒混合物的堆積理論不能準(zhǔn)確計(jì)算實(shí)際顆?；旌衔锏亩逊e率. 因此，他將顆粒的實(shí)際分布描述為具有同樣形式的分布，即“統(tǒng)計(jì)類似”，并提出了Andreason方程：

(2)

其中， Y(D)為小于粒徑D的累計(jì)篩余百分?jǐn)?shù)； D為顆粒粒徑； Dmax為顆粒體中最大顆粒粒徑； n為模型分布指數(shù).

式(2)描述了包含無限小的尺寸顆粒，這種情況雖然在實(shí)際的顆粒體系中不可能存在，但是Andreason[7]認(rèn)為，假如對于顆粒系統(tǒng)中最小顆粒的尺寸是有限小的或者是無限小的，則通過該方程計(jì)算出來的結(jié)果與實(shí)際值并沒有很大的差別. Andreason還通過實(shí)驗(yàn)結(jié)果得出了另一個結(jié)論，即對于顆粒各種分布的氣孔率隨著該方程中的顆粒分布指數(shù)n而減少，并認(rèn)為若使顆粒實(shí)際分布的氣孔率達(dá)到最小值， n的適宜取值范圍為0.33～0.50.

上述顆粒堆積的經(jīng)典理論，雖然存在一些不足，但對顆粒堆積模型的系統(tǒng)研究具有深遠(yuǎn)的理論影響和實(shí)用意義，如基于經(jīng)典堆積理論Andreason方程基礎(chǔ)而修正所得的Dinger-Funk方程[8]等.

2 顆粒堆積模型的研究現(xiàn)狀

2.1 顆粒堆積模型的發(fā)展

顆粒堆積模型是以數(shù)學(xué)方程為基礎(chǔ)，用來描述不同大小顆粒之間的幾何作用方式. 用模型來計(jì)算混凝土混合物的理論堆積密實(shí)度，是基于顆粒級配和顆粒群組的堆積密實(shí)度，且?guī)缀跛械念w粒堆積模型的數(shù)學(xué)方程本構(gòu)是一致的.

早在20世紀(jì)30年代，F(xiàn)urnas[9]首次用數(shù)學(xué)方程來描述顆粒堆積模型，提出了Furnas模型，但最初的Furnas模型是基于球形顆粒堆積理論，且只有在不考慮兩個粒徑顆粒間作用的前提下才能應(yīng)用該模型，即當(dāng)大顆粒占主導(dǎo)時，小顆粒能夠填充在大顆粒之間的空隙中，或以小顆粒為主時，大顆粒鑲嵌在小顆粒基體中. 然而，這與顆粒的實(shí)際堆積情況有較大差別. 在隨后的幾十年里，大量研究者基于Furnas理論和試驗(yàn)研究，對顆粒堆積模型進(jìn)行不斷的改善，提出了不同的顆粒堆積模型. 在這些模型中涉及到兩種效應(yīng)，即松動效應(yīng)和附壁效應(yīng)[9]. 松動效應(yīng)是指在以粗顆粒為主的情況下，細(xì)顆粒填充在粗顆粒堆積的空隙之中，并與空隙大小相比足夠大，如圖1；附壁效應(yīng)是指在以細(xì)顆粒為主的情況下，某些單個粗顆粒處于細(xì)顆粒群之中，而環(huán)繞粗顆粒接觸界面堆積的細(xì)顆粒會出現(xiàn)一些空隙量，如圖2. 這兩種效應(yīng)都隨著顆粒粒徑比值的增大而逐漸提高，但這兩種效應(yīng)的最大值均為1，即兩種顆粒材料粒徑相等且完全相互作用的情況(完全相互作用指的是某一尺寸顆粒的排列會受到鄰近的另一不同尺寸顆粒的影響).

圖1 松動效應(yīng)[9]Fig.1 Loosening effect[9]

圖2 附壁效應(yīng)[9]Fig.2 Wall effect[9]

1977年，Toufar等[10]提出了Toufar模型，并用該模型來預(yù)測直徑比在0.22～1.00的二元混合料的堆積密實(shí)度，它與A?m 等[11]在1968年提出的模型相同之處在于前者只考慮了二元混合料顆粒之間的附壁效應(yīng)，而忽視了顆粒間的松動效應(yīng). 1997年，Goltermann等[12]對Toufar模型進(jìn)行了修正，除了考慮顆粒之間的附壁效應(yīng)外，還包含了松動效應(yīng). 修正的Toufar模型對二元混合料堆積密實(shí)度的預(yù)測也較準(zhǔn)確，但是當(dāng)混合料中顆粒級別被劃分較多時，其預(yù)測堆積密實(shí)度的準(zhǔn)確性將會降低，甚至?xí)霈F(xiàn)較大的偏差[13].

1986年，Stovall等[14]通過將粉末顆粒堆積體系分為擠塞和非擠塞體系，提出了線性堆積密實(shí)度模型(linear packing density model，LPDM). 應(yīng)用該模型的條件是：假設(shè)在粉末顆粒堆積過程中，不考慮顆粒之間的靜電作用、范德華分子作用及顆粒團(tuán)聚作用等因素的影響，僅考慮顆粒間幾何作用對顆粒堆積密實(shí)度的影響. 此模型雖是現(xiàn)有粉末顆粒堆積模型中預(yù)測堆積密實(shí)度準(zhǔn)確度較高的模型，但由于其沒有考慮到替換機(jī)制[15]，即當(dāng)二元顆粒的粒度比大于某一臨界值時，其他顆粒的加入量對其堆積密實(shí)度無影響，故該模型的精度仍有待提高.

1998年提出的Dewar模型[16]通過利用正方形形狀顆粒的堆積方式來推導(dǎo)顆粒間的相互作用，這在一定程度上間接地考慮了顆粒的形狀和紋理. 由于該模型采用的是逐步回歸算法，且在模型中引入了兩個系數(shù)，黏性調(diào)整系數(shù)和經(jīng)驗(yàn)調(diào)整系數(shù)，從而局限了模型向顆粒形狀或者顆粒表面力作用效應(yīng)兩個方面去發(fā)展.

1999年，Larrard[17]提出了線性堆積密實(shí)度模型和可壓縮堆積模型(compressible packing model，CPM)，這兩種模型都需測量混凝土混合料的顆粒級配曲線. 線性堆積密實(shí)度模型是對Furnas模型的拓展和改進(jìn)，既可以預(yù)測多元混合料的堆積密實(shí)度，又考慮了混合料顆粒之間的幾何作用(附壁效應(yīng)和松動效應(yīng)). 而可壓縮堆積模型是基于線性堆積密實(shí)度模型和固體懸浮模型(solid suspension model, SSM)發(fā)展而成，它不僅考慮了混合料中不同尺寸顆粒間的相互作用，且涵蓋了混合料不同堆積形式對顆粒堆積密實(shí)度的影響，甚至通過大量實(shí)驗(yàn)對顆粒壓實(shí)指數(shù)K進(jìn)行了校正，并給出了K在不同堆積形式下的取值，如表1. 因此，CPM模型對預(yù)測任意分布及任意量的骨料顆粒組合而成的混合料堆積密實(shí)度具有較高的準(zhǔn)確度，并得到大多數(shù)研究者的認(rèn)可且被廣泛應(yīng)用.

表1 不同堆積形式下K的取值[17]

Kwan等[18-20]通過對顆粒堆積模型進(jìn)行進(jìn)一步研究，在原顆粒結(jié)構(gòu)堆積效應(yīng)(附壁效應(yīng)和松動效應(yīng))的基礎(chǔ)上，再引入了一個新的幾何作用效應(yīng)——楔入效應(yīng)，該效應(yīng)存在如圖3和圖4兩種情況. 圖3表明了當(dāng)混合料中組成部分主要是粗顆粒時，細(xì)顆粒不能填充在粗顆粒自身堆積之間的空隙，而是處于粗顆粒與粗顆粒之間，即將兩個相鄰的粗顆粒分離開；圖4表明了當(dāng)混合料中細(xì)顆粒占主要成分時，細(xì)顆粒環(huán)繞少量粗顆粒堆積，同樣使得兩個可以相接觸的粗顆粒分離. Kwan根據(jù)這3種效應(yīng)建立了三參數(shù)模型，并利用該模型對同一形狀(橢圓或菱角狀)顆粒混合物的堆積密實(shí)度進(jìn)行預(yù)測．結(jié)果表明，該三參數(shù)模型不但突破了線性堆積理論的局限性，同時具有較高的預(yù)測準(zhǔn)確度.

圖3 粗骨料占主導(dǎo)地位的楔入效應(yīng)[18]Fig.3 The wedging effect when the coarse particles are dominant[18]

圖4 細(xì)骨料占主導(dǎo)地位的楔入效應(yīng)[18]Fig.4 The wedging effect when the fine particles are dominant[18]

Knop等[21]基于Stovall和Larrard提出的線性堆積密實(shí)模型，對其部分進(jìn)行了簡單修正，并利用修正模型對水泥與不同尺寸的石粉顆粒任意混合組合的堆積密實(shí)度進(jìn)行預(yù)測．結(jié)果表明，其值與實(shí)際值的誤差在允許范圍之內(nèi).

2.2 離散單元法在顆粒體系中的應(yīng)用

近年來，隨著計(jì)算機(jī)技術(shù)的不斷發(fā)展，許多研究者已利用相應(yīng)的數(shù)學(xué)工具對顆粒體系在不同形式下的堆積過程進(jìn)行了模擬；尤其是離散單元法(discrete element model, DEM)在顆粒體系中的應(yīng)用已得到廣大研究者的認(rèn)可.

利用離散單元法進(jìn)行仿真分析，首先需熟悉其兩種顆粒接觸模型：硬球模型和軟球模型. 軟球模型[22]存在幾種主要的假設(shè)：① 模擬顆粒是一種形狀為圓形的剛體，且顆粒間的接觸均為點(diǎn)接觸；② 顆粒間接觸時可以允許有少量的重疊部分(見圖5中δij， vi和vj表示顆粒運(yùn)動方向， ri和rj為顆粒半徑， rij為兩個顆粒間的距離)；③ 當(dāng)某一顆粒受到擾動后，只對相鄰顆粒有一定的影響；④ 顆粒所受合力取決于顆粒間的相互作用. 而對于整個顆粒系統(tǒng)，既要滿足本構(gòu)方程，又需滿足平衡方程，本構(gòu)方程主要指顆粒接觸模型中顆粒間相互作用力與位移的關(guān)系，平衡方程則是指經(jīng)典力學(xué)中的牛頓第2運(yùn)動定律.

圖5 軟球模型中顆粒接觸碰撞[22]Fig.5 Particle contact collision in soft sphere model[22]

Wen等[22]通過利用DEM模擬了顆粒在不同形式下(滾筒和噴動床)的堆積過程，并分析了顆粒在這兩種形勢下的網(wǎng)格、色彩和位置等因素，評估了不同方法下顆粒攪拌指數(shù)的一些性能，研究表明，相鄰距離法對顆粒的攪拌過程具有更好的適應(yīng)性. Fan等[23]仿真了2維流化床上3種顆粒間無黏結(jié)力的彈射機(jī)制，并探討了外部施加的作用力對顆粒的影響，認(rèn)為顆粒之間的黏結(jié)力對顆粒的彈射和強(qiáng)度有影響. Cleary[24]在采用離散元模擬研磨機(jī)時，提出了一種對粉體或細(xì)顆粒多尺度模擬的方法，即微觀尺度模型采用適宜的剪切單元，嵌入到傳統(tǒng)的數(shù)字高程模型(digital elevation model, DEM)模型中；而宏觀尺度模型則為相應(yīng)的剪切單元給予剪切率和應(yīng)力，然后通過微觀尺度模型來預(yù)測粉體或細(xì)顆粒的研磨環(huán)境，并利用多尺度模型來仿真研磨的整個過程.

雖然這些模擬仿真顆粒堆積過程不是基于顆粒堆積模型，但有助于推動顆粒堆積模型的發(fā)展和應(yīng)用. 如時金廣[25]通過研究微小顆粒團(tuán)聚的堆積狀態(tài)，并通過利用EDEM軟件用數(shù)值模擬顆粒在自然狀態(tài)(圖6)和振動兩種方式下的堆積過程，并著重分析了振動時間和振動頻率對顆粒堆積密實(shí)度的影響，對可壓縮堆積模型進(jìn)行了優(yōu)化，能夠使優(yōu)化后的模型更加精確地預(yù)測混凝土混合料的堆積密實(shí)度.

圖6 模擬自然狀態(tài)下顆粒堆積[25]Fig.6 (Color online) Simulation of particle accumulation under the natural state[25]

3 顆粒堆積模型對比研究

近百年來，對顆粒堆積模型的研究一直處于發(fā)展階段，研究者提出的模型均有自身的優(yōu)點(diǎn)和局限性. 因此，不同的顆粒堆積模型對于混合料的堆積密實(shí)度的預(yù)測準(zhǔn)確度不完全一致，即在某些局限特定的條件下雖然相差不大，但在普遍適應(yīng)條件下還是存在很大差異.

Johansen等[26]通過將單一粒級的二元混合料進(jìn)行混合實(shí)驗(yàn)，得到了堆積密實(shí)度實(shí)驗(yàn)值，對A?m-Goff模型[12]和Toufar模型[11]進(jìn)行了比較. 結(jié)果表明，當(dāng)二元混合料的尺寸比小于0.22時，A?m-Goff模型較準(zhǔn)確；當(dāng)二元混合料的尺寸比大于0.22時，Toufar模型則更準(zhǔn)確.

Jones等[27-28]由實(shí)驗(yàn)所得的堆積密實(shí)度值，比較了文獻(xiàn)[12,14,16-17]分別提出的4種顆粒堆積模型，并指出這些模型與其他模型的準(zhǔn)確度不同；這些模型的預(yù)測準(zhǔn)確度隨著顆粒尺寸比的不同而變化，即當(dāng)兩個顆粒級別尺寸比大于0.4時，這4種模型對混合料的堆積密實(shí)度的預(yù)測準(zhǔn)確度較好且相差不大.

Fennis[29]對當(dāng)前使用較普遍的顆粒堆積模型進(jìn)行了闡述和比較，如Furnas模型、Toufar模型、Dewar模型、LPDM模型和CPM模型，認(rèn)為這些模型對于粗骨料在干堆積和濕堆積情況下的堆積密實(shí)度是一致的，但對于細(xì)顆粒體系，這些模型大都沒考慮用濕堆積的方法來計(jì)算混合料堆積密實(shí)度. 由于CPM模型考慮了顆粒壓實(shí)作用系數(shù)，被廣大研究者認(rèn)為是這些模型中預(yù)測堆積密實(shí)度最準(zhǔn)確的一種，并直接被用于計(jì)算粗細(xì)骨料混合物的干濕堆積密實(shí)度.

Chan等[13,18]通過試驗(yàn)研究比較了不同顆粒堆積模型的準(zhǔn)確性，根據(jù)模型中涉及到的顆粒間相互作用(松動效應(yīng)系數(shù)a、 附壁效應(yīng)系數(shù)b、 楔入效應(yīng)系數(shù)c和壓實(shí)指數(shù)k)， 把現(xiàn)有的顆粒堆積模型分為2參數(shù)模型(a, b)、 CPM模型(a,b, k)、 3參數(shù)模型(a,b, c)， 然后通過試驗(yàn)來比較這些模型的準(zhǔn)確性. 結(jié)果表明，這些模型的準(zhǔn)確性隨著二元混合料顆粒的尺寸比值和體積分?jǐn)?shù)的變化而變化，當(dāng)顆粒尺寸比值大于0.65時，所有顆粒堆積模型的預(yù)測精度都較為準(zhǔn)確；但當(dāng)比值小于0.65時，2參數(shù)模型的計(jì)算準(zhǔn)確度變得較低，甚至與實(shí)驗(yàn)值相差較大. 相對而言，3參數(shù)模型的計(jì)算精度是上述模型中最準(zhǔn)確，也是應(yīng)用最為廣泛的模型.

各模型所預(yù)測堆積密實(shí)度的準(zhǔn)確度對于材料的粒徑、類別要求差異較大，且各模型的適用范圍、特點(diǎn)也略有不同. 基于對上述各模型的深入研究，并結(jié)合模型的不斷優(yōu)化和演化，對其主要優(yōu)點(diǎn)、局限性、適用范圍做出相關(guān)總結(jié)，如表2.

表2 顆粒堆積模型對比

(續(xù)表1)

4 基于顆粒堆積模型的應(yīng)用研究

由于顆粒堆積密實(shí)度的研究具有重要的理論和實(shí)際意義，因而一直受到國內(nèi)外研究者的關(guān)注. 研究人員利用顆粒堆積模型不但研究了混凝土混合料的堆積問題，且在此基礎(chǔ)上研究了摻合料對整個混合料的堆積密實(shí)度的影響，進(jìn)一步探討了多尺寸顆粒堆積密實(shí)度對混凝土性能的影響[30-35].

牛全林等[36]根據(jù)Horsfield模型和Andrensen方程[7，37]探討了礦渣超細(xì)粉對膠凝材料(水泥)的密實(shí)填充效應(yīng). 結(jié)果表明，根據(jù)Horsfield模型確定的超細(xì)填充粉，具有良好的填充效應(yīng)，能有效改善硬化水泥漿體的孔徑分布，而基于Andrensen 方程，將礦渣超細(xì)粉摻入水泥后顆粒的粒徑分布狀態(tài)與最緊密堆積仍存在較大的差距，這由Andrensen 方程可知

(3)

其中， D為顆粒粒徑； Dmax為顆粒體中最大顆粒粒徑； n為模型分布指數(shù)，取值為1/3.

陳延信等[38]基于Andrensen方程和CPM模型及實(shí)驗(yàn)研究探討了提高粉體堆積密實(shí)的理論. 結(jié)果表明，用CPM模型來預(yù)測粉體顆粒的堆積密實(shí)度實(shí)驗(yàn)值為0.698 5，模型計(jì)算值為0.695 1，相對誤差僅為0.810 0%. 用該模型預(yù)測粉體顆粒體系的堆積密實(shí)度具有較高的可靠性，且經(jīng)過模型調(diào)質(zhì)后可以使混合體系的堆積效率提高6.0%.

龍廣成等[39]基于 A?m-Goff 模型定量分析了超細(xì)粉煤灰、磨細(xì)礦渣及硅粉3種礦物摻合料對新拌水泥凈漿堆積密實(shí)度的作用變化，研究了摻入這3種礦物料后漿體的流動性與相對密實(shí)度之間的關(guān)聯(lián)性，探討了礦物摻合料的密實(shí)填充作用對水泥膠砂抗壓強(qiáng)度的影響. 結(jié)果表明，在所做試驗(yàn)摻量范圍內(nèi)，這3種礦物摻合料能較大地改善新拌水泥漿體的密實(shí)度，其中硅灰對顆粒體系的作用效果最佳. 同時，新拌漿體相對密實(shí)度提高后，硬化后漿體的抗壓強(qiáng)度也有所提高.

Binici等[40-41]研究了摻入不同細(xì)度的摻合料(高爐礦渣和天然火山灰)對混合料堆積密實(shí)度的影響. 結(jié)果表明，高爐礦渣和天然火山灰細(xì)度越小，混合料的堆積密實(shí)度越大，對混凝土強(qiáng)度的提高越明顯.

Kwan等[42-43]基于三參數(shù)顆粒堆積模型研究了粉煤灰摻量為0～60%時對水泥基復(fù)合材料密實(shí)度、流動度和強(qiáng)度的影響. 結(jié)果表明，摻量為 40%的超細(xì)粉煤灰對水泥基材料的堆積密實(shí)度的貢獻(xiàn)最大，且超細(xì)粉煤灰有助于提高漿體黏度.

Ghoddousi等[44]研究了顆粒堆積密實(shí)度對自密實(shí)混凝土穩(wěn)定性和流變性能的影響. 研究表明，對于粗細(xì)砂混合，總是存在一個最佳值，使空隙率最低，松動效應(yīng)和附壁效應(yīng)減少. 研究同時得到了顆粒堆積密實(shí)度的一個優(yōu)化值，為0.804，比具有低靜態(tài)離析的最大堆積密實(shí)度值0.807略小. 研究者認(rèn)為利用濕堆積來計(jì)算自密實(shí)混凝土是一種可行的方法.

Nassim 等[45]通過比較利用可壓縮堆積模型、國內(nèi)規(guī)范(JGJT283—2012自密實(shí)混凝土應(yīng)用技術(shù)規(guī)程)和歐洲標(biāo)準(zhǔn)(BS EN 206-1)設(shè)計(jì)了自密實(shí)混凝土配合比. 結(jié)果表明，利用這3類標(biāo)準(zhǔn)進(jìn)行自密實(shí)混凝土配合比設(shè)計(jì)，可以得到28 d平均抗壓強(qiáng)度為(72.0±2.5)MPa的中等強(qiáng)度混凝土；同時，利用可壓縮堆積模型進(jìn)行設(shè)計(jì)的混凝土配合比中所使用膠凝材料用量相對較低，約為400 kg/m3. 因此，利用可壓縮堆積模型來制備自密實(shí)混凝土不但可以滿足工作性能和力學(xué)性能的需求，而且經(jīng)濟(jì)環(huán)保.

龔建清[46]論述了基于CPM模型計(jì)算粉體混合料堆積密實(shí)度時所存在的問題和相對應(yīng)的解決方法，并利用CPM模型分析了微粉在不同顆粒粒徑和摻入量的情況下，對顆粒混合料體系堆積密實(shí)度的影響. 實(shí)驗(yàn)表明，當(dāng)?shù)V物微粉在一定摻量范圍內(nèi)，顆粒體系的堆積密實(shí)度隨著摻量的增加而提高，其相對應(yīng)的抗壓強(qiáng)度也不斷提高；但當(dāng)摻量超過某一值時，顆粒體系的堆積密實(shí)度不再隨摻量的增加而提高. 同時，研究者從礦物細(xì)摻料的密實(shí)填充效應(yīng)角度，探討了礦物微粉的級配效應(yīng)對超高性能混凝土性能的影響，認(rèn)為當(dāng)超細(xì)粉煤灰和硅灰雙摻后顆粒體系比單摻超細(xì)粉煤灰后顆粒體系的顆粒堆積結(jié)構(gòu)更加密實(shí)，其顆粒級配更好，在水膠比較低(0.22)的情況下，密實(shí)填充效應(yīng)對顆粒體系堆積結(jié)構(gòu)更加密實(shí)的漿體作用更加明顯，更能充分發(fā)揮這些礦物微粉顆粒對顆粒體系結(jié)構(gòu)的密實(shí)填充效應(yīng)，對漿體抗壓強(qiáng)度的提高作用也就越明顯.

曹榮奎[47]利用CPM模型分析了礦物微粉(粉煤灰、礦渣和硅灰)以不同比例混合后，對其火山灰效應(yīng)和顆粒堆積密實(shí)度的影響；并重點(diǎn)對火山灰效應(yīng)和復(fù)合火山灰效應(yīng)進(jìn)行了模擬研究，通過數(shù)學(xué)工具擬合相關(guān)的參數(shù)(如堆積密實(shí)度作用值與火山灰效應(yīng)強(qiáng)度貢獻(xiàn)率間的關(guān)聯(lián)性)，提出了超高性能混凝土混合料配合比設(shè)計(jì)方法.

聶晶[48]根據(jù)CPM模型，對水泥基復(fù)合材料堆積密實(shí)度進(jìn)行了研究，如水泥-礦渣體系、水泥-粉煤灰體系等，并利用Andreasen方程計(jì)算了上述復(fù)合材料體系最緊密堆積狀態(tài)下的顆粒粒徑分布. 研究結(jié)果表明，Andreasen計(jì)算得到的顆粒粒徑分布下體系的堆積密實(shí)度與CPM模型計(jì)算結(jié)果一致；同時，研究者認(rèn)為礦物微粉對水泥顆粒堆積體系的密實(shí)效應(yīng)很大程度由礦物微粉顆粒的平均粒徑與水泥顆粒的平均粒徑之比決定.

李丙明[49]利用CPM模型對鋼渣、礦渣、超細(xì)粉煤灰(平均粒徑為5.51 μm)3種微粉摻入到水泥體系后混合料實(shí)際堆積密實(shí)度的變化進(jìn)行了研究. 結(jié)果表明，只有當(dāng)鋼渣、礦渣與超細(xì)粉煤灰復(fù)合摻入水泥后，粉體的堆積密實(shí)度才能提高；且粉煤灰的摻加比例越大，混合料的密實(shí)度越高. 研究者同時研究了上述3種復(fù)合微粉以不同摻量替代水泥后顆粒級配與堆積密實(shí)度之間的關(guān)系，及對混凝土性能的影響.

丑凱[50]在將CPM模型應(yīng)用于粉體顆粒體系(包括單摻硅灰體系和雙摻硅灰和粉煤灰體系)的堆積密實(shí)度計(jì)算時，發(fā)現(xiàn)堆積密實(shí)度與火山灰效應(yīng)在不同水膠比條件下存在一定關(guān)聯(lián)，提出了這兩種不同體系在28 d標(biāo)準(zhǔn)養(yǎng)護(hù)條件下，混凝土的堆積密實(shí)度貢獻(xiàn)率與火山灰效應(yīng)強(qiáng)度貢獻(xiàn)率之間的擬合通式為

1)單摻硅灰體系

(4)

其中， y為火山灰效應(yīng)貢獻(xiàn)率； x為堆積密實(shí)度貢獻(xiàn)率； A1、 A2、 x0和dx為擬合曲線參數(shù).

2)雙摻體系

(5)

其中， y為火山灰效應(yīng)貢獻(xiàn)率； x為堆積密實(shí)度貢獻(xiàn)率； y0、 xc、 A和w為擬合曲線參數(shù).

最后，丑凱[50]利用數(shù)學(xué)工具將堆積密實(shí)度與火山灰效應(yīng)進(jìn)行合理量化，并通過實(shí)例來驗(yàn)證兩者之間的關(guān)聯(lián)性，提出了一種有效設(shè)計(jì)超高性能混凝土配合比的方法.

何振偉[51]通過利用水化動力學(xué)理論模型，并基于CPM模型-粉煤灰混凝土強(qiáng)度預(yù)測模型和基于CPM 的彈性模量三重球模型的基礎(chǔ)上，提出了持續(xù)應(yīng)力下密封早齡期混凝土抗壓強(qiáng)度和彈性模量根據(jù)時間不斷變化而發(fā)展的理論模型，并根據(jù)實(shí)驗(yàn)對模型參數(shù)進(jìn)行擬合，對新模型相關(guān)參數(shù)進(jìn)行分析. 結(jié)果表明，所提新模型能很好地預(yù)測早期鋼管混凝土力學(xué)性能的發(fā)展變化. 這里提及的三重球模型是指將混凝土混合料看成由核心為漿體球體、中間為骨料、外層是漿體組成的一個完整球體，如圖7.

圖7 三重球模型基本單元[51]Fig.7 The basic unit of concrete based on the triple-sphere model[51]

Fennis[29]在CPM模型的基礎(chǔ)上，提出了與混凝土性能相關(guān)的兩個指數(shù)，即膠凝材料顆粒間隙指數(shù)和用水量富余指數(shù)，并將實(shí)際狀態(tài)下膠凝材料顆粒間隙指數(shù)(cement space factor， CSF)定義為

(6)

其中, Kb為膠凝材料顆粒的疊加壓實(shí)值； Vb為穩(wěn)定顆粒結(jié)構(gòu)狀態(tài)下固體顆粒的體積； Vw為混合料拌和時總的需水量.

(7)

其中， Vcem為單位體積混凝土中水泥含量； Va為單位體積混凝土中粗細(xì)骨料含量； Vm為單位體積混凝土中礦物摻合料含量； Bm為修正系數(shù)； at為混合料顆粒實(shí)際堆積密實(shí)度.

陳瑾祥[52]的實(shí)驗(yàn)研究表明，水灰比為0.45～0.57的純水泥混凝土，基于用水量富余指數(shù)所提出的工作性能預(yù)測模型與實(shí)測工作性能兩者之間的相關(guān)系數(shù)可達(dá)到0.96；而基于CSF指數(shù)所提出的標(biāo)準(zhǔn)混凝土強(qiáng)度預(yù)測模型與實(shí)測強(qiáng)度值之間的相關(guān)系數(shù)也較高，即3、7和28 d的相關(guān)系數(shù)分別為0.87、0.93和0.98.

5 結(jié) 語

綜上研究認(rèn)為，顆粒堆積模型是基于數(shù)學(xué)方程描述不同尺寸顆粒間幾何作用的計(jì)算方法，利用顆粒堆積模型可以計(jì)算混合料的堆積密實(shí)度并對其進(jìn)行優(yōu)化，使混合料達(dá)到物理最緊密的堆積狀態(tài)，從而減小顆粒體系的空隙率，減少水泥用量，最終實(shí)現(xiàn)綠色混凝土的配合比優(yōu)化設(shè)計(jì). 目前對各類型顆粒堆積模型的基礎(chǔ)及應(yīng)用研究尚可從以下入手完善.

1)對于固體顆粒體系，雖然大多數(shù)顆粒模型對其堆積密實(shí)度預(yù)測都較準(zhǔn)確，但在應(yīng)用這些模型預(yù)測顆粒堆積密實(shí)度時，均假設(shè)所有堆積顆粒的形狀是同一的，沒有考慮顆粒的形貌特征影響，故與堆積顆粒的真實(shí)情況不一致，使得預(yù)測理論值與實(shí)測值存在一定誤差. 建議引入形貌參數(shù)或形貌函數(shù)，以提高顆粒堆積模型的預(yù)測精度.

2)對于細(xì)顆粒體系，因其顆粒之間除了重力作用之外，還需考慮范德華力和靜電力等作用，從而使得顆粒之間的相互作用變得復(fù)雜. 因此有必要對顆粒堆積模型的作用機(jī)理進(jìn)行學(xué)入研究.

3)利用顆粒堆積模型，計(jì)算粗細(xì)骨料-膠凝材料混合堆積密實(shí)度的研究目前存在空白.

4)顆粒堆積模型在混凝土中的應(yīng)用，主要側(cè)重于計(jì)算混凝土混合料顆粒的物理堆積密實(shí)度，而膠凝材料水化后顆粒系統(tǒng)結(jié)構(gòu)會發(fā)生變化，這對顆粒堆積密實(shí)度會產(chǎn)生影響，而已有的顆粒堆積模型不再適用于準(zhǔn)確計(jì)算其堆積密實(shí)度. 如何考慮物理堆積和膠凝材料化學(xué)反應(yīng)所產(chǎn)生的對堆積密實(shí)度的協(xié)同作用效應(yīng)，還需更為系統(tǒng)深入地研究.

5)隨著高性能計(jì)算機(jī)的發(fā)展，計(jì)算仿真模擬在混凝土中的應(yīng)用已成為研究者的有效工具，但目前對于混凝土硬化全過程較全面的模擬仿真分析仍有待推動.

/ References：

[1] Gao Tianming, Shen Lei, Shen Ming, et al. Analysis on differences of carbon dioxide emission from cement production and their major determinants[J]. Journal of Cleaner Production, 2015, 103(23): 160-170.

[2] Suhendro B. Toward green concrete for better sustainable environment[J]. Procedia Engineering, 2014, 95(95): 305-320.

[3] Yoshioka K, Obata D, Nanjo H, et al. New ecological concrete that reduces CO2emissions below zero level new method for CO2capture and storage[J]. Energy Procedia, 2013, 37: 6018-6025.

[4] HiguchiI T, Morioka M, Yoshioka I, et al. Development of a new ecological concrete with CO2emissions below zero[J]. Construction and Building Materials, 2014, 67: 338-343.

[5] Yu R, Spesz P, Brouwers H J H. Development of an eco-friendly ultra-high performance concrete (UHPC) with efficient cement and mineral admixtures uses[J]. Cement and Concrete Composites, 2015, 55(1): 383-394.

[6] Westman A E R, Hugill H R. The packing of particles[J]. Journal of the American Ceramic Society ，1930, 13(10): 767-779.

[7] Andreason A H M. Uber die beziehung zwischen kornabstufung und zwischenraum in produckten aus losen kornern[J]. Kolloid Z, 1930, 50(3): 217-228.

[8] Dinger D R, Funk J E. Particle packing: review of packing theories[C]// Proceeding of Fine Particle Society 13th Annual Meeting. Chicago, USA: American Ceramic Society, 1982.

[9] Furnas C C. Grading aggregates I-mathematical relations for beds of broken solids of maximum density[J]. Industrial and Engineering Chemistry， 1931, 23(9): 1052-1028.

[10] Toufar W, Born E K M. Berechnung der packungsdichte von korngemischen[J]. Aufbereitungs-Technik, 1977, 18(11): 603-608.

[11] A?m R B, Goff P L. Effet de paroi dans les empilements désordonnés desphères et application à la porosité de mélanges binaries[J]. Powder Technology, 1968,1(5): 281-290.

[12] Goltermann P, Johansen V, Palbol L. Packing of aggregates: an alternative tool to determine the optimal aggregate mix[J]. American Concrete Institute Materials Journal, 1997, 94(5)：435-443.

[13] Chan K W, Kwan A K H. Evaluation of particle packing models by comparing with published test results[J]. Particuology, 2014, 16(5): 108-115.

[14] Stovall T, Bull M. Linear packing density model of grain mixtures[J]. Powder Technology, 1986, 48(1): 1-12.

[15] 王海兵,劉 詠,黃伯云, 等. 粉末顆粒線性堆積密實(shí)度模型的改進(jìn)[J].粉末冶金技術(shù), 2001, 19(4):208-211. Wang Haibing, Liu Yong, Huang Boyun, et al. Improvement of linear packing density model of powder particles[J]. Powder Metallurgy Technology, 2001, 19(4):208-211.(in Chinese)

[16] Dewar J D. Computer modelling of concrete mixtures[M]. London: Routledge, 1998.

[17] Larrard F. Concrete mixture proportioning: a scientific approach[M]. London: Routledge, 1999.

[18] Kwan A K H. A 3-parameter particle packing model incorporating the wedging effect [J]. Powder Technology, 2013, 237(3): 172-179.

[19] Wong V, Kwan A K H. A 3-parameter model for packing density prediction of ternary mixes of spherical particles[J]. Powder Technology, 2014, 238(1): 357-367.

[20] Kwan A K H, Wong V, Fung W W S. A 3-parameter packing density model for angular rock aggregate particles[J]. Powder Technology, 2015, 239(2): 154-162.

[21] Knop Y, Peled A. Packing density modeling of blended cement with limestone having different particle size[J]. Construction and Building Materials, 2016, 102: 44-50.

[22] Wen Yuanyun, Liu Malin, Liu Bing, et al. Comparative study on the characterization method of particle mixing index using DEM method[J]. Procedia Engineering, 2015, 102: 1630-1642.

[23] Fan Haojie, Mei Dengfei, Tian Fengguo, et al. DEM simulation of different particle ejection mechanisms in a fluidized bed with and without cohesive interparticle forces[J]. Powder Technology, 2016, 288: 228-240.

[24] Cleary P W. A multiscale method for including fine particle effects in DEM models of grinding mills[J]. Minerals Engineering, 2015, 84: 88-99.

[25] 時金廣.基于堆積模型的膠凝材料顆粒系統(tǒng)的離散單元法分析仿真研究[D]. 深圳：深圳大學(xué)，2016. Shi Jinguang. The research of discrete element method simulation for particle packing-based cementitious materials[D]. Shenzhen: Shenzhen University, 2016.(in Chinese)

[26] Johansen V, Andersen P J. Particle packing and concrete properties[J]. Materials Science of Concrete ， 1989, 2: 111-146.

[27] Jones M R, Zheng L, Newlands M D. Comparison of particle packing models for proportioning concrete constituents for minimum voids ratio[J]. Materials and Structures, 2002, 35(5): 301-309.

[28] Standish N, Borger D E. The porosity of particulate mixtures[J]. Powder Technology, 1979, 22(1): 121-125.

[29] Fennis S A A M. Design of ecological concrete by particle packing optimization[D]. Netherlands: Delft University of Technology, 2011.

[30] Kwan A K H, Wong H H C. Packing density of cementitious materials: packing and flow of OPC+PFA+CSF[J]. Materials and Structures, 2008, 41: 773-784.

[31] Gallias J L, Kara-Ali R, Bigas J P. The effect of fine mineral admixtures on water requirement of cement pastes[J]. Cement Concrete Composites, 2000, 30(10): 1543-1549.

[32] Peng Yanzhou, Hu Shuguang, Ding Qingjun. Dense packing properties of mineral admixtures in cementitious material[J]. Particuology, 2009, 7(5): 399-402.

[33] Knop Y, Peled A. Setting behavior of blended cement with limestone: influence of particle size and content[J].Materials and Structures, 2016, 49(1): 439-452.

[34] Felekogˇlu B. Effects of PSD and surface morphology of micro-aggregates on admixture requirement and mechanical performance of micro-concrete[J]. Cement Concrete Composites, 2007, 29(6): 481-489.

[35] Knop Y, Peled A, Cohen R. Influence of limestone particle size distributions and contents on blended cement properties[J]. Construction and Building Materials, 2014, 71: 26-34.

[36] 牛全林，馮乃謙，楊 靜.礦渣超細(xì)粉作用機(jī)理的探討[J].建筑材料學(xué)報, 2002, 5(1)：84-89. Niu Quanlin, Feng Naiqian, Yang Jing. Discussion on the action mechanism of superfine slag powder[J]. Journal of Building Materials, 2002, 5(1)：84-89.(in Chinese)

[37] Fayed M E, Otten L. Handbook of powder science and technology [M]. New York: Chapman and William Hall, 1997.

[38] 陳延信,吳 鋒,胡亞茹. 提高粉體堆積密實(shí)度的理論與實(shí)驗(yàn)研究[J]. 煤炭轉(zhuǎn)化, 2012(1): 40-43. Chen Yanxin, Wu Feng, Hu Yaru. Theoretical and experimental study on increasing the bulk density of powder[J]. Coal Conversion, 2012(1): 40-43.(in Chinese)

[39] 龍廣成，謝友均，王新友.礦物摻合料對新拌水泥漿體密實(shí)性能的影響[J].建筑材料學(xué)報, 2002, 5(1): 21-25. Long Guangcheng, Xie Youjun, Wang Xinyou. Effect of mineral admixtures on the compaction property of fresh cement paste[J]. Journal of Building Materials, 2002, 5(1)：21-25.(in Chinese)

[40] Binici H, Aksogan O, Cagatay I H, et al. The effect of particle size distribution on the properties of blended cements incorporating GGBFS and natural pozzolan (NP) [J]. Powder Technology, 2007, 177(3): 140-147.

[41] Givi A N, Rashid S A, Aziz F N A, et al. Assessment of the effects of rice husk ash particle size on strength, water permeability and workability of binary blended concrete[J]. Construction and Building Materials, 2008, 24(11): 2145-2150.

[42] Kwan A K H, Chen J J. Adding fly ash microsphere to improve packing density, flowability and strength of cement paste[J]. Powder Technology, 2013, 234:19-25.

[43] Kwan A K H, Li Y. Effects of fly ash microsphere on rheology, adhesiveness and strength of mortar[J]. Construction and Building Materials, 2013, 42: 137-45.

[44] Ghoddousi P, Javid A A S, Sobhani J. Effects of particle packing density on the stability and rheology of self-consolidating concrete containing mineral admixtures[J]. Construction and Materials, 2014, 53: 102-109.

[45] Nassim S, Mahfoud B. Composition of self compacting concrete (SCC) using the compressible packing model, the chinese method and the european standard[J]. Construction and Building Materials, 2013, 43: 382-328.

[46] 龔建清. 超高性能混凝土的級配效應(yīng)研究[D]. 長沙：湖南大學(xué), 2008. Gong Jianqing. Study on the gradation effect of ultra high performance concrete[D]. Changsha: Hunan University, 2008.(in Chinese)

[47] 曹榮奎. 基于CPM模型和比強(qiáng)度法的超高性能混凝土配合比設(shè)計(jì)研究[D]. 長沙：湖南大學(xué), 2008. Cao Rongkui. Study on mix proportion design of ultra-high performance concrete based on CPM model and specific strength[D]. Changsha: Hunan University, 2008.(in Chinese)

[48] 聶 晶. 基于可壓縮堆積模型的水泥基復(fù)合材料性能研究[D]. 長沙：湖南大學(xué), 2008. Nie Jing. Study on performance of cement-based composite materials based on compressible packing model[D]. Changsha: Hunan University, 2008.(in Chinese)

[49] 李丙明. 低堿度鋼渣配制復(fù)合膠凝材料研究[D]. 西安： 西安建筑科技大學(xué), 2009. Li Bingming. The study on preparing composite cementitious materials with low-alkalinity steel slag[D]. Xi’an: Xi'an University of Architecture and Technology, 2009.(in Chinese)

[50] 丑 凱. 超高性能混凝土堆積密實(shí)度和火山灰效應(yīng)量化研究[D]. 長沙：湖南大學(xué), 2010. Chou Kai. Quantitative analysis on packing density and pozzolanic effect of ultra-high performance concrete[D]. Changsha: Hunan University, 2010.(in Chinese)

[51] 何振偉. 持續(xù)應(yīng)力下密封早齡期混凝土力學(xué)性能依時發(fā)展規(guī)律研究[D]. 北京：北京交通大學(xué), 2013. He Zhenwei. Development of mechanical properties of sealed concrete in early age under sustained stress action[D]. Beijing: Beijing Jiaotong University, 2013.(in Chinese)

[52] 陳瑾祥. 基于可壓縮堆積模型的綠色混凝土材料性能研究[D]. 深圳: 深圳大學(xué), 2015. Chen Jinxiang. Compressible packing model based of green concrete[D]. Shenzhen: Shenzhen University, 2015.(in Chinese)

【中文責(zé)編：坪 梓；英文責(zé)編：之 聿】

Application of particle packing model in concrete

Long Wujian?, Zhou Bo, Liang Peijian, and Sun Roujia

Guangdong Provincial Key Laboratory of Durability for Marine Civil Engineering, College of Civil Engineering,Shenzhen University, Shenzhen 518060, Guangdong Province, P.R.China

The main advantages and limitations of particle packing models, as well as their applications in the field of composite particle systems are pointed out by comparing and analyzing the differences among various particle packing models. The analysis show that, in order to reduce the void volume and achieve the most compact accumulation state of the composite particle system, the particle packing model can be used to calculate the packing density, thus the volume of cementitious materials in the concrete can be reduced and the optimal design of low binder green concrete can be realized.

concrete; particle packing model; packing density; particle system; loosening effect; wall effect; void ratio; mix proportions

Received：2016-08-29；Accepted：2016-11-20

Foundation：National Natural Science Foundation of China (51278306，51578341); Shenzhen Science and Technology Project (JCYJ20140418095735540)

? Corresponding author：Associate professor Long Wujian. E-mail: longwj@szu.edu.cn

：Long Wujian, Zhou Bo, Liang Peijian, et al. Application of particle packing model in concrete[J]. Journal of Shenzhen University Science and Engineering, 2017, 34(1): 63-74.(in Chinese)

TU 528

A

10.3724/SP.J.1249.2017.01063

國家自然科學(xué)基金資助項(xiàng)目(51278306，51578341)；深圳市科技計(jì)劃資助項(xiàng)目(JCYJ20140418095735540)

龍武劍(1977—)，男，深圳大學(xué)副教授、博士. 研究方向：高性能混凝土材料設(shè)計(jì)及性能. E-mail： Longwj@szu.edu.cn

引 文：龍武劍，周 波，梁沛堅(jiān)，等. 顆粒堆積模型在混凝土中的應(yīng)用[J]. 深圳大學(xué)學(xué)報理工版，2017，34(1)：63-74.