地震作用下筒倉結(jié)構(gòu)貯料側(cè)壓力計算方法

周長東 張泳 邱意坤 梁立燦 阿斯哈

摘? ?要：為了對筒倉結(jié)構(gòu)進(jìn)行合理的抗震設(shè)計，得到實際工程中實用的貯料側(cè)壓力，引入“施衛(wèi)星”地震作用下貯料側(cè)壓力計算方法，以工程設(shè)計中能夠獲得的材料參數(shù)為基礎(chǔ)，確定貯料側(cè)壓力計算中所需要的修正參數(shù);然后利用ABAQUS建立筒倉結(jié)構(gòu)數(shù)值模型進(jìn)行計算分析，根據(jù)獲得的貯料側(cè)壓力變化曲線擬合得出相應(yīng)參數(shù)的修正值;綜合不同參數(shù)，歸納出考慮筒倉-貯料相互作用的貯料側(cè)壓力修正公式;最后，將貯料側(cè)壓力修正公式所得的計算值與既有試驗數(shù)據(jù)及規(guī)范設(shè)計值進(jìn)行對比分析. 結(jié)果表明，采用本文公式得到的計算值與試驗值吻合較好，能夠真實合理地反映實際地震作用下貯料側(cè)壓力分布;考慮強(qiáng)震作用下貯料側(cè)壓力按照規(guī)范設(shè)計值設(shè)計偏于危險，而按照本文計算公式設(shè)計更加經(jīng)濟(jì)可靠.

關(guān)鍵詞：鋼筋混凝土筒倉;地震作用;貯料側(cè)壓力;計算方法;抗震設(shè)計

中圖分類號：TU279.7? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼：A

Calculating Method on Storage Side Pressure

of Silo Structures under Earthquake

ZHOU Changdong?，ZHANG Yong，QIU Yikun，LIANG Lican，A Siha

（School of Civil Engineering，Beijing Jiaotong University，Beijing 100044，China）

Abstract：In order to get reasonable seismic design for silo structures and obtain the practical calculated value of storage side pressure， the “Shi Weixing” calculation method of storage side pressure under the earthquake is introduced at first. Based on the material parameters which can be obtained in engineering design， the revised parameters needed in the calculation of storage side pressure are determined. Then， the numerical model of silo structure is established with ABAQUS to calculate the side pressure variation curve and the correction values of the corresponding parameters are obtained by curve fitting. Modified formulas of side pressure considering silo-storage interaction are summarized by synthesizing different parameters. Finally， the calculated values obtained from the modified formula of storage side pressure are compared with the existing test data and the design values of the code. The results show that the calculated values are in good agreement with the experimental values and can reflect the real and reasonable storage side pressure. Compared with the results， the storage side pressure under the action of strong earthquake is dangerous according to the design value of the code， which means that the proposed calculation formula in this paper is more economical and reliable.

Key words：reinforced concrete silo;seismic action;storage side pressure;calculation method;seismic design

在鋼筋混凝土筒倉的抗震設(shè)計中，通常將筒倉當(dāng)作懸臂結(jié)構(gòu)來計算整個結(jié)構(gòu)根部的內(nèi)力，并采用截面強(qiáng)度設(shè)計，忽略了貯料側(cè)壓力對筒倉壁的影響. 實際上，該側(cè)壓力變化將在倉壁的圓周方向上引起彎曲應(yīng)力，這對于筒倉這類薄壁結(jié)構(gòu)是不利的.

近年來國內(nèi)外學(xué)者對于筒倉貯料側(cè)壓力進(jìn)行了廣泛的研究. 施衛(wèi)星等[1-2]對2個1 ∶ 10鋼筋混凝土圓筒倉結(jié)構(gòu)模型進(jìn)行了振動臺試驗，研究了地震作用下筒倉的動力響應(yīng)和破壞形式，根據(jù)試驗結(jié)果提出了儲煤對倉壁側(cè)壓力計算方法; Holler等[3]對動態(tài)激發(fā)下的存儲顆粒材料的筒倉進(jìn)行了有限元數(shù)值模擬，并將采用該模型獲得的數(shù)值模擬結(jié)果與在法國Saclay的振動臺試驗結(jié)果進(jìn)行了比較分析;張翀等[4-5]分析了各種卸料工況下倉壁應(yīng)變的變化規(guī)律;趙松[6]進(jìn)行了鋼筒倉側(cè)壓力模型試驗，得到了筒倉結(jié)構(gòu)靜止及卸料狀態(tài)下的貯料側(cè)壓力分布曲線;周長東等[7]針對筒倉靜力作用提出了一種散粒體的亞塑性本構(gòu)模型，分析了貯料對倉壁的靜止側(cè)壓力;龐照昆等[8]通過模型試驗方法，研究了筒倉在側(cè)壁卸料與中心卸料的動態(tài)超壓現(xiàn)象.

