会员注册 | 登录 | 微信快捷登录 支付宝快捷登录 QQ登录 微博登录 | 帮助中心 人人文库renrendoc.com美如初恋!
站内搜索 百度文库

热门搜索: 直缝焊接机 矿井提升机 循环球式转向器图纸 机器人手爪发展史 管道机器人dwg 动平衡试验台设计

水利工程论文-坡面降雨入渗产流规律的数值模拟研究.doc水利工程论文-坡面降雨入渗产流规律的数值模拟研究.doc -- 2 元

宽屏显示 收藏 分享

资源预览需要最新版本的Flash Player支持。
您尚未安装或版本过低,建议您

水利工程论文坡面降雨入渗产流规律的数值模拟研究摘要本文采用运动波理论和两次改进后的GreenAmpt入渗模型建立了坡面降雨入渗产流的动力学模型,并得到了实验资料的良好验证。运用该模型研究了简单坡面上降雨入渗产流的动力学规律,分析了雨强、土壤初始含水量、渗透系数、坡面阻力,以及坡长、坡度等因素对坡面产流过程的影响规律,得出了一些有益的结论。关键词入渗产流坡面动力学1概述雨水降落在坡面上将产生雨水的聚集并形成坡面水流。坡面水流是土壤水蚀过程的主要动力,搞清产流的动力学特点是进一步研究侵蚀过程规律的基础。坡面水流不同于一般明渠流动,其水深极浅一般只有几毫米,沿程不断有质量源和动量源加入,使其随时间和空间有较大的变化。且坡面流的坡度较一般河渠陡得多,边界条件也更为复杂。这些特点使得对坡面水流的研究有相当的难度。坡面产流研究已有很长历史,但对它的数学求解还只有三十多年。60年代后期Woolhiser和Ligget1967将运动波模型引入坡面水流研究,大大简化了计算工作,促进了研究的发展。运动波模型是从一维圣维南方程简化而来,其基本假设是水流的能坡和底坡相等,并借助Chezy阻力公式得到流量和水深的关系。Woolhiser和Ligget的研究结果表明在运动波波数k10时,运动波模型可以很好地描述坡面水流运动。而实际坡面流的运动波波数一般远大于10沈冰等,1996。因此,运动波近似是一种较好的数学描述方式。其后,又有对运动波理论的修正Ponce,1978,Govindaraju,1988,保留了水深的沿程变化项,相当于压力梯度,被称为扩散波模型。该模型扩展了适用的参数范围,但并无实质性改进,因此实际应用仍以运动波为主。也有使用完整圣维南方程求解实际问题的戚隆溪,1997。土壤入渗过程的研究也有很长历史,从1911年提出概念明确形式简单的GreenAmpt积水入渗模型开始,相继有Horton1940,Philip1957等模型出现,但GA模型仍以其简单的形式,明晰的物理概念,良好的扩展性和可信的应用效果受到广泛重视,特别是经过MeinLarson1973和Chu1978的两次改进,使其可应用于不均匀降雨的入渗计算,更使它成为最有效和应用最广泛的模型。在国内,GA模型尚未受到重视,Horton模型曾得到相当广泛的运用,但其参数的物理意义明显不如GA模型明晰。也有研究者使用更基本的土壤水分运动微分方程,但所需的参数更加难于获取,计算也更为复杂。本文工作旨在建立物理概念明晰的降雨入渗产流综合计算模式,并用以研究简单坡面的产流过程,分析各主要因素的影响和各主要因素的影响和各主要参量的变化规律。以期对坡面产流的动力学规律有清楚的认识。2计算模式坡面流运动十分复杂,目前主要采用运动波理论、扩散波或完整圣维南方程进行描述。正如前文所述,运动波近似理论在大多数情况下可以很好地描述坡面流运动过程,且计算简单。因此本文仍采用一维运动波理论,即坡面流基本方程为1此处第二式直接使用了水力学中熟知的Chezy公式和Manning公式。其中,x为沿坡面向下的坐标,t为时间s,h为水深m,q为单宽流量m2/s,p为降雨强度m/s,此处假设降雨方向垂直向下,i为入渗率m/s,S0为坡面坡度,S0sinθ,θ为坡面倾角,n为Manning糙率系数。土壤的入渗过程对坡面流的形成和流动过程影响很大,本文采用形式简单、物理概念明晰的GA入渗模型,其计算方程为idI/dtK1θSθiS/IIKtSθSθiln1I/SθSθi2其中K为土壤饱和导水率渗透系数m/s,θS为土壤饱和含水率,即有效孔隙率,θi为土壤初始含水率,S为土壤吸力m,I为累积入渗量m。经典的GreenAmpt模型是干土积水入渗模型,其前提是在整个入渗过程中地表始终有积水。MeinLarson1973年将其推广应用至降雨入渗的情况。设有稳定的雨强p,只有p大于土壤的入渗能力时,地表才能形成积水。而在降雨的初始阶段,全部降雨都渗入地下。由GA模型知,入渗率是随累积入渗量的增加而减小的。设想当累积入渗量达到某一值时,ip,此时开始积水,称此累积入渗量为Ip。因此由GA模型入渗公式可以导出开始积水时的Ip值IpθSθiS/p/K13开始积水时间由tpIp/p给出。因此整个过程的入渗率可表示为ipt≤tpiK1θSθiS/Ittp4式中的I为积水开始后的累积入渗量包含未积水时段的入渗量在内。由于不是由t0开始积水,I的计算须采用修正后的公式KttptSISθSθiln1I/SθSθittp5tS表示假设由t0开始积水,到入渗量IIp或ip时所需时间,可理解为一个虚拟时间,可计算如下KtSIpSθSθiln1Ip/SθSθi6改进的主要思想是将整个过程假设为从一开始就是积水入渗,这样该曲线在积水后部分相对于实际入渗曲线将向左平移tptS,将这条曲线向右平移tptS,再加上积水前的入渗强度等于降雨强度的关系,就得到真实的入渗过程。但稳定降雨在实际应用中远不能满足需要,Chu1978将MeinLarson1973改进的GA模型再作推广,使其可应用于变化的降雨过程。基本作法是,对每个计算时段将地表状态分为四种情况1.开始无积水,结束无积水2.开始无积水,结束有积水3.开始有积水,结束有积水4.开始有积水,结束无积水在每一时段开始,已知降雨总量与入渗总量,剩余总量。根据两个因子判断时段结束时是否有积水。若时段开始无积水,使用因子cu,若时段开始有积水,使用因子cp,其表达式为cuPtnRtn1KSM/iKcu0时段结束将积水,cu0时段结束仍有积水,cu0积水消失7其中M代表θSθi,Ptn代表tn时刻降雨总量,Rtn1代表tn1时刻剩余总量。可以证明,时段结束时积水与否与此两因子的正负等价。当iK时,始终无积水,不用此两因子判断。3模型求解模型的求解也分为两部分。运动波方程是一个非线性的对流型方程,求解采用了简单的一阶显式迎风格式。经过比较,一阶显式迎风格式对此问题的计算结果在幅值、相位、守恒性几个方面的综合效果较其他一些格式为优。我们曾使用过Preissmann四点偏心格式,它在计算此类坡面薄层水流运动波方程时耗散过大,表现为达到平衡产流的时间偏快。我们也曾试过中心差分格式和高阶迎风格式,结果发现标准的中心差分格式仍然和在普通的线性对流型方程中的表现一样
编号:201312122148225229    大小:21.52KB    格式:DOC    上传时间:2013-12-12
  【编辑】
2
关 键 词:
专业文献 建筑水利 精品文档 水利工程
温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 人人文库网仅提供交流平台,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
  人人文库网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

暂无评论,赶快抢占沙发吧。

当前资源信息

4.0
 
(2人评价)
浏览:4次
奋斗不息上传于2013-12-12

官方联系方式

客服手机:13961746681   
2:不支持迅雷下载,请使用浏览器下载   
3:不支持QQ浏览器下载,请用其他浏览器   
4:下载后的文档和图纸-无水印   
5:文档经过压缩,下载后原文更清晰   

相关资源

相关资源

相关搜索

专业文献   建筑水利   精品文档   水利工程  
关于我们 - 网站声明 - 网站地图 - 友情链接 - 网站客服客服 - 联系我们
copyright@ 2015-2017 人人文库网网站版权所有
苏ICP备12009002号-5