高皓宸　李偉　張旭東

摘要：為了實(shí)現(xiàn)混凝土抗壓強(qiáng)度及彈性模量的精細(xì)化數(shù)值計(jì)算，基于雙尺度理論建立了混凝土雙尺度數(shù)值模型。該理論以有限元方法為載體，通過(guò)正交試驗(yàn)對(duì)混凝土抗壓強(qiáng)度和彈性模量的影響因素進(jìn)行了分析，并進(jìn)行了力學(xué)試驗(yàn)驗(yàn)證。研究結(jié)果表明：1）雙尺度理論在混凝土數(shù)值模擬中的應(yīng)用具有可行性，運(yùn)用該方法的數(shù)值結(jié)果與試驗(yàn)結(jié)果差值在4.22%之內(nèi);2）混凝土的抗壓強(qiáng)度受細(xì)觀參數(shù)影響高度顯著，其中砂漿性能對(duì)整體抗壓強(qiáng)度的影響比例可達(dá)24.35%;3）混凝土的彈性模量受細(xì)觀參數(shù)影響顯著，其中粗骨料性能對(duì)整體彈性模量的影響比例可達(dá)6.55%。研究結(jié)果揭示了混凝土內(nèi)部的粗骨料、砂漿、界面過(guò)渡區(qū)與孔隙對(duì)整體性能的影響效應(yīng)，為優(yōu)化和設(shè)計(jì)新型混凝土材料提供了新思路。

關(guān)鍵詞：非金屬建筑材料;雙尺度理論;有限元方法;混凝土數(shù)值模擬;正交試驗(yàn)

中圖分類(lèi)號(hào)：TU5281文獻(xiàn)標(biāo)識(shí)碼：A

DOI： 10.7535/hbgykj.2022yx02009

Study on optimizing mechanical properties of concrete based on two-scale approach

GAO Haochen LI Wei ZHANG Xudong

（1School of Urban Planning and Municipal Engineering，Xi′an Polytechnic University，Xi′an，Shaanxi 710600，China;2Network Information Center，Xi′an Jiaotong University，Xi′an，Shaanxi? 710049，China）

Abstract：In order to realize the refined calculation of concrete compressive strength and elastic modulus，a two-scale numerical model of concrete was established based on the two-scale approachThe finite element method was taken as the carrier of this approach，the influencing factors of concrete compressive strength and elastic modulus were analyzed through orthogonal experiment，and the result was verified by mechanical testThe results show that：1） the application of two-scale approach in numerical simulation of concrete materials is feasible，and the difference between the numerical results and the experimental results is within 4.22%;2） the compressive strength of concrete is significantly affected by meso parameters，and the influence proportion of mortar on the overall compressive strength is up to 24.35%;3） the elastic modulus of concrete is significantly affected by meso parameters，and the influence proportion of coarse aggregate on the overall elastic modulus is up to 6.55%The effects of coarse aggregate，mortar，interface transition zone and pores in concrete on the overall performance is revealed through the research results，which provide new methods for the optimization and design of new concrete materials

Keywords： non-metallic building material;two-scale approach;finite element method;numerical simulation of concrete;orthogonal experiment

混凝土是土木工程中最常用的材料之一。在宏觀尺度下，混凝土被假定為一種各向同性的均質(zhì)材料;在細(xì)觀尺度下，混凝土表現(xiàn)為由粗骨料、硬化的水泥砂漿、兩者之間的界面過(guò)渡區(qū)（interfacial transition zone，ITZ）以及孔隙等構(gòu)成的多相復(fù)合材料。