目前的研究大多集中于靜力狀態(tài)以及卸料狀態(tài)下貯料側(cè)壓力的變化[9-13]，主要通過理論分析法、有限單元法以及離散元法進(jìn)行分析，而少有對于地震作用下貯料側(cè)壓力分布的研究. 考慮地震作用時，離散元法雖然可以較好地模擬顆粒之間的相互作用，但輸入地震波的困難讓該方法難以應(yīng)用;有限元法數(shù)值計算的精度及穩(wěn)定性高，但是由于其對特定結(jié)構(gòu)進(jìn)行數(shù)值分析，并不能得出一個普遍的貯料側(cè)壓力計算方法;理論分析法在地震作用下受多種因素的限制，并不能單獨通過公式推導(dǎo)得出所需的計算參數(shù)[14].

本文首先對地震作用下考慮筒倉-貯料相互作用的貯料側(cè)壓力計算方法進(jìn)行理論分析，確定貯料側(cè)壓力計算中所需要的修正參數(shù);然后利用ABAQUS建立筒倉結(jié)構(gòu)數(shù)值模型進(jìn)行計算，獲取與影響參數(shù)相應(yīng)的貯料側(cè)壓力變化曲線;隨后對貯料側(cè)壓力變化曲線進(jìn)行擬合，得到相應(yīng)參數(shù)在貯料側(cè)壓力計算中的修正值;之后綜合不同參數(shù)擬合結(jié)果，歸納出考慮筒倉-貯料相互作用的貯料側(cè)壓力修正公式;最后，將根據(jù)貯料側(cè)壓力修正公式所得的計算值與既有試驗數(shù)據(jù)及規(guī)范設(shè)計值進(jìn)行對比分析，驗證計算公式的合理性.

1? ?地震作用下貯料側(cè)壓力計算

目前針對地震作用下貯料側(cè)壓力計算方法的研究中，僅有施衛(wèi)星等通過試驗與理論相結(jié)合的形式，推導(dǎo)出符合相應(yīng)試驗?zāi)Ｐ偷牡卣鹱饔孟沦A料側(cè)壓力計算方法[2]. 因此，本文結(jié)合數(shù)值模擬的方法，在“施衛(wèi)星”計算方法的基礎(chǔ)上進(jìn)行修正，得到可在實際工程應(yīng)用的地震作用下貯料側(cè)壓力計算方法.

1.1? ?“施衛(wèi)星”計算方法

施衛(wèi)星等[2]提出以計算筒倉貯料質(zhì)心處的煤側(cè)壓力為主，再通過各種修正系數(shù)去計算貯料上層及底層的側(cè)壓力的計算方法.

當(dāng)?shù)卣鹱饔糜谕矀}時，在其貯料質(zhì)心處的絕對加速度為Sa，對稱壓力Pce不可能使貯料產(chǎn)生加速度，而偏心壓力Pee（θ）使貯料產(chǎn)生了Sa加速度，假設(shè) Pee（θ）的分布滿足：

考慮到貯料質(zhì)心加速度是貯料側(cè)壓力及貯料層間摩擦力共同作用的結(jié)果，因此使得計算所得側(cè)壓力與實際值有所差別，施衛(wèi)星對該公式進(jìn)行了修正：

式中：Ce由試驗結(jié)果確定，與測點位置、臺面加速度、貯料內(nèi)摩擦角和濕度等因素有關(guān).

為求解出工程設(shè)計時考慮的最大貯料側(cè)壓力Pmax，需先計算最小貯料側(cè)壓力Pmin：

式中：h為貯料高度;Cp為超壓系數(shù);Ph為Janssen靜壓力.

綜上，考慮筒倉-貯料相互作用的貯料側(cè)壓力計算方法需要確定6個參數(shù)：筒倉半徑、貯料質(zhì)心處加速度、計算截面位置、貯料容重、內(nèi)摩擦角、貯料與倉壁摩擦因數(shù).

1.2? ?簡化計算方法

應(yīng)用公式（3）進(jìn)行地震作用下貯料側(cè)壓力計算時，需要輸入貯料質(zhì)心處加速度. 而工程設(shè)計時，無法直接得到貯料質(zhì)心加速度，可通過綜合考慮測點位置、貯料內(nèi)摩擦角及摩擦因數(shù)等因素，對臺面加速度進(jìn)行相應(yīng)修正后[2]，得到：

式中：Cg為考慮地震強(qiáng)度的修正系數(shù);Ch為考慮貯料高度的修正系數(shù);C?為考慮貯料內(nèi)摩擦角的修正系數(shù);Cμ為考慮摩擦因數(shù)的修正系數(shù).

