数学建模论文符号

1、一摘要 第二问中，本文将贮料的分布情况分别在整体流动和管状流动（漏斗状流动）下进行考虑。对于仓内贮料的分布情况，可通过贮料容量随时间的变化关系进行描述，并根据贮料容量与出入料速度之间的等量关系建立颗粒流微分方程模型。但由于出料速度在两个状态时存在很大差异故需分开求解，在整流状态下可结合能量守恒定律，得出筒仓内物料上表面的高度与出料速度的关系；在管状流动（漏斗状流动）状态下，由于贮料上表面曲面不规则且较为复杂，通过简化流速公式求解出料速度。最后将出料速度代入微分方程可求得贮料容量随时间的表达式，即可表征物料分布情况。二问题分析 在探讨筒仓仓内产品的分布与堆积问题时，为了描述分布与堆积状况，考虑到

2、不同状态时仓筒贮料容量会存在差异可以反映堆积特征，故可用容量来进行描述。首先先对煤炭筒仓内物料流动形态进行分析，查阅文献可知贮料流动形态大致分为两阶段，一种为整体流动，即卸料时贮料整体向下运动；一种为管状流动（漏斗状流动），即卸料时贮料从其内部的流动腔中流动。卸料时，筒仓上部颗粒先进行整体流动，当贮料到达一定高度后转变为管状流动状态。接着进一步描述筒仓贮料的分布情况，考虑到贮料容量由入放料速度之差决定，可分别建立基于能量守恒的颗粒流微分方程模型与基于公式的颗粒流微分方程模型，对贮料流动平衡前后两种状态进行深入描述。3 符号说明符号表示符号意义筒仓中贮料实际容量入料速度出料口速度同时入放料开始后

3、时间物体动能物体势能筒仓各部分实际质量筒仓各部分贮料实际高度筒仓各部分横截面积流动系数重力加速度水力半径四模型假设1.物料与仓壁、物料与物料之间的摩擦系数为0。2.煤质可视为流态的颗粒流。3.不考虑卸料口成拱堵塞等特殊情况。4.假设产品入料的速度为可人为控制的定值。五模型建立与求解5.1.模型分析 卸料开始时，靠近仓壁的贮料与中心线的贮料运动位移大致相同，使得贮料水平截面保持水平，此时为整体流动状态；当颗粒到达一定高度时，靠近仓壁贮料位移比中心线附近贮料运动的慢，这时流动状态为管道流动，如图一所示。图一：贮料流动形式 接着进一步描述流动过程中贮料的分布情况，由容量与入放料速率的关系,可引入颗粒

4、流微分方程模型。考虑两种状态下放料速率的影响因素不同，故应分开求解。对于整体流动，其运动形式较为规则，可采用能量守恒定律求解放料口速度。反之，管状流动运动形式不规则，运用能量守恒分析非常复杂且与实际情况相差较大，故这里采用公式简化计算放料口速度。5.2.颗粒流微分方程建立与求解 根据筒仓容量变化的物理规律可知，贮料容量变化速率等于入料速率与出料速率之差，查阅文可知贮料投放速率为一人为可控定值，因此只需考虑筒仓出料速率，而出料速度既与贮料容量和流动方式有关，因此需分开考虑。在此建立以下颗粒流微分方程来总体考虑容量变化情况： 设为时间，为仓内贮料总量，为筒仓入料速率，为筒仓出料速率，可建立如下关系

5、式： 对于整体流动与管状流动，出料速率的求解方式不同，故建立两种模型分别求解。5.2.1能量守恒定律求解整流状态下出料速率计算的基本依据是能量守恒定，筒仓贮料流动前后的机械能总量相等，即式中：贮料流动前动能； 贮料流动前势能； 贮料流动后动能； 贮料流动后势能如图2所示，贮料流动前动能，只有势能；如图3贮料流动后情况，贮料向下流动了一段距离，此时贮料既有动能，又有势能。故该能量守恒方程可化为： 即 图2：仓斗贮料流动情况 图3：仓斗贮料流动后情况对于贮料流动前势能，包括仓斗上部圆柱体贮料势能与仓斗下部椎体贮料势能，取仓斗中流出物料的下表面为零势能面，如图3所示线；对于贮料流动后势能，包括仓斗上