WITTMANN等[1]首次提出了“混凝土的細(xì)觀結(jié)構(gòu)”，將基于細(xì)觀結(jié)構(gòu)的研究思路和方法引入到混凝土材料宏觀性能的研究中，該研究為探索混凝土力學(xué)性能提供了一種精細(xì)化方法。劉光廷等[2]提出的隨機(jī)骨料模型已成為混凝土細(xì)觀結(jié)構(gòu)建模使用較為廣泛的方法之一。馬懷發(fā)等[3]采用隨機(jī)骨料法對(duì)鋼筋混凝土梁進(jìn)行了抗彎實(shí)驗(yàn)分析。趙吉坤[4]采用對(duì)砂漿性能弱化的方法來(lái)模擬孔隙和微裂紋對(duì)混凝土宏觀力學(xué)性能及細(xì)觀破壞的影響。宋來(lái)忠等[5]將復(fù)合材料的增強(qiáng)顆粒形狀構(gòu)建為星形，提供了創(chuàng)建更加接近真實(shí)試件的幾何模型方法。田夢(mèng)云等[6]研究了各種因素對(duì)混凝土宏觀應(yīng)力-應(yīng)變關(guān)系和最終破壞模式的影響。

上述研究中，針對(duì)混凝土細(xì)觀結(jié)構(gòu)中的界面過(guò)渡區(qū)以及孔隙實(shí)體對(duì)混凝土整體性能和破壞模式的影響研究還不夠深入，且研究的數(shù)值模擬邊界條件也存在著與實(shí)際混凝土力學(xué)試驗(yàn)不符的問(wèn)題?；诖?，本研究將雙尺度理論作為理論基礎(chǔ)，根據(jù)自編程序模擬四相混凝土的幾何結(jié)構(gòu)，通過(guò)模擬不同工況下的數(shù)值實(shí)現(xiàn)混凝土宏觀等效力學(xué)參數(shù)的計(jì)算。

1混凝土數(shù)值模擬方法

1.1細(xì)觀-宏觀雙尺度理論

混凝土是一種典型的多相復(fù)合材料，根據(jù)LI等[7]的復(fù)合材料雙尺度理論，對(duì)具有兩相或多相復(fù)合材料建立坐標(biāo)系，如圖1所示[7]。

基于復(fù)合材料力學(xué)量的漸進(jìn)展開(kāi)式和攝動(dòng)理論，將雙尺度位移場(chǎng)分解為宏觀和細(xì)觀2個(gè)尺度。

1.2混凝土數(shù)值模型

1.2.1混凝土雙尺度模型

本研究的數(shù)值模擬方法以雙尺度理論為基礎(chǔ)，而雙尺度公式的形式為有限元微分方程組，因此有限元法相對(duì)于離散元法等更適合作為數(shù)值模擬的載體。筆者將包括粗骨料、水泥砂漿、界面過(guò)渡區(qū)與孔隙在內(nèi)的四相模型用于混凝土抗壓強(qiáng)度及彈性模量的有限元數(shù)值模擬?；炷链止橇习葱螤羁煞譃樗槭c卵石。王曉琴等[8]提出不規(guī)則粗骨料顆粒比圓形粗骨料顆粒更能反映真實(shí)骨料的特點(diǎn)，計(jì)算結(jié)果也更接近于宏觀試驗(yàn)。故模型中的碎石骨料被模擬為隨機(jī)多邊形，單個(gè)粗骨料示例如圖2 a）所示。其中，所有粗骨料在整體模型中的分布采用了隨機(jī)均勻算法。同時(shí)，粗骨料的級(jí)配滿(mǎn)足二維模型常用的Walraven表達(dá)式，如式（11）所示，且基于文獻(xiàn)[9]的算法保證骨料之間不會(huì)發(fā)生相交情況。

混凝土內(nèi)部常常含有孔隙或微裂縫，它們對(duì)混凝土力學(xué)性能和失效破壞模式會(huì)產(chǎn)生影響。陳青青等[10]提出圓形、橢圓形、多邊形孔隙對(duì)計(jì)算精度影響有限，而圓形孔隙的計(jì)算效率明顯高于其他兩種。因此模型中的孔隙被模擬為圓形，孔隙在混凝土模型中滿(mǎn)足隨機(jī)均勻分布。

