金立兵 梁新亞 霍承鼎 王振清 王珍

摘要：地下糧倉可充分利用地下空間，具有較好的防火、防毒、防爆等性能，并能利用淺層地能實現(xiàn)準低溫儲糧，具有節(jié)能、低損、保證糧食品質的優(yōu)點，對保障糧食安全及可持續(xù)發(fā)展意義重大。采用工程足尺試驗與數(shù)值分析相結合的研究方法，對大型地下混凝土筒倉竣工驗收前倉壁的力學性能進行研究，通過對比試驗結果與數(shù)值模擬結果，驗證數(shù)值分析方法的合理性與有效性，進而對倉壁在最不利荷載工況下的力學性能進行數(shù)值分析。結果表明：地下混凝土筒倉倉壁內力的數(shù)值模擬結果與試驗結果吻合較好;由于倉壁較厚且與倉底和倉頂間接觸部位剛度較大，導致倉壁內外側徑向應力表現(xiàn)出相反的變化規(guī)律;倉壁徑向應力在倉壁底部位置處最高，環(huán)向應力在距倉壁頂部約2/3位置處最高;在最不利荷載工況下，倉壁徑向應力和環(huán)向應力隨深度表現(xiàn)出與實際工況下相似的變化規(guī)律，且相同深度下應力較大。

關鍵詞：地下筒倉;足尺試驗;數(shù)值分析;倉壁;力學性能

中圖分類號：TU926 文獻標志碼：A 文章編號：20966717（2020）03004006

Abstract：

The underground silos can take full advantages of the underground space， which is of better fireproof， anti-toxicity， explosionproof and other merits. In the meantime， the underground silos can use the shallow geothermal energy to realize quasilow temperature storage. With the advantages of energysaving， lowconsumption and grain qualityensuring， it plays a significant role in ensuring the grain safety and sustainable development of China. The mechanical properties of the silo wall before completion acceptance of large underground concrete silos were studied by the combination of engineering fullscale test and numerical analysis. By comparing the test results with the engineeringscale testing results，， the rationality and effectiveness of the numerical analysis method were verified. Furthermore， the mechanical properties of the silo wall under the most unfavorable load conditions were numerically analyzed. The results show that the internal forces of the numerical simulation results of the silo wall were in good agreement with that of experimental results. Due to the relatively large thickness of the silo wall， and the rigidity of the contact part between the silo wall and the silo bottom and the silo roof， the radial stress on the inside and outside of the silo wall shows the opposite law. The radial stress of the silo wall is highest at the bottom of the silo wall， and the hoop stress is highest at about 2/3 of the position from the top of the silo wall. Under the most unfavorable load conditions， the radial stress and hoop stress of the silo wall show similar changes with the burial depth as in the actual working condition， and the stress is greater at the same depth.

Keywords：underground silo; full scale test; numerical analysis; silo wall; mechanical properties

糧食安全是關系國民經濟發(fā)展、社會穩(wěn)定和國家自立的全局性重大戰(zhàn)略問題。地下糧倉可以利用淺層地能實現(xiàn)糧食準低溫儲藏，具有保證糧食品質、節(jié)地、節(jié)能、無污染、全壽命周期成本低等突出優(yōu)點[12]。中國從仰韶文化的原始社會就已經開始應用地下糧倉（窖）進行糧食儲存，在漢、隋、唐、明、清等朝代也均有地下倉窖的遺跡被發(fā)現(xiàn)[3]。至今，中國多數(shù)的?。ㄖ陛犑小⒆灾螀^(qū)）都建有地下糧倉。但是，已建成的地下糧倉存在倉容較小、機械化程度較低、交通不便、建設受地域限制等缺點，不便于日常生產管理和使用[4]。

其他國家對地下空間的開發(fā)多用于儲備石油、天然氣等戰(zhàn)略物資以及解決城區(qū)不斷擴張、人口不斷增加所造成的土地使用競爭激烈的問題[58]。同時，對地下工程用高強度混凝土材料、既有地下工程耐久性能評估與維護改造等方面進行了大量的研究[912]。中國對地下糧倉的研究主要集中在新倉型的探究、倉壁的受力性能分析與結構優(yōu)化以及維護結構傳熱規(guī)律等方面 [1317]。研究多通過模型試驗與數(shù)值分析完成，結果缺少相應工程足尺試驗以及其他理論方法和實踐經驗支持。由于數(shù)值分析結果受人為因素影響較大，故取得結果的有效性有待驗證。