在求解出地震作用下貯料側(cè)壓力變化值Pe后，需要確定地震作用下最小貯料側(cè)壓力Pmin來確定最大貯料側(cè)壓力Pmax. 根據(jù)文獻(xiàn)[15]中提供的筒倉在靜力狀態(tài)以及地震作用下的貯料側(cè)壓力數(shù)值模擬分析結(jié)果可看出，筒倉結(jié)構(gòu)在地震作用下的最小貯料側(cè)壓力Pmin < 靜力狀態(tài)下貯料側(cè)壓力，與施衛(wèi)星定義下的Pmin計算公式相悖. 在工程設(shè)計時通常將最大貯料側(cè)壓力Pmax作為主要設(shè)計影響因素，因此對地震作用下最大貯料側(cè)壓力Pmax進(jìn)行簡化計算，綜合考慮式（3）與式（7）對最大貯料側(cè)壓力Pmax進(jìn)行修正，即：

式中：Cg，max為考慮地震強(qiáng)度的最大貯料側(cè)壓力修正系數(shù);Ch，max為考慮貯料高度的最大貯料側(cè)壓力修正系數(shù); C?，max為考慮貯料內(nèi)摩擦角的最大貯料側(cè)壓力修正系數(shù);Cμ，max為考慮摩擦因數(shù)的最大貯料側(cè)壓力修正系數(shù).

筒倉結(jié)構(gòu)設(shè)計時，根據(jù)結(jié)構(gòu)所在地區(qū)的環(huán)境參數(shù)和設(shè)定的材料參數(shù)，建立有限元分析模型，選取合理的地震波進(jìn)行筒倉結(jié)構(gòu)的地震響應(yīng)分析，獲取公式（9）中的修正系數(shù)，即可計算得到設(shè)計所需的最大貯料側(cè)壓力.

2? ?筒倉結(jié)構(gòu)有限元分析模型

2.1? ?工程概況

以浙江某糧庫的筒倉結(jié)構(gòu)為工程背景，將鋼筋混凝土筒倉群中的單一筒倉作為設(shè)計模擬的原型（如圖1所示）. 建筑物總高度35.0 m，糧食裝載高度27.0 m，筒倉單倉內(nèi)徑12.0 m，筒倉壁厚0.22 m;環(huán)梁尺寸為0.4 m × 0.8 m，門洞尺寸為1.8 m × 2.7 m，窗洞尺寸為0.9 m × 1.5 m. 抗震設(shè)防烈度為7度，設(shè)計地震組為第一組. 抗震設(shè)防類別為乙類，框架抗震設(shè)防等級為二級，場地類別為Ⅲ類[16].

2.2? ?ABAQUS有限元模型

采用ABAQUS軟件對所選取的單倉結(jié)構(gòu)建立分析模型（見圖1）. 選用混凝土損傷塑性本構(gòu)模擬倉壁的混凝土材料（見表1），混凝土的單軸本構(gòu)曲線采用《混凝土結(jié)構(gòu)設(shè)計規(guī)范》（GB 50010—2010）附表C.2中的推薦公式;選用理想彈塑性本構(gòu)模擬鋼筋（見表2）;選用摩爾庫倫本構(gòu)模擬貯料（見表3）.

利用ABAQUS中S4單元模擬筒倉倉壁，利用C3D8R單元模擬貯料;利用ABAQUS自帶的“接觸對”來模擬倉壁與貯料的接觸關(guān)系，其中選取倉壁內(nèi)表面為主接觸面，兩者之間的摩擦因數(shù)為0.4.

2.3? ?模型驗證

由模型試驗[14]知，滿倉狀態(tài)下，在地震烈度為7度多遇和7度基本時，上海人工波SHW2對單倉y方向影響較大，因此在筒倉數(shù)值模型底部施加SHW2波進(jìn)行動力時程分析. 為驗證單倉模型的合理性，提取單倉模型y方向單元最大加速度以及最大相對位移、自振頻率值，與相應(yīng)的筒倉結(jié)構(gòu)模型試驗結(jié)果進(jìn)行對比，對比結(jié)果詳見圖2和圖3.

通過有限元模擬分析，得到滿倉狀態(tài)下數(shù)值模型的自振頻率為2.617 Hz. 振動臺試驗獲取單倉一階頻率后，通過相似系數(shù)返回原型結(jié)構(gòu)得到的原型結(jié)構(gòu)自振頻率為2.631 Hz. 如圖2和圖3所示，將有限元模擬提取的最大加速度和最大相對位移值與單倉試驗數(shù)據(jù)進(jìn)行比較，偏差均在5%之內(nèi). 因此，本文建立的單倉數(shù)值模型與振動臺試驗?zāi)Ｐ臀呛陷^好，證明了數(shù)值模型的合理性.