水泥砂漿在混凝土中起著膠結(jié)作用?；炷恋钠茐耐l(fā)生于砂漿內(nèi)部以及砂漿與骨料的交界處，砂漿強(qiáng)度性能對(duì)混凝土力學(xué)特性影響高度顯著[11]。

RAO等[12]認(rèn)為在砂漿與骨料的交界處存在ITZ，這是由于壁效應(yīng)以及骨料顆粒周?chē)嗨a(chǎn)生的孔隙而形成的。ITZ的材質(zhì)與砂漿差別不大，但是其內(nèi)部疏松多孔，因此ITZ容易成為混凝土整體的薄弱部位。ITZ幾何由已建成的粗骨料模型進(jìn)行縮放而成，其厚度為100 μm，如圖2 b）所示。

1.2.2邊界條件與網(wǎng)格劃分

雙尺度模型處于單軸受壓狀態(tài)時(shí)的邊界條件應(yīng)設(shè)置為底邊中點(diǎn)的橫向與縱向位移約束、底邊其余點(diǎn)的縱向位移約束以及頂邊的縱向位移荷載[13]，如圖3所示。

網(wǎng)格劃分及細(xì)部結(jié)構(gòu)如圖3所示，網(wǎng)格單元為四邊形為主的線(xiàn)性平面應(yīng)變單元。其中，對(duì)于橫截面尺寸沿縱軸不變、受有平行于橫截面而且不沿縱軸變化的面力或約束的物體，平面應(yīng)變模型可以將其簡(jiǎn)化為二維模型進(jìn)行分析，相較于三維模型提高了計(jì)算效率。

1.3材料參數(shù)

基于文獻(xiàn)[14]—[16]中的相關(guān)規(guī)范參數(shù)賦予混凝土各組成相力學(xué)參數(shù)。一般情況下，粗骨料的抗壓強(qiáng)度遠(yuǎn)高于水泥砂漿和界面過(guò)渡區(qū)，數(shù)值模擬時(shí)，粗骨料被假定為線(xiàn)彈性材料[17]，粗骨料的參數(shù)如表1所示。水泥砂漿屬于脆性材料，其抗拉強(qiáng)度遠(yuǎn)小于其抗壓強(qiáng)度，因此砂漿的塑性損傷本構(gòu)關(guān)系應(yīng)分為抗壓與抗拉參數(shù)[5-6]，水泥砂漿的強(qiáng)度等級(jí)如表2所示。界面過(guò)渡區(qū)只存在于骨料與砂漿表面的交界處，其力學(xué)性能一般為砂漿的60%～80%[18-19]。界面過(guò)渡區(qū)的力學(xué)性能可以等價(jià)為水泥砂漿力學(xué)參數(shù)按一定比例的“弱化”，筆者將水泥砂漿與界面過(guò)渡區(qū)之間的弱化系數(shù)設(shè)置為08，如表3所示。

1.4數(shù)值模擬方法總流程

基于雙尺度理論和有限元方法的混凝土雙尺度模型數(shù)值模擬總流程圖如圖4所示，總體可分為以下4步。

步驟1：輸入模型的細(xì)觀組成相參數(shù)，進(jìn)而生成各組成相的幾何模型。

步驟2：將生成的各組成相幾何模型逐一隨機(jī)投放至指定區(qū)域，同時(shí)進(jìn)行碰撞檢測(cè)判斷，直至各組成相的占比達(dá)到輸入的預(yù)期值。

步驟3：完成幾何結(jié)構(gòu)建模，并對(duì)整個(gè)模型進(jìn)行網(wǎng)格劃分，賦予各組成相力學(xué)參數(shù)及設(shè)置邊界條件等步驟。

步驟4：將模型提交計(jì)算，基于雙尺度方法進(jìn)行后處理，輸出所期望的混凝土宏觀力學(xué)參數(shù)。

2混凝土數(shù)值模擬正交試驗(yàn)設(shè)計(jì)