筆者基于工程足尺試驗，采用試驗研究與數(shù)值分析相結合的方法，對地下大型混凝土筒倉在竣工驗收前倉壁的受力情況進行分析研究，從而提出一種能有效分析地下混凝土筒倉受力的數(shù)值分析模型。

1工程試驗

1.1工程概況

試驗倉位于河南省鄭州市中牟縣萬三路與萬洪路交叉口的金地集團糧食物流園區(qū)內。工程場地的地質情況如表1所示。結構選用強度等級為C40的自防水混凝土，抗?jié)B等級為P12，倉底下部選用C20的素混凝土墊層，倉壁厚度為300 mm。鋼筋選用HRB400E級鋼筋。倉體外徑25.6 m，倉底深19.56 m，倉容3 500 t，倉頂覆土深度1.5 m，所處場地的地下水類型為潛水，歷史最高水位-1.0 m[18]。

1.2倉壁內力檢測

1.2.1測點布置

在試驗倉倉壁內部鋼筋上沿深度方向每隔2 m位置對稱安裝鋼筋計，用于測量壁內鋼筋的徑向應力和環(huán)向應力;在倉壁外部沿深度方向每隔2 m位置安裝壓力計，用于測量倉壁所受到的來自外部的水土復合壓力。鋼筋計和壓力計分布如圖1所示。

1.2.2試驗數(shù)據(jù)采集及結果

試驗選用基康儀器股份有限公司和金土木工程儀器股份有限公司生產的鋼筋計和壓力計，型號、量程以及精度見表2。

2基于實際工況的倉壁數(shù)值分析

2.1模型建立與邊界條件

利用有限元軟件ABAQUS對地下混凝土筒倉結構建模，由于倉壁的厚高比遠小于1/10，符合殼單元的構造要求。所以，模型的倉頂蓋、倉壁、倉底、倉內核心支撐筒均選用S4R殼單元建模，倉頂群梁選用梁單元建模。由于倉體的施工采用滑模形式，因此，建立模型時將倉底與倉壁之間、倉壁與倉頂蓋之間、倉頂蓋與中間支撐筒之間、倉底與中間支撐筒之間均設置為剛接。并對地下混凝土筒倉抗浮懸臂部位設置z向約束，模型如圖4所示。

2.2參數(shù)選擇與施加荷載

混凝土強度等級C40，密度取2 400 kg/m3，根據(jù)《混凝土結構設計規(guī)范》（GB 50010—2015），軸心抗壓強度設計值取19.1 N/mm2，混凝土抗拉強度設計值取1.71 N/mm2，鋼筋混凝土彈性模量取3.25×104 N/mm2，泊松比取0.2。

根據(jù)現(xiàn)場測試結果顯示的地下水位約位于-5.88 m，因此，造成了以地下水位為分界點，倉壁所受側壓力隨深度變化曲線的斜率發(fā)生較大變化。為方便對地下倉模型施加荷載，取深度為-5.76 m設置的測點為分界點，對倉壁壓力計實測側壓力的數(shù)據(jù)分兩段進行線性擬合，擬合公式為

2.3結果對比分析

通過有限元軟件模擬計算得到倉壁的內力云圖，結果如圖5所示。提取倉壁沿z向對應鋼筋計位置的環(huán)向和徑向應力與試驗數(shù)據(jù)進行對比分析，結果如圖6所示。

采用擬合側壓力作為荷載進行倉壁內力計算，由于擬合壓力值與實測值在-7.760 m測點處偏差相對其他測點較大，因此，該位置處的環(huán)向應力實測值與計算值存在相對較大的偏差，其余3個測點的實測值與理論值都吻合較好，驗證了模型的有效性。

根據(jù)倉體結構和受力形式，倉底與倉壁、倉壁與倉頂接觸部位可視為剛接，同時，在這兩個連接部位設有環(huán)梁，剛度較大。而且倉壁厚度相對較大，造成倉壁徑向應力呈現(xiàn)出外側與內側規(guī)律相反的情況;而環(huán)向則由于不受彎矩影響或影響較小，兩側的環(huán)向應力總體呈現(xiàn)出變化規(guī)律基本一致的狀態(tài)。

通過分析得到，徑向應力在倉壁外側的頂部位置最先表現(xiàn)為壓應力，隨著深度的增加，壓應力逐漸減小，約在距倉壁頂部1/6位置處轉變?yōu)槔瓚?，并在距倉壁頂部2/3處拉應力達到最大，后隨深度增加，再次轉變?yōu)閴簯Γ译S深度增加而增大;徑向應力在倉壁內側頂部至底部所顯現(xiàn)的規(guī)律與外側呈相反規(guī)律。徑向應力在倉壁底部位置最高。

