基于拉格朗日插值的挡潮排涝涵闸调洪演算.pdfVIP

第O2期 (总第 381期) 吉 林 水 利 2014年02月 [文章编号]1009—2846(2014)02—0042—04 基于拉格朗日插值的挡潮排涝涵闸调洪演算 夏明月 (辽宁江河水利水电新技术设计研究院，辽宁 沈 阳 110003) [摘要]为了探究涵闸进 口水位一流量关系未知且外江潮水对 内江洪水具有顶托作用的调洪演算过程 ，在传统 的列表试算方法基础上提 出了单时段 内根据涵 闸运行原则且基于拉格 朗 日il1次插值多项式的泄流前及泄流后容 积水深插值调洪算法．为了计算方便，采用C++编写拉格朗日插值程序。将改进算法的计算结果代入传统的列 表试算法中进行验证．计算结果表明只有个别时段水深值具有微小偏差．改进后的算法具有较高的精度。改进 后的算法不仅可以避免传统算法的大量试算．为传统算法提供准确的试算初始点，而且对于类似实际工程的调 洪演算具有一定的参考价值 [关键词]拉格朗日插值；潮水顶托；列表试算；泄流过程线 [中图分类号]Tv663 [文献标识码]B 当涵 闸受到外江潮水顶托影响时．在某些工 ()= +口1戈+…+ 况下涵闸会关闸挡水．此时来流洪水会在库区滞 =yolo()+y111(戈)+…+yl() 留．使得水位增加。如果 闸顶高程的设计没有考虑 潮水的顶托作用 ．会对工程的安全造成一定的威 = ∑ () (1) 胁。因此 ．探讨涵 闸在受到外江潮水顶托的影响下 式中：lo(x)，…，fn()为 n次插值基函数，其表 如何进行调洪演算是很有必要的 达式为 ： 水库调洪演算方法 。是水库防洪计算的基础[1]。 调洪演算方法大致有列表试算法、半图解法和图 (2) 解法。随着计算机的普及．为了使计算程序化．出 其中 现了数值一解析解法圆、牛顿迭代法f31、人工神经网 ∞，1()=( 0)(—I)…( ) (3) 络方法4[1等 本文是在列表试算法的基础上引入拉 格朗 日插值水深容积 ．并按照涵 闸的运行原则及 ∞叶1()=(·0)…(— 1)(Xi--~．1)…(i— ) (4) 上下游水深与涵洞高度的关系进行泄水流量 的计 算 ．从而得到泄水流量过程线 2 调洪演算的基本原理 l 拉格朗 日型 1／,次插值 多项式 调洪演算的本质思想即计算时段 内进入库区 的平均洪水量与平均下泄量之差等于上游库容的 若函数 )在 n+1个不 同点 。，…， 上的函 变化量 。进行调洪计算主要是为了确定水库的调 数值分别为yo，…，yn，则唯一存在一个次数不超过 洪库容 、最大泄量 q 以及设计／校核洪水位 。其 n的多项式 Pn()去描述函数 ()，且 (兢) ，