发明名称 一种叶片式抛送装置出料管结构形状优化设计方法
摘要 本发明公开了一种叶片式抛送装置出料管结构形状优化设计方法,其特征在于:步骤如下:1)、通过高速摄像试验观察出料管中物料的运动规律,结合理论分析建立出料管中物料运动的数学模型;2)、采用欧拉‑拉格朗日法对出料管中的气固两相流动进行数值仿真,其中物料运动规律采用用户自定义函数,为上述已建物料运动数学模型,并将数值仿真结果与高速摄像试验结果比较验证其可靠性;3)、在相同工况情况下,选取出料管结构形状及参数,对其内部气固两相流动进行仿真,利用物料离开出料管的运动速度的仿真结果,确定出料管结构形状的优化结果。本发明缩短了开发周期,而且提高了设计效率、设计质量和一次开发的成功率。
申请公布号 CN103617331B 申请公布日期 2016.08.17
申请号 CN201310675254.9 申请日期 2013.12.13
申请人 内蒙古工业大学 发明人 翟之平;周蕊;周岗;朱明新;杨忠义
分类号 G06F17/50(2006.01)I;G06F9/455(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京国坤专利代理事务所(普通合伙) 11491 代理人 郭伟红
主权项 一种叶片式抛送装置出料管结构形状优化设计方法,其特征在于:步骤如下:1)、通过高速摄像试验观察出料管中物料的运动规律,结合理论分析建立出料管中物料运动的数学模型;2)、采用计算流体力学软件Fluent中的欧拉‑拉格朗日法对出料管中的气固两相流动进行数值仿真,其中物料运动规律采用用户自定义函数,为上述已建物料运动数学模型,并将数值仿真结果与高速摄像试验结果比较验证其可靠性;3)、在相同工况情况下,选取出料管结构形状及参数,对其内部气固两相流动进行仿真,利用物料离开出料管的运动速度的仿真结果,确定出料管结构形状的优化结果;所述步骤1)中建立抛送装置出料管内物料运动的数学模型:由高速摄像试验可观察到物料颗粒离开抛送叶片后,以分散的方式沿出料直管垂直向上运动,直到与出料弯管上壁碰撞,然后以密集的方式沿弯管上壁运动直至被抛离出料管;(1)物料颗粒刚离开抛送叶片时,受到的离心力远大于气流所受到的离心力,因此物料速度υ<sub>p</sub>大于气流平均速度υ<sub>a</sub>,气流对物料的作用力P<sub>a</sub>为气动阻力,方向与物料运动方向相反;物料在重力、摩擦碰撞阻力及气动阻力作用下以分散的方式沿出料直管垂直向上运动,其速度逐渐下降,当下降到低于气流速度时,气流对物料产生向上的气动力,维持物料向上的速度;出料直管内物料运动的方程式为:<maths num="0001"><math><![CDATA[<mrow><msub><mi>m</mi><mi>p</mi></msub><mfrac><mrow><msup><mi>d</mi><mn>2</mn></msup><mi>l</mi></mrow><mrow><msup><mi>dt</mi><mn>2</mn></msup></mrow></mfrac><mo>=</mo><msub><mi>P</mi><mi>a</mi></msub><mo>-</mo><msub><mi>m</mi><mi>p</mi></msub><mi>g</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000968953630000011.GIF" wi="879" he="126" /></maths>方程式中,mp为物料颗粒质量,单位为kg;l为颗粒流运动路径,单位为m;t为颗粒运动时间,单位为s;g为重力加速度,单位为m/s2;Pa为气动力,单位为N;<img file="FDA0000968953630000021.GIF" wi="590" he="167" />其中υs为物料颗粒悬浮速度,单位为m/s;(2)物料进入出料弯管后与其外管壁碰撞,设碰撞为非弹性碰撞,碰撞后的平均速度为υ<sub>2</sub>,则<maths num="0002"><math><![CDATA[<mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>m</mi><mi>p</mi></msub><msubsup><mi>&upsi;</mi><mn>2</mn><mn>2</mn></msubsup><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>m</mi><mi>p</mi></msub><msubsup><mi>&upsi;</mi><mn>1</mn><mn>2</mn></msubsup><mo>-</mo><mn>1</mn><mo>/</mo><mn>2</mn><msub><mi>m</mi><mi>p</mi></msub><msubsup><mi>&upsi;</mi><mn>1</mn><mn>2</mn></msubsup><msup><mi>sin</mi><mn>2</mn></msup><mi>&alpha;</mi><mo>+</mo><msub><mi>m</mi><mi>p</mi></msub><msub><mi>&upsi;</mi><mn>1</mn></msub><msub><mi>sin&alpha;f&upsi;</mi><mi>m</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000968953630000022.GIF" wi="1553" he="119" /></maths>式中,υ<sub>1</sub>为物料与管壁碰撞前的平均速度,单位为m/s;α为碰撞角;f为物料与管壁的摩擦系数;υ<sub>m</sub>为碰撞过程中的平均速度,<img file="FDA0000968953630000023.GIF" wi="910" he="111" />单位为m/s;(3)物料与管壁碰撞后在离心力的作用下,与气流分离,在重力m<sub>p</sub>g、离心力<img file="FDA0000968953630000024.GIF" wi="192" he="119" />及摩擦力<img file="FDA0000968953630000025.GIF" wi="526" he="164" />作用下,以密集的方式沿弯管上壁运动,物料运动方程为:<maths num="0003"><math><![CDATA[<mrow><msub><mi>m</mi><mi>p</mi></msub><mfrac><mrow><msup><mi>d</mi><mn>2</mn></msup><mi>l</mi></mrow><mrow><msup><mi>dt</mi><mn>2</mn></msup></mrow></mfrac><mo>+</mo><msub><mi>m</mi><mi>p</mi></msub><mi>g</mi><mi> </mi><mi>c</mi><mi>o</mi><mi>s</mi><mfrac><mi>l</mi><mi>r</mi></mfrac><mo>+</mo><mi>f</mi><mo>&lsqb;</mo><mfrac><msub><mi>m</mi><mi>p</mi></msub><mi>r</mi></mfrac><msup><mrow><mo>(</mo><mfrac><mrow><mi>d</mi><mi>l</mi></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msub><mi>m</mi><mi>p</mi></msub><mi>g</mi><mi> </mi><mi>s</mi><mi>i</mi><mi>n</mi><mfrac><mi>l</mi><mi>r</mi></mfrac><mo>&rsqb;</mo><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000968953630000026.GIF" wi="1446" he="166" /></maths>式中,r为出料弯管外壁曲率半径,单位为m;所述步骤2)中基于欧拉‑拉格朗日法建立出料管中气固两相流动数值仿真模型:选取整个叶片式抛送装置的流动空间作为计算区域,在CFD的前处理软件GAMBIT中生成计算区域的实体模型,并对计算区域进行网格划分;计算区域的进Ι口设为速度进口边界条件,出口设为压力出口,给定标准大气条件作为压力边界值;采用欧拉‑拉格朗日法模拟出料管内的气固两相流动,即在欧拉坐标系下考察气流的运动,在拉格朗日坐标系下研究固体颗粒的运动;选取Fluent软件中的DPM(Discrete Phase Model)模型,考虑物料流和气流的相互作用,采用耦合计算;其中出料管中物料运动规律采用用户自定义UDF(User‑Defined Functions)函数,为上述已建物料运动数学模型;将数值模拟结果与出料管内物料运动的高速摄像试验结果比较,验证数学模型及数值计算结果的可靠性。
地址 010051 内蒙古自治区呼和浩特市新城区爱民街49号内蒙古工业大学机械学院