CN105160096A - 一种变转速喷头三维动态水量分布模型建立方法 - Google Patents
一种变转速喷头三维动态水量分布模型建立方法 Download PDFInfo
- Publication number
- CN105160096A CN105160096A CN201510543947.1A CN201510543947A CN105160096A CN 105160096 A CN105160096 A CN 105160096A CN 201510543947 A CN201510543947 A CN 201510543947A CN 105160096 A CN105160096 A CN 105160096A
- Authority
- CN
- China
- Prior art keywords
- water distribution
- shower nozzle
- factor
- establishing
- model
- 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
Landscapes
- Nozzles (AREA)
Abstract
本发明提供了一种变转速喷头三维动态水量分布模型建立方法,属于节水灌溉喷灌系统中水量分布模型的建立方法。具体为:在确定喷头的型号、喷洒试验收集方法及布置方式的情况下设计在不同影响因素下单喷头的水量分布的试验。利用单喷头的水量分布实验数据,在确定喷头组合方式和组合间距后,计算出组合喷灌均匀度。分析不同因素下组合喷灌均匀度数据,确定主要影响因素,次要影响因素及各个因子直接的互作关系。最后建立三维动态水量分布模型。该建立方法为多因素下水量分布模型的建立提供理论依据。
Description
技术领域
本发明涉及节水灌溉喷灌系统中水量分布模型的建立方法,特别是一种变转速喷头三维动态水量分布模型建立方法。
背景技术
节水灌溉不仅是我国国民经济和社会可持续发展所要求的,也是我国农业资源,尤其是水资源短缺、水土资源配置失衡等严峻形势所决定的。喷灌作为一种先进的节水灌溉技术,能提高灌溉水利用率,具有广阔的发展空间。喷灌均匀度是指在喷灌面积上水量分布的均匀程度,它是喷灌系统的重要性能,是衡量喷灌质量的重要指标和喷灌系统规划设计中的重要参数。喷洒的水量分布是决定喷洒均匀性重要因素之一,影响作物的品质,而影响喷洒水量分布的因素众多,主要包括:喷头的流量、工作压力、喷头的转动均匀性、喷洒时间、风速风向、喷嘴结构形式等。现有喷灌水量分布模型多采用线性回归的方法,在距喷头的距离与喷灌强度的二维空间建立。然而,水量分布模型应结合不同影响因素,建立相应的径向、周向与喷灌强度随时间变化的三维水量分布模型,会更接近不同类型喷头的实际喷洒情况,更准确地预测出灌水效率等参数,使模型更具有普遍适用性。因此提出一种变转速喷头三维动态水量分布模型确定方法具有重要意义。
经检索,目前还没有相关的申报专利。
发明内容
针对现有技术中存在不足,本发明提供了一种变转速喷头三维动态水量分布模型建立方法。
本发明是通过以下技术手段实现上述技术目的的。
一种变转速喷头三维动态水量分布模型建立方法,包括以下步骤,
1.一种变转速喷头三维动态水量分布模型建立方法,其特征在于,设计方法按下列步骤进行:
步骤一,确定喷头型号:
选定喷头型号,确定选定喷头的工作压力范围,以及喷头旋转一圈转速的变化次数;
步骤二,确定喷洒试验收集方法及布置方式:
使用通用的雨量筒进行水量收集,并采用径向布置的方式排列雨量筒;
步骤三,确定单喷头的水量分布的影响因素并进行试验:
通过径向分布的雨量筒法收集不同影响因素下单喷头的喷灌强度hi、雨滴直径D、雨滴累积频率f的试验数据,其中水量分布由喷灌强度hi得到,水滴分布由雨滴直径D与雨滴累积频率f得到;
步骤四,确定喷头组合方式和组合间距:
采用矩形组合,喷头与喷头间横向距离为1m,纵向距离为1m;
步骤五,利用单喷头的水量分布实验数据,计算出组合喷灌均匀度:
使用三次样条插值法处理单喷头组合的数据,然后利用克里斯琴森均匀度计算公式进行计算得到组合喷灌均匀度;
三次样条差值法的计算模型如图1所示,三次样条插值法的计算方法为:
计算喷灌区域内任意一点P的水深,通过径向插值分别计算出各辐射线上与网格点P距喷头相同远处Ai的降水深Di(i=1、2、...、6),利用同一圆周上的数据采用与径向插值方法相同的插值方法计算出所求网格点P的降水深Dp;
进行径向插值,通过Li(i=1、2、…、6)上测得的雨量筒数据dij(i=1、2、…、6;j=1、2、…、n),其中n为每条辐射线上雨量筒的个数,选用一种插值方法,分别求出Ai(i=1、2、…、6)点的降水深Di:Di=Fi(dij,r)(i=1、2、…、6;j=1、2、…、n),式中:Fi—径向插值函数;
进行周向插值,通过上一步得到的Ai点的降水深Di(i=1、2、…、6),对极角α采用与径向相同的插值方法进行周向插值,求出P点的降水深DP:D=F(Dj,α)(i=1、2、…、6),式中,F—径向插值函数;
划分区间[a,b],其中,Δ:a=x0<x1<…<xn=b,同时建立在区间[a,b]上的函数f(x),若函数分别满足下列条件:
一致通过n+1个型值点(xi,yi),即s(xi)=f(xi)=yi(i=0,1,2,...,n);
二阶导数连续,即s(x)∈c[a,b];
三次分段,即在每一个小区间[xi-1,xi]上均为三次多项式,
这样的s(x)称为[a,b]上以xi(i=0,1,2,...,n)为节点的三次插值样条函数,三次插值样条函数的几何图形称为三次样条曲线,s(x)的公式如下:
将处理后的数据代入克里斯琴森均匀度计算公式,
式中Cu—喷灌均匀系数
hi—第i个测点的降水深
h—喷洒面积上各测点平均降水深
n—测点数目。
步骤六,建立影响因素的数值关系式:
分析不同压力p,喷嘴直径d,流量系数μ,转速r、喷头仰角α,安装高度H因素下组合喷灌均匀度数据,使用非线性工具、多维尺度分析,进行非线性回归分析,依据组间均方差与组内均方差之比最大的原则来进行判别,确定主要影响因素,次要影响因素及各个因子直接的互作关系,建立各因素的线性函数X(p,r),其中:p,r是相互影响的,关系式为:F(x)=f(X(p,r),d,μ,α,h);
步骤七,建立三维动态水量分布模型:
在不同因素下,单喷头水量分布实验,通过三次样条插值法计算出矩形组合下的喷灌均匀强度,利用SPSS软件分析主要、次要影响因素及互助关系建立水量分布模型,之后考虑相应的径向、周向与喷灌强度随时间变化,继续使用SPSS软件,拟合出从三维水量分布模型,F(x)=f(X(p,r),d,μ,α,h,t),并最终确定各因素的影响大小。
本发明的技术效果是:1.本发明涉及一种变转速喷头三维动态水量分布模型建立方法,能结合不同影响因素,建立相应的径向、周向与喷灌强度随时间变化的三维水量分布模型,会更接近不同类型喷头的实际喷洒情况,使模型更具有普遍适用性;2.根据模型建立方法可以计算出不同类型喷头的三维动态水量分布,有利于喷灌前沿科学合理的选择合适的喷头。
附图说明
图1为三次三次样条插值法的计算模型图。
图2为喷头喷洒试验雨量筒径向布置示意图。
图3为试验系统示意图。
图4为喷头矩形组合方式示意图。
图5为模型建立流程图。
附图标记说明如下:1.储水槽2.离心泵3.输水管4.流量计5.压力表6.喷头7.雨量筒。
具体实施方式
下面结合附图以及具体实施例对本发明作进一步的说明,但本发明的保护范围并不限于此。
本发明是一种变转速喷头三维动态水量分布模型建立方法,以尼尔森R33喷头建立三维动态水量分布模型为例,该方法的步骤包括:
步骤一:确定喷头型号
使用HR33喷头,该喷头的特点是工作压力范围为150kPa-350kPa,,每旋转一圈内转速变化大于6次。
步骤二:确定喷洒试验收集方法及布置方式
使用雨量筒收集,雨量筒的布置方式为径向布置。雨量筒与雨量筒的间距为1m。如图1
步骤三:设计在不同影响因素下单喷头的水量分布的试验。
如图2所示,通过径向分布的雨量筒法收集不同影响因素下单喷头的喷灌强度hi的实验数据如下表1,表2。水量分布由喷灌强度hi得到。
表1不同影响因素下单喷头的喷灌强度hi
表2不同影响因素下单喷头的喷灌强度hi
表3压力与转速试验关系
压力(kPa) | 150 | 175 | 200 | 250 |
转速(r/min) | 0.85 | 0.95 | 1.05 | 1.25 |
步骤四:确定喷头组合方式和组合间距。
采用矩形组合,喷头与喷头间横向距离为1m,纵向距离为1m。如图3。
步骤五:利用单喷头的水量分布实验数据,计算出组合喷灌均匀度
使用三次样条插值法处理单喷头组合的数据,然后利用克里斯琴森均匀度计算公式(1)进行计算如下表4。
三次样条插值法计算模型为:
计算喷灌区域内任意一点P的将水深,分两步进行。如图2-5所示,首先通过径向插值分别计算出各辐射线上与网格点P距喷头相同远处Ai的降水深Di(i=1、2、...、6),再利用同一圆周上的数据采用与径向插值方法相同的插值方法计算出所求网格点P的降水深Dp。
第一步径向插值:通过Li(i=1、2、…、6)上测得的雨量筒数据dij(i=1、2、…、6;j=1、2、…、n),其中n为每条辐射线上雨量筒的个数,选用一种插值方法,分别求出Ai(i=1、2、…、6)点的降水深Di:
Di=Fi(dij,r)(i=1、2、…、6;j=1、2、…、n)
式中:Fi——径向插值函数
第二步周向插值:通过上一步得到的Ai点的降水深Di(i=1、2、…、6),对极角α采用与径向相同的插值方法进行周向插值,求出P点的降水深DP:
D=F(Dj,α)(i=1、2、…、6)
式中:F——径向插值函数
给定区间[a,b]的一个划分:
Δ:a=x0<x1<…<xn=b
和区间[a,b]上的一个函数f(x)。
若函数满足下列条件:
1.一致通过n+1个型值点(xi,yi),即
s(xi)=f(xi)=yi(i=0,1,2,...,n)
2.二阶导数连续,即s(x)∈c[a,b]
3.三次分段,即在每一个小区间[xi-1,xi]上均为三次多项式。
这样的s(x)称为[a,b]上以xi(i=0,1,2,...,n)为节点的三次插值样条函数,三次插值样条函数的几何图形称为三次样条曲线。s(x)的公式如下:
克里斯琴森均匀度计算公式:
式中Cu—喷灌均匀系数
hi—第i个测点的降水深
h—喷洒面积上各测点平均降水深
n—测点数目
表4不同因素下组合喷灌均匀度
步骤六:分析不同压力p,喷嘴直径d,流量系数μ,转速r、喷头仰角α,安装高度H因素下组合喷灌均匀度数据,使用非线性工具,进行非线性回归分析试验数据,借用一元方差分析的思想,即依据组间均方差与组内均方差之比最大的原则来进行判别,确定了主要影响因素为压力p,喷嘴直径d,流量系数μ,转速r。次要影响因素为喷头仰角α,安装高度H;转速与压力具有互作用的关系。同时借方差分析的工具构造线性判别函数,建立它们的关系式:
F(x)=1.23p2/3α0.4d5/2μ0.7H1.1
其中,p,r是相互影响的关系式为X(p,r)。
步骤七:建立三维动态水量分布模型。
在不同因素下,单喷头水量分布实验,使用三次样条插值法计算出矩形组合下,喷灌均匀强度,利用SPSS软件,分析主要,次要影响因素及互助关系建立水量分布模型,之后考虑相应的径向、周向与喷灌强度随时间变化,继续使用SPSS软件,拟合出从三维水量分布模型,技术路线如图4。
建立它们的关系式为:
F(x)=1.23p2/3α0.4d5/2μ0.7H1.1t0.8
该模型能结合不同影响因素(压力p,喷嘴直径d,流量系数μ,转速r、喷头仰角α,安装高度H因素),确定压力p,喷嘴直径d,流量系数μ,转速r为主要影响因素,喷头仰角α,安装高度H为次要影响因素,建立相应的径向、周向与喷灌强度随时间变化的三维水量分布模型。该模型由于明确了主次影响因素,并加入了时间变量,所以会更接近的实际喷洒情况。同时,根据模型计算HR33喷头的三维动态水量分布,有利于判断在实际喷灌需求中该喷头是否适合。
所述实施例为本发明的优选的实施方式,但本发明并不限于上述实施方式,在不背离本发明的实质内容的情况下,本领域技术人员能够做出的任何显而易见的改进、替换或变型均属于本发明的保护范围。
Claims (1)
1.一种变转速喷头三维动态水量分布模型建立方法,其特征在于,设计方法按下列步骤进行:
步骤一,确定喷头型号:
选定喷头型号,确定选定喷头的工作压力范围,以及喷头旋转一圈转速的变化次数;
步骤二,确定喷洒试验收集方法及布置方式:
使用通用的雨量筒进行水量收集,并采用径向布置的方式排列雨量筒;
步骤三,确定单喷头的水量分布的影响因素并进行试验:
通过径向分布的雨量筒法收集不同影响因素下单喷头的喷灌强度hi、雨滴直径D、雨滴累积频率f的试验数据,其中水量分布由喷灌强度hi得到,水滴分布由雨滴直径D与雨滴累积频率f得到;
步骤四,确定喷头组合方式和组合间距:
采用矩形组合,喷头与喷头间横向距离为1m,纵向距离为1m;
步骤五,利用单喷头的水量分布实验数据,计算出组合喷灌均匀度:
使用三次样条插值法处理单喷头组合的数据,然后利用克里斯琴森均匀度计算公式进行计算得到组合喷灌均匀度;
步骤六,建立影响因素的数值关系式:
分析不同压力p,喷嘴直径d,流量系数μ,转速r、喷头仰角α,安装高度H因素下组合喷灌均匀度数据,使用非线性工具、多维尺度分析,进行非线性回归分析,依据组间均方差与组内均方差之比最大的原则来进行判别,确定主要影响因素,次要影响因素及各个因子直接的互作关系,建立各因素的线性函数X(p,r),其中:p,r是相互影响的,关系式为:F(x)=f(X(p,r),d,μ,α,h);
步骤七,建立三维动态水量分布模型:
在不同因素下,单喷头水量分布实验,通过三次样条插值法计算出矩形组合下的喷灌均匀强度,利用SPSS软件分析主要、次要影响因素及互助关系建立水量分布模型,之后考虑相应的径向、周向与喷灌强度随时间变化,继续使用SPSS软件,拟合出从三维水量分布模型,F(x)=f(X(p,r),d,μ,α,h,t),并最终确定各因素的影响大小。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510543947.1A CN105160096B (zh) | 2015-08-31 | 2015-08-31 | 一种变转速喷头三维动态水量分布模型建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510543947.1A CN105160096B (zh) | 2015-08-31 | 2015-08-31 | 一种变转速喷头三维动态水量分布模型建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105160096A true CN105160096A (zh) | 2015-12-16 |
CN105160096B CN105160096B (zh) | 2018-06-26 |
Family
ID=54800952
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510543947.1A Active CN105160096B (zh) | 2015-08-31 | 2015-08-31 | 一种变转速喷头三维动态水量分布模型建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105160096B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105678079A (zh) * | 2016-01-07 | 2016-06-15 | 江苏大学 | 一种变转速喷头外流场水滴累积频率的计算方法及其改进方法 |
CN105760611A (zh) * | 2016-02-25 | 2016-07-13 | 江苏大学 | 一种低压均匀喷洒喷头空间流道的优化设计方法 |
CN106527130A (zh) * | 2016-11-01 | 2017-03-22 | 江苏大学 | 一种瞬时喷灌强度的计算方法 |
CN106954531A (zh) * | 2017-04-24 | 2017-07-18 | 中国水利水电科学研究院 | 圆形喷灌机喷头配置方法及灌溉方法 |
CN106982711A (zh) * | 2017-03-02 | 2017-07-28 | 江苏大学镇江流体工程装备技术研究院 | 一种直行式喷灌机不同组合间距平均喷灌强度的优化方法 |
CN107577888A (zh) * | 2017-09-18 | 2018-01-12 | 巧夺天宫(深圳)科技有限公司 | 消防喷淋头设计系统、设计方法、电子设备及计算机程序产品 |
CN115879695A (zh) * | 2022-11-14 | 2023-03-31 | 四川大学锦江学院 | 一种基于K-means算法的灌溉喷头布局方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1908615A (zh) * | 2006-08-11 | 2007-02-07 | 西北农林科技大学 | 一种基于连续系统的喷灌均匀度测量方法 |
CN1908617A (zh) * | 2006-08-15 | 2007-02-07 | 西北农林科技大学 | 一种基于标准差的喷灌分布均匀系数测定方法 |
CN103697931A (zh) * | 2013-12-19 | 2014-04-02 | 江苏大学 | 一种变量喷洒单喷头均匀性的评价方法 |
-
2015
- 2015-08-31 CN CN201510543947.1A patent/CN105160096B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1908615A (zh) * | 2006-08-11 | 2007-02-07 | 西北农林科技大学 | 一种基于连续系统的喷灌均匀度测量方法 |
CN1908617A (zh) * | 2006-08-15 | 2007-02-07 | 西北农林科技大学 | 一种基于标准差的喷灌分布均匀系数测定方法 |
CN103697931A (zh) * | 2013-12-19 | 2014-04-02 | 江苏大学 | 一种变量喷洒单喷头均匀性的评价方法 |
Non-Patent Citations (3)
Title |
---|
李德镇等: "基于正交试验的喷嘴喷射特性模拟分析", 《车用发动机》 * |
韩文霆等: "《中国农业工程学会2007年学术会论文集》", 1 August 2007 * |
韩文霆等: "正三角形组合喷灌均匀度计算方法", 《农业机械学报》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105678079A (zh) * | 2016-01-07 | 2016-06-15 | 江苏大学 | 一种变转速喷头外流场水滴累积频率的计算方法及其改进方法 |
CN105678079B (zh) * | 2016-01-07 | 2018-02-27 | 江苏大学 | 一种变转速喷头外流场水滴累积频率的计算方法及其改进方法 |
CN105760611B (zh) * | 2016-02-25 | 2019-08-27 | 江苏大学 | 一种低压均匀喷洒喷头空间流道的优化设计方法 |
CN105760611A (zh) * | 2016-02-25 | 2016-07-13 | 江苏大学 | 一种低压均匀喷洒喷头空间流道的优化设计方法 |
WO2017143694A1 (zh) * | 2016-02-25 | 2017-08-31 | 江苏大学 | 一种低压均匀喷洒喷头空间流道的优化设计方法 |
US10984149B2 (en) | 2016-02-25 | 2021-04-20 | Jiangsu University | Optimization design method for spatial flow passage of low-pressure even spray nozzle |
CN106527130A (zh) * | 2016-11-01 | 2017-03-22 | 江苏大学 | 一种瞬时喷灌强度的计算方法 |
US11504730B2 (en) | 2016-11-01 | 2022-11-22 | Jiangsu University | Method for calculating instantaneous sprinkler strength |
WO2018082129A1 (zh) * | 2016-11-01 | 2018-05-11 | 江苏大学 | 一种瞬时喷灌强度的计算方法 |
CN106527130B (zh) * | 2016-11-01 | 2019-04-02 | 江苏大学 | 一种瞬时喷灌强度的计算方法 |
CN106982711A (zh) * | 2017-03-02 | 2017-07-28 | 江苏大学镇江流体工程装备技术研究院 | 一种直行式喷灌机不同组合间距平均喷灌强度的优化方法 |
CN106982711B (zh) * | 2017-03-02 | 2019-11-05 | 江苏大学镇江流体工程装备技术研究院 | 一种直行式喷灌机不同组合间距平均喷灌强度的优化方法 |
CN106954531B (zh) * | 2017-04-24 | 2020-07-17 | 中国水利水电科学研究院 | 圆形喷灌机喷头配置方法及灌溉方法 |
CN106954531A (zh) * | 2017-04-24 | 2017-07-18 | 中国水利水电科学研究院 | 圆形喷灌机喷头配置方法及灌溉方法 |
CN107577888A (zh) * | 2017-09-18 | 2018-01-12 | 巧夺天宫(深圳)科技有限公司 | 消防喷淋头设计系统、设计方法、电子设备及计算机程序产品 |
CN115879695A (zh) * | 2022-11-14 | 2023-03-31 | 四川大学锦江学院 | 一种基于K-means算法的灌溉喷头布局方法 |
CN115879695B (zh) * | 2022-11-14 | 2024-02-09 | 四川大学锦江学院 | 一种基于K-means算法的灌溉喷头布局方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105160096B (zh) | 2018-06-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105160096A (zh) | 一种变转速喷头三维动态水量分布模型建立方法 | |
CN103697931B (zh) | 一种变量喷洒单喷头均匀性的评价方法 | |
CN102402185A (zh) | 基于模糊控制的非充分灌溉控制方法 | |
Ge et al. | Analysis of kinetic energy distribution of big gun sprinkler applied to continuous moving hose-drawn traveler | |
CN106951674B (zh) | 一种无历史资料小流域洪水重现期推求方法 | |
CN101720654A (zh) | 草坪节水灌溉预报方法及系统 | |
Dwomoh et al. | Field performance characteristics of fluidic sprinkler | |
Alabas | Evaluation the hydraulic performance of drip irrigation system with multi cases | |
Smith et al. | A decision support model for travelling gun irrigation machines | |
CN103714563A (zh) | 一种曲线形农田作业区域边界建模方法 | |
CN107918398B (zh) | 一种基于多重优化的集群无人机协同定位方法 | |
CN110110457B (zh) | 一种复杂地形风电场的噪声分布预测方法 | |
CN103955872A (zh) | 一种结合蚁群算法和压力相似性分析的管网压力分区方法 | |
Stambouli et al. | Performance of new agricultural impact sprinkler fitted with plastic nozzles | |
Fordjour et al. | Effect of riser height on rotation uniformity and application rate of the dynamic fluidic sprinkler | |
CN110264101A (zh) | 一种盐渍农田的作物水分生产力的优化方法 | |
CN109559017A (zh) | 青海气区气井泡沫排水采气实施前技术指标界限及其确定 | |
CN105468819B (zh) | 一种利用响应面法优化垂直摇臂式喷头工作参数的方法 | |
CN102621075B (zh) | 一种水稻抽穗期自动检测的方法 | |
CN107622145A (zh) | 一种均匀热成形模具型面水管流速的方法 | |
CN108470015A (zh) | 一种平移式喷灌机组合均匀度的计算方法 | |
CN107515996B (zh) | 一种达里厄水轮机导流罩型线的优化设计方法 | |
CN103506237B (zh) | 一种正方形喷洒域喷头的布置方法 | |
CN106570321B (zh) | 一种平移式喷灌机不同工况喷洒水深的计算方法 | |
CN106527130B (zh) | 一种瞬时喷灌强度的计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |