江 浩，李榮建, ，閆 蕊，劉軍定，張 媛

（1. 西安理工大學(xué) 巖土工程研究所，陜西 西安 710048；2. 長安大學(xué) 國土資源部巖土工程開放研究實(shí)驗(yàn)室，陜西 西安 710054）

1 引 言

模型試驗(yàn)研究是一種比較直觀的科學(xué)研究方法，可根據(jù)實(shí)際需要進(jìn)行試驗(yàn)設(shè)計(jì)，通過模型試驗(yàn)對(duì)各種復(fù)雜的巖土工程結(jié)構(gòu)的工作性能進(jìn)行試驗(yàn)研究，因此，巖土工程中通過各類模型試驗(yàn)探討土-結(jié)構(gòu)系統(tǒng)力學(xué)響應(yīng)機(jī)制逐漸成為試驗(yàn)研究的重要的發(fā)展方向之一，土工模型試驗(yàn)越來越受到研究者的廣泛關(guān)注[1－4]，并取得了豐富的研究成果。

模型試驗(yàn)要求模型與原型之間遵循一定的相似準(zhǔn)則，使得模型可以再現(xiàn)原型的特性，但模型試驗(yàn)研究的可靠性在基于幾何相似的基礎(chǔ)上才取決于模型相似材料的相似性以及數(shù)據(jù)采集的可靠性。幾何相似比尺、相似材料力學(xué)指標(biāo)與實(shí)際材料符合程度直接關(guān)系到模型試驗(yàn)數(shù)據(jù)的準(zhǔn)確程度，為確保在正確的尺寸縮尺模型中材料和原型材料性狀相似、應(yīng)力和應(yīng)變相等，關(guān)鍵的問題是正確設(shè)計(jì)模型試驗(yàn)相似比尺，選擇和研制更加嚴(yán)格滿足相似比尺的模型材料，才能夠較好地模擬實(shí)際情況[5－7]。

在隧道模型研究方面，陳正發(fā)[8]、周小文等[9]對(duì)隧洞結(jié)構(gòu)在彎曲和軸向荷載作用下的相似條件和計(jì)算依據(jù)進(jìn)行了分析，指出了隧道襯砌結(jié)構(gòu)原型和模型量值之間必需要有一套合理的比例縮放關(guān)系。李榮建等[10－11]指出了模型在力學(xué)響應(yīng)模擬和破壞模擬中具體的模型理論比尺關(guān)系和關(guān)鍵問題，分析了不同層次的相似比尺要求下模型結(jié)構(gòu)材料的選擇與設(shè)計(jì)方法。熊田芳等[12]以西安地鐵二號(hào)線為研究對(duì)象，設(shè)計(jì)了幾何相似比尺1:50的隧道物理模型試驗(yàn)，采用高強(qiáng)度替代材料制作了隧道模型，研究了隧道穿越地裂縫時(shí)襯砌與土體相互作用的力學(xué)響應(yīng)。但是，模型試驗(yàn)中幾何相似比尺的設(shè)計(jì)存在一定問題，可能使得通過隧道模型反演的原型隧道應(yīng)力和彎矩等方面產(chǎn)生較大的誤差，不能合理地模擬實(shí)際情況。

為了評(píng)價(jià)幾何相似差異性對(duì)隧道模型和原型之間應(yīng)力與彎矩的影響程度，本文首先討論隧道結(jié)構(gòu)模型試驗(yàn)的相似比尺關(guān)系，分析模型試驗(yàn)中相似材料的選用原則和要求。然后基于兩種不同的幾何相似程度，分析討論幾何相似程度對(duì)隧道模型和隧道原型結(jié)構(gòu)之間應(yīng)力與彎矩結(jié)果的影響，藉此論證模型試驗(yàn)中幾何相似程度決定了模型相似比尺設(shè)計(jì)的合理性，是模型試驗(yàn)?zāi)芊窈侠怼⒏呔鹊胤从澈湍M原型結(jié)構(gòu)的關(guān)鍵性制約因素之一。

2 隧道模型試驗(yàn)相似比尺與相似材料

2.1 模型試驗(yàn)相似理論

模型試驗(yàn)的理論基礎(chǔ)是相似理論，相似理論既是指導(dǎo)模型試驗(yàn)設(shè)計(jì)的指南，又是分析原型力學(xué)響應(yīng)的依據(jù)。

根據(jù)相似理論確定相似關(guān)系的原理是：若一個(gè)物理過程中共有n個(gè)物理量，k個(gè)基本量綱，則可以組成n-k個(gè)獨(dú)立的無量綱參數(shù)組合，無量綱參數(shù)組合簡稱“π數(shù)”。設(shè)一個(gè)物理過程為式（1），含有n個(gè)參數(shù)X1, X2, …, Xn：

