地震作用下筒仓结构贮料侧压力计算方法

周长东 张泳 邱意坤 梁立灿 阿斯哈

摘   要：为了对筒仓结构进行合理的抗震设计，得到实际工程中实用的贮料侧压力，引入“施卫星”地震作用下贮料侧压力计算方法，以工程设计中能够获得的材料参数为基础，确定贮料侧压力计算中所需要的修正参数;然后利用ABAQUS建立筒仓结构数值模型进行计算分析，根据获得的贮料侧压力变化曲线拟合得出相应参数的修正值;综合不同参数，归纳出考虑筒仓-贮料相互作用的贮料侧压力修正公式;最后，将贮料侧压力修正公式所得的计算值与既有试验数据及规范设计值进行对比分析. 结果表明，采用本文公式得到的计算值与试验值吻合较好，能够真实合理地反映实际地震作用下贮料侧压力分布;考虑强震作用下贮料侧压力按照规范设计值设计偏于危险，而按照本文计算公式设计更加经济可靠.

关键词：钢筋混凝土筒仓;地震作用;贮料侧压力;计算方法;抗震设计

中图分类号：TU279.7                              文献标志码：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

在鋼筋混凝土筒仓的抗震设计中，通常将筒仓当作悬臂结构来计算整个结构根部的内力，并采用截面强度设计，忽略了贮料侧压力对筒仓壁的影响. 实际上，该侧压力变化将在仓壁的圆周方向上引起弯曲应力，这对于筒仓这类薄壁结构是不利的.

近年来国内外学者对于筒仓贮料侧压力进行了广泛的研究. 施卫星等[1-2]对2个1 ∶ 10钢筋混凝土圆筒仓结构模型进行了振动台试验，研究了地震作用下筒仓的动力响应和破坏形式，根据试验结果提出了储煤对仓壁侧压力计算方法; Holler等[3]对动态激发下的存储颗粒材料的筒仓进行了有限元数值模拟，并将采用该模型获得的数值模拟结果与在法国Saclay的振动台试验结果进行了比较分析;张翀等[4-5]分析了各种卸料工况下仓壁应变的变化规律;赵松[6]进行了钢筒仓侧压力模型试验，得到了筒仓结构静止及卸料状态下的贮料侧压力分布曲线;周长东等[7]针对筒仓静力作用提出了一种散粒体的亚塑性本构模型，分析了贮料对仓壁的静止侧压力;庞照昆等[8]通过模型试验方法，研究了筒仓在侧壁卸料与中心卸料的动态超压现象.

目前的研究大多集中于静力状态以及卸料状态下贮料侧压力的变化[9-13]，主要通过理论分析法、有限单元法以及离散元法进行分析，而少有对于地震作用下贮料侧压力分布的研究. 考虑地震作用时，离散元法虽然可以较好地模拟颗粒之间的相互作用，但输入地震波的困难让该方法难以应用;有限元法数值计算的精度及稳定性高，但是由于其对特定结构进行数值分析，并不能得出一个普遍的贮料侧压力计算方法;理论分析法在地震作用下受多种因素的限制，并不能单独通过公式推导得出所需的计算参数[14].

本文首先对地震作用下考虑筒仓-贮料相互作用的贮料侧压力计算方法进行理论分析，确定贮料侧压力计算中所需要的修正参数;然后利用ABAQUS建立筒仓结构数值模型进行计算，获取与影响参数相应的贮料侧压力变化曲线;随后对贮料侧压力变化曲线进行拟合，得到相应参数在贮料侧压力计算中的修正值;之后综合不同参数拟合结果，归纳出考虑筒仓-贮料相互作用的贮料侧压力修正公式;最后，将根据贮料侧压力修正公式所得的计算值与既有试验数据及规范设计值进行对比分析，验证计算公式的合理性.

1   地震作用下贮料侧压力計算

目前针对地震作用下贮料侧压力计算方法的研究中，仅有施卫星等通过试验与理论相结合的形式，推导出符合相应试验模型的地震作用下贮料侧压力计算方法[2]. 因此，本文结合数值模拟的方法，在“施卫星”计算方法的基础上进行修正，得到可在实际工程应用的地震作用下贮料侧压力计算方法.

1.1   “施卫星”计算方法

施卫星等[2]提出以计算筒仓贮料质心处的煤侧压力为主，再通过各种修正系数去计算贮料上层及底层的侧压力的计算方法.

当地震作用于筒仓时，在其贮料质心处的绝对加速度为Sa，对称压力Pce不可能使贮料产生加速度，而偏心压力Pee（θ）使贮料产生了Sa加速度，假设 Pee（θ）的分布满足：

