CN108228963A - 一种烟叶烤机内部流场分析与调整方法 - Google Patents
一种烟叶烤机内部流场分析与调整方法 Download PDFInfo
- Publication number
- CN108228963A CN108228963A CN201711220794.2A CN201711220794A CN108228963A CN 108228963 A CN108228963 A CN 108228963A CN 201711220794 A CN201711220794 A CN 201711220794A CN 108228963 A CN108228963 A CN 108228963A
- Authority
- CN
- China
- Prior art keywords
- flow field
- machine
- orifice plate
- formula
- coefficient
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- A—HUMAN NECESSITIES
- A24—TOBACCO; CIGARS; CIGARETTES; SIMULATED SMOKING DEVICES; SMOKERS' REQUISITES
- A24B—MANUFACTURE OR PREPARATION OF TOBACCO FOR SMOKING OR CHEWING; TOBACCO; SNUFF
- A24B3/00—Preparing tobacco in the factory
- A24B3/10—Roasting or cooling tobacco
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Manufacturing Of Cigar And Cigarette Tobacco (AREA)
- Manufacture Of Tobacco Products (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种烟叶烤机内部流场分析与调整方法,通过简化拷机内部设置的风机和孔板的计算模型,建立烟草烤机内部空气流动区域的三维模型,对其进行流体动力学计算,获得拷机内部流场信息,以不均匀系数为评价标准,统计流场内关键区域的流场参数,对烤机内的流场均匀性进行评价,不均匀系数越小表明气流分布越均匀;以原始方案中烤机下方孔板的开孔率和倾斜角度为因素变量,设计正交实验,得到最优因素水平组合;本发明对烟草烤机内部空气的速度场、静压场进行分析和评价,为提高流场均匀性提供参考,可为改进烤机设计提供指导。
Description
技术领域
本发明涉及烟草生产领域,具体的说是一种烟叶烤机内部流场分析与调整方法。
背景技术
在对烟叶进行处理的过程中需要使烟叶湿度达到要求的范围。烟草烤机先将空气加热,使用风机将温度高、湿度低的空气经过下层孔板和网袋吹向烟叶,带走烟叶内水分。烟叶平铺在网袋和上层孔板之间,空气流场的均匀程度直接影响到烟叶湿度的均匀性。传统的实验测量方法时间和经济成本较高,而计算流体力学模拟的方法能够以低成本和较快的速度模拟烤机内流场,为改进烤机设计提供指导。
发明内容
本发明的目的是提供一种烟叶烤机内部流场分析与调整方法,对烟草烤机内部空气的速度场、静压场进行分析和评价,为提高流场均匀性提供参考。
为解决上述技术问题,本发明采用的技术方案为:
一种烟叶烤机内部流场分析与调整方法,包括以下步骤:
步骤一:对烟草烤机内部空气流动区域建立三维模型;
步骤二:以流体力学中的二维风机模型为基础,通过压差跳跃的方法模拟流场中的风机,保留风机壳及出口影响空气流动的结构特征;
步骤三:采用二维风机模型,在计算模型中赋予平面一个跳跃的压差,这个压差可以为常数,也可以为速度的函数,压差与速度的关系用公式(1)表示:
(1)
式中:为风机两侧压差;为压差的多项式的系数;为垂直于风扇的速度;
步骤四:建立孔板的计算流体力学模型,以不同速度为计算的输入,以孔板两侧对应的压降为输出,得到一组自变量和因变量;多孔介质的参数依据公式(2)确定;
(2)
式中:为多孔介质两侧压差;为多孔介质厚度;为多孔介质的渗透率系数;为多孔介质的阻力系数;为空气密度;为空气粘度;
针对多孔介质特征参数计算结果的正负不一致问题,将孔板以及孔板前后区域视为一个整体,将其压降和厚度代入公式(2)进行计算,得到多孔介质的特征参数;
步骤五:以不均匀系数为评价标准,统计流场内关键区域的流场参数,对烤机内的流场均匀性进行评价;
流场参数不均匀系数计算方法如下:
在工作区取n个测点,测得其流场参数,求算数平均值:
(3)
计算均方根偏差:
(4)
不均匀系数为:
(5)
不均匀系数越小表明气流分布越均匀;
步骤六:以原始方案中烤机下方孔板的开孔率和倾斜角度为因素变量,依据正交实验设计方法,设计正交实验表;
步骤七:对正交实验表内每个方案进行数值模拟实验,对得到的流场计算速度和静压不均匀系数,并对正交实验结果进行极差分析,得到最优因素水平组合;
步骤八:依据最优因素水平组合,再次进行数值模拟实验,将结果与原始方案和正交实验方案进行对比,得到最优调整方案。
公式(2)描述了压降和速度的关系;对速度和压降进行二次曲线的拟合,得到公式(2)的系数,依据公式(2)计算,得到多孔介质的特征参数。
设置入口为压力入口,湍流模型为模型,流体介质为理想气体。
采用正交实验与不均匀系数评价标准相结合的方法,选取最优实验因素水平组合为最终装置的优化方案;其中,不均匀系数计算中涉及的参数均由数值实验获得。
本发明的有益效果:
本发明提供的烟叶烤机内部流场分析与调整方法,通过建立烟草烤机内部空气流动区域的三维模型,对其进行流体动力学计算,获得拷机内部流场信息,以不均匀系数为评价标准,统计流场内关键区域的流场参数,对烤机内的流场均匀性进行评价,不均匀系数越小表明气流分布越均匀;以原始方案中烤机下方孔板的开孔率和倾斜角度为因素变量,设计正交实验,得到最优因素水平组合;本发明对烟草烤机内部空气的速度场、静压场进行分析和评价,为提高流场均匀性提供参考,可为改进烤机设计提供指导。
附图说明
图1 本发明烟草烤机结构示意图;
图2 本发明烟草烤机气流组织方案模拟优化流程图;
附图标记:1、进风口,2、换热器,3、风机,4、电机,5、出风口,6、上孔板,7、网袋,8、倾斜孔板,9、整个工作区域。
具体实施方式
下面结合具体实施方式对本发明做进一步的阐述。
一种烟叶烤机内部流场分析与调整方法,包括以下步骤:
步骤一:对烟草烤机内的空气流域、孔板、网袋、风机进行三维建模,并做适当的简化;
步骤二:以流体力学中的二维风机模型为基础,通过压差跳跃的方法模拟流场中的风机,保留风机壳及出口影响空气流动的结构特征;在计算模型中直接使用三维模型模拟风机的叶片旋转会增加模型前处理和模拟消耗的时间,适当简化风机模型可以加快计算速度,同时不影响计算结果准确性;设置入口为压力入口,湍流模型为模型,流体介质为理想气体,初始化之后开始计算;
步骤三:采用二维风机模型,在计算模型中赋予平面一个跳跃的压差,这个压差可以为常数,也可以为速度的函数,压差与速度的关系用公式(1)表示:
(1)
式中:为风机两侧压差;为压差的多项式的系数;为垂直于风扇的速度;
步骤四:建立孔板的计算流体力学模型,以不同速度为计算的输入,以孔板两侧对应的压降为输出,得到一组自变量和因变量;多孔介质的参数依据公式(2)确定;公式2描述了压降和速度的关系,对速度和压降进行二次曲线的拟合,得到公式2的系数,依据公式2计算,得到多孔介质的特征参数;
孔板上分布数千个小孔,造成计算的预处理、计算过程、结果后处理等过程繁琐、缓慢。为减小计算量,加快计算速度,同时不影响模拟结果的准确性,将孔板视为二维多孔介质,通过制定多孔介质的渗透率、阻力系数、厚度等特征参数,使计算过程中多孔介质能够取代孔板;
(2)
式中:为多孔介质两侧压差;为多孔介质厚度;为多孔介质的渗透率系数;为多孔介质的阻力系数;为空气密度;为空气粘度;
针对多孔介质特征参数计算结果的正负不一致问题,将孔板以及孔板前后区域视为一个整体,将其压降和厚度代入公式(2)进行计算,得到多孔介质的特征参数;
步骤五:以不均匀系数为评价标准,统计流场内关键区域的流场参数,对烤机内的流场均匀性进行评价;
流场参数不均匀系数计算方法如下:
在工作区取n个测点,测得其流场参数,求算数平均值:
(3)
计算均方根偏差:
(4)
不均匀系数为:
(5)
不均匀系数越小表明气流分布越均匀;如果此区域局部存在较大的参数梯度,而对整个区域只求一个不均匀系数,该局部流场的不均匀性就会被忽略,影响最终的烟叶烘干效果。因此将该长方体工作区域依据坐标平均分为九块,对每一块取样点中的速度、压力,计算速度、压力的不均匀系数;
步骤六:以原始方案中烤机下方三块孔板的开孔率和倾斜角度为因素变量,依据四因素三水平的正交实验设计方法,设计正交实验表;正交实验表中的因素水平依据原始方案中的因素水平加减百分之二十确定;
步骤七:对正交实验表内每个方案进行数值模拟实验,依据公式(3)~(5),对每一个实验的模拟结果计算速度和静压不均匀系数,并对每一组实验的速度和静压不均匀系数分别进行极差分析,得到两组最优因素水平组合;
步骤八:依据两组最优水平组合再次进行数值模拟实验,将实验结果与原始方案进行对比,以不均匀系数为评价标准,分别得到速度与静压均匀性最优的实验方案。
依据速度、压力不均匀系数,采取向复烤机内流场添加导流板、改变孔板孔径和间距等方法,使各块不均匀系数减小并且接近,改善流场分布均匀性。
本发明通过建立烟草烤机内部空气流动区域的三维模型,以不均匀系数为评价标准,统计空气流场内关键区域的流场参数,对烤机内的流场均匀性进行评价,不均匀系数越小表明气流分布越均匀;以原始方案中烤机下方孔板的开孔率和倾斜角度为因素变量,设计正交实验,得到最优因素水平组合;本发明对烟草烤机内部空气的速度场、静压场进行分析和评价,为提高流场均匀性提供参考,可为改进烤机设计提供指导。
Claims (4)
1.一种烟叶烤机内部流场分析与调整方法,其特征在于,包括以下步骤:
步骤一:对烟草烤机内部空气流动区域建立三维模型;
步骤二:以流体力学中的二维风机模型为基础,通过压差跳跃的方法模拟流场中的风机,保留风机壳及出口影响空气流动的结构特征;
步骤三:采用二维风机模型,在计算模型中赋予平面一个跳跃的压差,这个压差可以为常数,也可以为速度的函数,压差与速度的关系用公式(1)表示:
(1)
式中:为风机两侧压差;为压差的多项式的系数;为垂直于风扇的速度;
步骤四:建立孔板的计算流体力学模型,以不同速度为计算的输入,以孔板两侧对应的压降为输出,得到一组自变量和因变量;多孔介质的参数依据公式(2)确定;
(2)
式中:为多孔介质两侧压差;为多孔介质厚度;为多孔介质的渗透率系数;为多孔介质的阻力系数;为空气密度;为空气粘度;
针对多孔介质特征参数计算结果的正负不一致问题,将孔板以及孔板前后区域视为一个整体,将其压降和厚度代入公式(2)进行计算,得到多孔介质的特征参数;
步骤五:以不均匀系数为评价标准,统计流场内关键区域的流场参数,对烤机内的流场均匀性进行评价;
流场参数不均匀系数计算方法如下:
在工作区取n个测点,测得其流场参数,求算数平均值:
(3)
计算均方根偏差:
(4)
不均匀系数为:
(5)
不均匀系数越小表明气流分布越均匀;
步骤六:以原始方案中烤机下方孔板的开孔率和倾斜角度为因素变量,依据正交实验设计方法,设计正交实验表;
步骤七:对正交实验表内每个方案进行数值模拟实验,对得到的流场计算速度和静压不均匀系数,并对正交实验结果进行极差分析,得到最优因素水平组合;
步骤八:依据最优因素水平组合,再次进行数值模拟实验,将结果与原始方案和正交实验方案进行对比,得到最优调整方案。
2.如权利要求1所述的烟叶烤机内部流场分析与调整方法,其特征在于:拷机下方布置的均流孔板简化为二维多孔介质模型,公式(2)描述了孔板两侧压降和速度的关系;对速度和压降进行二次曲线的拟合,得到公式(2)的系数,依据公式(2)计算,得到多孔介质的特征参数。
3.如权利要求1所述的烟叶烤机内部流场分析与调整方法,其特征在于:设置入口为压力入口,湍流模型为模型,流体介质为理想气体。
4.如权利要求1所述的烟叶烤机内部流场分析与调整方法,其特征在于:采用正交实验与不均匀系数评价标准相结合的方法,选取最优实验因素水平组合为最终装置的优化方案;其中,不均匀系数计算中涉及的参数均由数值实验获得。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711220794.2A CN108228963B (zh) | 2017-11-29 | 2017-11-29 | 一种烟叶烤机内部流场分析与调整方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711220794.2A CN108228963B (zh) | 2017-11-29 | 2017-11-29 | 一种烟叶烤机内部流场分析与调整方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108228963A true CN108228963A (zh) | 2018-06-29 |
CN108228963B CN108228963B (zh) | 2021-04-09 |
Family
ID=62653042
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711220794.2A Active CN108228963B (zh) | 2017-11-29 | 2017-11-29 | 一种烟叶烤机内部流场分析与调整方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108228963B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109448767A (zh) * | 2018-12-21 | 2019-03-08 | 郑州云海信息技术有限公司 | 一种用于存储设备检测的模拟风道 |
CN115859508A (zh) * | 2022-11-23 | 2023-03-28 | 北京百度网讯科技有限公司 | 流场分析方法、元件模型生成方法、训练方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2548095Y (zh) * | 2002-02-27 | 2003-04-30 | 歙县农业机械化管理局 | 热风烘干装置 |
CN107391845A (zh) * | 2017-07-23 | 2017-11-24 | 南京理工大学 | 一种烟草除杂装置流场的数值模拟方法 |
-
2017
- 2017-11-29 CN CN201711220794.2A patent/CN108228963B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2548095Y (zh) * | 2002-02-27 | 2003-04-30 | 歙县农业机械化管理局 | 热风烘干装置 |
CN107391845A (zh) * | 2017-07-23 | 2017-11-24 | 南京理工大学 | 一种烟草除杂装置流场的数值模拟方法 |
Non-Patent Citations (5)
Title |
---|
B.Y. GUO等: "Numerical modelling of the gas flow through perforated plates", 《CHEMICAL ENGINEERING RESEARCH AND DESIGN》 * |
庄江婷等: "烟草烘箱内部气流组织的优化", 《能源技术》 * |
杨圣印: "全地面起重机冷却系统空气流场仿真分析", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
王龙飞等: "基于CFD的新型密集烤房结构优化研究", 《江苏农业科学》 * |
金莎等: "基于CFD 方法的孔板送风气流组织优化研究", 《低温工程》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109448767A (zh) * | 2018-12-21 | 2019-03-08 | 郑州云海信息技术有限公司 | 一种用于存储设备检测的模拟风道 |
CN109448767B (zh) * | 2018-12-21 | 2024-04-26 | 郑州云海信息技术有限公司 | 一种用于存储设备检测的模拟风道 |
CN115859508A (zh) * | 2022-11-23 | 2023-03-28 | 北京百度网讯科技有限公司 | 流场分析方法、元件模型生成方法、训练方法及装置 |
CN115859508B (zh) * | 2022-11-23 | 2024-01-02 | 北京百度网讯科技有限公司 | 流场分析方法、元件模型生成方法、训练方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN108228963B (zh) | 2021-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106354977B (zh) | 一种利用cfd模型对室内热环境的分析方法 | |
CN107223195A (zh) | 用于hvac系统的可变风量建模 | |
CN101976276A (zh) | 一种空调器风道结构优化设计方法 | |
CN111551339B (zh) | 一种翼身融合飞机通气模型畸形流道专用测压耙设计方法 | |
CN108228963A (zh) | 一种烟叶烤机内部流场分析与调整方法 | |
CN109753716A (zh) | 基于流场仿真的核/火电汽轮机组流体激励数值计算方法及系统 | |
CN108536922A (zh) | 飞机及发动机的内外流场一体化数值模拟方法 | |
CN105550467B (zh) | 高速风洞斜孔壁壁面流动模型的建立方法 | |
Gan | Interaction between wind and buoyancy effects in natural ventilation of buildings | |
CN102607658B (zh) | 一种基于浓度法的复杂结构通道内气体流量测量方法 | |
CN108256262A (zh) | 一种轴对称射流稳压腔参数设计的数值模拟方法 | |
CN103557183A (zh) | 偏心的风扇外壳 | |
CN109186815A (zh) | 一种低温高马赫数测试用探针温度标定装置 | |
CN107179206A (zh) | 一种横流湿式冷却塔阻力特性试验装置及测试方法 | |
CN116467817A (zh) | 一种基于上进风燃烧器的风道设计方法 | |
CN106225158B (zh) | 送风风向控制方法及装置 | |
Wu et al. | Numerical simulation of harmful gas distribution in a range hood with an improved flow channel | |
CN108446452B (zh) | 一种混流泵叶轮鲁棒优化设计方法 | |
CN204188642U (zh) | 一种机箱散热风量测试装置 | |
US20150169797A1 (en) | Providing a model of impact of a cooling infrastructure | |
CN111157896B (zh) | 基于恒转矩控制的测试方法 | |
CN106126803B (zh) | 制冷系统模拟方法和装置 | |
CN116282299B (zh) | 热法海水淡化组合系统大数据采集系统及系统优化方法 | |
CN109829184A (zh) | 3d打印系统吸气口管道结构及其优化设计方法、装置 | |
CN219233465U (zh) | 一种实验室多通道导流风道 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |