基于模型的直升机航空施药飘移预测方法技术

本发明专利技术涉及一种基于模型的直升机航空施药飘移预测方法，属于农用航空施药技术领域。该方法模拟农用直升机和喷嘴的雾化面，确定预测域和预测域的边界条件，并运行Fluent，得出施药直升机行进及旋翼产生的气流受到机身阻挡后在雾化面形成的风场，通过迭代运算，模拟边界条件下预测域各网格内的药液浓度和药液密度变化，获得药液雾滴的沉积分布，完成直升机的航空施药飘移预测。本发明专利技术将CFD思想引入到农用航空施药技术领域，可以快速有效地模拟即将进行的航空施药药液沉积分布情况及飘移范围等，从而为准确、安全、便捷评估农业航空施药作业创造条件。

【技术实现步骤摘要】

本专利技术涉及一种航空施药飘移预测方法，属于农用航空施药

技术介绍
农业航空在上世纪得到了长足的发展，而随着人们对生态环境的越来越重视，进入本世纪以来，不再一味地追求航空施药的高效性，而同时考虑到航空施药的安全性和对环境的污染问题，而航空施药的农药飘移是最为突出的问题。通过对药液雾滴的收集，检测航空施药的农药飘移量，确定安全施药带，在农业航空施药的时候能够有效地避免对非施药区域的环境污染以及人身安全的危害。传统方法是在试验或作业场地空中，设置方向与飞行方向垂直架设聚乙烯线(筒形)、地面按照一定的行间距大量放置采样卡或水敏纸等，进行药液雾滴或示踪剂的收集，再通过专业设备仪器·如荧光光度计等对所收集的药液或示踪剂进行沉积量分析，从而确定药液雾滴的飘移距离和浓度，这样的方法简单有效，但是需要消耗大量的人力物力，且设备价格昂贵，操作专业性强。另外，上述方法都仅局限于在航空施药作业之后，对药液雾滴的沉积分布量数据进行采集，再进一步分析农药飘移范围，工作量大，且缺乏时效性和预测性。目前亟需找到一种操作简易、成本低廉，而且准确度较高的方法，能够在航空施药作业之前对药液雾滴飘移进行预测，评估农业航空施药作业的安全性。我国农业航空发展正处于起步且蓬勃发展阶段，在航空施药安全意识和评估手段上与国外发达国家相比差距很大。据申请人了解，我国国内至今没有见到相关航空施药飘移模型的文献。因此找到一种适合我国国情的农用航空施药飘移预测方法十分重要。
技术实现思路
本专利技术解决的技术问题是提出一种，该方法可以快速有效地模拟即将进行的航空施药药液沉积分布情况及飘移范围等，从而为准确、安全、便捷评估农业航空施药作业创造条件。目前，美国 ANSYS 集团的 CFD 商业软件 FLUENT ( Aavid Thermal Technologies,Inc.开发)用来模拟从不可压缩到高度可压缩范围内的复杂流动。由于采用了多种求解方法和多重网格加速收敛技术，因而FLUENT能达到足够的收敛速度和求解精度。灵活的非结构化网格和基于解的自适应网格技术及成熟的物理模型，使FLUENT在转换与湍流、传热与相变、化学反应与燃烧、多相流、旋转机械、动/变形网格、噪声、材料加工、燃料电池等方面已得到广泛应用。为了解决上述技术问题，本专利技术提出的技术方案是一种，在载有Fluent软件包的计算机中，按以下步骤进行预测 第一步、模拟直升机一实测施药直升机机身曲面参数，输入计算机，进行网格化处理,生成机身网格 第二步、模拟喷嘴雾化面一根据安装在施药直升机机身上的喷嘴个数以及喷嘴之间的距离，按下式确定喷嘴的喷量L=I S = I^vsi 上式中L为喷量(单位为升)，I.为每亩所需喷量(单位为升/亩)，S为作业面积(单位为亩)V是直升机作业飞行速度(单位为米/秒)，s是直升机作业有效喷幅(单位为米)， 是时间(早位为秒)； 并测得雾滴粒径的尺寸，然后将喷嘴的喷量和雾滴粒径的尺寸输入计算机，生成预定分隔精度的喷嘴雾化面网格 第三步、确定预测域——以施药直升机旋翼轴中心线与机身上平面交点为基准点，以距所述上平面5±0. 5和8±0. 5倍旋翼直径的平面分别为上边界和下边界，以距所述中心线15± I倍的旋翼直径圆周面为周向边界，建立圆柱形预测域，并生成预定精度的预测域网格 第四步、输入预测域边界条件一将实测外界风速、温度、压力以及喷嘴内的药液浓度和药液密度输入计算机，作为边界条件； 第五步、预测航空施药飘移——运行Fluent，得出施药直升机行进及旋翼产生的气流受到机身阻挡后在雾化面形成的风场，通过迭代运算(例如高斯-赛德尔迭代法)，模拟边界条件下预测域各网格内的药液浓度和药液密度变化，获得药液雾滴的沉积分布，完成直升机的航空施药飘移预测。在本专利技术中，喷嘴的开口底面为喷射面，雾化面为垂直于喷射轴线的一定距离上雾滴组成的平面。本专利技术创造性的将CFD (计算流体力学)思想引入到农用航空施药
，对实际的直升机尺寸进行线性化和拟合，结合实时的作业参数，实现了对航空施药作业的药液飘移情况的预测。 对上述技术方案的改进是在第一步、第二步及第三步中均使用贴体坐标法生成机身、雾化面以及预测域的网格图。对上述技术方案的进一步改进是在第三步中确定预测域的方法为，以直升机旋翼轴中心线与机身上平面交点为基准点，所述预测域的上边界距直升机5倍旋翼直径距离，预测域的下边界面距直升机8倍旋翼直径距离，预测域的周向边界面距离直升机15倍旋翼直径距离。预测域选取的目的是取到受影响的流场区域，预测域选取如果过小，则仿真出的最终结果不准确，不能够完全体现药液雾滴飘移的距离；预测域选取如果过大，会造成运算的时间过长，对资源造成浪费。对上述技术方案的再进一步改进是在第五步中采用高斯-赛德尔算法对预测域进行迭代计算，在网格化后的预测域上使用有限体积法将所述控制方程离散成线性方程组，并求解离散得到的线性方程组。对上述技术方案的更进一步改进是在执行第一步前，对实测的直升机的机身生成曲面，具体方法如下 A、对机身用非均匀有理B样条(NURBS)曲线进行过度； B、对应边界进行相容性处理； C、生成NURBS的u向和V向直纹面；D、生成基于NURBS的张量积曲面； E、NURBS边界曲面生成。本专利技术的有益效果是本专利技术根据我国旋翼机的使用现状，引入计算流体力学(CFD)方法，对直升机进行了模型化处理，并引入实际作业的参数和条件，构建了直升机施药模型，通过结合实际作业参数和环境因子进行CFD计算，对整个作业过程动态模拟，来分析复杂的农业航空施药中药液雾滴的沉积、飘移规律，准确地预测了航空施药药液的飘移范围。本专利技术将飞机尾流、翼尖涡流、直升机旋翼下旋气流和机身周边空气扰动纳入到对雾滴的影响因素，通过模拟仿真的手段来预测雾滴的运动和地面沉积情况，从而实现低成本、快速、有效地预测出即将作业的效果。附图说明下面结合附图对本专利技术的基于模型的作进一步说明。·图I是本专利技术实施例的直升机机身的网格图。图2是本专利技术实施例的雾化面网格图。图3是本专利技术实施例的预测域网格化图。图4是本专利技术实施例预测出来的空气中等浓度线分布图。图5是本专利技术实施例预测出来的地面上等浓度线分布图。具体实施例方式本实施例中直升机的翼型剖面段长度L=1200mm ;桨叶弦长C=140mm ;桨叶安装角(r=0. 2m处)amin=2. 4。，amax=12. I。，总矩可调范围0%—99%，工作总矩为50% ;翼型剖面段翼根段扭转角为9. 05°，翼型剖面段翼尖段扭转角为3.05°。本实施例的，包括以下步骤 第一步、模拟直升机一实测施药直升机机身曲面参数，输入计算机，进行网格化处理,生成机身网格图。本实施例中首先对本专利技术在执行第一步前对实测的直升机的机身生成曲面，具体方法如下 A、对机身用非均匀有理B样条(NURBS)曲线进行过度； B、对应边界进行相容性处理； C、生成NURBS的u向和V向直纹面； D、生成基于NURBS的张量积曲面； E、NURBS边界曲面生成。本实施例使用贴体坐标法生成直升机的机身曲面的网格图,直升机机身的网格图如图I所示。贴体坐标法的具体算法可以参考任玉新、陈海昕编著的《计算流体力学基础》，清华大学出版本文档来自技高网...