3? ?最大貯料側(cè)壓力計算參數(shù)分析

采用增量動力分析（Incremental Dynamic Analysis，IDA）方法，針對不同地震波的峰值地面加速度（Peak Ground Acceleration，PGA）進(jìn)行調(diào)幅并在數(shù)值模型底部輸入，得到時程分析結(jié)果. 為獲取式（9）中的相應(yīng)參數(shù)，需要擬合不同工況下貯料側(cè)壓力求解修正參數(shù). 在進(jìn)行數(shù)值模擬計算時，單一模型輸入的參數(shù)包括了地震動強(qiáng)度、貯料高度、內(nèi)摩擦角以及摩擦因數(shù)，因此獲得的貯料側(cè)壓力分布曲線不能直接通過擬合得到單一的修正參數(shù). 本文假定當(dāng)PGA = 0.1g、? = 28°、μ = 0.40條件下的筒倉貯料側(cè)壓力修正系數(shù)為1，即Cg = 0.1，max = 1、 C? = 28°，max = 1、Cμ = 0.4，max = 1，在此條件下，通過變換不同參數(shù)進(jìn)行數(shù)值模擬，將所得貯料側(cè)壓力分布曲線進(jìn)行擬合.

3.1? ?地震波的選取

根據(jù)原型筒倉的工程概況，按《建筑抗震設(shè)計規(guī)范》[17]的要求選取兩條天然波和一條人工波進(jìn)行時程分析，分別選用El-Centro波、NGA#178號地震波以及上海人工地震波SHW2波輸入數(shù)值模型底部進(jìn)行研究. 圖4所示為所選地震波的平均反應(yīng)譜，與建筑抗震設(shè)計規(guī)范中的設(shè)計反應(yīng)譜吻合較好.

圖5、圖6、圖7分別為考慮地震動強(qiáng)度、內(nèi)摩擦角以及摩擦因數(shù)影響下的最大貯料側(cè)壓力分布曲線，可發(fā)現(xiàn)貯料高度在不同參數(shù)中對最大貯料側(cè)壓力的影響最為明顯. 因此，首先考慮以PGA = 0.1g時筒倉最大貯料側(cè)壓力曲線為基準(zhǔn)，對考慮貯料高度的最大貯料側(cè)壓力參數(shù)Ch，max進(jìn)行擬合;隨后對不同參數(shù)影響下最大貯料側(cè)壓力分布曲線進(jìn)行擬合，確定不同參數(shù)影響下的最大貯料側(cè)壓力修正系數(shù).

3.2? ?考慮貯料高度的最大貯料側(cè)壓力參數(shù)

采用PGA = 0.1g、貯料參數(shù)為? = 28°、 μ = 0.40的筒倉模型進(jìn)行數(shù)值模擬計算，選用此參數(shù)可排除地震動強(qiáng)度、內(nèi)摩擦角及貯料與倉壁的摩擦因數(shù)的影響. 為求解考慮高度的最大貯料側(cè)壓力參數(shù)，將式（8）、式（9）進(jìn)行變形，得：

式中：Pd，max是倉壁受到的最大動壓力.

圖8所示為筒倉結(jié)構(gòu)在PGA = 0.1g地震作用下沿高度方向修正系數(shù)分布曲線. 可以看出，筒倉沿高度方向的修正系數(shù)曲線可分為四段，分別為漏斗處、筒倉倉壁相對高度0 ～ 0.25、相對高度0.25 ～ 0.75、相對高度0.75 ～ 1.0. 選取筒倉倉壁相對高度為0、0.25、0.5、0.75、1處的修正系數(shù)進(jìn)行分段擬合，如圖9所示.

3.3? ?考慮地震動強(qiáng)度的最大貯料側(cè)壓力參數(shù)

考慮貯料高度對筒倉側(cè)壓力計算的影響后， 采用貯料參數(shù)為? = 28°、 μ = 0.40的筒倉模型進(jìn)行數(shù)值模擬計算，選用此參數(shù)可以排除內(nèi)摩擦角、貯料與倉壁的摩擦因數(shù)的影響. 為求解考慮地震動強(qiáng)度的最大貯料側(cè)壓力參數(shù)，將式（8）（9）進(jìn)行變形，得：