基于前述數(shù)值模擬方法，混凝土內(nèi)部的多種因素都有可能對(duì)混凝土抗壓強(qiáng)度、彈性模量產(chǎn)生影響。正交試驗(yàn)可以高效分析多因素多水平對(duì)混凝土抗壓強(qiáng)度、彈性模量的影響效應(yīng)，正交試驗(yàn)中設(shè)置各因素的原因如下。

1）粗骨料的面積分?jǐn)?shù)（因素A）可以反映出實(shí)際混凝土配合比的改變對(duì)整體強(qiáng)度性能的影響，粗骨料占比的改變也就意味著砂漿及界面過(guò)渡區(qū)面積分?jǐn)?shù)的相對(duì)改變。

2）混凝土內(nèi)部的孔隙率（因素B）可以反映出實(shí)際混凝土中存在的初始缺陷對(duì)整體性能的影響，孔隙率衡量著混凝土中是否存在初始缺陷及初始缺陷大小。

3）粗骨料的種類(lèi)（因素C）可以反映出實(shí)際混凝土中選用不同種類(lèi)的粗骨料對(duì)整體力學(xué)性能的影響。

4）水泥砂漿的性能（因素D）可以反映出實(shí)際混凝土中的水泥砂漿采用不同的強(qiáng)度等級(jí)對(duì)整體力學(xué)性能的影響。

正交試驗(yàn)采用了L（4）正交表，各因素與水平的設(shè)置如表4所示。其中，因素C和D涉及了不同水平的材料屬性，具體如表1—3所示。

3結(jié)果與討論

3.1混凝土數(shù)值模擬結(jié)果與試驗(yàn)驗(yàn)證

基于正交試驗(yàn)完成數(shù)值模擬的全流程后，通過(guò)后處理得到的各工況下混凝土的應(yīng)力-應(yīng)變曲線(xiàn)如圖5所示。

為了驗(yàn)證數(shù)值模擬方法的準(zhǔn)確性，嚴(yán)格按照模擬參數(shù)設(shè)計(jì)配合比進(jìn)行混凝土立方體抗壓強(qiáng)度試驗(yàn)。試驗(yàn)所用的電液伺服壓力試驗(yàn)機(jī)型號(hào)為YNS-Y2000，試驗(yàn)中的工況2，3，4代表試塊的破壞情況如圖6所示。相同工況下的力學(xué)試驗(yàn)結(jié)果與數(shù)值模擬結(jié)果的對(duì)比情況如圖7所示。

試塊的破壞機(jī)理可以根據(jù)第二強(qiáng)度理論進(jìn)行分析，任意應(yīng)力狀態(tài)下，只要ε達(dá)到極限值ε/E，材料就發(fā)生斷裂，如式（12）所示。

混凝土單軸受壓時(shí)，若試驗(yàn)機(jī)的承臺(tái)與試塊的接觸面潤(rùn)滑良好，試塊將沿垂直于壓力的方向開(kāi)裂，這正是ε的方向，這與圖6 b）至圖6 d）中試塊的縱向或斜縱向開(kāi)裂情況較為一致。斜縱向的開(kāi)裂模式也與混凝土內(nèi)部的粗骨料分布有關(guān)，當(dāng)裂縫發(fā)展的路徑中存在粗骨料時(shí)，裂縫會(huì)繞開(kāi)粗骨料，并沿著粗骨料表面的界面過(guò)渡區(qū)繼續(xù)開(kāi)展，若界面過(guò)渡區(qū)縱向的鄰近位置存在孔隙，則裂縫會(huì)優(yōu)先向著容易產(chǎn)生應(yīng)力集中的孔隙處開(kāi)展。