環(huán)向應力在倉壁外側的頂部位置最先表現(xiàn)為較小的壓應力，隨著深度的增加，壓應力逐漸變大，大約在距倉頂2/3位置處，壓應力達到最大，后隨深度增大持續(xù)減小;環(huán)向應力在倉壁內側頂部至底部應力變化規(guī)律同外側基本一致，不同點為倉壁內側的頂部以及底部應力表現(xiàn)為拉應力。環(huán)向應力在距倉頂2/3位置最高。

3最不利荷載工況下倉壁數(shù)值分析

3.1施加荷載

以未裝糧的試驗倉在場地地下水達到歷史最高水位下所受到水土復合壓力作為最不利荷載工況[19]。在此荷載工況下對倉壁內力進行研究。根據(jù)地質情況，對倉壁所受水土復合壓力進行計算。根據(jù)倉壁側壓力取值的相關研究，靜止土壓力與主動土壓力大小相差不大，而被動土壓力與二者相差較大[11]。因此，對靜止土壓力和被動土壓力進行計算，其中，σ1為靜止土壓力，σ2為被動土壓力，如圖7所示。

3.2模擬結果

通過有限元模擬得到工況二倉壁的內力云圖如圖8所示。倉壁沿z向對應鋼筋計位置的環(huán)向應力與徑向應力結果如圖9所示。

通過分析數(shù)值模擬結果可以看出，在最不利荷載工況下，倉壁徑向應力和環(huán)向應力隨深度表現(xiàn)出與實際工況下基本一致的變化規(guī)律，且相同深度下應力較大。倉壁徑向應力最大值在倉壁底部位置，環(huán)向應力最大位置在距倉壁頂部約2/3處。

4結論

通過工程足尺試驗與數(shù)值分析結合的方法對地下混凝土筒倉竣工驗收前倉壁力學性能進行研究，并對最不利荷載工況下倉壁的受力進行數(shù)值分析，得到以下結論：

1）通過分析對比模擬結果與工程試驗結果，二者能夠較好吻合，驗證了所建立模型的合理性。

2）分別對不同工況下倉壁的力學性能進行數(shù)值分析，結果表明，倉壁的徑向應力在倉壁底部位置最高，環(huán)向應力在距倉壁頂部約2/3位置處最高。

3）考慮到大型地下混凝土筒倉在施工以及使用過程中容易受到各施工、運輸機械等荷載作用造成地面超載，從而對地下糧倉結構安全性能造成不利影響，建議在設計時采用被動土壓力作為設計用土壓力。參考文獻：

[1] 鄭培， 王振清， 余漢華. 我國地下倉的簡史及特點[J]. 糧油加工， 2009（1）： 9094.

ZHENG P， WANG Z Q， YU H H. Brief history and characteristics of underground granaries in China[J]. Cereals and Oils Processing， 2009（1）： 9094. （in Chinese）

[2] ?PEI Y B， FENG Y. The underground grain storage environment regionalization [J]. Advanced Materials Research， 2013， 807/808/809： 20432047.

[3] 陳桂香， 岳龍飛， 王振清， 等. 鋼筋混凝土地下糧倉準靜態(tài)溫度場數(shù)值模擬[J]. 中國糧油學報， 2014， 29（3）： 7983.

CHEN G X， YUE L F， WANG Z Q， et al.Quasistatic temperature field simulation of reinforced concrete underground warehouse grain [J]. Journal of the Chinese Cereals and Oils Association， 2014， 29（3）： 7983. （in Chinese）

[4] 余漢華， 王錄民， 王振清， 等. 我國地下糧倉應用的現(xiàn)狀及前瞻[J]. 河南工業(yè)大學學報（自然科學版）， 2008， 29（6）： 7981.

YU H H， WANG L M， WANG Z Q， et al.Present condition and outlook for underground silos [J]. Journal of Henan University of Technology（Natural Science Edition）， 2008， 29（6）： 7981. （in Chinese）

[5] ?ZEVGOLIS I， MAVRIKOS A， KALIAMPAKOS D. Construction， storage capacity and economics of an underground warehousing – logistics center in Athens， Greece [J]. Tunnelling and Underground Space Technology， 2004， 19（2）： 165173.