任何一種可以用數(shù)學(xué)方程定義的物理現(xiàn)象都可以用無量綱數(shù)π來定義。相似的物理過程其相應(yīng)的π數(shù)必然相等，因而可根據(jù)上述規(guī)律來確定相似條件，這就是Bockingham π定理。

2.2 模型試驗(yàn)中地鐵隧道相似比尺的理論分析

基于上述Bockingham π定理，確定相似條件的具體方法有方程式分析法和量綱分析法。采用方程式分析法建立相似關(guān)系方便明確，但必須在模型試驗(yàn)設(shè)計(jì)之前對(duì)所研究的物理過程中各物理量之間的函數(shù)關(guān)系，亦即對(duì)試驗(yàn)結(jié)果與試驗(yàn)條件之間的關(guān)系提出明確的數(shù)學(xué)方程式。量綱分析法是根據(jù)方程式量綱均衡性原理進(jìn)行的，僅需明確影響所研究物理現(xiàn)象的物理量以及明確量測這些物理量所需要的量綱即可。

為確保在尺寸減小的模型試驗(yàn)中各種材料和原型材料的性狀相似，試驗(yàn)測定的數(shù)據(jù)有明確的對(duì)應(yīng)原型，原型和模型之間必須有一套相適應(yīng)的比例縮放關(guān)系。因此，相似比的設(shè)計(jì)是整個(gè)模型試驗(yàn)中的一個(gè)關(guān)鍵環(huán)節(jié)。

將隧道襯砌視為薄板結(jié)構(gòu)，用薄板彎曲的控制方程來推導(dǎo)相似比。根據(jù)彈塑性力學(xué)中薄板彎曲的微分方程：

設(shè)板的厚度為h、x、y為2個(gè)互相垂直的坐標(biāo)，橫向均布力為q，薄板的撓曲為w，彎曲剛度為Q，則方程變?yōu)?/p>

根據(jù)相似常數(shù)λ的定義，假定模型與原型有以下關(guān)系：

式中：z為板厚度方向的坐標(biāo)。采用上述相同的方法推導(dǎo)，可以得到原型應(yīng)力：

即為襯砌彎曲作用中地下結(jié)構(gòu)模型與原型應(yīng)力的換算公式。將模型試驗(yàn)所得的彎曲應(yīng)力mσ代入式（20），可推得地下結(jié)構(gòu)原型的彎曲應(yīng)力pσ。

2.3 相似材料的選用

模型試驗(yàn)是采用與原型力學(xué)性質(zhì)相似的材料，按一定關(guān)系制成模型，應(yīng)該具有原型結(jié)構(gòu)的全部或主要特征。在普通重力場內(nèi)，模型試驗(yàn)規(guī)模是有限的，模型中各點(diǎn)的應(yīng)力水平仍然比原型低得多，原型材料的應(yīng)力-應(yīng)變關(guān)系特別是非線性關(guān)系在模型中不能得到真實(shí)體現(xiàn)。通常在通過離心場或氣囊壓力實(shí)施后，這一問題可以得到一定程度的解決。模型試驗(yàn)成果主要依賴于模型相似材料的相似性以及數(shù)據(jù)采集的可靠性，因此，相似材料力學(xué)指標(biāo)與實(shí)際材料符合程度直接關(guān)系到模型試驗(yàn)數(shù)據(jù)的準(zhǔn)確程度。

進(jìn)行相似模型試驗(yàn)必須遵循相似條件，即要求模型不僅幾何形狀相似，而且在模型試驗(yàn)過程中所包含的各物理量或主要物理量應(yīng)與原型相似，由于不同種類的材料各有其特殊性質(zhì)，同一種材料，在應(yīng)力-應(yīng)變關(guān)系曲線中的不同階段，所表現(xiàn)出的力學(xué)性質(zhì)也不盡相同，包括受力條件在內(nèi)的各種外界條件都可以導(dǎo)致材料性質(zhì)的多變性。模型試驗(yàn)前必須對(duì)所模擬的對(duì)象的有關(guān)問題進(jìn)行具體分析，找出合適的相似材料，制作與原型相似的模型。

