明渠恒定非均匀流.ppt

1 第六章明渠恒定非均匀流人工渠道或天然河道中的水流绝大多数是非均匀流 明渠非均匀流的特点是明渠的底坡线 水面线 总水头线彼此互不平行 明渠非均匀流分为明渠非均匀渐变流和明渠非均匀急变流 本章着重研究明渠中恒定非均匀渐变流的基本特性及其水力要素 主要是水深 沿程变化的规律 具体地说 就是要分析水面线的变化及其计算 以便确定明渠边墙高度 以及回水淹没的范围等 通常把明渠均匀流的水深称为正常水深h0 2 6 1明渠水流的三种流态 扰动 在流场的某一点或者某一个区域 由于某种原因 使流动参数发生变化 这种变化叫做扰动 波 扰动区域与未扰动区域的分界面 3 微弱扰动的一维传播 非定常过程 4 6 1明渠水流的三种流态 注意 波速与流体质点速度的区别 5 6 1明渠水流的三种流态 当v 0时 水流静止 干扰波能向四周以一定的速度传播 在t 0 1 2 3 4s 分别有水滴滴入o点 研究t 4s的流动图象静水中传播的微波速度vw c 称为相对波速 6 当v vw时 水流为缓流 干扰波能向上游和下游传播 在t 0 1 2 3 4s 分别有水滴滴入o点 研究t 4s的流动图象 7 当v vw时 水流为急流 干扰波不能向上游传播 只能向下游传播 马赫椎内 在t 0 1 2 3 4s 分别有水滴滴入o点 研究t 4s的流动图象 马赫角 马赫锥的半顶角 即圆锥的母线与来流速度方向之间的夹角 8 当v vw时 水流为临界流 在t 0 1 2 3 4s 分别有水滴滴入o点 研究t 4s的流动图象 9 6 1明渠水流的三种流态明渠水流有和大气接触的自由表面 与有压流不同 具有独特的水流流态 即缓流 临界流和急流三种 静水中传播的微波速度vw称为相对波速 当v 0时 水流静止 干扰波能向四周以一定的速度传播 当v vw时 水流为缓流 干扰波能向上游和下游传播 当v vw时 水流为临界流 当v vw时 水流为急流 干扰波不能向上游传播 只能向下游传播 马赫椎内 10 微波波速的计算 以一竖直平板在平底矩形棱柱体明渠中激起一个干扰微波 观察者随波前行 对上述的运动坐标系水流作恒定非均匀流动 不计摩擦力对1 1和2 2断面建立连续性和能量方程 11 联解上两式 并令 消去得令 则微波波速 明渠断面为任意形状时 式中 为断面平均水深 A为断面面积 B为水面宽度 12 实际工程中微波传播的绝对速度对临界流断面平均流速恰好等于微波相对波速对临界流有佛汝德数 佛汝德数的 物理意义是 流速与相对波速之比 13 显然 当Fr 1 水流为缓流 当Fr 1 水流为临界流 当Fr 1 水流为急流 佛汝德数的 物理意义是 过水断面单位重量液体平均动能与平均势能之比的二倍开平方 14 佛汝德数的 物理意义 即佛汝德数的力学意义是 代表水流的惯性力和重力两种作用的对比关系 15 6 2断面比能与临界水深 明渠中水流的流态也可从能量的角度来分析 一 断面比能 比能曲线如图所示渐变流 若以0 0为基准面 则过水断面上单位重量液体所具有的总能量为 16 如果我们把参考基准面选在渠底这一特殊位置 把对通过渠底的水平面0 0 所计算得到的单位能量称为断面比能 并以来表示 则在实用上 因一般明渠底坡较小 可认为故常采用 17 当流量Q和过水断面的形状及尺寸一定时 断面比能仅仅是水深的函数 即Es f h 以图表示则称为 比能曲线 18 为什么 19 因在过水断面上 代入上式有若取 则有因而对断面比能曲线有 20 二 临界水深相应于断面单位能量最小值的水深称为临界水深 以hk表示 由临界流方程当流量和过水断面形状及尺寸给定时 利用上式即可求解临界水深 注以脚标表示临界水深时的水力要素 21 1 矩形断面明渠临界水深的计算上式中为单宽流量 22 23 2 断面为任意形状时 临界水深的计算 24 1 试算法当给定流量Q及明渠断面形状 尺寸后 6 15 式的左端为一定值 该式的右端乃仅仅是水深的函数 于是可以假定若干个水深h 从而可算出若干个与之对应的值 当某一值刚好与相等时 其相应的水深即为所求的临界水深hK 25 2 图解法图解法的实质和试算法相同 当假定不同的水深h时 可得出若干相应的值 然后将这些值点绘成h 关系曲线图 见图 在该图的轴上 量取其值为的长度 由此引铅垂线与曲线相交于C点 C点所对应的h值即为所求hK 26 3 等腰梯形断面临界水深计算若明渠过水断面为梯形 且两侧边坡相同 在这种情况下 可应用一种简便图解法 现将其原理简述如下 对于等腰梯形断面有 代入 6 15 式可得 令 1 将上式两端同除以后开立方则得 27 6 20 上式中 b为梯形断面的底宽 上式左端实际上表示一个与梯形断面底宽相等的矩形断面的临界水深 为了与欲求的梯形断面的临界水深相区别将其以来表示 即令 6 21 28 若将 6 20 式两端同乘以可得 6 22 上式移项后可得 6 23 29 30 31 求解梯形断面临界水深的方法 1 求出与梯形断面底宽相等的矩形断面的临界水深 2然后根据梯形断面已知m b值算出 3再由关系曲线上查出相应的值 从而可算出梯形断面的值 32 当h 时 Fr1 为急流 根据所给流量及断面尺寸 应用上述方法求出临界水深以后 也可用来判断 33 例6 1一矩形断面明渠 流量Q 30m3 s 底宽b 8m 要求 1 用计算及图解法求渠中临界水深 2 计算渠中实际水深h 3m时 水流的弗劳德数 微波波速 并据此以不同的角度来判别水流的流态 解 1 求临界水深或者 由附图 右下角hK q关系图上可查出当q 3 75m3 s m时 hK 1 13m 34 2 当渠中水深h 3m时渠中流速弗劳德数微波波速临界流速 35 从水深看 因h hK 故渠中水流为缓流 以Fr为标准 因Fr 1 水流为缓流 以微波波速与实际水流流速作比较 因 微波可以向上游传播 故水流为缓流 以临界流速与实际水流流速作比较 因 故水流为缓流 36 例6 2一梯形断面渠道 底宽b为5m 边坡系数m为1 要求 计算通过流量分别为Q1为10m3 s Q2为15m3 s Q3为20m3 s时的临界水深 解 1 绘制h 关系曲线因对梯形断面先假定若干h 计算相应的值 计算成果见表6 1 37 表6 1根据表中数值 绘制h 关系曲线 如图6 8所示 2 计算各级流量下的值 并由图中查读临界水深 38 当时 由图查得当时 由图查得当时 由图查得 39 例6 3已知梯形断面渠道 b为45m m为2 0 Q为500m3 s 要求 用图解法求临界水深解 因查附图III右下角曲线 得计算 根据 查同一图左上角曲线得为0 97则临界水深 40 6 3临界底坡 缓坡与陡坡 到目前我们知道了三种水深 均匀流正常水深非均匀流水深临界水深 若已知明渠断面形状及尺寸 当流量一定时 均匀流情况下可将底坡与渠中正常水深的关系 p213 为 明渠均匀流的正常水深h0恰好与临界水深hk相等时 此坡度定义为临界底坡 当底坡i增大时 正常水深h0将减小 反之当i减小时正常水深h0将增大 41 在临界底坡上作均匀流时 满足临界流的条件式另一方面又要同时满足均匀流的基本方程式联立可得临界底坡的计算式为 42 一个坡度为i 0 的明渠 与其相应 即同流量 同断面尺寸 同糙率 的临界底坡相比较可能有三种情况 为缓坡 为陡坡 为临界坡 43 在明渠均匀流的情况下 若则若则若则用底坡的类型就可以判别水流的流态 即在缓坡上水流为缓流 在陡坡上水流为急流 在临界坡上水流为临界流 这种判别只能适用于均匀流的情况 在非均匀流时就不一定了 44 45 46 47 48 49 6 4临界水深的一些实例 临界水深对分析明渠流具有重要的意义 1 河道或渠道中如知道临界水深并量取该断面的尺寸 其流量就能简便地估算出来 2 可将发生临界水深断面作为控制断面据此来推求上下游水面曲线 50 6 4临界水深的一些实例 一 当渠道底坡自陡坡变为缓坡时此时水流会产生一种水面突然跃起的特殊水力现象叫水跃 水跃自水深小于临界水深跃入大于临界水深 其间必经过临界水深 K K线与底坡i无关 51 二 当渠道底坡自缓坡变为陡坡时渠道中均匀流由缓流变为急流时 水流会产生水面降落现象 叫做水跌 52 三 当缓坡渠道末端自由跌落时 53 三 当缓坡渠道末端自由跌落时 54 四 当水流自水库进入陡坡渠道时水库中水流为缓流 而陡坡渠道中均匀流为急流 水流由缓流过渡到急流时 必经过临界水深 55 6 5明渠恒定非均匀渐变流的微分方程式在底坡为i的明渠渐变流中 沿水流方向任取一微分流段ds 对微分流段上 下游断面建立能量方程如下 56 又因化简可得非均匀流沿程水头损失尚无精确的计算方法 仍然近似地采用均匀流公式计算 局部水头损失用代入可得 57 若明渠底坡i值小于1 10 实用上一般都采用 并用铅垂水深代替垂直于槽底水深上式即为明渠恒定非均匀渐变流的基本微分方程式 58 6 5明渠恒定非均匀渐变流的微分方程式下面研究明渠水深沿流程的变化规律 一 水深沿程变化的微分方程式式中一般情况下 A f h s 对非棱柱体明渠非均匀渐变流水深沿流程变化的微分方程式为 59 棱柱体明渠 60 二 水位沿流程变化的微分方程式在天然河道中 常用水位的变化来反映非均匀渐变流变化规律 因又有所以即 61 二 水位沿流程变化的微分方程式将 35 代入 27 得非均匀渐变流的水位沿流程变化微分方程为上式对棱柱体及非棱柱体均适用 主要用于探讨天然河道水流的水位变化规律 62 6 6棱柱体明渠中恒定非均匀渐变流水面曲线分析棱柱体明渠非均匀渐变流微分方程式为 上式表明水深h沿流程s的变化是和渠道底坡i及实际水流的流态有关 按底坡性质分为正坡 平坡和逆坡 对正坡明渠可分为缓坡 陡坡和临界坡三种情况 63 正坡明渠中 水流有可能做均匀流动 因而存在正常水深h0 另一方面它也存在临界水深 64 在平底及逆坡棱柱体明渠中 因不可能有均匀流 不存在正常水深h0 仅存在临界水深 所以只能画出与渠底相平行的临界水深线K K 下图乃是平底和逆坡棱柱体明渠中K K线情况 65 1 a区凡实际水深h既大于 又大于即凡是在K K线和N N线二者之上的范围称为a区 66 2 b区凡是实际水深h介于和之间的范围称为b区 b区可能有两种情况 K K线在N N线之下 缓坡明渠 或K K线在N N线之上 陡坡明渠 无论那种情况都属于b区 对于平底和逆坡棱柱体明渠 因不存在N N线 或者可以设想N N线在无限远处 所以只存在b区与c区 67 3 c区凡是实际水深h既小于又小于的区域 即在N N线及K K线二者之下的区域 对于平底和逆坡棱柱体明渠 因不存在N N线 所以只存在b区与c区 K K线之上为b区 K K线之下为c区 68 对于正底坡 存在均匀流由式 5 8 得 69 对于缓坡i ik a1型壅水曲线如右图所示 70 对于缓坡i ik a1型壅水曲线如右图所示 71 b1型降水曲线如左上图所示 72 b1型降水曲线如左上图所示 73 c1型壅水曲线如左下图所示 74 缓坡水面线的型式及实例见下图 75 陡坡水面线a2 h hk h0 K K0 Fr0b2 hk h h0 K K0 Fr 1 故dh dsh0 h K1 故dh ds 0 76 临界坡水面线a3 h hk h0 K K0 Fr0c3 hk h0 h K1 故dh ds 0 77 平坡和逆坡水面线b0和b Fr1 故dh ds 0 78 79 例6 5试讨论分析下图所示两段断面尺寸及糙率相同的长直棱柱体明渠 由于底坡变化所引起渠中非均匀流水面变化形式 已知上游及下游渠道底坡均为缓坡 但i2 i1 80 解 为分析渠中水面变化 首先分别画出上 下游渠道的K K线及N N线 由于上 下游渠道断面尺寸相同 故两段渠道的临界水深均相等 而上 下游渠道底坡不等 故正常水深则不等 因i1h02 下游渠道的N N线低于上游渠道的N N线 因渠道很长 在上游无限远处应为均匀流 其水深为正常水深h01 下游无限远处亦为均匀流 其水深为正常水深h02 81 由上游较大的水深h01要转变到下游较小的水深h02 中间必经历一段降落的过程 水面降落有三种可能 1 上游渠中不降 全在下游渠中降落 2 完全在上游渠中降落 下游渠中不降落 3 在上 下游渠中分别都降落一部分 82 在上述三种可能情况中 若按照第一种或第三种方式降落 那么必然会出现下游渠道中a区发生降水曲线的情况 前面已经论证 缓坡a区只能存在的是壅水曲线 所以第一 第三两种降落方式不能成立 惟一合理的方式是第二种 即降水曲线全部发生在上游渠道中 由上游很远处趋近于h01的地方 逐渐下降至分界断面处水深达到h02 而下游渠道保持为h02的均匀流 所以上游渠道水面曲线为b1型降水曲线 1 上游渠中不降 全在下游渠中降落 2 完全在上游渠中降落 下游渠中不降落 3 在上 下游渠中分别都降落一部分 83 84 85 i1 ik i2 0 i3 ik 从转折截面左右分析 86 i4 ik i3 ik i2 ik i1 0 87 88 从急流过渡到缓流要发生水跃 水跃的位置有三种 h01的共轭水深用h01 来表示 若h02 h01 水跃在变坡处发生 称为临界水跃 89 若h02 h01 水跃在下游渠道发生 称为远驱式水跃 90 若h02 h01 水跃在上游渠道中发生 称为淹没水跃 91 从急流过渡到缓流要发生水跃 水跃的位置有三种 h01的共轭水深用h01 来表示 若h02 h01 水跃在变坡处发生 称为临界水跃 若h02h01 水跃在上游渠道中发生 称为淹没水跃 92 6 7明渠恒定非均匀渐变流水面线曲的计算 逐段试算法一 基本计算公式明渠恒定非均匀渐变流的基本微分方程式为 略去局部损失 得 流段的平均水力坡度一般采用 93 平均值或可用以下三种方法之一计算 或者 94 二 计算方法逐段推算实际计算可能有两种类型 1 已知流段两端的水深 求流段的距离 对于棱柱体明渠 可直接计算 对于非棱柱体明渠 必须采用试算法 2 已知流段一端的水深和流段长 求另一端断面水深 须采用试算法 95 例6 6一长直棱柱体明渠 底宽b为10m m为1 5 n为0 022 i为0 0009 当通过流量Q为45m3 s时 渠道末端水深h为3 4m 求 计算渠道中的水面曲线 解 一 由于渠道底坡大于零 应首先判别渠道是缓坡或是陡坡 水面曲线属于哪种类型 本题条件与例6 4相同 由例6 4计算已知 再计算均匀流水深 因 96 由附图 查得 所以 因 故渠道属于缓坡 又因下游渠道末端水深大于正常水深 所以水面线一定在a区 水面线为型壅水曲线 型水面曲线上游端以正常水深线为渐近线 取曲线上游端水深比正常水深稍大一点 即 97 二 计算水面曲线首先列出各计算公式 式中 98 今以 求两断面间之距离 将有关已知数值代入上列公式中 分别求得 99 100 1 374 10 4 101 其余各流段的计算完全相同 采用列表法进行 102 三 根据表6 2的数值 绘制水面曲线见图6 23 103 例6 7某一边墙成直线收缩的矩形渠道 如图6 24 渠长60m 进口宽b1为8m 出口宽b2为4m 渠底为反坡 i为 0 001 粗糙系数n为0 014 当Q为18m3 s时 进口水深h1为2m 要求计算中间断面及出口断面水深 104 解 渠道宽度逐渐收缩 故为非棱柱体明渠 求指定断面的水深 必须采用试算法仍引用公式 6 40 来计算 105 一 计算中间断面的水深已知中间断面宽度b为6m 今假定其水深h为1 8m 按下列各式计算有关水力要素 106 107 将以上各值列于表6 3中 又因进口断面宽度及水深已知 按以上公式计算进口断面的各水流要素 将计算结果列于表6 3中 根据表6 3中有关数值 代入 6 40 式中 算出为 108 109 计算得到为93 4m 与实际长度30米相差甚远 说明前面所假设之水深1 8米与实际不符合 必须重新假设 故已假设中间断面水深为1 9米 按以上程序计算 得到为29 58米 与实际长度非常接近 所以可认为中间断面水深为1 9米 二 出口断面水深的计算与前面的计算方法完全一样 不再赘述 从表6 3看出 出口水深应为1 5米 110 6 8河渠恒定非均匀流的流量与糙率 一流量的计算若已知河道或明渠的水面线以及断面的尺寸 底坡等 需要估算河渠中流量 而河渠中水流常属于非均匀流 应按照非均匀流情况来处理 例如根据洪水痕迹估算洪水流量乃属于这种类型 非均匀渐变流的水位沿流程变化微分方程为 111 二粗糙系数的计算对于具有显著非均匀特征的水流 必须按照非均匀流来处理 将满宁公式代入下式得 112 例6 8某河流经过洪水调查找到历史上某次洪水的若干痕迹 测得有关资料列于表6 4 河段上 下游断面间距 s为183m 水位落差 z为2 59m 113 例6 9岷江上游干流某河段长128m 曾进行水力要素及河道断面的实测 当河道中流量为1620m3 s时 河段有关资料列于表6 5 若不计局部损失 试计算该河段粗糙系数 114 6 9河道水面曲线的计算在计算河道水面线之前 先要收集有关水文 泥沙及河道地形等资料 如河道粗糙系数 河道纵横剖面图等 由于河道断面形状极不规则 在计算河道水面曲线时 划分计算流段应注意以下几个方面 1 要求每个计算流段内 过水断面形状 尺寸以及粗糙系数 底坡等变化都不太大 115 2 在一个计算流段内 上 下游断面水位差不能过大 一般对平原河流取0 2 1 0米 山区河流取1 0 3 0米 3 每个计算流段内没有支流流入或流出 若河道有支流存在 必须把支流放在计算流段的入口或出口 对加入的支流最好放在流段的进口附近 流出的支流放在流段的出口 116 117 一 试算法由得式中得到 118 并把方程中同一断面的水力要素列在等式的同一端 得到 方程的两端各是上游水位和下游水位函数 119 试算法的步骤 1 若已知下游断面的水位Zd 按式求出函数值 2 假定若干上游断面水位Zu 按式 6 51 的第一式可算出相应的若干f Zu 值 并绘制Zu f Zu 关系曲线 如图 120 在图的横坐标上截取一段向曲线作铅垂线交曲线于A点 A点之纵坐标值即为所求之上游断面水位 该水位对另一个流段来讲 又是下游水位 循着上述步骤可依次求出上游各断面的水位 从而得到全河道的水面曲线 121 二 图解法河道水面曲线图解法很多 所有图解法归纳起来可分为两大类型 一种类型是考虑了流速水头和局部损失的图解法 另一种类型是忽略流速水头和局部损失的图解法 从本质上看 后者是在前者基础上的简化 只要掌握了前一种方法 第二种情况就迎刃而解 122 二 计入流速水头和局部损失的图解法 123 图解前先作好各计算流段的断面函数曲线即和关系曲线 以水位为纵坐标 以下游及上游断面团数为横坐标并分别布置在水位