考虑到贮料质心加速度是贮料侧压力及贮料层间摩擦力共同作用的结果，因此使得计算所得侧压力与实际值有所差别，施卫星对该公式进行了修正：

式中：Ce由试验结果确定，与测点位置、台面加速度、贮料内摩擦角和湿度等因素有关.

为求解出工程设计时考虑的最大贮料侧压力Pmax，需先计算最小贮料侧压力Pmin：

式中：h为贮料高度;Cp为超压系数;Ph为Janssen静压力.

综上，考虑筒仓-贮料相互作用的贮料侧压力计算方法需要确定6个参数：筒仓半径、贮料质心处加速度、计算截面位置、贮料容重、内摩擦角、贮料与仓壁摩擦因数.

1.2   简化计算方法

应用公式（3）进行地震作用下贮料侧压力计算时，需要输入贮料质心处加速度. 而工程设计时，无法直接得到贮料质心加速度，可通过综合考虑测点位置、贮料内摩擦角及摩擦因数等因素，对台面加速度进行相应修正后[2]，得到：

式中：Cg为考虑地震强度的修正系数;Ch为考虑贮料高度的修正系数;C?为考虑贮料内摩擦角的修正系数;Cμ为考虑摩擦因数的修正系数.

在求解出地震作用下贮料侧压力变化值Pe后，需要确定地震作用下最小贮料侧压力Pmin来确定最大贮料侧压力Pmax. 根据文献[15]中提供的筒仓在静力状态以及地震作用下的贮料侧压力数值模拟分析结果可看出，筒仓结构在地震作用下的最小贮料侧压力Pmin < 静力状态下贮料侧压力，与施卫星定义下的Pmin计算公式相悖. 在工程设计时通常将最大贮料侧压力Pmax作为主要设计影响因素，因此对地震作用下最大贮料侧压力Pmax进行简化计算，综合考虑式（3）与式（7）对最大贮料侧压力Pmax进行修正，即：

式中：Cg，max为考虑地震强度的最大贮料侧压力修正系数;Ch，max为考虑贮料高度的最大贮料侧压力修正系数; C?，max为考虑贮料内摩擦角的最大贮料侧压力修正系数;Cμ，max为考虑摩擦因数的最大贮料侧压力修正系数.

筒仓结构设计时，根据结构所在地区的环境参数和设定的材料参数，建立有限元分析模型，选取合理的地震波进行筒仓结构的地震响应分析，获取公式（9）中的修正系数，即可计算得到设计所需的最大贮料侧压力.

2   筒仓结构有限元分析模型

2.1   工程概况

以浙江某粮库的筒仓结构为工程背景，将钢筋混凝土筒仓群中的单一筒仓作为设计模拟的原型（如图1所示）. 建筑物总高度35.0 m，粮食装载高度27.0 m，筒仓单仓内径12.0 m，筒仓壁厚0.22 m;环梁尺寸为0.4 m × 0.8 m，门洞尺寸为1.8 m × 2.7 m，窗洞尺寸为0.9 m × 1.5 m. 抗震设防烈度为7度，设计地震组为第一组. 抗震设防类别为乙类，框架抗震设防等级为二级，场地类别为Ⅲ类[16].

2.2   ABAQUS有限元模型

采用ABAQUS软件对所选取的单仓结构建立分析模型（见图1）. 选用混凝土损伤塑性本构模拟仓壁的混凝土材料（见表1），混凝土的单轴本构曲线采用《混凝土结构设计规范》（GB 50010—2010）附表C.2中的推荐公式;选用理想弹塑性本构模拟钢筋（见表2）;选用摩尔库伦本构模拟贮料（见表3）.

利用ABAQUS中S4单元模拟筒仓仓壁，利用C3D8R单元模拟贮料;利用ABAQUS自带的“接触对”来模拟仓壁与贮料的接触关系，其中选取仓壁内表面为主接触面，两者之间的摩擦因数为0.4.

2.3   模型验证

由模型试验[14]知，满仓状态下，在地震烈度为7度多遇和7度基本时，上海人工波SHW2对单仓y方向影响较大，因此在筒仓数值模型底部施加SHW2波进行动力时程分析. 为验证单仓模型的合理性，提取单仓模型y方向单元最大加速度以及最大相对位移、自振频率值，与相应的筒仓结构模型试验结果进行对比，对比结果详见图2和图3.