6、部圆柱体贮料势能、仓斗下部椎体贮料势能与仓斗流出贮料势能，则有： 式中：仓斗上部圆柱体贮料质量； 仓斗中流出贮料的质量； 仓斗上部圆柱体贮料长度； 仓斗下部椎体贮料长度； 仓斗中流出物料的长度； 上部圆柱体横截面积； 下部椎体出口处面积而由图2图3可知下部椎体贮料势能在流动前后无变化，即，故把式代入式得：以下是分步对动能的计算过程： 此处由三部分组成，分别是仓斗上部圆柱体物料的动能 ，仓斗下部锥体物料的动能和从仓斗中流出物料的动能 ，计算图例见图三。再设上部分圆柱体物料的流速为,从仓斗中流出物料的流速为.可得如下关系式：说明：此处上仓都中流出的物料是不均匀的，因此在实际当中上面物料流动慢而下面

7、流动快，但由于假设很小导致上下流速相差不大，故可认为上下流速相等。在实际中仓斗上部和下部流动不均匀，但此处假设流动是均匀的。 对于计算要相对复杂一些，需用微积分思想求解：图4：下部椎体示意图 图4为仓斗下部锥体，可在图中距离锥体顶部为的地方取微圆 柱体，取圆柱体横截面积为，另设物料比重为，则则圆柱体质量为：圆柱体速度为：因此圆柱体动能为：同时对进行积分可得：将关系式代入公式，可得表达式： 但是若在实际中使用以上公式则比较繁琐，考虑实际情况一般,故，同时有，故可得到以下简化表达式： 此即为整流状态下出料速率公式。5.2.2 Janssen公式求解管状流状态下出料速率 根据Janssen公式，在漏

8、斗、管状流型下，卸料孔作中心布置时，其排料速度可按下列公式计算 式中：流动系数，决定于物料的流动性及其粒度。对于干燥易流的散体物料，可取。（这里将产品视为干燥易流的散体物料） 重力加速度（）； 水力半径（）。 卸料孔处的水利半径按如下式中：圆形排料孔的直径物料的最大典型物料尺寸综上，产品的流速的计算公式如下： 此即为管道流状态下出料速率公式。由于在标准形式下，产品的流速与产品的堆积高度无关，与卸料口的水力半径有关。根据题意，讨论对象为只有一个进口，一个开口的的理想筒仓，水力半径为一定值，故产品的流速为一定值。5.2.3.颗粒流微分方程模型此时即可得到两种状态下仓底出料速率，即为将代入中，输入m

9、atlab进行求解可得： 式中：此处表达式中函数是欧米加函数（或乘数对数），其满足的关系，常出现在微分方程符号解中，当其中各个符号赋予具体值时可以求得其为一个常数。则此式可表示贮料总量随时间的变化情况。 此即为最终颗粒流微分方程模型。6 模型评价与改进 6.1模型优缺点本文的颗粒流微分方程模型是建立在对贮料流态的定性分析之上，考虑了整体流与管状流两种形态，并分别对其进行了定量分析，给出两种形态下具体的表达式，准确描述出容量的分布堆积情况。定性与定量的结合，可以较为全面的对贮料分布堆积情况进行描述。但此微分方程模型仍存在不足，对于贮料流态的分析上不够深入，并未在整体流与管状流之下深入分析贮料分布变化情况；对于微分方程建立上，不能够求得统一的模型对所有流态进行描述，只做到分别描述流态，因此存在缺陷。 6.2模型改进在求解整体流状态下的出料速率中，由于假定了物料与仓壁，物料与物料的摩擦系数为0，所以求得速率比实际值大，故此处应加入一个修正系数,则可将式子变为。其中：K为与物料性质有关的参数。7 参考文献翟振威 原国平 张峰涛，筒仓贮料流态的颗粒流数值模拟 ，山西建筑 ，第34卷35期，200