試驗(yàn)室內(nèi)可利用相似材料，按相似原理制作與原型相似的模型，借助試驗(yàn)儀器、儀表觀測模型的力學(xué)參數(shù)、變形狀態(tài)以及應(yīng)力分布規(guī)律，根據(jù)在模型上研究的成果，推斷原型中可能發(fā)生的力學(xué)現(xiàn)象，達(dá)到利用相似模型的研究來解決工程中的實(shí)際問題的目的。因此，進(jìn)行相似模型研究，相似材料必須滿足下面要求：（1）相似材料的主要力學(xué)性質(zhì)應(yīng)與模擬的原型的結(jié)構(gòu)相似。（2）材料的力學(xué)性能穩(wěn)定，不受外界的影響。（3）容易成型，方便制作。（4）材料來源廣，成本低廉。

3 隧道模型幾何相似程度分析與比較

本文以西安地鐵二號(hào)線區(qū)間隧道為例，作為模型試驗(yàn)中幾何相似程度分析的案例。地鐵隧道直徑為930 cm，內(nèi)襯采用鋼筋混凝土結(jié)構(gòu)，襯砌厚50 cm，模擬的地鐵隧道總長度為90 m。

3.1 縮尺模型的設(shè)計(jì)

文獻(xiàn)[12]中設(shè)計(jì)的模型試驗(yàn)中圍巖試驗(yàn)土樣取自地鐵隧道施工現(xiàn)場，隧道模型襯砌采用有機(jī)玻璃制作，模型襯砌厚度為1 cm，制作模型襯砌的有機(jī)玻璃材料的泊松比為0.372。

值得注意的是：有機(jī)玻璃材料的彈性模量約為混凝土襯砌的 1/10，混凝土襯砌材料泊松比約為0.17，因此，彈性模量和泊松比的相似常數(shù)不為1。

3.2 隧道模型尺寸確定依據(jù)與檢驗(yàn)

為了模型制作和試驗(yàn)方便，隧道模型采用圓形橫斷面，如圖1所示。

圖1 隧道模型示意圖Fig.1 Sketches of tunnel model

文獻(xiàn)[12]中根據(jù)已確定的幾何相似常數(shù)lλ=1/50，確定有機(jī)玻璃隧道模型襯砌結(jié)構(gòu)尺寸。則隧道模型長度為

根據(jù)嚴(yán)格的相似比，采用本文推導(dǎo)的相似比尺公式所計(jì)算的隧道模型襯砌的壁厚為2.07 cm。但文獻(xiàn)[12]則忽略了相差較大的原型襯砌材料與模型襯砌材料的彈性模量和泊松比，并且認(rèn)為這對(duì)襯砌結(jié)構(gòu)應(yīng)力狀態(tài)誤差的影響很小，并求得隧道模型襯砌壁厚為1.0 cm。計(jì)算分析表明，這兩種結(jié)果差異明顯較大，說明了文獻(xiàn)[12]中幾何比尺的設(shè)計(jì)存在問題。

4 隧道模型與原型力學(xué)響應(yīng)差異分析

4.1 隧道模型和隧道原型襯砌的應(yīng)力比較

文獻(xiàn)[12]在模型試驗(yàn)中通過對(duì)土體施加附加壓力的途徑來實(shí)現(xiàn)圍巖應(yīng)力條件的相似比等于1.0的條件，并且認(rèn)為由于試驗(yàn)所模擬的圍巖應(yīng)力場與原型近似一致，則模型與原型的隧道襯砌結(jié)構(gòu)的應(yīng)力一致，模型中襯砌結(jié)構(gòu)的應(yīng)力變化規(guī)律可直接反映到原型，即

根據(jù)模型試驗(yàn)理論，對(duì)隧道原型和通過對(duì)土體施加附加壓力的隧道模型進(jìn)行受力分析。圖2為隧道原型在自重應(yīng)力下的受力分布圖，圖3為隧道模型表面施加附加應(yīng)力的受力分布圖。

圖2 隧道原型的自重應(yīng)力分布Fig.2 Vertical stress distribution of prototype tunnel

圖3 模型表面施加附加應(yīng)力的應(yīng)力分布Fig.3 Equivalence stress distribution in tunnel model

為了驗(yàn)證試驗(yàn)的差異性，參考了文獻(xiàn)[12]某一測點(diǎn)測試值（見圖4），當(dāng)正交地裂縫發(fā)生25 cm錯(cuò)動(dòng)位移時(shí)，在距離地裂縫0.5 m處該測點(diǎn)上盤內(nèi)模型襯砌結(jié)構(gòu)應(yīng)力為-2.006 4 MPa，文獻(xiàn)[12]認(rèn)為該點(diǎn)對(duì)應(yīng)的原型隧道襯砌結(jié)構(gòu)的應(yīng)力pσ= -2.006 4 MPa。

圖4 隧道模型襯砌結(jié)構(gòu)測點(diǎn)布置Fig.4 Arrangement of measuring points inlining of tunnel model