圖10所示是筒倉結(jié)構(gòu)在不同地震動強(qiáng)度下的 Cg，max分布曲線. 由圖10可知，筒倉結(jié)構(gòu)在不同高度處Cg，max的變化規(guī)律不同. 圖11分別對位于筒倉漏斗處、筒倉倉壁相對高度為0、0.25、0.5、0.72、1處不同地震動強(qiáng)度的Cg，max進(jìn)行擬合，可以發(fā)現(xiàn)：在相對高度>0的筒倉倉壁部位，PGA = 0.1g ～ 0.5g時筒倉地震強(qiáng)度修正系數(shù)呈線性分布，PGA = 0.6g ～ 1.0g時筒倉地震強(qiáng)度修正系數(shù)的分布形式為二次函數(shù);在筒倉結(jié)構(gòu)漏斗處，PGA = 0.1g ～ 0.4g時筒倉結(jié)構(gòu)地震強(qiáng)度修正系數(shù)呈線性分布，PGA = 0.5g ～ 1.0g時修正系數(shù)分布形式為斜率不同的直線. 因此，分段討論分布函數(shù)公式.

3.4? ?考慮內(nèi)摩擦角的最大貯料側(cè)壓力參數(shù)

為研究貯料之間的不同摩擦因數(shù)對貯料側(cè)壓力的影響，考慮選用? = 25°、28°、30°、33°、35°五種內(nèi)摩擦角進(jìn)行計算. 本節(jié)采用貯料參數(shù)為PGA = 0.1g、 μ = 0.40的筒倉模型進(jìn)行數(shù)值模擬計算，選用此參數(shù)可以排除地震動強(qiáng)度、貯料與倉壁的摩擦因數(shù)的影響. 為求解考慮內(nèi)摩擦角的最大貯料側(cè)壓力參數(shù)，將式（8）、式（9）進(jìn)行變形，得：

圖12是筒倉結(jié)構(gòu)在PGA = 0.1g的地震作用下不同摩擦角修正系數(shù)分布曲線. 由圖12可知，隨著貯料內(nèi)摩擦角的增大，考慮內(nèi)摩擦角的最大貯料側(cè)壓力參數(shù)逐漸減小;當(dāng)位于筒倉結(jié)構(gòu)漏斗處時，隨著內(nèi)摩擦角的變化，考慮內(nèi)摩擦角的最大貯料側(cè)壓力參數(shù)變化量基本相同;當(dāng)位于筒倉結(jié)構(gòu)中上部時，相同高度處隨著內(nèi)摩擦角的變化，考慮內(nèi)摩擦角的最大貯料側(cè)壓力參數(shù)變化量有較大差異. 因此，對該參數(shù)進(jìn)行分段考慮，分別選取不同內(nèi)摩擦角的漏斗處修正系數(shù)及筒倉倉壁相對高度為0、0.25、0.5、0.72、1處修正系數(shù)平均值進(jìn)行擬合，如圖13所示.

3.5? ?考慮摩擦因數(shù)的最大貯料側(cè)壓力參數(shù)

為研究貯料與倉壁之間不同摩擦因數(shù)對貯料側(cè)壓力的影響，考慮采用μ = 0.30、0.35、0.40、0.45、0.50五種摩擦因數(shù)進(jìn)行計算. 采用貯料參數(shù)為PGA = 0.1 g、? = 28°的筒倉模型進(jìn)行數(shù)值模擬計算，選用此參數(shù)可以排除地震動強(qiáng)度、貯料內(nèi)摩擦角的影響. 為求解考慮摩擦因數(shù)的最大貯料側(cè)壓力參數(shù)，將式（8）、式（9）進(jìn)行變形，得：

圖14所示是筒倉結(jié)構(gòu)在PGA = 0.1g地震作用下不同摩擦因數(shù)的修正系數(shù)分布曲線. 由圖14可知，隨著筒倉倉壁與貯料摩擦因數(shù)的增大，Cμ，max逐漸減小;當(dāng)位于筒倉結(jié)構(gòu)漏斗處時，隨著摩擦因數(shù)的變化，Cμ，max的變化量基本相同;當(dāng)位于筒倉結(jié)構(gòu)中下部時，摩擦因數(shù)引起的筒倉結(jié)構(gòu)Cμ，max的變化量隨著高度的增加基本不變;當(dāng)位于筒倉結(jié)構(gòu)上部時，摩擦因數(shù)引起的筒倉結(jié)構(gòu)Cμ，max的變化量隨著高度的增加而減小. 因此，考慮將高度與摩擦因數(shù)的影響結(jié)合，進(jìn)行分段考慮.

4? ?計算值與試驗數(shù)據(jù)對比

王錄民等[18]通過模擬地震振動臺試驗對鋼筋混凝土單倉及群倉模型進(jìn)行研究，觀察了貯料在地震作用下的側(cè)壓力變化規(guī)律. 根據(jù)第3節(jié)中得到的地震作用下修正的貯料側(cè)壓力計算公式，結(jié)合試驗中的數(shù)據(jù)進(jìn)行對比，對比分析結(jié)果如圖17所示.