通过有限元模拟分析，得到满仓状态下数值模型的自振频率为2.617 Hz. 振动台试验获取单仓一阶频率后，通过相似系数返回原型结构得到的原型结构自振频率为2.631 Hz. 如图2和图3所示，将有限元模拟提取的最大加速度和最大相对位移值与单仓试验数据进行比较，偏差均在5%之内. 因此，本文建立的单仓数值模型与振动台试验模型吻合较好，证明了数值模型的合理性.

3   最大贮料侧压力计算参数分析

采用增量动力分析（Incremental Dynamic Analysis，IDA）方法，针对不同地震波的峰值地面加速度（Peak Ground Acceleration，PGA）进行调幅并在数值模型底部输入，得到时程分析结果. 为获取式（9）中的相应参数，需要拟合不同工况下贮料侧压力求解修正参数. 在进行数值模拟计算时，单一模型输入的参数包括了地震动强度、贮料高度、内摩擦角以及摩擦因数，因此获得的贮料侧压力分布曲线不能直接通过拟合得到单一的修正参数. 本文假定当PGA = 0.1g、? = 28°、μ = 0.40条件下的筒仓贮料侧压力修正系数为1，即Cg = 0.1，max = 1、 C? = 28°，max = 1、Cμ = 0.4，max = 1，在此条件下，通过变换不同参数进行数值模拟，将所得贮料侧压力分布曲线进行拟合.

3.1   地震波的选取

根据原型筒仓的工程概况，按《建筑抗震设计规范》[17]的要求选取两条天然波和一条人工波进行时程分析，分别选用El-Centro波、NGA#178号地震波以及上海人工地震波SHW2波输入数值模型底部进行研究. 图4所示为所选地震波的平均反应谱，与建筑抗震设计规范中的设计反应谱吻合较好.

图5、图6、图7分别为考虑地震动强度、内摩擦角以及摩擦因数影响下的最大贮料侧压力分布曲线，可发现贮料高度在不同参数中对最大贮料侧压力的影响最为明显. 因此，首先考虑以PGA = 0.1g时筒仓最大贮料侧压力曲线为基准，对考虑贮料高度的最大贮料侧压力参数Ch，max進行拟合;随后对不同参数影响下最大贮料侧压力分布曲线进行拟合，确定不同参数影响下的最大贮料侧压力修正系数.

3.2   考虑贮料高度的最大贮料侧压力参数

采用PGA = 0.1g、贮料参数为? = 28°、 μ = 0.40的筒仓模型进行数值模拟计算，选用此参数可排除地震动强度、内摩擦角及贮料与仓壁的摩擦因数的影响. 为求解考虑高度的最大贮料侧压力参数，将式（8）、式（9）进行变形，得：

式中：Pd，max是仓壁受到的最大动压力.

图8所示为筒仓结构在PGA = 0.1g地震作用下沿高度方向修正系数分布曲线. 可以看出，筒仓沿高度方向的修正系数曲线可分为四段，分别为漏斗处、筒仓仓壁相对高度0 ～ 0.25、相对高度0.25 ～ 0.75、相对高度0.75 ～ 1.0. 选取筒仓仓壁相对高度为0、0.25、0.5、0.75、1处的修正系数进行分段拟合，如图9所示.

3.3   考虑地震动强度的最大贮料侧压力参数

考虑贮料高度对筒仓侧压力计算的影响后， 采用贮料参数为? = 28°、 μ = 0.40的筒仓模型进行数值模拟计算，选用此参数可以排除内摩擦角、贮料与仓壁的摩擦因数的影响. 为求解考虑地震动强度的最大贮料侧压力参数，将式（8）（9）进行变形，得：

图10所示是筒仓结构在不同地震动强度下的 Cg，max分布曲线. 由图10可知，筒仓结构在不同高度处Cg，max的变化规律不同. 图11分别对位于筒仓漏斗处、筒仓仓壁相对高度为0、0.25、0.5、0.72、1处不同地震动强度的Cg，max进行拟合，可以发现：在相对高度>0的筒仓仓壁部位，PGA = 0.1g ～ 0.5g时筒仓地震强度修正系数呈线性分布，PGA = 0.6g ～ 1.0g时筒仓地震强度修正系数的分布形式为二次函数;在筒仓结构漏斗处，PGA = 0.1g ～ 0.4g时筒仓结构地震强度修正系数呈线性分布，PGA = 0.5g ～ 1.0g时修正系数分布形式为斜率不同的直线. 因此，分段讨论分布函数公式.

3.4   考虑内摩擦角的最大贮料侧压力参数