經(jīng)過本文分析，在嚴(yán)格的控制相似比尺關(guān)系的條件下，根據(jù)式（25）計(jì)算得到該點(diǎn)對(duì)應(yīng)的隧道襯砌的原型應(yīng)力pσ′=-2.261 2 MPa。pσ與pσ′這兩者之間的差值為0.254 8 MPa，這個(gè)差值還是比較大的，其分析誤差超過了12%。同理，如果按照式（26）進(jìn)行計(jì)算，結(jié)果差異將會(huì)更加顯著。因此，本文認(rèn)為文獻(xiàn)[12]模型中襯砌結(jié)構(gòu)的應(yīng)力變化規(guī)律可直接反映到原型是不合理的，與實(shí)際情況有較大差異。

4.2 隧道模型和原型襯砌結(jié)構(gòu)彎矩比較

材料力學(xué)彎曲正應(yīng)力的計(jì)算公式為

圖5 隧道模型縱向正應(yīng)力與彎矩分布Fig.5 Distribution of longitudinal normal stress and bending moment of tunnel model

由式（30）可得隧道模型襯砌彎矩 Mm和隧道原型Mp的比例關(guān)系式為

本文考慮到原型襯砌材料與模型襯砌材料的彈性模量和泊松比相差較大時(shí)，在嚴(yán)格控制相似比條件下，對(duì)式（31）進(jìn)行推導(dǎo)以便得出模型和原型的彎矩反演關(guān)系。首先，可令

由此可以得到隧道襯砌原型彎矩 Mp和隧道模型彎矩Mm之間的關(guān)系式，即

對(duì)比式（39）和式（40）可知，本文和文獻(xiàn)[12]在兩種不同幾何相似程度情況下得到的隧道襯砌原型彎矩Mp和隧道模型彎矩Mm之間的關(guān)系式差異很大，其分析誤差達(dá)到了142.1%。

本文提出的隧道襯砌原型彎矩 Mp和隧道模型彎矩 Mm的關(guān)系式具有嚴(yán)格的相似比尺的理論基礎(chǔ)，可以通過隧道模型彎矩正確地反演原型隧道相應(yīng)的彎矩，可以合理地模擬原型隧道的實(shí)際情況。

5 結(jié) 論

（1）為確保在尺寸減小的隧道襯砌模型中各種材料和原型材料的性狀相似，試驗(yàn)測定的數(shù)據(jù)具有明確對(duì)應(yīng)的原型，討論了隧道襯砌結(jié)構(gòu)原型和模型之間量值之間的比例縮放關(guān)系和模型材料的選擇。

（2）分析了兩種幾何縮尺模型，第一種是沒有按照嚴(yán)格相似比進(jìn)行隧道模型幾何尺寸的確定，第二種是嚴(yán)格控制相似比進(jìn)行模型幾何尺寸的確定，這兩種情況計(jì)算的隧道模型襯砌厚度差值較大。分析表明幾何相似是必要條件，不能輕易放松這個(gè)條件，否則數(shù)據(jù)分析的偏差較大而不能較好地模擬實(shí)際原型。

（3）推導(dǎo)了兩種縮尺隧道模型條件下應(yīng)力和彎矩的關(guān)系式，在沒有嚴(yán)格控制相似比和嚴(yán)格控制相似比尺兩種不同條件下得到的隧道原型的應(yīng)力與彎矩差異較大，應(yīng)力的誤差超過12%，彎矩的誤差超過 140%。表明只有嚴(yán)格控制相似比條件下隧道模型才能合理、高精度地反映隧道原型。

隧道襯砌結(jié)構(gòu)的原型材料與模型材料的彈性模量和泊松比相差較大時(shí)，隧道襯砌結(jié)構(gòu)的厚度必須按照嚴(yán)格的相似比來計(jì)算，才能夠較好地模擬實(shí)際情況。

[1] 陳仁朋, 王作洲, 蔣紅光, 等. I型軌道-路基系統(tǒng)動(dòng)力荷載放大系數(shù)模型試驗(yàn)研究[J]. 巖土力學(xué), 2013, 34(4):1045－1052.CHEN Ren-peng, WANG Zuo-zhou, JIANG Hong-guang,et al. Experimental study of dynamic load magnification factor for type I track-subgrade system[J]. Rock and Soil Mechanics, 2013, 34(4): 1045－1052.

