采取excel 求解气垫式调压室水位涌波过程.pdfVIP

第４０卷 第１９期 人　民　长　江 Ｖｏｌ．４０，Ｎｏ．１９ ２００９年 １０月　 　 Ｙａｎｇｔｚｅ　Ｒｉｖｅｒ　 　Ｏｃｔ．， ２００９ 　　文章编号：１００１－４１７９（２００９）１９－００２１－０３ 采用 Ｅｘｃｅｌ求解气垫式调压室水位涌波过程 欧 传 奇 （河海大学水利水电工程学院，江苏 南京２１００９８） 摘要：因气垫调压室水位涌波方程中包含具有绝对值的平方阻尼项，难以获得其解析计算公式。而借助于水力 过渡过程计算软件的电算法虽然计算精度高，却要依赖专业的计算程序，且数据处理及建模的工作量也大。为 此，结合数值计算方法的特点及常用软件Ｅｘｃｅｌ的优势，设计了采用Ｅｘｃｅｌ求解水位涌波过程的便捷方法。此 法具有求解思路清晰、操作简单等优点，是一种既实用又便于推广的数值解算方法。 关　键　词：涌波；Ｅｘｃｅｌ；龙格—库塔法；气垫式调压室 中图分类号：ＴＶ７３２．５　　　文献标识码：Ａ 　　目前，气垫式调压室的最高、最低涌浪多半是采用水力过渡 业计算软件（Ｍａｔｌａｂ、Ｍａｔｈｅｍａｔｉｃａ等），具有一定的优势。不过， ［１，２］ 过程计算软件或程序经电算 得到的，虽然计算精度高，但针 数值计算方法对设计人员提出了较高的要求，且需依赖一些专 对的是确定的输水系统，不仅要依赖于专业的计算软件或程序， 业软件，应用和推广存在困难，因而一直不被重视。实际上，若 还需要处理大量数据（如机组特性曲线的处理）。为此，通过简 能寻求计算简单、思路清晰，则仅需常用软件就能求解水位涌波 化机组过流特性以推导形式简单的气垫式调压室涌波方程，并 方程的数值计算方法，其实用价值就能得到充分体现。 设法获得水位涌波计算的解析计算公式或实用便捷的求解方 与Ｗｏｒｄ一样，Ｅｘｃｅｌ也属微软ｏｆｆｉｃｅ中的常用办公软件，使 法，在工程运用中意义重大。 用非常普遍，且绝大多数设计技术人员都很熟悉。其特点是集 若假定机组引用流量在瞬间发生改变，则由动量方程、连续 数据处理、绘图和扩展编程为一体。由于数值计算的实质就是 方程、气体状态方程以及气垫式调压室的体型结构可以导出气 建立各计算时步之间的递推关系并逐步积分求解，而Ｅｘｃｅｌ软 垫式调压室水位涌波所满足的微分方程为： 件的自动填充功能可方便地描述此类递推关系，据此本文设计 ２ －ｍ ｋ·Ｌ ·ｚ″＝ＲＱ －ｚ－Ｐ［（１－ｋｚ／Ｖ） －１］－ 了采用Ｅｘｃｅｌ求解水位涌波方程的便捷方法。 １ ｇｆ ｔ １ １ １ １ 　　　　Ｒ｜ｋｚ′＋Ｑ ｜（ｋｚ′＋Ｑ）－Ｒ｜ｋｚ′｜ｋｚ′ ｔ １ ２ １ ２ ｓ １ １ 　 （１） １　水位涌波方程的数值解算模型 { Ｑ －Ｑ Ｌ ｚ（０）＝０，ｚ′（０）＝ １ ２，Ｌ ＝ ｋ ｇｆ ∑ｇｆ 常微分方程的数值解法有很多成熟的计算方法，如欧拉法、 １ 式中ｚ为以调压室初始水位为基准水位涌波，ｍ；Ｑ、Ｑ 分别为 龙格—库塔法等。为了保证计算精度以及方便求解，拟采用定 １ ２ ３ 步长四阶Ｒｕｎｇｅ－Ｋｕｔｔａ法求解水位涌波方程，并对此模型进行