圖7中，試驗(yàn)測(cè)得的工況2，3，4下混凝土的抗壓強(qiáng)度分別為17.96，20.14，20.97 MPa，其與相同工況下數(shù)值模擬結(jié)果的差值分別為2.39%，4.20%，4.22%;試驗(yàn)測(cè)得的工況2，3，4下混凝土的彈性模量分別為26 632，29 217，31 180 MPa，其與相同工況下數(shù)值模擬結(jié)果的差值分別為1.78%，1.58%，1.40%。試驗(yàn)與數(shù)值模擬得出的單軸受壓應(yīng)力-應(yīng)變曲線(xiàn)在彈性階段擬合較好，在工程精度下，可認(rèn)為該數(shù)值模擬方法具有較高的準(zhǔn)確性。此外，該數(shù)值模擬的有效性也可通過(guò)與傳統(tǒng)有限元方法進(jìn)行對(duì)比得到驗(yàn)證?？紤]到該數(shù)值模型設(shè)置了多種細(xì)觀參數(shù)，在計(jì)算之前是未知其宏觀性能的，但是傳統(tǒng)有限元方法計(jì)算時(shí)需要已知混凝土的宏觀性能，兩者較難通過(guò)控制變量進(jìn)行對(duì)比。

工況2試塊在單軸受壓數(shù)值模擬中的應(yīng)力-應(yīng)變曲線(xiàn)如圖8所示。由于混凝土的抗拉強(qiáng)度遠(yuǎn)小于其抗壓強(qiáng)度，其破壞模式往往表現(xiàn)為受拉開(kāi)裂破壞，[JP4]圖8曲線(xiàn)中不同階段的特征點(diǎn)處的受拉損傷情況如圖9所示。

圖8和圖9中，混凝土所受壓應(yīng)力達(dá)到臨界點(diǎn)（a）時(shí)，損傷首先產(chǎn)生在界面過(guò)渡區(qū)表面，而后沿著鄰近的孔隙貫通展開(kāi)，孔隙處產(chǎn)生了明顯的應(yīng)力集中現(xiàn)象?；炷了軌簯?yīng)力達(dá)到峰值點(diǎn)（b）時(shí)，損傷的路徑已發(fā)展于骨料之間的砂漿上，并且在兩側(cè)邊附近以及骨料分布密集的部位出現(xiàn)了較密集的損傷。隨著加載繼續(xù)，損傷達(dá)到限值并逐漸發(fā)展成混凝土內(nèi)部的微裂縫。混凝土所受壓應(yīng)力達(dá)到拐點(diǎn)（c）時(shí)，裂縫繼續(xù)擴(kuò)展、貫通，并逐漸發(fā)展為密集的宏觀裂縫，內(nèi)部結(jié)構(gòu)賴(lài)以傳遞荷載的路徑不斷減少?；炷了軌簯?yīng)力達(dá)到收斂點(diǎn)（d）時(shí)，宏觀裂縫分布密集且已貫穿了混凝土整體，其中較為明顯的是開(kāi)展在大粒徑粗骨料兩側(cè)的“倒V型”裂縫，同時(shí)在模型的兩側(cè)向也產(chǎn)生了密集的縱向裂縫，此時(shí)混凝土已失去繼續(xù)承壓能力。

圖9 d）收斂點(diǎn)（d）處的損傷情況與圖6 b）中試驗(yàn)的破壞情況對(duì)比顯示：兩者均在中部出現(xiàn)了明顯的縱向斜裂縫，同時(shí)兩側(cè)邊附近產(chǎn)生了密集的縱向裂縫，可見(jiàn)數(shù)值模擬損傷情況與試驗(yàn)開(kāi)裂破壞情況具有較好的一致性。

3.2結(jié)果分析

3.2.1正交試驗(yàn)結(jié)果

基于正交試驗(yàn)對(duì)混凝土雙尺度模型進(jìn)行不同工況的數(shù)值模擬后，所得結(jié)果如表5所示。

3.2.2極差分析

基于正交試驗(yàn)的極差分析結(jié)果如表6與表7所示。

由表6可知，4種因素對(duì)混凝土的抗壓強(qiáng)度影響程度從大到小依次為水泥砂漿的性能（因素D）、混凝土內(nèi)部的孔隙率（因素B）、粗骨料的面積分?jǐn)?shù)（因素A）、粗骨料的種類(lèi)（因素C）。