为研究贮料之间的不同摩擦因数对贮料侧压力的影响，考虑选用? = 25°、28°、30°、33°、35°五种内摩擦角进行计算. 本节采用贮料参数为PGA = 0.1g、 μ = 0.40的筒仓模型进行数值模拟计算，选用此参数可以排除地震动强度、贮料与仓壁的摩擦因数的影响. 为求解考虑内摩擦角的最大贮料侧压力参数，将式（8）、式（9）进行变形，得：

图12是筒仓结构在PGA = 0.1g的地震作用下不同摩擦角修正系数分布曲线. 由图12可知，随着贮料内摩擦角的增大，考虑内摩擦角的最大贮料侧压力参数逐渐减小;当位于筒仓结构漏斗处时，随着内摩擦角的变化，考虑内摩擦角的最大贮料侧压力参数变化量基本相同;当位于筒仓结构中上部时，相同高度处随着内摩擦角的变化，考虑内摩擦角的最大贮料侧压力参数变化量有较大差异. 因此，对该参数进行分段考虑，分别选取不同内摩擦角的漏斗处修正系数及筒仓仓壁相对高度为0、0.25、0.5、0.72、1处修正系数平均值进行拟合，如图13所示.

3.5   考虑摩擦因数的最大贮料侧压力参数

为研究贮料与仓壁之间不同摩擦因数对贮料侧压力的影响，考虑采用μ = 0.30、0.35、0.40、0.45、0.50五种摩擦因数进行计算. 采用贮料参数为PGA = 0.1 g、? = 28°的筒仓模型进行数值模拟计算，选用此参数可以排除地震动强度、贮料内摩擦角的影响. 为求解考虑摩擦因数的最大贮料侧压力参数，将式（8）、式（9）进行变形，得：

图14所示是筒仓结构在PGA = 0.1g地震作用下不同摩擦因数的修正系数分布曲线. 由图14可知，随着筒仓仓壁与贮料摩擦因数的增大，Cμ，max逐渐减小;当位于筒仓结构漏斗处时，随着摩擦因数的变化，Cμ，max的变化量基本相同;当位于筒仓结构中下部时，摩擦因数引起的筒仓结构Cμ，max的变化量随着高度的增加基本不变;当位于筒仓结构上部时，摩擦因数引起的筒仓结构Cμ，max的变化量随着高度的增加而减小. 因此，考虑将高度与摩擦因数的影响结合，进行分段考虑.

4   计算值与试验数据对比

王录民等[18]通过模拟地震振动台试验对钢筋混凝土单仓及群仓模型进行研究，观察了贮料在地震作用下的侧压力变化规律. 根据第3节中得到的地震作用下修正的贮料侧压力计算公式，结合试验中的数据进行对比，对比分析结果如图17所示.

如图17所示，7度罕遇和8度罕遇地震作用下，根据简化的最大贮料侧压力计算公式得到的模型结构贮料侧压力计算值与试验数据吻合较好;7度基本地震作用下筒仓结构中部计算值与试验值有所差异，其原因在于模型试验时，在地震强度较小时，中部贮料由于散粒体的固结作用未能充分运动，因此试验值相对于计算值偏小;在7度基本、7度罕遇、8度罕遇地震作用下贮料底部及中上部计算值相对于试验数据稍微偏大，因此使用简化的最大贮料侧压力计算公式能合理反映实际地震作用下侧压力分布，保证了采用该公式设计时结构的安全性.

5   计算值与规范设计值对比

目前，比较经典的筒仓压力分布理论主要有Janssen理论、Rankine理论、Airy理论以及修正的Coulomb理论等. 我国在进行深仓贮料压力计算[19]时，以Janssen公式为理论基础，同时在设计时考虑相应的放大系数：

式中：Ph为高度h处仓壁受到的法向侧压力;γ为贮料的重力密度，kN/m3;μ为筒仓仓壁与贮料的摩擦因数;ρ为筒仓结构的水力半径，m;k为主动侧压力系数;Ch為深仓贮料水平压力修正系数.

在考虑筒仓侧压力静力计算时，取Ch = 1进行计算;在进行筒仓结构设计时，需进行深仓贮料水平侧压力修正，见表4[19].

修正系数Ch 备注（参数） 仓壁 1. 当筒仓贮料高度hn与筒仓内径dn的比大于3.0时，Ch应乘以1.1;

2. 流动性较差的贮料，Ch宜乘以0.9;

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

