第36卷第22期农业工程学报 V ol.36 No.22
158 2020年11月Transactions of the Chinese Society of Agricultural Engineering Nov. 2020
基于MATLAB与COMSOL联合仿真的梯形迷宫滴头 胡宇祥，彭军志，殷飞，刘喜峰，李娜
精氨酸
（吉林农业科技学院，吉林 132101）
摘要：传统的滴灌灌水器流道优化设计方法是通过正交试验设计，进行有限元分析，并配合试验验证，试验过程强度大、造价高。该研究将智能优化算法与有限元分析相结合，联合MATLAB与COMSOL仿真软件对滴灌灌水器流道参数优化设计。首先，通过MATLAB读取梯形迷宫流道灌水器在COMSOL软件中的仿真计算结果后，传递给遗传算法，并以流态指数为优化目标，流道齿参差值、齿高、单元数、流道转角、齿间距为设计变量，求解出常压（100 kPa工作压力）下的灌水器流道参数。利用SPSS软件回归分析法，建立了5个变量与流态指数之间的多元线性回归数学模型。分析结果表明，5个变量对流态指数影响程度排序为齿参差值、齿高、单元数、流道转角、齿间距。进一步分析灌水器流道内水力特性可知，滴头流道内低流速区域流速达到0.1～1.2 m/s，过流面积占流道截面的75.1%，有堵塞风险，高流速 区域流速达到2.8～3.8 m/s，中流速区域流速达到1.5～2.5 m/s，流道内压力沿流道长度呈线性变化。实物样机测试结果表明，滴头流量模拟计算值与实测值之间的平均误差为6.1%，优化结果精度高，可为梯形迷宫流道灌水器流道参数优化设计提供参考依据。
关键词：优化；设计；滴灌灌水器；结构参数；COMSOL；MATLAB
doi：10.11975/j.issn.1002-6819.2020.22.017
中图分类号：S275.6 文献标志码：A 文章编号：1002-6819(2020)-22-0158-07
胡宇祥，彭军志，殷飞，等. 基于MATLAB与COMSOL联合仿真的梯形迷宫滴头流道优化[J]. 农业工程学报，2020，36(22)：158-164. doi：10.11975/j.issn.1002-6819.2020.22.017
Hu Yuxiang, Peng Junzhi, Yin Fei, et al. Optimization of trapezoidal labyrinth emitter channel based on MATLAB and COMSOL co-simulation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(22): 158-164. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2020.22.017
0 引 言
滴灌灌水器是滴灌系统中最为关键的元件，其使用过程中蒸发损失小，不产生地面径流，几乎没有深层渗漏。随着农业用水形势逐年严峻以及国家“节水增粮行动”项目的大力推广，滴灌技术成为最有效的灌溉方式，是保障增产增收的重要手段[1]。滴灌灌水器结构尺寸对滴灌系统灌水的均匀性、抗堵塞能力影响显著。过大的流道结构参数将造成灌溉均匀度下降，过小的流道结构参数将造成灌水器流道堵塞，严重时会使整个系统无法正常工作[2-3]。因此，研究流道结构参数与滴头水力性能之间的关系一直是节水灌溉领域的研究热点。
目前，国内外学者已在滴灌灌水器流道结构参数优化方面取得了一系列的研究成果，颜廷熠等利用正交试验法，研究了流道结构参数对滴头水力性能的影响[4-6]；王文娥利用FLUENT 软件对内镶齿形片状滴头的水力性能进行了数值模拟，获得了滴头压力-流量关系式并分析了滴头内速度、压力分布[7]；李永欣等建立了迷宫滴头的CFD数值模型，并对滴头的压力流量关系、流道内部
收稿日期：2020-08-16 修订日期：2020-11-06
基金项目：吉林省教育厅“十三五”科学研究规划项目（JJKH20180739KJ）作者简介：胡宇祥，博士。Email：*********************的压力和流速分布进行了数值模拟计算，利用原型滴头和滴头放大模型实测值对模型和模拟计算结果进行了试验验证[8]；陈爱平利用 GAMBIT软件建立了输送管的仿真模型，通过 FLUENT 软件的数值模拟计算，得到了不同滴灌装置输送管结构参数下的流场分布特性[9]；
Souza 等通过试验方式，设计了一种灌溉滴头，测定了滴头流量系数和流态指数[10]；Li等利用多元回归分析，建立了流量与流道几何参数之间的线性回归方程[11]；Wei等利用CFD对滴头水利特性进行了仿真研究[12]。目前，上述研究多是利用正交试验方法，将流道结构参数设计成多因素多水平的试验，通过CFD 方法或者试验法获得滴头的宏观水力特性，即压力与流量间的关系，分析流道结构参数与流态指数之间关系[13]。但是通过正交试验设计进行数据分析后得到的优选结果，只是因素间不同水平的相互组合，难于确定数据变化规律。水平设计的少，则自变量定义域不完整，试验设计具有随机性，结果不够精确[14-15]；如果水平设计的多，则试验数据体量大，试验次数多，尤其是费用昂贵的试验，弊端更加凸显。因此，需要一种能够快速准确求解灌水器流道尺寸的计算方法。
鉴于此，本文提出联合MATLAB与COMSOL仿真软件进行滴灌灌水器流道参数优化设计，采用COMSOL 软件对梯形迷宫流道灌水器进行仿真，求解最佳流态指
第22期 胡宇祥等：基于MATLAB 与COMSOL 联合仿真的梯形迷宫滴头流道优化
159
数下的滴灌灌水器结构参数，研究流道齿参差值、齿高、单元数、流道转角、齿间距这5个变量对灌水器流道参数的影响，并在此基础上建立数学预测模型，以期为快速确定灌水器流道结构参数提供理论依据以及科学参考。
1 研究方法
1.1 滴灌灌水器流道数学模型 1.1.1 设计变量
迷宫形滴头相比贴片形式滴头抗堵性能好，灌水均匀，滴头整体性强。Netafim 公司作为滴灌技术的发明者，一直在引领滴灌领域技术变革。本文以该公司迷宫滴头流道参数的取值范围作为参考依据，其流道结构如图1所示。
注：l 为齿间距，mm ；w 为流道宽度，mm ；θ 为流道转角，（°）；h 为齿高，mm ；j 为齿参差值，mm 。
Note:l is inter-tooth space, mm; w is width of labyrinth, mm; θ is angle of labyrinth, (°); h is height of tooth, mm; j is tooth diversification index, mm.
图1 灌水器流道结构及参数示意图
Fig.1 Schematic and parameter diagram of emitter channel
本文选取梯形迷宫流道滴头进行研究，影响其水力性能的关键结构参数包括流道宽度 w 、流道转角 θ
、齿高 h 、齿间距 l 、齿参差值j ，灌水器流道单元数n 。θ、h 、l 、j 、n 这5个参数可以确定梯形齿道单元平面形状，流道宽度 w 按式（1）计算[16-17]。
22sin(arctan 2cot w h j
l h θθ
=+⨯-+ （1） 算例中，流道转角 θ、齿高 h 、齿间距 l 、齿参差值j 和流道单元数n 是设计变量。根据已有的研究文献，各变量取值范围如表1 所示[13]。 表1 灌水器的结构参数取值范围
Table 1 Range of structural parameters of emitter
参数Parameters
值Values
流道转角Angle of labyrinth θ/(°) 120～140 齿高Height of tooth h/mm 0.5～1.5 齿间距Inter-tooth space l/mm 1～2 流道单元数Units number of labyrinth n 10～20 齿参差值Tooth diversification index j/mm
-~22
h h
1.1.2 目标函数
滴灌滴头的流态指数x 表示流量对压力发生改变的敏感程度，流态指数越小，滴头的水力性能越好[18]。x =1表示流量变化与压力以相同的百分率变化，x =0表
示流量不随压力变化而变化。压力补偿式灌水器，流态指数趋近于0；非压力补偿式灌水器，流态指数介于0~1之间。非压力补偿式灌水器x 值越趋近于0.5，流量对压力敏感程度就越小，水力性能越优异[19]。因此，建立目标函数如式（2）所示，要求非压力补偿式灌水器流态指数趋近于0.5，保证灌水器水力性能，提高出流均匀性。 F =min|x -0.5| （2） 式中F 为目标函数；x 为流态指数，按照SL/T67.1—94要求，可由流量与压力关系[20]计算得出： Q =K d H x （3） 式中Q 为滴头流量，L/h ；H 为工作压力， kPa ；x 为流
近浅海观测网态指数；K d 为流量系数。
从0.5～1.5倍工作压力（100 kPa ）范围内，均匀选取9个工作压力值作为自变量，输入COMSOL 软件作为滴头进口压力，通过仿真计算出滴头出口流量Q 作为因变量，利用MATLAB 进行最小二乘法计
算，建立基于乘幂函数形式的流量与压力关系，得到回归系数K d 和x 。 1.1.3 边界条件
滴头流道进水口设为工作压力，出口设为当地大气压，计算出口流速。初始状态设置为滴头内充满水，滴头以外为空气，空气压力为大气压[7]。本文采用精度高、计量小的标准壁面函数，减小梯形滴头流道的壁面对紊流的影响[21-22]。 1.2 遗传算法原理
与传统的优化算法相比，遗传算法生成初始体开始搜索，按照基因的表达方式，目标函数值转换得到适应数值的信息，而不需要借助其他辅助信息以及参数。遗传算法有3个基本操作：选择算子、交叉算子和变异算子，各个算子之间即相互配合又相互竞争[23]。遗传算法作为一种全局优化算法，相对于粒子算法、蚁算法和鱼算法等智能算法，其基本原理简单易懂，初始条件少，广泛应用于不同领域。本文利用MATLAB 自带的遗传算法工具箱进行求解。
1.3 COMSOL 与MATLAB 联合仿真在滴灌灌水器流道优化中应用
COMSOL （COMSOL Multiphysics ）是多物理场仿真软件，可用于温度、电磁、流体等多物理场之间的相互耦合用。 COMSOL 中包含livelink for MATLAB 模块，能够实现与MATLAB 软件之间数据的相互传递、相互交流，这样可以利用 MATLAB 及其工具箱进行预处理、模型设定和后处理，进行多物理场模型仿真和遗传算法优化[24]。
COMSOL 与MATLAB 联合仿真优化滴头流道结构流程是，首先运行MATLAB 遗传算法，生成设计变量、设定遗传优化的主要参数，然后在COMSOL 中完成滴头流道结构的参数化建模与计算，最后使用MATLAB 读取COMSOL 仿真结果，计算目标函数，完成遗传算法寻优化。COMSOL 与MATLAB 两者之间流转交换步骤如下：
农业工程学报（ ） 2020年
160
1）运行MATLAB 遗传算法，生成流道转角 θ、齿高h 、齿间距 l 、齿参差值j 和流道单元数n 设计变量。
2）将设计变量和已知变量输入COMSOL ，进行滴头流道结构的参数化建模与仿真计算。
3）将仿真计算结果传给MATLAB ，拟合出Q ～H 公式。
4）设定遗传算法中选择、交叉、变异参数，开始计算目标函数值。
5）判断是否满足收敛条件，如果是，转向步骤6）；否则，转向步骤1），生成下一代设计变量。
6）求出最优解，计算结束。 1.4 线性回归分析原理
回归分析（Regression analysis ）是确定2种或2种以上变量间相互依赖定量关系的一种统计分析方法。线性回归分析是最为常用的方法之一，其目的是在因变量和一个或多个自变量之间建立一种关系。该方法可用于因素分析，统计预测，调整因素等，通过求解标准化回归系数，可以分析自变量对因变量的影响程度[25]。求解过程包括多重共线性检验和计算标准化回归系数。
方差膨胀系数（Variance Inflation Factor ，VIF ）是衡量多元线性回归模型中多重共线性严重程度的一种度量。它表示回归系数估计量的方差与假设自变量间不线性相关时方差相比的比值。当出现严重共线性问题时，会导致分析结果不准确，因而需要及时进行处理。VIF 计算如式（4）所示。
VIF 2i 11-R = （4） 式中VIF 为方差膨胀系数；R i 为自变量X i 对其余自变量作回归分析的负相关系数。VIF>5说明存在共线性问题，VIF>10说明存在严重的多重共线性问题，模型构建效果不理想，需要进行处理。
当分析通过多重共线性检验后，需要计算标准化回归系数，用于比较自变量对因变量的影响程度。标准化回归系数绝对值越大说明该自变量对因变量的影响越大。若回归方程的形式如式（5）所示。 Y =b 0+b 1X 1+...+b j X j +...+b m X m （5） 式中Y 是估计值；m 为变量个数，取5；多项式系数b j 是非标准化系数，通过最小二乘法求得。标准化回归系数为
Beta =b j j
Y σσ （6）
式中Beta 为标准化回归系数；σj 为X j 的标准差；σY 为Y 的标准差。 1.5 数值模拟
本文利用COMSOL 软件建立梯形迷宫滴头流道模型，基于标准 k -ε 湍流模型进行仿真，计算中对梯
形滴头流道划分网格。梯形流道结构壁薄且有较多转角，在对流道划分网格时，这些结构对网格质量的影响十分显著。为了避免转角结构给网格划分带来困难，采用平面
计算模型来研究流道的水力特性，进行平面模型网格划分[22]。采用的网格大小为0.01～l /20 mm （l 为齿间距），因为过密划分网格对于提高计算精度并不明显，还会增加计算时间。本文采用渐变方式划分网格，对于规则、流动变化缓慢的区域，采用较稀疏的网格；对于不规则、流动变化大的区域采用细密网格，疏密网格间平滑过渡，滴头进口设置为100 kPa 工作压力。
图2 迷宫流道的网格划分
Fig.2 Grid generation of labyrinth channel
1.6 试验方法
为了验证联合仿真计算方法的可行性，评价梯形迷宫滴头流道灌水器实际的工作性能，利用CNC 精雕机制作亚克力材料滴头测试样本。按照SL/T67.1—94规范要求，在0.5～1.5倍工作压力（100 kPa ）范围内，排列25个滴头按流量由小到大排列编号，取第3、12、13、23号4个灌水器为试样[26-27]。由小到大调节滴头的工作压力，在测试压力范围内均匀分布9个压力点，用量筒测量4个灌水器每一个压力点的出水量，时间不少于3 min ，试验重复 3 次，3次所测水量之差不得大于2%，取平均值计算
流量，进而得到试样灌水器的压力与流量关系曲线。将试样灌水器每一组压力下的平均流量进行回归，绘制灌水器压力-流量关系曲线。
2 结果与讨论
2.1 优化结果与验证
在COMSOL 建模基础上，利用MATLAB 遗传算法工具箱进行优化，开展梯形迷宫滴头流道优化联合仿真计算。遗传算法参数选择如下：初始种规模为50，遗传终止代数MAXGEN= 50，选择算子采用赌方式，选择概率为1/3；交叉算子采用两点交叉算子，交叉概率为0.75；变异算子采用高斯变异，变异概率为0.02。灌水器的结构参数优化结果如表2所示。
表2 灌水器结构参数优化结果
Table 2 Optimization results of emitter structure parameters
参数Parameters
值Values 流道转角Angle of labyrinth θ/(°) 125 齿高 Height of tooth h/mm 1.81 齿间距 Inter-tooth space l/mm 1.11 流道单元数 Units number of labyrinth n 13 齿参差值 Tooth diversification index j/mm
0.3
根据表2，通过回归分析得到灌水器流量与压力关系
为Q =0.466 9H 0.499 9
，相关系数R 2为0.99，拟合较好。根据试验结果绘制压力-流量关系曲线如图3所示，由实测数据得到滴头试样流量与压力的函数关系为
第22期胡宇祥等：基于MATLAB与COMSOL联合仿真的梯形迷宫滴头流道优化161
Q=0.332 3H0.557 6，R2 = 0.968。对比仿真优化结果，在常压工作压力下，实测平均流量普遍偏小。利用SPSS软件，对仿真优化结果与试样实测结果进行相关性分析，相关性系数为0.994，根据显著性检验方法计算得到P 值为0.000 1，小于0.001，表明结果有较强相关性。根据GB/T 17187—2009规定，滴头流量模拟计算值与实测值之间的平均误差为6.1%，小于7%，设计的滴头满足国标要求[28-29]。
图3 优化结果与实测结果压力-流量曲线对比Fig.3 Comparison of pressure - flow curve between optimized
results and measured results
2.2 滴头流道流速
在100 kPa工作条件下，优化参数的梯形迷宫流道灌水器流道的流速矢量分布和速度局部放大图如图4所示。从滴头流道流速分布图可知，各单元之间流速分布一致。根据滴头流道速度局部放大图，可以把梯形滴头流道内部分为3个区域：低流速区域A，处于齿道转角处，速度0.1～1.2 m/s；高流速区域B，处于轴心位置，速度2.8～3.8 m/s；中流速区域C，处于A、B区域之间的过渡区域，速度1.5～2.5 m/s。由低流速区域形状可知，梯形齿尖已经趋近于三角形。基于梯形迷宫滴头流道建立优化模型，得到趋于三角形迷宫滴头流道优化结果，可见三角形滴头相比梯形滴头优势明显[18,30]。
注：A表示低流速区域，B表示高流速区域，C表示中流速区域。
Note: A represents low velocity area, B represents high velocity region, and C represents medium velocity region.
图4 滴头流道流速分布
Fig.4 Velocity distribution in emitter channel
从流道内的流速分布图可知，流道进口段流速较低，是长流道滴头容易堵塞位置，与试验研究较为符合[31]。由图4可知，齿道转角部位A区域流速近似为0，该区域流速慢，颗粒易发生沉积，从定性角度来说，更易发生堵塞。以单个齿道为例，A区域面积占截面面积75.1%，从灌水器流道的结构设计出发，下一步可将滴头流道水力性能与抗堵性能作为优化目标，将流道形式设计改为圆弧形边壁，减小
低流速区域面积，从而提高灌水器的抗堵塞性能，优化流道设计。
2.3 滴头流道压力
最优方案的梯形迷宫流道灌水器流道压力分布如图 5所示，由于整个流道为全封闭状态，没有与外界的能量交换，流道内部流体流动过程中产生了沿程水头损失和局部水头损失，压力由进口到出口逐级减小，从100 kPa下降为0，13个单元中每个单元段水头损失大约为7.7 kPa，压降基本相同，压力沿流道长度呈线性变化。
健康无忧网图5 滴头流道压力分布
Fig.5 Pressure distribution in emitter channel
2.4 流道结构参数对滴头水力性能的影响
为分析流道转角、齿高、齿间距、单元数、齿参差值5个设计变量对流态指数x的影响，以仿真优化过程中历代遗传算法计算结果数据为基础，利用SPSS软件进行线性回归分析，结果如表3所示。
表3 流道结构参数逐步回归分析结果
Table 3 Results of stepwise regression analysis of flow channel
structure parameters
设计变量
Design variables
非标准化系数
Coefficient of
non-
standardization
标准化回归
系数
Standardization
coefficient
容差
Tolerance
方差膨胀
系数
Variance
Inflation
Factor VIF 常数项0.479 2 - - -
流道转角
Angle of labyrinth θ/(°)
0.000 3 0.109 0.977 1.023
齿间距
Inter-tooth space l/mm
0.009 2 0.089 0.998 1.002
齿高
Height of tooth h/mm
-0.034 8 -0.459 0.926 1.080 齿参差值
Tooth diversification
index j/mm
0.051 5 0.596 0.395 2.531
单元数
Units number of
大象洗完澡为什么要往身上弄泥巴labyrinth n
-0.001 0 -0.326 0.407 2.458 由表3得到5个变量与流态指数之间的多元线性回归数学模型为
x=0.479 2+0.000 3θ+0.009 2l-0.034 8h
+0.051 5j-0.001n （7）从表3可知，最大的VIF为2.531<5，说明5个变量不存在共线性问题，说明此模型是可靠的。该线形模型的相关系数R 2=0.914，显著性P值为0.001，说明5个变量与流态指数之间线性关系很显著。通过比较5个设计变量标准化回归系数，得到设计变量对流态指数x影响程度由大到小的排序为：齿参差值、齿高、单元数、流道转角、齿间距。齿参差值、齿高2个参数对流态系数大小起主要影响作用。齿参差值、流道转角、齿间距3个参数标准化回归系数（Beta）为正值，说明增加齿道的齿参差值、流道转角大小可以增大流态指数；齿高、单元数2个参数标准化回归系数（Beta）为负值，说明减少齿高、单元数可以增大流态指数。由分析结果可知，
农业工程学报（） 2020年162
100 kPa工作条件下，流道结构参数中的齿参差值对滴头的水力性能影响较大。
3 结 论
本文通过MATLAB与COMSOL联合仿真计算求解出常压下最优梯形迷宫流道灌水器流道结构参数，流道内低流速区域流速达到0.1～1.2 m/s，过流面积占流道截面的75.1%，有堵塞风险；高流速区域流速达到 2.8～3.8 m/s，中流速区域流速达到1.5～2.5 m/s，流道内压力沿流道长度呈线性变化；仿真结果与实测结果的流量误差为6.1%；得到5个自变量对因变量影响程度由大到小的排序为：齿参差值、齿高、单元数、流道转角、齿间距，建立了流态指数与5个参数之间的多元线性回归数学模型，具有较好的代表性。
本文研究了常压（100 kPa）工作条件下的流道结构尺寸，未涉及低压、高压条件，研究结论有一定的局限性。以流态指数作为判断滴头水力性能的指标外，还需要对滴头抗堵塞性能做进一步研究。受滴灌系统工程造价、作物品种等因素影响，应该开展不同工作压力、不同流量条件下，灌水器流道结构参数优化研究，为开发新产品、推广应用提供科学依据。
[参 考 文 献]
[1] 寇宝峰，丁林. 低压灌溉滴头流量及毛管铺设对玉米
生长及产量的影响[J]. 甘肃农业科技，2020，533(6)：
29-33.
Kou Baofeng, Ding Lin. Effects of drop flow rate and capillary laying on growth and yield of corn under low pressure irrigation[J]. Gansu Agricultural Science and Technology, 2020, 533(6): 29-33. (in Chinese with English
abstract)
[2] 杨彬，吴勇，王建东，等. 齿型迷宫灌水器抗堵塞性能分
析与结构优化模拟[J]. 灌溉排水学报，2019，38(8)：84-89，
128.
Yang Bin, Wu Yong, Wang Jiandong, et al. Anti-clogging
performance analysis and structure optimization simulation
for dental labyrinth emitter[J]. Journal of Irrigation and Drainage, 2019, 38(8): 84-89, 128. (in Chinese with English
abstract)
[3] 王浩翔，张新燕，牛文全，等. 低压滴灌施肥条件下温度
对滴头堵塞的影响[J]. 灌溉排水学报，2020，39(3)：63-71.
Wang Haoxiang, Zhang Xinyan, Niu Wenquan, et al. The
effect of temperature on emitter clogging in low-pressure drip
fertilization system, the effect of temperature on emitter clogging in low-pressure drip fertilization system[J]. Journal
of Irrigation and Drainage, 2020, 39(3): 63-71. (in Chinese
with English abstract)
[4] 颜廷熠，周芸，仵峰. 迷宫滴头CFD模拟精度影响因素的
正交试验分析[J]. 节水灌溉，2009(6)：25-28.巴黎公约
Yan Tingyi, Zhong Yun, Wu Feng. Orthogonal experimental
analysis of influencing factors of cfd simulation precision of
labyrinth emitter[J]. Water Saving Irrigation, 2009(6): 25-28.
(in Chinese with English abstract)
[5] 王建东，龚时宏，李光永，等. 低压下流道结构参数对锯
齿型滴头水力性能影响的试验研究[J]. 水利学报，2014，
45(1)：72-78.
Wang Jiandong, Gong Shihong, Li Guangyong, et al. The
influence of geometrical parameters of dental flow passage of
labyrinth emitter on the hydraulic performance under low
working pressure[J]. Journal of Hydraulic Engineering, 2014,
45(1): 72-78. (in Chinese with English abstract)
[6] 王建东，李光永，邱象玉，等. 流道结构形式对滴头水力
性能影响的试验研究[J]. 农业工程学报，2005(增刊1)：
100-103.
Wang Jiandong, Li Guangyong, Qiu Xiangyu. Effects of flow
passage forms on hydraulic performance of emitters[J].
Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2005(Supp.1):
100-103. (in Chinese with English abstract)
[7] 王文娥，王福军. 迷宫滴头水力特性非定常数值模拟研究[J].
水利学报，2010，41(3)：332-337.
Wang Wener, Wang Fujun. Numerical simulation of
unsteady flow in labyrinth emitters of drip irrigation
system[J]. Journal of Hydraulic Engineering, 2010, 41(3):
332-337. (in Chinese with English abstract)
[8] 李永欣，李光永，邱象玉，等. 迷宫滴头水力特性的计算
流体动力学模拟[J]. 农业工程学报，2005，21(3)：12-16.
Li Yongxin, Li Guangyong, Qiu Xiangyu, et al. Modeling of
hydraulic characteristics through labyrinth emitter in dripirrigation using computational fluid dynamics[J].
Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2005, 21(3): 12-16.
(in Chinese with English abstract)
[9] 陈爱平. 农田滴灌装置输送管流场稳定性和结构优化设计[J].
农机化研究，2017(9)：75-79.
Chen Aiping. The flow field stability and structure optimization of drip irrigation device[J]. Journal of Agricultural Mechanization Research, 2017(9): 75-79. (in
Chinese with English abstract)
[10] Souza W D J , Sinobas L R , Sánchez, Raúl, et al. Prototype
emitter for use in subsurface drip irrigation: manufacturing,
hydraulic evaluation and experimental analyses[J].
小说月报2012
Biosystems Engineering, 2014(128):41-51.
[11] Li Y K, Yang P L, Xu T W, et al. Hydraulic property and
flow characteristics of three labyrinth flow paths of drip
irrigation emitters under micro-pressure[J]. Transactions of
the Asabe, 2009, 52(4): 1129-1138.
[12] Wei Qingsong, Shi Yusheng, Dong Wenchu, et al. Study on
hydraulic performance of drip emitters by computational
fluid dynamics[J]. Agricultural Water Management, 2006,
84(1): 130-136.
[13] 郭霖，白丹，程鹏，等. 滴灌灌水器三角形迷宫流道优化
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议