如圖17所示，7度罕遇和8度罕遇地震作用下，根據(jù)簡化的最大貯料側(cè)壓力計算公式得到的模型結(jié)構(gòu)貯料側(cè)壓力計算值與試驗數(shù)據(jù)吻合較好;7度基本地震作用下筒倉結(jié)構(gòu)中部計算值與試驗值有所差異，其原因在于模型試驗時，在地震強(qiáng)度較小時，中部貯料由于散粒體的固結(jié)作用未能充分運動，因此試驗值相對于計算值偏小;在7度基本、7度罕遇、8度罕遇地震作用下貯料底部及中上部計算值相對于試驗數(shù)據(jù)稍微偏大，因此使用簡化的最大貯料側(cè)壓力計算公式能合理反映實際地震作用下側(cè)壓力分布，保證了采用該公式設(shè)計時結(jié)構(gòu)的安全性.

5? ?計算值與規(guī)范設(shè)計值對比

目前，比較經(jīng)典的筒倉壓力分布理論主要有Janssen理論、Rankine理論、Airy理論以及修正的Coulomb理論等. 我國在進(jìn)行深倉貯料壓力計算[19]時，以Janssen公式為理論基礎(chǔ)，同時在設(shè)計時考慮相應(yīng)的放大系數(shù)：

式中：Ph為高度h處倉壁受到的法向側(cè)壓力;γ為貯料的重力密度，kN/m3;μ為筒倉倉壁與貯料的摩擦因數(shù);ρ為筒倉結(jié)構(gòu)的水力半徑，m;k為主動側(cè)壓力系數(shù);Ch為深倉貯料水平壓力修正系數(shù).

在考慮筒倉側(cè)壓力靜力計算時，取Ch = 1進(jìn)行計算;在進(jìn)行筒倉結(jié)構(gòu)設(shè)計時，需進(jìn)行深倉貯料水平側(cè)壓力修正，見表4[19].

修正系數(shù)Ch 備注（參數(shù)） 倉壁 1. 當(dāng)筒倉貯料高度hn與筒倉內(nèi)徑dn的比大于3.0時，Ch應(yīng)乘以1.1;

2. 流動性較差的貯料，Ch宜乘以0.9;

3. 具有特殊促流裝置的筒倉，Ch宜乘以1.1～1.4. ]

將經(jīng)過第3節(jié)相關(guān)系數(shù)修正后得到的地震作用下貯料側(cè)壓力簡化計算公式所得的計算結(jié)果與我國規(guī)范對筒倉貯料壓力設(shè)計值進(jìn)行對比，如圖18所示. 可以發(fā)現(xiàn)，貯料壓力設(shè)計值隨高度增大，側(cè)壓力逐漸減小. 當(dāng)高度h小于1 m時，貯料側(cè)壓力簡化計算值隨高度增大而減小;當(dāng)1 m ≤ h ≤ 1.5 m時，計算值隨高度增大而增大;當(dāng)h > 1.5 m時，計算值隨高度增大而減小. 這是由于規(guī)范設(shè)計值是通過放大系數(shù)對靜力作用下貯料側(cè)壓力進(jìn)行修正;而地震作用下，由于筒倉結(jié)構(gòu)中下部填充密實，貯料與筒倉倉壁相互作用較小;筒倉結(jié)構(gòu)中上部貯料則由于填充密實度較差，因此動力作用下與筒倉倉壁相互作用較大，造成筒倉結(jié)構(gòu)中上部出現(xiàn)貯料側(cè)壓力增大的情況. 筒倉結(jié)構(gòu)頂部則由于貯料較少，因此在h > 1.5 m時，側(cè)壓力隨高度增大而減小.

如圖18所示，在7度基本烈度地震作用下，筒倉結(jié)構(gòu)中下部動態(tài)側(cè)壓力均小于規(guī)范設(shè)計值，可見在7度抗震設(shè)防地區(qū)按照規(guī)范設(shè)計值對倉壁進(jìn)行承載力計算偏于保守;而筒倉中上部動態(tài)側(cè)壓力計算值均大于規(guī)范設(shè)計值，按照規(guī)范設(shè)計筒倉中上部偏于危險;在7度罕遇烈度及8度罕遇烈度下，筒倉動態(tài)側(cè)壓力均大于規(guī)范設(shè)計值，可見在考慮強(qiáng)震作用下按照規(guī)范設(shè)計值偏于危險. 因此采用本文修正的地震作用下貯料側(cè)壓力簡化計算公式進(jìn)行筒倉結(jié)構(gòu)設(shè)計，能夠保證筒倉結(jié)構(gòu)抗震設(shè)計的安全性.