由表7可知，4種因素對(duì)混凝土的彈性模量影響程度從大到小依次為水泥砂漿的性能（因素D）、粗骨料的面積分?jǐn)?shù)（因素C）、粗骨料的種類(lèi)（因素A）、混凝土內(nèi)部的孔隙率（因素B）。

粗骨料的面積分?jǐn)?shù)由30%提高至35%，混凝土的抗壓強(qiáng)度提高了5.42%，彈性模量提高了2.78%;粗骨料的面積分?jǐn)?shù)由35%提高至40%，混凝土的抗壓強(qiáng)度提高了0.36%，彈性模量提高了6.37%;粗骨料的面積分?jǐn)?shù)由40%提高至45%，混凝土的抗壓強(qiáng)度提高了3.47%，彈性模量提高了3.02%?？梢?jiàn)，粗骨料在混凝土中起骨架作用，其具有較大的剛度與強(qiáng)度，因此粗骨料在混凝土內(nèi)部被分配了較大的應(yīng)力，從而改善了其余組成相的受力情況。

孔隙率由0%提高至1%，混凝土的抗壓強(qiáng)度減少了7.05%，混凝土的彈性模量減少了0.65%;孔隙率由1%提高至2%，混凝土的抗壓強(qiáng)度減少了5.96%，混凝土的彈性模量減少了1.75%;孔隙率由2%提高至3%，混凝土的抗壓強(qiáng)度減少了5.02%，混凝土的彈性模量減少了0.28%?？梢?jiàn)，隨著混凝土孔隙率的減少，混凝土的密實(shí)度提高，混凝土的抗壓強(qiáng)度也隨之提高。同時(shí)，孔隙率的提高會(huì)對(duì)混凝土的剛度產(chǎn)生輕微的削弱效果。

粗骨料種類(lèi)由大理巖變?yōu)槭規(guī)r，混凝土的抗壓強(qiáng)度減少了2.39%，混凝土的彈性模量增加了6.55%;粗骨料種類(lèi)由石灰?guī)r變?yōu)榛◢弾r，混凝土的抗壓強(qiáng)度減少了0.66%，混凝土的彈性模量增加了4.59%;粗骨料種類(lèi)由花崗巖變?yōu)樾鋷r，混凝土的抗壓強(qiáng)度減少了3.03%，混凝土的彈性模量增加了3.54%?？梢?jiàn)，混凝土的性能受到粗骨料種類(lèi)的影響，且粗骨料種類(lèi)對(duì)彈性模量的影響較為明顯。

水泥砂漿的性能由M15提高至M20，混凝土的抗壓強(qiáng)度提高了24.35%，彈性模量提高了7.75%;水泥砂漿的性能由M20提高至M25，混凝土的抗壓強(qiáng)度提高了15.10%，彈性模量提高了4.52%;水泥砂漿的性能由M25提高至M30，混凝土的抗壓強(qiáng)度提高了14.26%，彈性模量提高了4.60%?？梢?jiàn)，水泥砂漿是混凝土的基體，其力學(xué)性能極大程度上決定了混凝土整體的力學(xué)性能。此因素可看作水泥砂漿與ITZ對(duì)混凝土力學(xué)性能的共同影響。

效應(yīng)曲線(xiàn)圖是對(duì)極差分析表的直觀描述與總結(jié)，如圖10所示。

3.2.3方差分析

基于正交試驗(yàn)的方差分析結(jié)果如表8與表9所示。

如表8所示，對(duì)于抗壓強(qiáng)度指標(biāo)：水泥砂漿性能達(dá)到了高度顯著水平（α=0.01），孔隙率達(dá)到了顯著水平（α=0.05），粗骨料面積分?jǐn)?shù)和粗骨料種類(lèi)沒(méi)有達(dá)到顯著水平。