[2] 李天斌, 田曉麗, 韓文喜, 等. 預(yù)加固高填方邊坡滑動(dòng)破壞的離心模型試驗(yàn)研究[J]. 巖土力學(xué), 2013, 34(11):3061－3070.LI Tian-bin，TIAN Xiao-li， HAN Wen-xi, et al.Centrifugal model tests on sliding failure of a pile-stabilized high fill slope[J]. Rock and Soil Mechanics, 2013, 34(11): 3061－3070.

[3] 馬慶宏, 朱大勇, 雷先順, 等. 無黏性土斜坡地基承載力模型試驗(yàn)研究[J]. 巖土工程學(xué)報(bào), 2014, 36(7): 1271－1280.MA Qing-hong, ZHU Da-yong, LEI Xian-shun, et al.Model tests on bearing capacity of footing on sand slopes[J]. Chinese Journal of Geotechnical Engineering,2014, 36(7): 1271－1280.

[4] 周健, 周韻鴻, 李飛, 等. 高填黃土明洞土拱效應(yīng)及土壓力減載模型分析[J]. 巖土工程學(xué)報(bào), 2014, 36(3): 555－561.ZHOU Jian, ZHOU Xun-hong, LI Fei, et al. Model analysis of earth pressure load high fill open cut tunnel in loess area[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(3): 555－561.

[5] 范鶴, 劉斌, 范澤, 等. 高填土涵洞相似材料模型試驗(yàn)與數(shù)值模擬[J]. 吉林大學(xué)學(xué)報(bào), 2008, 38(2): 400－403.FAN He, LIU Bin, FAN Ze, et al. Model test with similar material and numerical simulation of culvert with high fills[J]. Journal of Jilin University(Engineering and Technology Edition), 2008,38(2): 400－403.

[6] 羅先啟, 程圣國, 張振華, 等. 地質(zhì)力學(xué)電磁場模型試驗(yàn)相似理論研究[J]. 巖土力學(xué), 2011, 32(4): 1036－1039.LUO Xian-qi, CHENG Sheng-guo, ZHANG Zhen-hua,et al. Study of similarity theory of geomechanical model test in electromagnetic field[J]. Rock and Soil Mechanics, 2011, 32(4): 1036－1039.

[7] 李榮建, 于玉貞, 呂禾, 等. 動(dòng)力離心模型試驗(yàn)中微混凝土抗滑模型樁的設(shè)計(jì)[J]. 清華大學(xué)學(xué)報(bào)(自然科學(xué)版), 2009, 49(9): 66－70.LI Rong-jian, YU Yu-zhen, LU He, et al. Design of model concrete piles to reinforced slopes for dynamic centrifuge tests[J]. Journal of Tsinghua University(Science and Technology), 2009, 49(9): 66－70.

[8] 陳正發(fā). 黏性土地基中地鐵隧道動(dòng)力離心模型試驗(yàn)系統(tǒng)開發(fā)[D]. 北京：清華大學(xué)，2005.

[9] 周小文, 濮家騮. 隧洞結(jié)構(gòu)受力及變形特征的離心模型試驗(yàn)研究[J]. 清華大學(xué)學(xué)報(bào)(自然科學(xué)版), 2001,41(8): 110－116.ZHOU Xiao-wen, PU Jia-liu. Centrifuge model test study of the earth pressure and deformation of tunnel lining[J].Journal of Tsinghua University(Science and Technology), 2001, 41(8): 110－116.

[10] 李榮建, 閆蕊, 劉軍定, 等. 土工動(dòng)力離心模型試驗(yàn)中物理破壞模擬比較及關(guān)鍵問題[J]. 地震工程學(xué)報(bào), 2014,36(3): 718－726.LI Rong-jian, YAN Rui, LIU Jun-ding. A comparison of dynamic response and failure simulations in geotechnical dynamic centrifugal model tests and the key issues[J].China Earthquake Engineering Journal, 2014, 36(3):718－726.

[11] 李榮建, 劉軍定, 李海濤, 等. 土工離心模型試驗(yàn)中動(dòng)力破壞模擬能力分析[J]. 科技導(dǎo)報(bào), 2013, 31(27): 43－48.LI Rong-jian, LIU Jun-ding , LI Hai-tao, et al. Analysis on modeling capability of dynamic failure in centrifugal modeling test[J]. Science & Technology Review, 2013,31(27): 43－48.

[12] 熊田芳, 邵生俊, 王天明, 等. 西安地鐵正交地裂縫隧道的模型試驗(yàn)研究[J]. 巖土力學(xué), 2010, 31(1): 179－186.XIONG Tian-fang，SHAO Sheng-jun，WANG Tian-ming,et al. Research on physical model experiment of Xi’an metro tunnel orthogonally crossing ground fissure[J].Rock and Soil Mechanics, 2010, 31(1): 179－186.