6? ?結(jié)? ?論

本文首先對筒倉結(jié)構(gòu)在地震作用下的貯料側(cè)壓力計算理論進(jìn)行分析，在考慮工程實際后，對施衛(wèi)星提出的計算方法進(jìn)行了相關(guān)影響參數(shù)的修正研究，得到以下結(jié)論：

1）考慮貯料高度、地震動強(qiáng)度、貯料內(nèi)摩擦角和貯料與倉壁之間的摩擦因數(shù)對貯料側(cè)壓力的影響，進(jìn)行了相關(guān)系數(shù)的修正，給出了地震作用下貯料側(cè)壓力的簡化計算公式.

2）根據(jù)本文提出的地震作用下貯料側(cè)壓力的簡化公式計算得到的模型結(jié)構(gòu)貯料側(cè)壓力與試驗數(shù)據(jù)進(jìn)行對比，計算值比試驗數(shù)據(jù)稍微偏大，本文推薦公式能夠真實合理地反映地震作用下貯料側(cè)壓力分布，并且也保證了根據(jù)該公式設(shè)計時的筒倉結(jié)構(gòu)的安全性.

3）根據(jù)本文提出的地震作用下貯料側(cè)壓力的簡化計算公式計算得到的模型結(jié)構(gòu)動態(tài)側(cè)壓力與我國規(guī)范對筒倉貯料壓力設(shè)計值進(jìn)行對比，在7度抗震設(shè)防地區(qū)按照規(guī)范設(shè)計時，筒倉下部偏于保守，中上部偏于危險;在考慮強(qiáng)震作用時按照規(guī)范設(shè)計筒倉結(jié)構(gòu)均偏于危險. 采用本文推薦公式計算地震作用下的貯料側(cè)壓力更能保證筒倉結(jié)構(gòu)設(shè)計的安全.

參考文獻(xiàn)

[1]? ? 施衛(wèi)星，朱伯龍. 鋼筋混凝土圓形筒倉地震反應(yīng)試驗研究[J]. 特種結(jié)構(gòu)，1994，11（4）：55—58.

SHI W X，ZHU B L. The experimental study of earthquake response of RC cylinder silos [J]. Special Structures，1994，11（4）：55—58. （In Chinese）

[2]? ? 施衛(wèi)星，朱伯龍. 地震動下煤倉貯料側(cè)壓力研究[J]. 特種結(jié)構(gòu)，1993，10（1）：16—18.

SHI W X，ZHU B L. Study on lateral pressure of coal silo under earthquake [J]. Special Structures，1993，10（1）：16—18. （In Chinese）

[3]? ? HOLLER S，MESKOURIS K. Granular material silos under dynamic excitation：numerical simulation and experimental validation [J]. Journal of Structural Engineering，2006，132（10）：1573—1579.

[4]? ? 張翀，舒贛平. 落地式鋼筒倉卸料的模型試驗研究[J]. 東南大學(xué)學(xué)報（自然科學(xué)版），2009，39（3）：531—535.

ZHANG C，SHU G P. Experimental research of ground-supported steel silo models discharge [J]. Journal of Southeast University （Natural Science Edition），2009，39（3）：531—535. （In Chinese）

[5]? ? 張翀，舒贛平. 顆粒形狀對顆粒流模擬雙軸壓縮試驗的影響研究[J]. 巖土工程學(xué)報，2009，31（8）：1281—1286.

ZHANG C，SHU G P. Effect of particle shape on biaxial tests simulated by particle flow code [J]. Chinese Journal of Geotechnical Engineering，2009，31（8）：1281—1286. （In Chinese）

[6]? ? 趙松. 筒倉貯料壓力分析及其應(yīng)用[D]. 武漢：武漢理工大學(xué)土木工程與建筑學(xué)院，2013：19—23.

ZHAO S. Analysis and application for silo storage pressure [D]. Wuhan：College of Civil Engineering and Architecture，Wuhan University of Technology，2013：19—23. （In Chinese）

[7]? ? 周長東，郭坤鵬，孟令凱，等. 鋼筋混凝土筒倉散料的靜力相互作用分析[J]. 同濟(jì)大學(xué)學(xué)報（自然科學(xué)版），2015，43（11）：1656—1661.

ZHOU C D，GUO K P，MENG L K，et al. Static interaction analysis of RC silo structure-granular material [J]. Journal of Tongji University （Natural Science），2015，43（11）：1656—1661. （In Chinese）

[8]? ? 龐照昆，原方，李丙濤，等. 筒倉側(cè)壁與中心卸料的兩種不同卸料方式的動態(tài)超壓研究[J]. 河南科技，2018，33（22）：113—114.

