㊀第25卷第2期洁净煤技术Vol．25㊀No．2㊀㊀2019年3月Clean Coal TechnologyMar.㊀2019㊀下行移动床颗粒流动特性冷模试验研究杨㊀帆,郑㊀雄,李建波,王泉海,薛景文,陆㊀凯,卢啸风(重庆大学低品位能源利用技术及系统教育部重点实验室,重庆㊀400044)摘㊀要:针对目前大宗固体颗粒余热回收常存在换热不充分㊁气固流动阻力大㊁回收得到的余热质量较低等问题,提出了一种气固交叉流动移动床高温颗粒冷却技术方案,并在自行设计的工业级试验装置上利用CFB 锅炉炉渣完成了冷态条件下的颗粒流动特性试验研究㊂通过取样与高速摄影仪拍摄相结合的方式对颗粒流动特性进行测量,分析讨论了颗粒流通截面尺寸对颗粒流动特性的影响和错流段空截面风速对颗粒流动稳定性的影响㊂结果表明:直流段截面中部颗粒下行速度基本一致,在忽略边壁影响的条件下,计算得到直流段正面和侧面流动指数M F1均大于0．3,认为直流段颗粒流动为整体流状态;通过对不同位置颗粒取样,发现各取样点颗粒的质量分数最大相差约14．4%,在实际工程中可认为扩大段颗粒整体流动均匀,但不同位置的颗粒粒径分布存在差异;随着横向风速的增大,颗粒携带速率增大,流过下行移动床的水平风速不能超过1．45m /s ㊂关键词:移动床;颗粒流;流动特性;冷模试验中图分类号:TQ051．5㊀㊀㊀文献标志码:A㊀㊀㊀文章编号:1006-6772(2019)02-0089-07移动阅读收稿日期:2019-01-02;责任编辑:白娅娜㊀㊀DOI :10.13226/j.issn.1006-6772.19010201基金项目:国家重点研发计划资助项目(2016YFB0600201)作者简介:杨㊀帆(1990 ),男,湖北恩施人,硕士研究生,研究方向为颗粒余热回收㊂E -mail :20120734@ ㊂通讯作者:卢啸风,男,教授,主要从事循环流化床锅炉㊁洁净煤燃烧技术相关研究㊂E -mail :xfluke@引用格式:杨帆,郑雄,李建波,等.下行移动床颗粒流动特性冷模试验研究[J].洁净煤技术,2019,25(2):89-95.YANG Fan,ZHENG Xiong,LI Jianbo,et al.Experimental study on the particle flow characteristics in a downward moving bed[J].Clean Coal Technology,2019,25(2):89-95.Experimental study on the particle flow characteristics in adownward moving bedYANG Fan,ZHENG Xiong,LI Jianbo,WANG Quanhai,XUE Jingwen,LU Kai,LU Xiaofeng(Key Laboratory of Low -Grade Energy Utilization Technology &System ,Chongqing University ,Chongqing ㊀400044,China )Abstract :In terms of the problems that exsit in the waste heat recovery of bulk solid particles,such as insufficient heat transfer,high gas -solid flow resistance and low quality of the waste heat recovery,a cooling technology for high temperature particle using gas -solid cross -flow moving bed was proposed in this study.The experimental study on particle flow characteristics under cold operating conditions was car-ried out with CFB boiler ash via a self -designed industrial experimental apparatus.Particle flow characteristics were measured by sampling and high -speed photography.The effects of the cross -sectional size of the particle flow on particle flow characteristics and cross -sectional gas speed on particle flow stability were analyzed and discussed.The results show that the downward velocity of particles in the central part of the straight flow section is basically constant.Neglecting the influence of the side walls,the calculated flow index M F1of the front and side section of straight flow section exceed 0.3,which show that the particle flow was considered as a bulk flow under the experimental conditions.By sampling the particles at different locations,it is found that the maximum difference of particle mass ratio is about 14.4%,which could be considered that the particle flow in the enlarged section was uniform in practical engineering,while the particle size distri-bution at different locations is different.With increasing the transverse air velocity,the particle carrying rate increases.The air velocity in the downward moving bed can not exceed 1.45m /s.Key words :moving bed;granular flow;flow characteristic;experimental study0㊀引㊀㊀言大宗固体颗粒余热回收问题广泛存在于建材[1-2]㊁电力[3-4]㊁冶金[5-6]㊁化工[7-8]等行业,如干法粒化后的高炉渣颗粒温度在1000ħ以上[5],陶粒生产企业生产的陶粒温度可达800ħ[9],大型循环982019年第2期洁净煤技术第25卷流化床(CFB)锅炉从炉膛排出的灰渣温度在850~ 900ħ[10]㊂为提高能源利用效率,需对这些高温颗粒物料进行冷却,同时进行余热回收㊂高温固体颗粒的冷却方式主要采用间接式冷却和接触式冷却[11]㊂间接式冷却通常以水作为冷却剂,受水温升高限制以及换热效率与换热管磨损的矛盾,很难置换出高品质(高温)余热,而低品质的余热难以有效利用;接触式冷却通常采用气体与高温固体颗粒直接接触冷却㊂具体的气固接触方式[11]有移动床换热(如水泥生产工艺中的篦冷机[1])㊁流化床换热(如循环流化床锅炉的流化床冷渣器[3-4])以及气力输送冷却㊂选择冷却方式时需考虑高温颗粒的物理特性(粒径㊁比重等)以及要求的冷却参数(高温固体颗粒冷却速率㊁颗粒或冷却介质的最终温度以及冷却介质耗量等)㊂实践表明,采用气体(空气㊁氮气或烟气)冷却高温固体颗粒,既可获得较好的高温颗粒冷却效果,也可获得较高品质的余热㊂3种气固换热技术[1-4]广泛用于各行业的高温固体颗粒冷却与余热回收过程,但都很难获得与被冷却物料的最高温度接近的换热介质终温㊂学者针对这些高温固体颗粒冷却设备提出了改进方法㊂如为提高篦冷机的热回收效率,邵卫[1]采取优化篦冷机配风方式降低风机功耗,李心宁[12]提出在篦冷机尾部增加一个余热取风口降低直接排出的低温风风温,提高了篦冷机效率,但由于篦冷机采用分段冷却,篦冷机优化后直接排出的低温风和400~500ħ中温风仍存在㊂为提高流化床冷渣器的余热回收效率,甘露[3]提出的双喷动床式冷渣器沿用了复合式流化床冷渣器[4]处理底渣时采用的粗细颗粒分选分离并分级冷却的思路,通过合理配风能较好控制冷渣器的风冷比例,但换热后的热风温度仅为300~400ħ,为底渣温度(850ħ左右)的一半㊂以上多从优化传统换热器结构出发提高余热回收效率,受传统换热器的结构原理限制,优化后空气耗量大㊁空气温升小㊁余热难以全部利用等问题仍存在㊂为解决传统高温固体颗粒气固直接接触换热时的颗粒流动问题,满足大宗固体颗粒余热回收要求换热量大㊁换热充分㊁气固流动阻力小㊁回收的余热质量较高等要求,本文提出了一种气固交叉流动移动床高温颗粒冷却技术方案,并在自行设计的工业级试验装置上完成了冷态条件下的颗粒流动特性试验研究㊂1㊀试㊀㊀验1．1㊀试验装置及试验物料水平气流作用下移动床冷态试验台如图1所示,试验台本体由钢板和角钢焊接而成,图1中A㊁B㊁C处安装有机玻璃板,便于观察试验台内的物料流动情况㊂颗粒由截面较小的入口段(300mmˑ150mm)经扩大段后进入颗粒截面较大的直流段(300mmˑ500mm),试验台内气体2次水平穿过颗粒层,错流段采用百叶窗式挡板将气固通道分开㊂本文宽度和高度方向分别为x㊁z,深度方向y垂直于xz平面向外㊂图1㊀水平气流作用下移动床冷态试验台Fig．1㊀Moving bed cold experimental apparatus underhorizontal air flow试验时,颗粒从料斗加入,并沿颗粒通道向下流动㊂颗粒的下行流动速率通过调节试验台底部星型卸料器的转速来控制㊂风速由热线风速仪在线测量,本文中风速v均表示错流段空床截面风速㊂为接近实际,试验所用物料为某电厂CFB锅炉炉渣,主要物理性质见表1,粒径分布如图2所示㊂表1㊀锅炉灰渣物理性质Table1㊀Physical properties of boiler ash真密度/(kg㊃m-3)散密度/(kg㊃m-3)平均粒径/mm最大粒径d max/mm自然堆积空隙率ε堆积角φ/(ʎ) 2283．61146．05．84300．530．5 1．2㊀试验内容1．2．1㊀颗粒流通截面尺寸对颗粒流动特性影响为保证气固换热均匀,颗粒在直流段中的流动状态应呈现整体流流动状态㊂受颗粒通道入口段及直流段截面尺寸差异的影响,需研究颗粒通道扩大段内颗粒的流动规律㊂1．2．2㊀错流段空截面风速对颗粒流动稳定性影响气体以不同气流速率流过颗粒料层时,会影响09杨㊀帆等:下行移动床颗粒流动特性冷模试验研究2019年第2期图2㊀床料筛分特性曲线Fig．2㊀Screening characteristic curves of bed material 颗粒整体流动和局部流动特性,可能导致换热不均匀,因此研究不同风速下颗粒携带速率v b(错流段单位面积单位时间内吹入风室的颗粒质量)及其粒径分布,从而判断风速对床层稳定性的影响㊂1．3㊀测量及分析方法1．3．1㊀颗粒均匀流动特性的测量颗粒在经过颗粒通道扩大段后沿深度方向的粒径分布是移动床整体颗粒流动均匀性的重要参数,本文通过取样与高速摄影仪拍摄相结合的方式,对颗粒流动特性进行测量㊂拍摄点位于图1中A㊁B㊁C处,拍摄后的图像利用Matlab软件处理,得到颗粒下移轨迹和颗粒通道扩大段倾斜面附近的空腔形状及空腔面积A c(图4(b)空腔图片中黑色区域面积),操作步骤如图3所示㊂按照文献[13-14]的数据处理方式,处理后图像如图4所示㊂取样方式如图5所示,试验运行稳定后利用取样盒在取样点M1㊁M2㊁M3㊁M4取一定时间内的下落颗粒,对取样用不同筛分度的标准筛进行筛分并精确称量(电子秤精度为0．1g),得到颗粒整体下移速率v d ㊂图3㊀图像处理的操作步骤Fig．3㊀Operational steps of image processing 1．3．2㊀颗粒流动均匀性分析方法一般颗粒在床内的流动主要有2种形式:整体流和漏斗流[15]㊂整体流是全部颗粒以相对均匀的速率向出口移动,下料速率稳定,下料密度均匀,下料顺序为先进先出;而漏斗流的特点是床内颗粒速率梯度大,中心线处的颗粒速率大,边壁处颗粒速率小,颗粒有先进后出的现象[16]㊂在本装置中,直流图4㊀图像处理后效果Fig．4㊀Result of imageprocessing图5㊀颗粒取样方式Fig．5㊀Particle sampling method段呈现整体流对颗粒均匀冷却有重要意义㊂定义流动指数M F1为壁面附近颗粒的平均速率v w与中心处的颗粒平均速率v c的比值,即M F1=v w/v c(1)㊀㊀根据Johanson和Jenike[17]研究可知,M F1=0．3是整体流向漏斗流转变的分界线,M F1>0．3时颗粒流动为整体流,M F1<0．3时颗粒流动为漏斗流[18-19]㊂本文考虑到粗糙壁面对颗粒下移造成阻塞,取离壁面0．045m处颗粒的平均速率作为壁面附近颗粒的平均速率v w㊂1．3．3㊀颗粒层流动稳定性评价方法由于气固直接接触,当风速增大到一定值后,颗粒会被风携带至风室,颗粒层稳定性被破坏㊂由于没有相关计算式,本文暂用临界流化风速作为破坏颗粒层稳定性风速的参考值㊂利用浙江大学[11]提出的准则式对临界流化风速进行描述,即Re mf=0．0882Ar0．528(2)式中,Ar=2ˑ104~7ˑ104;Re mf为临界流化风速对应雷诺数,计算式为Re mf=u mf d p/υg(3)式中,u mf为临界流化速度;υg为气体的运动黏度;d p 为颗粒平均粒径,计算式为d p=ΦpðX i d i(4)式中,Φp为颗粒的球形度;X i为某一筛上的残留质量分数;d i为某一残留质量所在的筛孔直径与上一筛孔直径的平均值,计算式为192019年第2期洁净煤技术第25卷d i =d i d i +1(5)㊀㊀对式(2)㊁(3)变形后可得到u mf =0．294d 0．584pυ0．056g ρp -ρg ρg æèöø0．528(6)其中,ρp 为颗粒真密度;ρg 为气体密度㊂通过式(6)可得到颗粒层对应的临界流化风速㊂1．4㊀试验工况颗粒整体下移速率v d 分别为4㊁8mm /s,其对应的风速v 为0．65㊁0．83㊁1．05㊁1．25㊁1．43㊁1．70m /s㊂2㊀试验结果与讨论2．1㊀直流段颗粒流动特性v d =4mm /s 时,得到直流段正面和侧面不同时刻颗粒的下移位置如图6所示㊂可知,直流段同一时刻中部颗粒的位置基本在一条水平线上,呈整体流状态;而受边壁的影响,边壁处颗粒下移速率慢,甚至不移动;在局部区域发现颗粒下移路径并非为一条垂直直线,即存在颗粒水平移动的现象㊂通过计算v d =4㊁8mm /s 时颗粒在t =0~13s 内的位移得到正面和侧面的颗粒下移速率与水平速率,如图7所示㊂可知,颗粒水平速率很小,在0附近;下移速率呈中间颗粒速率大㊁两边速率小的趋势,与蔡九菊等[20]研究结果一致㊂进一步分析下移速率,发现颗粒位置比(颗粒在x 或y 方向坐标值的绝对值与该方向颗粒通道截面尺寸的比值)在0．15~0．85时,颗粒下移速率基本相同,靠近边壁颗粒下移速率为0,这与试验台壁面不平整有关㊂图6㊀直流段不同时刻颗粒位置Fig．6㊀Particle position mapat different time in straightflow section图7㊀直流段初始截面颗粒速率Fig．7㊀Particle rate at initial cross section in straight flow section㊀㊀忽略边壁的影响,分别计算直流段正面和侧面在v d =4mm /s 和v d =8mm /s 时的流动指数M F1,得到M F1,Z (4)=0．48㊁M F1,C (4)=0．44㊁M F1,Z (8)=0．60㊁M F1,C (8)=0．52(下标Z 和C 分别代表正面和29杨㊀帆等:下行移动床颗粒流动特性冷模试验研究2019年第2期侧面),均大于0．3,故认为在试验条件下,直流段颗粒流动为整体流状态㊂2．2㊀风速对颗粒层稳定性的影响图8为不同风速条件下进入二风室和四风室的颗粒携带速率v b 以及携带颗粒粒径分布㊂由图8(a)可知,随着风速的增大,二风室和四风室的颗粒携带速率均增大,其中四风室的颗粒携带速率明显大于二风室㊂由颗粒力链理论[21-22]可知,在上部载荷作用下,错流段内部颗粒被挤压形成力链网络,小颗粒是力链的组成部分或内嵌于力链网络,使其更难逃离力链网络㊂错流段上部载荷来自错流段上部颗粒重力,由于上部错流段所受上部载荷要小于下部错流段所受上部载荷,故四风室携带速率大于二风室携带速率㊂由图8(b)㊁(c)可知,随着风速的增大,二风室和四风室中细颗粒所占份额逐渐减小,粗颗粒所占份额逐渐增大㊂由图8可知,与其他工况相比,风速v =1．77m /s 时,四风室颗粒携带速率v b 明显增大,且粗颗粒份额也明显增大,可认为大量颗粒进入风室,即颗粒层稳定性被破坏㊂利用式(6)计算平均粒径所对应的流化风速为1．45m /s(图9),可认为风速超过1．45m /s 时,颗粒层稳定性被破坏㊂因此对于试验物料,为保证颗粒层稳定,风速不能超过1．45m /s㊂图8㊀风室颗粒携带特性Fig．8㊀Particle carrier characteristics inwindbox图9㊀颗粒平均粒径对应的临界流化风速Fig．9㊀Critical fluidization air speed corresponding to averageparticle diameter2．3㊀扩大段中下行颗粒的流动行为颗粒整体下移速率v d =4mm /s 时,扩大段颗粒运动动态如图10所示㊂通过空腔面积A c 来判断颗粒的流动㊂由图10可知,空腔面积A c 增大与减小交替出现,这是由于随着颗粒层整体向下移动,倾斜面附近空腔面积增大,由于倾斜壁面倾斜角41ʎ大于颗粒堆积角30ʎ,紧挨倾斜面的颗粒向右下移动填充空腔,导致空腔面积减小,颗粒整体下移与颗粒填充基本同时进行,空腔面积A c 也会随之增大或减小㊂当空腔中某点被颗粒占据时,空腔中会出现断点,如图中a 点所示,且断点会沿空腔方向向下移动㊂空腔面积A c 始终在20cm 2附近波动㊂图10㊀扩大段周期动态图Fig．10㊀Periodic dynamics diagram in enlarged section颗粒下移速率v d =4mm /s 时,扩大段颗粒不同时刻位置如图11所示㊂可知,颗粒整体向右下迁移,越靠近倾斜面,颗粒水平迁移距离越大,下移距离也越大,而越靠近竖直面,颗粒迁移距离越小㊂通过计算t =4s 与t =0间的颗粒位移得到颗粒下移速率与水平速率,如图12所示㊂由图12可知,随着颗粒位置由竖直壁面向倾斜壁面变化,颗粒水平速率与下移速率均增大㊂靠近竖直壁面下移速率与水平速率均较小;靠近倾斜壁面颗粒下移速率与水平速率较大,甚至出现颗粒水平速率大于颗粒下移速率的现象,处于中间位置颗粒下移速率与水平392019年第2期洁净煤技术第25卷图11㊀扩大段颗粒不同时刻位置Fig．11㊀Particle position map at different time inenlarged section速率比较平稳,呈现整体向右下移动,且颗粒速率较大和较小的区域占比较小,约为30%㊂图12㊀扩大段初始截面颗粒速率Fig．12㊀Particle rate at initial cross section in enlargedsection综上可知,虽然颗粒在扩大段局部出现速率较大或较小的流动,但颗粒在扩大段大部分区域呈现整体流;虽然倾斜面间断性出现20cm 2左右的空腔,但颗粒基本能填充整个扩大段㊂实际运行时设置一定长度的入口段,可有效弥补扩大段因空腔对气固流动造成的影响㊂2．4㊀扩大段颗粒粒径分布特性通过称量取样点M 1㊁M 2㊁M 3㊁M 4的取样质量得到各取样点颗粒质量分数(与取样总质量之比),如图13所示㊂在各取样点颗粒质量分数中,处于中间位置的占比略高于两边,即在颗粒下移流动过程中,沿深度方向中间位置的颗粒速率略高于两边的颗粒速率㊂各取样点颗粒质量分数均在25%左右,最大值与最小值相差14．4%,在工程实践中,可认为扩大段内颗粒整体流动是均匀的㊂图13㊀各取样点质量占比Fig．13㊀Particle mass ratio at each location取样点M 1㊁M 2㊁M 3㊁M 4的颗粒粒径分布如图14所示㊂由图14可知,颗粒粒径分布整体趋势基本相同;但各取样点粒径分布有差别,中位径d 50(M 4)>d 50(M 3)>d 50(M 1)>d 50(M 2),说明颗粒经入口段进入扩大段,大颗粒会沿着倾斜面流动至倾斜面末端,而细颗粒因体积小,在沿倾斜面移动过程中从大颗粒间的缝隙下移,水平迁移距离短,故细颗粒在取样点M 2累计最大㊂图14㊀各取样点颗粒粒径分布Fig．14㊀Particle size distribution at each location综上可知,扩大段内颗粒整体流动是均匀的,但不同的位置颗粒粒径分布存在差异,在实际应用过程中应注意此类问题,并对其进行优化㊂3㊀结㊀㊀论1)提出一种针对大宗固体颗粒余热回收的下行移动床技术方案,并在冷态试验台上利用CFB 锅炉炉渣进行了下行颗粒流动特性试验㊂颗粒移动过程中,直流段截面中部颗粒下行速度基本一致,呈整体流状态;而受边壁的影响,边壁处颗粒下移速率较慢,甚至不移动;在局部区域发现颗粒存在水平移动现象;忽略边壁的影响,计算的直流段正面和侧面流动指数M F1均大于0．3,故认为在试验条件下,直流段颗粒流动为整体流状态㊂2)随着横向风速的增大,二风室和四风室的颗49杨㊀帆等:下行移动床颗粒流动特性冷模试验研究2019年第2期粒携带速率均增大,且细颗粒占比逐渐减小,粗颗粒占比逐渐增大,其中四风室的颗粒携带速率明显大于二风室的颗粒携带速率;在风速大于1．45m/s 后,四风室携带速率明显增大,粗颗粒占比明显增大,判断此时颗粒层稳定性被破坏,即对于试验物料,流过下行移动床的水平风速不能超过1．45m/s㊂3)在扩大段,颗粒填充整个扩大段,但倾斜壁面会出现20cm2左右的空腔;扩大段内的大部分区域呈整体流,靠近竖直壁面流动速率较小,靠近倾斜壁面流动速率较大;通过对不同位置颗粒取样,发现取样点颗粒质量占比最大相差约14．4%,在实际工程中可认为扩大段颗粒整体流动均匀;分析取样点颗粒粒径分布,发现粒径分布整体趋势基本相同㊂参考文献(References):[1]㊀邵卫.篦冷机内传热机理与实验研究[D].青岛:山东大学,2017.[2]㊀王美琪.水泥篦冷机高温渗流换热规律及温度预测模型研究[D].秦皇岛:燕山大学,2016.[3]㊀甘露.循环流化床锅炉双喷动床式冷渣器研究[D].重庆:重庆大学,2014.[4]㊀曾兵.循环流化床锅炉选择性排渣冷却系统研究[D].重庆:重庆大学,2012.[5]㊀李朋.高炉渣余热回收及碳资源协同减排应用基础研究[D].沈阳:东北大学,2013.[6]㊀刘军祥.高炉渣余热回收装置传热特性实验研究[D].沈阳:东北大学,2009.[7]㊀郑斌,刘永启,李瑞阳,等.高温煅烧石油焦排料过程余热回收[J].化工进展,2015,34(6):1539-1543.ZHENG Bin,LIU Yongqi,LI Ruiyang,et al.Experimental investi-gation on waste heat reutilization of calcined petroleum coke[J].Chemical Industry and Engineering Progress,2015,34(6):1539-1543.[8]㊀陈光辉,李升大,陶少辉,等.焦炉余热综合利用研究进展[J].化工进展,2018,37(10):3799-3805.CHEN Guanghui,LI Shengda,TAO Shaohui,et al.Application and research of process of comprehensive utilization of coke oven waste heat[J].Chemical Industry and Engineering Progress,2018,37(10):3799-3805.[9]㊀范锦忠.提高陶粒回转窑热效率的有效措施[J].墙材革新与建筑节能,2006(11):29-32.FAN Jinzhong.Measure for improving the thermal efficiency of ro-tary ceramisite kiln[J].Wall Material Innovation and Energy Sav-ing in Buildings,2006(11):29-32.[10]㊀邬万竹.660MW超超临界CFB锅炉冷渣器选型技术经济性研究[J].洁净煤技术,2017,23(2):103-107.WU Wanzhu.Technical economic study on type selection ofslag cooler for a660MW ultra-supercritical CFB boiler[J].Clean Coal Technology,2017,23(2):103-107. [11]㊀岑可法.循环流化床锅炉理论㊁设计与运行[M].北京:中国电力出版社,1998.[12]㊀李心宁.篦冷机冷却系统的数值模拟与实验研究[D].青岛:山东大学,2016.[13]㊀朱建新.基于图像法的流化床内颗粒测试系统及管式流化床内颗粒混合研究[D].杭州:浙江大学,2004.[14]㊀孙苏皖.气固流化床内宽粒径分布颗粒的流动特性研究[D].重庆:重庆大学,2016.[15]㊀陶珍东,郑少华.粉体工程与设备[M].北京:化学工业出版社,2010.[16]㊀陶贺,金保昇,钟文琪.不同物性对椭球形颗粒在移动床中流动特性影响的模拟研究[J].中国电机工程学报,2011,31(5):68-75.TAO He,JIN Baosheng,ZHONG Wenqi.Effect of particle proper-ties on the flow behaviors of ellipsoidal particles in the movingbed[J].Proceedings of the CSEE,2011,31(5):68-75. [17]㊀JOHANSON J R,JENIKE A W.Stress and velocity fields in gravi-ty flow of bulk solids[M].Salt Lake City:University of Utah,1962:161-162.[18]㊀BROWN R L,RICHARDS J C.Kinematics of the flow of dry pow-ders and bulk solids[J].Rheologica Acta,1965,4(3):153-165.[19]㊀POLDERMAN H G,BOOM J,HILSTER E D,et al.Solids flowvelocity profiles in mass flow hoppers[J].Chemical Engineer-ing Science,1987,42(4):737-744.[20]㊀蔡九菊,董辉,傅巍.竖式移动床颗粒流动的研究[J].东北大学学报(自然科学版),2007,28(11):1599-1603.CAI Jiuju,DONG Hui,FU Wei.Study on particle flow inshaft moving beds[J].Journal of Northeastern University(Natu-ral Science),2007,28(11):1599-1603.[21]㊀孙其诚,辛海丽,刘建国,等.颗粒体系中的骨架及力链网络[J].岩土力学,2009,30(S1):83-87.SUN Qicheng,XIN Haili,LIU Jianguo,et al.Skeleton andforce chain network in granular systems[J].Rock and soil Me-chanics,2009,30(S1):83-87.[22]㊀孙其诚,王光谦.颗粒物质力学导论[M].北京:科学出版社,2009.59。