【技术保护点】
一种基于模型的直升机航空施药飘移预测方法，在载有Fluent软件包的计算机中，按以下步骤进行预测：第一步、模拟直升机——实测施药直升机机身曲面参数，输入计算机，进行网格化处理，生成机身网格图；第二步、模拟喷嘴雾化面——根据安装在施药直升机机身上的喷嘴个数以及喷嘴之间的距离，按下式确定喷嘴的喷量：上式中为喷量（单位为升），为每亩所需喷量（单位为升/亩），为作业面积（单位为亩），是直升机作业飞行速度（单位为米/秒），是直升机作业有效喷幅（单位为米），是时间（单位为秒）；并测得雾滴粒径的尺寸，然后将喷嘴的喷量和雾滴粒径的尺寸输入计算机，生成预定分隔精度的喷嘴雾化面网格图；第三步、确定预测域——以施药直升机旋翼轴中心线与机身上平面交点为基准点，以距所述上平面5±0.5和8±0.5倍旋翼直径的平面分别为上边界和下边界，以距所述中心线15±1倍的旋翼直径圆周面为周向边界，建立圆柱形预测域，并生成预定精度的预测域网格图；第四步、输入预测域边界条件——将实测外界风速、温度、压力以及喷嘴内的药液浓度和药液密度输入计算机，作为边界条件；第五步、预测航空施药飘移——运行Fluent，得出施药直升机行进及旋翼产生的气流受到机身阻挡后在雾化面形成的风场，通过迭代运算，模拟边界条件下预测域各网格内的药液浓度和药液密度变化，获得药液雾滴的沉积分布，完成直升机的航空施药飘移预测。2012104110158100001dest_path_image001.jpg,11761dest_path_image002.jpg,928901dest_path_image004.jpg,2012104110158100001dest_path_image005.jpg,646322dest_path_image006.jpg,2012104110158100001dest_path_image007.jpg,2012104110158100001dest_path_image009.jpg...

【技术特征摘要】

【专利技术属性】
技术研发人员：薛新宇，张宋超，孙竹，常春，梁建，秦维彩，周立新，孔伟，周良富，蔡晨，王宝坤，
申请(专利权)人：农业部南京农业机械化研究所，
类型：发明
国别省市：