如表9所示，對(duì)于彈性模量指標(biāo)：水泥砂漿性能達(dá)到了高度顯著水平（α=0.01），粗骨料種類(lèi)達(dá)到了高度顯著水平（α=0.01），粗骨料面積分?jǐn)?shù)達(dá)到了高度顯著水平（α=0.01），孔隙率沒(méi)有達(dá)到顯著水平。

4結(jié)語(yǔ)

將雙尺度理論作為理論基礎(chǔ)，通過(guò)有限元方法完成了混凝土雙尺度模型的數(shù)值模擬全過(guò)程，達(dá)到了對(duì)混凝土宏觀力學(xué)性能精細(xì)化計(jì)算與優(yōu)化的目的，并基于數(shù)值模擬參數(shù)進(jìn)行了混凝土立方體抗壓強(qiáng)度試驗(yàn)來(lái)驗(yàn)證模擬的準(zhǔn)確性，結(jié)論如下。

1）基于混凝土雙尺度模型進(jìn)行的數(shù)值模擬可以預(yù)測(cè)混凝土的宏觀力學(xué)性能。該方法已實(shí)現(xiàn)參數(shù)化、程序化，其參數(shù)可根據(jù)研究的需要進(jìn)行靈活調(diào)整。

2）雙尺度理論和有限元方法首次被結(jié)合用于混凝土等效力學(xué)性能的精細(xì)化計(jì)算中，力學(xué)試驗(yàn)結(jié)果證明該路徑可行。

3）混凝土的宏觀力學(xué)性能極大程度上受到細(xì)觀組成相的結(jié)構(gòu)與性能的影響。粗骨料的作用更多在于增大混凝土剛度，其對(duì)抗壓強(qiáng)度的影響有限;砂漿的強(qiáng)度等級(jí)與混凝土的抗壓強(qiáng)度直接相關(guān)，砂漿性能決定了混凝土力學(xué)性能的上限;界面過(guò)渡區(qū)與孔隙是損傷與裂縫容易產(chǎn)生的部位，兩者決定了混凝土力學(xué)性能的下限。

目前混凝土中粗骨料的參數(shù)采用了線(xiàn)彈性本構(gòu)關(guān)系，但實(shí)際中粗骨料是一種彈脆性材料，未來(lái)的精細(xì)化計(jì)算應(yīng)該將粗骨料參數(shù)模擬為更加真實(shí)的本構(gòu)關(guān)系。另外，目前雙尺度理論在混凝土數(shù)值模擬中的應(yīng)用僅限于材料階段，未來(lái)雙尺度理論在混凝土數(shù)值模擬中的應(yīng)用應(yīng)向著結(jié)構(gòu)階段發(fā)展。

參考文獻(xiàn)/References：

[1]WITTMANN F H，ROELFSTRA P E，SADOUKI HSimulation and analysis of composite structures[J].Materials Science and Engineering，1985，68（2）：239-248.

[2]劉光廷，王宗敏用隨機(jī)骨料模型數(shù)值模擬混凝土材料的斷裂[J].清華大學(xué)學(xué)報(bào)（自然科學(xué)版），1996，36（1）：84-89LIU Guangting，WANG ZongminNumerical simulation study of fracture of concrete materials using random aggregate model[J]. Journal of Tsinghua University （Science and Technology），1996，36（1）：84-89.

[3]馬懷發(fā)，陳厚群，黎保琨混凝土試件細(xì)觀結(jié)構(gòu)的數(shù)值模擬[J].水利學(xué)報(bào)，2004，35（10）：27-35MA Huaifa，CHEN Houqun，LI BaokunMeso-structure numerical simulation of concrete specimens[J].Journal of Hydraulic Engineering，2004，35（10）：27-35.

[4]趙吉坤混凝土四相復(fù)合模型的三維細(xì)觀破壞模擬[J].土木建筑與環(huán)境工程，2009，31（4）：37-43ZHAO Jikun3D meso-scale failure simulation of four-phase composite concrete[J].Journal of Civil，Architectural & Environmental Engineering，2009，31（4）：37-43.