将经过第3节相关系数修正后得到的地震作用下贮料侧压力简化计算公式所得的计算结果与我国规范对筒仓贮料压力设计值进行对比，如图18所示. 可以发现，贮料压力设计值随高度增大，侧压力逐渐减小. 当高度h小于1 m时，贮料侧压力简化计算值随高度增大而减小;当1 m ≤ h ≤ 1.5 m时，计算值随高度增大而增大;当h > 1.5 m时，计算值随高度增大而减小. 这是由于规范设计值是通过放大系数对静力作用下贮料侧压力进行修正;而地震作用下，由于筒仓结构中下部填充密实，贮料与筒仓仓壁相互作用较小;筒仓结构中上部贮料则由于填充密实度较差，因此动力作用下与筒仓仓壁相互作用较大，造成筒仓结构中上部出现贮料侧压力增大的情况. 筒仓结构顶部则由于贮料较少，因此在h > 1.5 m时，侧压力随高度增大而减小.

如图18所示，在7度基本烈度地震作用下，筒仓结构中下部动态侧压力均小于规范设计值，可见在7度抗震设防地区按照规范设计值对仓壁进行承载力计算偏于保守;而筒仓中上部动态侧压力计算值均大于规范设计值，按照规范设计筒仓中上部偏于危险;在7度罕遇烈度及8度罕遇烈度下，筒仓动态侧压力均大于规范设计值，可见在考虑强震作用下按照规范设计值偏于危险. 因此采用本文修正的地震作用下贮料侧压力简化计算公式进行筒仓结构设计，能够保证筒仓结构抗震设计的安全性.

6   结   论

本文首先对筒仓结构在地震作用下的贮料侧压力计算理论进行分析，在考虑工程实际后，对施卫星提出的计算方法进行了相关影响参数的修正研究，得到以下结论：

1）考虑贮料高度、地震动强度、贮料内摩擦角和贮料与仓壁之间的摩擦因数对贮料侧压力的影响，进行了相关系数的修正，给出了地震作用下贮料侧压力的简化计算公式.

2）根据本文提出的地震作用下贮料侧压力的简化公式计算得到的模型结构贮料侧压力与试验数据进行对比，计算值比试验数据稍微偏大，本文推荐公式能够真实合理地反映地震作用下贮料侧压力分布，并且也保证了根据该公式设计时的筒仓结构的安全性.

3）根据本文提出的地震作用下贮料侧压力的简化计算公式计算得到的模型结构动态侧压力与我国规范对筒仓贮料压力设计值进行对比，在7度抗震设防地区按照规范设计时，筒仓下部偏于保守，中上部偏于危险;在考虑强震作用时按照规范设计筒仓结构均偏于危险. 采用本文推荐公式计算地震作用下的贮料侧压力更能保证筒仓结构设计的安全.

参考文献

[1]    施卫星，朱伯龙. 钢筋混凝土圆形筒仓地震反应试验研究[J]. 特种结构，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]    施卫星，朱伯龙. 地震动下煤仓贮料侧压力研究[J]. 特种结构，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]. 东南大学学报（自然科学版），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]. 岩土工程学报，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]    赵松. 筒仓贮料压力分析及其应用[D]. 武汉：武汉理工大学土木工程与建筑学院，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]. 同济大学学报（自然科学版），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]    庞照昆，原方，李丙涛，等. 筒仓侧壁与中心卸料的两种不同卸料方式的动态超压研究[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]    李智峰，彭政，蒋亦民. 粮仓内颗粒压力的测量：Janssen 行为及其偏差[J]. 物理学报，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]  韩阳，李东桥，陈家豪，等. 筒仓静态储粮的边界压力及仓壁摩擦力试验研究[J]. 农业工程学报，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]. 土木工程与管理学报，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]  王培军，朱立，王庆华，等. 库侧卸料筒仓内物料流动及侧压力分布实验研究[J]. 中国粉体技术，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]. 徐州：中国矿业大学力学与土木工程学院，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]  張逯见. 筒承式立筒群仓结构模型模拟地震振动台试验研究[D]. 郑州：河南工业大学土木建筑学院，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]  郭坤鹏. 筒仓贮料压力分布机制及地震易损性研究[D]. 北京：北京交通大学土木建筑工程学院，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]  许启铿，王录民，丁永刚. 筒承式粮食立筒仓结构模态及地震反应分析[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建筑抗震设计规范[S]. 北京：中国计划出版社，2010：19—20.

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

[18]  王录民，张昊，张逯见，等. 地震作用下立筒群仓贮料侧压力试验研究[J]. 四川建筑科学研究，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 钢筋混凝土筒仓设计标准[S]. 北京：中国计划出版社，2017：7—11.

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