[6] ?GENTILE R J. Geology and utilization of underground space in metropolitan Kansas City area， USA [J]. Environmental Geology， 1997， 29（1/2）： 1116.

[7] ?DELMASTRO C， LAVAGNO E， SCHRANZ L. Energy and underground [J]. Tunnelling and Underground Space Technology， 2016， 55： 96102.

[8] ?YU KULIKOVA E. Defects of urban underground structure and their prediction [J]. IOP Conference Series： Materials Science and Engineering， 2018， 451（1）： 012108.

[9] ?MIN K H， JUNG H C， YANG J M， et al. Shrinkage characteristics of highstrength concrete for large underground space structures [J]. Tunnelling and Underground Space Technology， 2010， 25（2）： 108113.

[10] ?WINKLER B， HOFSTETTER G， LEHAR H. Application of a constitutive model for concrete to the analysis of a precast segmental tunnel lining [J]. International Journal for Numerical and Analytical Methods in Geomechanics， 2004， 28（78）： 797819.

[11] ?CURIELESPARZA J， CANTOPERELLO J， CALVO M A. Establishing sustainable strategies in urban underground engineering [J]. Science and Engineering Ethics， 2004， 10（3）： 523530.

[12] ?GHAFARI N. Corrosion control in underground concrete structures using double waterproofing shield system （DWS） [J]. International Journal of Mining Science and Technology， 2013， 23（4）： 603611.

[13] 鄭培， 劉銀來， 王振清. 我國地下糧倉新倉型初探及其有限元分析[J]. 河南工業(yè)大學學報（自然科學版）， 2009， 30（3）： 9194.

ZHENG P， LIU Y L， WANG Z Q. Brief introduction and finite element analysis of new underground warehouse in China [J]. Journal of Henan University of Technology（Natural Science Edition），2009， 30（3）： 9194. （in Chinese）

[14] ?王海濤， 李向前. 基于Galerkin法的地下糧倉圍護結構傳熱特性數(shù)值模擬[J]. 農業(yè)工程學報， 2014， 30（12）： 241248.

WANG H T， LI X Q.Numerical simulation of heat transfer characteristics of building envelopes of underground warehouse based on elementfree Galerkin method [J]. Transactions of the CSAE， 2014， 30（12）： 241248. （in Chinese）

[15] 賈騰南， 王振清， 丁永剛. 施工階段鋼筋砼地下糧食筒倉倉壁受力分析[J]. 河南科學， 2014， 32（2）： 216219.

JIA T N， WANG Z Q， DING Y G. Wall stress analysis of reinforced concrete underground grain silos in construction stage [J]. Henan Science， 2014， 32（2）： 216219. （in Chinese）

[16] 熊曉莉， 金立兵， 史笑珂， 等. 鋼筋混凝土圓形地下糧倉倉壁承載力計算方法研究[J]. 河南工業(yè)大學學報（自然科學版）， 2015， 36（1）： 9599.

XIONG X L， JIN L B， SHI X K， et al. Study on wall bearing capacity calculation method of reinforced concrete underground granary [J]. Journal of Henan University of Technology（Natural Science Edition）， 2015， 36（1）： 9599. （in Chinese）

[17] ?熊曉莉， 金立兵， 王振清. 鋼筋混凝土地下糧倉倉壁土壓力取值方法及倉壁結構受力分析[J]. 應用基礎與工程科學學報， 2016， 24（1）： 103114.

XIONG X L， JIN L B， WANG Z Q. Earth pressure and bearing capacity analysis on the wall of reinforced concrete underground granary [J]. Journal of Basic Science and Engineering， 2016， 24（1）： 103114. （in Chinese）

[18] 金立兵， 胡潁， 祁繼鵬. 基于遺傳算法的地下鋼筋混凝土筒倉倉壁的結構優(yōu)化[J]. 河南工業(yè)大學學報（自然科學版）， 2017， 38（3）： 99102， 122.

JIN L B， HU Y， QI J P. Structural optimization of underground reinforced concrete silo wall based on genetic algorithm [J]. Journal of Henan University of Technology（Natural Science Edition）， 2017， 38（3）： 99102， 122. （in Chinese）

[19] 涂成順. 地下大直徑鋼筋混凝土筒倉倉壁結構分析及比較[D]. 鄭州： 河南工業(yè)大學， 2011.

TU C S. Wall structure analysis and comparison for largediameter reinforced concrete underground silo [D]. Zhengzhou： Henan University of Technology， 2011.（in Chinese）

（編輯王秀玲）