[5]宋來(lái)忠，廖大乾混凝土細(xì)觀結(jié)構(gòu)的數(shù)值模擬研究[J].三峽大學(xué)學(xué)報(bào)（自然科學(xué)版），2018，40（6）：56-60SONG Laizhong，LIAO DaqianResearch on numerical simulation of meso-structure of concrete[J].Journal of China Three Gorges University（Natural Sciences），2018，40（6）：56-60.

[6]田夢(mèng)云，張恩，曹瑞東，等基于細(xì)觀尺度的混凝土單軸力學(xué)性能仿真計(jì)算分析[J].應(yīng)用力學(xué)學(xué)報(bào)，2020，37（3）：975-981TIAN Mengyun，ZHANG En，CAO Ruidong，et alMeso-scale simulation analysis of uniaxial mechanical behavior of concrete[J].Chinese Journal of Applied Mechanics，2020，37（3）：975-981.

[7]LI W，JIA Y J，LI L XA two-scale approach to numerically predict the strength and degradation of composites[J].Applied Composite Materials，2018，25（1）：85-103.

[8]王曉琴，楊名超，肖明，等自密實(shí)混凝土單軸壓縮試驗(yàn)細(xì)觀模擬[J]中國(guó)科技論文，2020，15（12）：1410-1416WANG Xiaoqin，YANG Mingchao，XIAO Ming，et alMeso-scale simulation on uniaxial compression test of self-compacting concrete[J].China Sciencepaper，2020，15（12）：1410-1416.

[9]ERICSON C實(shí)時(shí)碰撞檢測(cè)算法技術(shù)[M].劉天慧譯北京：清華大學(xué)出版社，2010.

[10]陳青青，張煜航，張杰，等含孔隙混凝土二維細(xì)觀建模方法研究[J].應(yīng)用數(shù)學(xué)和力學(xué)，2020，41（2）：182-194CHEN Qingqing，ZHANG Yuhang，ZHANG Jie，et al Study on a 2D mesoscopic modeling method for concrete with voids[J].Applied Mathematics and Mechanics，2020，41（2）：182-194.

[11]CHANDRA S，BERNTSSON LLightweight Aggregate Concrete Microstructure[M]. Norwich：William Andrew Publishing，2002.

[12]RAO G A，PRASAD B K PInfluence of type of aggregate and surface roughness on the interface fracture properties[J].Materials and Structures，2004，37（5）：328-334.

[13]CHEN Hongbing，XU Bin，WANG Jiang，et alParametric analysis on compressive strain rate effect of concrete usingmesoscale modeling approach[J].Construction and Building Materials，2020，246：118375.

[14]《工程地質(zhì)手冊(cè)》編委會(huì)工程地質(zhì)手冊(cè)[M].北京：中國(guó)建筑工業(yè)出版社，2007.

[15]杜修力，王陽(yáng)，路德春水泥砂漿材料的靜動(dòng)態(tài)應(yīng)力應(yīng)變關(guān)系研究[J].土木工程學(xué)報(bào)，2010，43（sup2）：119-126DU Xiuli，WANG Yang，LU DechunStudy on static and dynamic stress-strain relationship of cement mortar material[J].China Civil Engineering Journal，2010，43（sup2）：119-126.

[16]GB 175—2007通用硅酸鹽水泥[S].

[17]WRIGGERS P，MOFTAH S OMesoscale models for concrete：Homogenisation and damage behaviour[J].Finite Elements in Analysis and Design，2006，42（7）：623-636.

[18]ASBRIDGE A H，PAGE C L，PAGE M MEffects of metakaolin，water/binder ratio and interfacial transition zones on the microhardness of cement mortars[J].Cement and Concrete Research，2002，32（9）：1365-1369.

[19]LIAO Kuoyu，CHANG Pingkun，PENGA Y N，et alA study on characteristics of interfacial transition zone in concrete[J].Cement and Concrete Research，2004，34（6）：977-989.