PANG Z K，YUAN F，LI B T，et al. Research on dynamic overpressure of two different unloading ways for side wall and central unloading of silos [J]. Henan Science and Technology，2018，33（22）：113—114. （In Chinese）

[9]? ? 李智峰，彭政，蔣亦民. 糧倉內(nèi)顆粒壓力的測量：Janssen 行為及其偏差[J]. 物理學(xué)報，2014，63（10）：104503-1-8.

LI Z F，PENG Z，JIANG Y M. Measurements of granular pressure in silo：Janssen behaviour and deviation [J]. Acta Physica Sinica，2014，63（10）：104503-1-8. （In Chinese）

[10]? 韓陽，李東橋，陳家豪，等. 筒倉靜態(tài)儲糧的邊界壓力及倉壁摩擦力試驗研究[J]. 農(nóng)業(yè)工程學(xué)報，2018，34（13）：296—302.

HAN Y，LI D Q，CHEN J H，et al. Experimental study on boundary pressure and wall friction under static grain storage in silo [J]. Transactions of the Chinese Society of Agricultural Engineering，2018，34（13）：296—302. （In Chinese）

[11]? 程奇鵬，孫巍巍，盧賽. 基于PFC3D的淺圓倉偏心卸料離散元研究[J]. 土木工程與管理學(xué)報，2016，33（6）：43—47.

CHENG Q P，SUN W W，LU S. Discrete element analysis of squat silo under eccentric discharge by PFC3D [J]. Journal of Civil Engineering and Management，2016，33（6）：43—47. （In Chinese）

[12]? 王培軍，朱立，王慶華，等. 庫側(cè)卸料筒倉內(nèi)物料流動及側(cè)壓力分布實驗研究[J]. 中國粉體技術(shù)，2016，22（2）：1—6.

WANG P J，ZHU L，WANG Q H，et al. Experimental studies on flow pattern and wall pressure distribution in wall discharged silo [J]. China Powder and Technology，2016，22（2）：1—6. （In Chinese）

[13]? 梁斌昌. 基于離散元的筒倉卸料模擬研究[D]. 徐州：中國礦業(yè)大學(xué)力學(xué)與土木工程學(xué)院，2018：43—70.

LIANG B C. Study on silo discharges simulation base on DEM [D]. Xuzhou：School of Mechanics and Civil Engineering，China University of Mining and Technology，2018：43—70. （In Chinese）

[14]? 張逯見. 筒承式立筒群倉結(jié)構(gòu)模型模擬地震振動臺試驗研究[D]. 鄭州：河南工業(yè)大學(xué)土木建筑學(xué)院，2010：14—72.

ZHANG L J. Study of shaking table tests on the model of cylindrical-supporting group silos structures [D]. Zhengzhou：College of Civil Engineering and Architecture，Henan University of Technology，2010：14—72. （In Chinese）

[15]? 郭坤鵬. 筒倉貯料壓力分布機(jī)制及地震易損性研究[D]. 北京：北京交通大學(xué)土木建筑工程學(xué)院，2016：19—32.

GUO K P. Lateral pressure distribution and seismic vulnerability assessment of circular reinforced concrete silo[D]. Beijing：School of Civil Engineering，Beijing Jiaotong University，2016：19—32. （In Chinese）

[16]? 許啟鏗，王錄民，丁永剛. 筒承式糧食立筒倉結(jié)構(gòu)模態(tài)及地震反應(yīng)分析[J]. 世界地震工程，2015，31（3）：55—61.

XU Q K，WANG L M，DING Y G. Modal and earthquake response analysis of cylinder-supported grain silo structures [J]. World Earthquake Engineering，2015，31（3）：55—61. （In Chinese）

[17]? GB 50011—2010建筑抗震設(shè)計規(guī)范[S]. 北京：中國計劃出版社，2010：19—20.

GB 50011—2010 Code for seismic design of buildings [S]. Beijing：China Planning Press，2010：19—20. （In Chinese）

[18]? 王錄民，張昊，張逯見，等. 地震作用下立筒群倉貯料側(cè)壓力試驗研究[J]. 四川建筑科學(xué)研究，2012，38（2）：175—178.

WANG L M，ZHANG H，ZHANG L J，et al. Experimental study on lateral pressures acting on the wall of a group of granular silos induced by earthquakes [J]. Sichuan Building Science，2012，38（2）：175—178. （In Chinese）

[19]? GB 50077—2017 鋼筋混凝土筒倉設(shè)計標(biāo)準(zhǔn)[S]. 北京：中國計劃出版社，2017：7—11.

GB 50077—2017 Standard for design of reinforced concrete silos [J]. Beijing：China Planning Press，2017：7—11. （In Chinese）