基于CFD-DEM的排肥用波纹管结构优化设计与试验_张一帆.pdf
《基于CFD-DEM的排肥用波纹管结构优化设计与试验_张一帆.pdf》由会员分享,可在线阅读,更多相关《基于CFD-DEM的排肥用波纹管结构优化设计与试验_张一帆.pdf(11页珍藏版)》请在咨信网上搜索。
1、浙江农业学报 Acta Agriculturae Zhejiangensis,2023,35(1):191 201http:/www zjnyxb cn张一帆,何瑞银,段庆飞,等 基于 CFD-DEM 的排肥用波纹管结构优化设计与试验J 浙江农业学报,2023,35(1):191 201DOI:10.3969/j issn 1004-1524.2023.01.21收稿日期:2021-09-24基金项目:国家重点研发计划(2016YFD0300908);江苏省农业科技自主创新资金(CX(19)2012)作者简介:张一帆(1996),男,江苏南京人,硕士研究生,主要从事农机现代设计理论及方法研究。
2、E-mail:njauzhangyifan163 com*通信作者,何瑞银,E-mail:ryhe_njau163 com基于 CFD-DEM 的排肥用波纹管结构优化设计与试验张一帆,何瑞银*,段庆飞,徐勇(南京农业大学 工学院,江苏 南京 210031)摘要:为探究波纹管结构参数与输送均匀性的关系,获得最优参数组合以提高气力排肥均匀性。本文基于CFD-DEM 耦合仿真技术构建了以排肥均匀性变异系数为响应值的二次响应面模型,在单因素试验的基础上采用 BDD 法研究了波纹管幅宽、波纹间距和波长之间交互作用对排肥均匀性的影响,并获得了波纹管的最佳参数组合,最后进行台架试验对仿真模型与优化结果进行验
3、证。仿真试验结果表明,排肥均匀性变异系数随波长、幅宽、波纹间距的增加先减小后增加,波纹管结构参数对排肥均匀性的影响显著性顺序为幅宽 波长 波纹间距。波纹管最优结构参数为波长 18.103 mm,波纹间距 12.158 mm,幅宽 8.863 mm,此时排肥均匀性变异系数为 7.5%。单因素台架试验中尿素颗粒仿真预测值与台架实际值变化趋势一致,数值基本重合,最优参数组合波纹管排肥均匀性变异系数与仿真试验预测值相对误差的均值为 5.84%,说明仿真模型精确可靠。复合肥颗粒变化趋势与尿素颗粒基本一致,说明仿真模型具有普适性。本研究可为波纹管的使用与优化提供参考。关键词:波纹管;响应面;仿真试验;结构
4、优化;排肥均匀性中图分类号:S220.3文献标志码:A文章编号:1004-1524(2023)01-0191-11Numerical analysis of flow characteristics and structural optimization of bellows based onCFD-DEMZHANG Yifan,HE Ruiyin*,DUAN Qingfei,XU Yong(College of Engineering,Nanjing Agricultural University,Nanjing 210031,China)Abstract:In order to explo
5、re the relationship between bellows structural parameters and conveying uniformity,theoptimal parameter combination was obtained to improve the uniformity of pneumatic fertilizer discharge Based onCFD-DEM coupling simulation technology,this paper constructed a quadratic response surface model with t
6、he varia-tion coefficient of fertilizer discharge uniformity as the response value On the basis of single factor test,BDD methodwas used to study the influence of the interaction between bellows width,bellows spacing and wavelength on fertilizerdischarge uniformity,and the best parameter combination
7、 of bellows was obtained,Finally,the bench test was car-ried out to verify the simulation model and optimization results The simulation results showed that the variation coef-ficient of fertilizer discharge uniformity firstly decreased and then increased with the increase of wavelength,widthand ripp
8、le spacing The order of significance of the influence of bellows structural parameters on fertilizer discharge u-niformity was width wavelength ripple spacing The optimal structural parameters of bellows were wavelength 18.103mm,ripple spacing 12.158 mm and width 8.863 mm At this time,the variation
9、coefficient of fertilizer discharge uni-formity was 7.5%In the single factor bench test,the change trend of urea particle simulation predicted value wasconsistent with the actual value of the bench,and the values basically coincided The average value of the relative errorbetween the variation coeffi
10、cient of corrugated pipe fertilizer discharge uniformity of the optimal parameter combinationand the predicted value of the simulation test was 5.84%,which showed that the simulation model was accurate and relia-ble The change trend of compound fertilizer particles was basically consistent with that
11、 of urea particles,which showedthat the simulation model was universal This study can provide a reference for the use and optimization of bellowsKey words:corrugated pipe;response surface;simulation test;structural optimization;fertilizer uniformity肥料是作物的“粮食”,合理使用肥料可以改良土壤肥力、提高单位面积产量1 3。然而随着肥料使用量的增大与
12、不科学施肥方式的增多,国内施肥行业产生了肥料用量大、利用率低、面源污染严重等问题4 6。传统粗放的施肥方式已不适合新时代绿色农业的发展要求,急需新型施肥技术提高肥料利用率,减少使用量7 9。气力输肥系统由于排肥均匀性好,输送效率高近年来在施肥领域得到了广泛应用10 11。波纹管具有良好的理化特性,许多专家学者在气力输肥系统中使用波纹管作为关键输送部件并开展了相关研究。杨庆璐等12 研究了波纹管直径对气固两相流动性的影响,通过颗粒碰撞次数与颗粒分布图判断均匀性。张晓辉等13 以气流场内速度分布均匀性作为评价指标,通过对比直管与褶皱管内部速度场分布判断输送性能。戴亿政等14 通过空气流的单相仿真分
13、析了波纹管结构对内部气体速度场的影响。李中华等15 研究了波纹管内部单相压力场,得出了波纹管结构导致气流压力平缓降低。常金丽16 设计了一种褶皱管,用来增压提升输送均匀性。以上学者大都只针对波纹管某一结构参数变量进行定性研究,鲜有研究不同结构参数变量之间交互作用对输送均匀性的影响。在判断输送均匀性时大多依据非标准化的流场分布图,没有参数化的指标进行均匀性判断。本文旨在通过 CFD-DEM 耦合仿真分析技术,在单因素试验的基础上研究波纹管幅宽、波纹间距和波长之间交互作用对输送均匀性的影响,通过仿真试验明确波纹管结构参数与输送均匀性的关系。构建以排肥均匀性变异系数为响应值的二次响应面模型,并利用响
14、应面分析法获得波纹管的最佳参数组合以提高气力排肥均匀性。最后通过台架试验对仿真模型与最优参数组合进行验证,确定仿真模型的可靠性。1基于CFD-DEM 的耦合仿真模型建立1.1几何模型建立参考戴亿政等14 的研究,波纹管几何结构如图 1-a 所示,图中 d 表示波纹管内径,A 表示幅宽,L 表示波长,S 表示波纹间距,H 表示总长。三维建模效果如图 1-b 所示。1.2Fluent 参数设置1.2.1数学模型本次研究模拟工作介质为不可压缩牛顿流体空气,忽略流体流动时黏性耗散作用,假设波纹管壁面为刚性17。连续性方程、动量方程、能量方程分别表示如下18:连续性方程:u=0。(1)动量方程:ut+u
15、()u=f p+2u。(2)能量方程:cpTt+u()T=2T+q。(3)a,波纹管几何结构图;b,三维建模效果图。a,Bellows geometric structure diagram;b,3D modeling rendering图 1波纹管结构示意图Fig 1Structural diagram of bellows291浙江农业学报第 35 卷第 1 期式中:u 为速度矢量,m s1;t 为时间,s;T 为温度,K;为密度,kgm3;f 为单位质量上的质量力,N kg1;p 为压强,Pa;为动力黏度,Pa s;cp为定压比热容,J kg1 K1;为导热系数,W m1 K1。1.2.
16、2网格划分与模拟方法利用前处理软件 Meshing 进行网格划分,由于波纹管本身曲率变化较大,对整体进行网格划分时采用 Tet/Hybrid 非结构性网格19。划分完成后通过软件自带的 Statistic 功能进行质量检验20。划 分 后 网 格 数 量 233 634,节 点 数 量80 911,平均偏度(skewness)为 0.25,正交质量(orthogonal quality)为 0.823,说明网格质量较好可以用来模拟,网格划分效果如图 2 所示。考虑到实际使用情况,将进出口适当延长以消除仿真过程中物理边界对进出口边界影响。模拟方法根据流动状态雷诺数判断,雷诺数计算公式如下:Re=
17、L。(4)式中:Re表示雷诺数;表示密度,kg m3;表示速度,ms1;L 表示特征长度,m;表示动力黏度,Pa s1。试验模拟工作介质为标准条件下空气,密度为 =1.29 kg m3,动力黏度 =1.79 105Pa s,特征长度 L=0.06 m。参考文献 21,模拟风速 =25 m s1,将空气物性参数与风速带入公式(4)可得波纹管入口处雷诺数 Re=108 136 4 000,因此本次模拟需要使用湍流模型。湍流模型采用标准 K-epsilon 模型,算法选择 Simple 二阶迎风差分格式离散。入口边界条件选择速度入口,出口边界条件选择压力出口,各项变量的收敛残差均为 106。1.3E
18、DEM 参数设置参考文献 22,仿真试验所用肥料颗粒为尿素,泊松比 0.51,弹性模量 8.9 107Pa,剪切模量 3.56 107Pa,密度 1 337 kgm3,粒径 3.77mm。在 EDEM 中接触模型均选择 Hertz-Mindlin(No Slip)碰撞接触模型,颗粒与颗粒之间碰撞恢复系数为 0.28,静摩系数 0.36,滚动摩擦系数0.15,颗粒与设备之间碰撞恢复系数为 0.36,静摩系数 0.41,滚动摩擦系数 0.0423。参考高观保24、杨庆璐等12 的研究,初步设定本研究中起始施肥速率 350 g s1。1.4耦合模型参数设置采用基于 Lagrangian 法的多相流模
19、型的耦合算法,设置动量亚松弛、体积亚松弛因子为0.3 25。由于气流速度相对于模型尺度较大,为捕捉到足够数据点,设置 Fluent 时间步长0.000 1 s,每100 步记录1 次,Edem 时间步长5 106s,记录间隔0.01 s。2单因素仿真试验2.1排肥均匀性指标建立利用虚拟网格划分技术,在波纹管靠近出口处设置网格组,该网格组以管道中心为基准均分为 6 个扇形网格,在仿真时长内对通过各个网格的颗粒数量进行累加计数,并由 EDEM 计算得出通过该网格颗粒的总质量,网格划分与计数效果如图 3 所示。再通过公式(5)、公式(6)、公式(7)计算排肥均匀性变异系数:图 2网格划分结果示意图F
20、ig 2Schematic diagram of meshing results391张一帆,等 基于 CFD-DEM 的排肥用波纹管结构优化设计与试验图 3划分网格箱体并对各个箱体内颗粒累加计数Fig 3Dividing grid box and accumulating countingthe number of particles in each boxaj=njx=1axnj;(5)Sj=njx=1(ax aj)2nj1;(6)Vj=Sjaj100%。(7)式中,aj为各网格累计颗粒质量平均值,单位 g;ax为第 x 个网格内累计颗粒质量,单位 g;Sj为各网格累计颗粒质量标准差;Vj
21、为排肥均匀性变异系数;nj为测定网格个数。2.2单因素试验设计与分析参考杨庆璐等12、戴亿政等14 的研究,模拟波纹管内径固定为60 mm,总长500 mm,以波长、波纹间距、幅宽 3 个最能反映波纹构造的因素作为试验因素,不同因素水平选取采用等距排列,每组试验重复 3 次取平均值,因素水平表如表 1所示。2.2.1幅宽对排肥均匀性的影响如图 4 所示,排肥均匀性变异系数随幅宽增加先减小后增加,在幅宽 A 6 mm 时变化幅度较大,幅宽 A 6 mm 后变化幅度较小,变异系数总体变化幅度 44%,推测幅宽处于 6 10 mm,有最小值点。综上,选取幅宽 A 为 6、8、10 mm 作为响应面试
22、验的 3 个水平。2.2.2波长对排肥均匀性的影响如图 5 所示,排肥均匀性变异系数随波长增加先减小后增加,波长 L 10 mm 时变化幅度较大,波长 L 10 mm 后变化幅度较小,变异系数总表 1单因素试验因素水平表Table 1Single factor test factor level tablemm水平Level因素 Factor波长Wavelength波纹间距Ripple spacing幅宽Width of cloth100225543101064151585202010图 4不同幅宽下排肥均匀性变异系数变化Fig 4Variation coefficient of fertil
23、izer uniformity underdifferent width of cloth图 5不同波长下排肥均匀性变异系数变化Fig 5Variation coefficient of fertilizer uniformity underdifferent wavelengths491浙江农业学报第 35 卷第 1 期体变化幅度 25%,推测波长处于 10 20 mm,有最小值点。综上,选取波长 L 分别为 10、15、20mm 作为响应面试验的 3 个水平。2.2.3波纹间距对排肥均匀性的影响如图 6 所示,排肥均匀性变异系数随波纹间距增加先平缓减少再缓慢增加,变异系数总体变化幅度 5%
24、。综上,选取波纹间距 S 分别为 0、10、20 mm 作为响应面试验的 3 个水平。图 6不同波纹间距下排肥均匀性变异系数变化Fig 6Variation coefficient of fertilizer uniformity underdifferent ripple spacing3基于 BDD 法的响应面优化试验3.1响应面模型建立为进一步研究不同波纹管结构参数之间交互关系对输送均匀性能的影响,进行 BDD 法响应面分析。选取幅宽、波纹间距、波长作为影响因素,排肥均匀性变异系数作为响应值。参考上一节分析结论,实验组因素编码及水平如表 2 所示。利用 Design Expert 10
25、软件计算仿真数据,得到响应面试验结果如表 3 所示。表 2BDD 法实验组因素编码及水平Table 2Factors and levels of experiment based onBDD method编码Code影响因素Influence factor水平 Level10+1A幅宽 Width of cloth/mm6810B波纹间距 Ripple spacing/mm01020C波长 Wavelength/mm101520表 3实验组 BDD 法响应面试验设计与结果Table 3Orthogonal test table based on BDD method序号Serialnumber
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 CFD DEM 排肥用 波纹管 结构 优化 设计 试验 张一帆
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。