中小型天然河道水面曲线可视化计算系统的开发研究

1、中小型天然河道水面曲线可视化计算系统的开发研究任翔 刘川顺 （水资源与水利水电工程国家重点实验室，武汉大学，武汉430072）摘要 在天然河道治理和堤防除险加固中，都要通过水力计算得出水面曲线，从而确定河道的水力要素。本文探讨了中小型天然河道水力要素计算的实用方法，并采用VB6.0语言编制出天然河道水面曲线可视化计算系统，应用该系统可简化工程设计人员的工作难度，大大提高设计效率。关键词 中小型天然河道；水力计算；水面曲线；软件；VB6.0；CAD1 概述目前计算中小型天然河道水面曲线的常用方法是试算法、图解法和迭代法，试算法计算工作量相当大，多用于计算流量不变和流段较少的情况，对于流量变化和流

2、段较多的情况，采用图解法可以节省一定的工作量。图解法大致分为两大类，一类是考虑流速水头和局部水头损失的图解法，该法过程繁琐；另一类是忽略流速水头和局部水头损失的图解法，如艾斯考夫法，该法为了免去计算的难度，略去某些次要水力因素的影响，必将带来很大的累计误差，然而采用二分迭代法能很好地避免工作量大、精确度不高的缺点。同时考虑到现有的中小型天然河道水面曲线计算程序不能完全顾及河道断面、流量的沿程变化和水工建筑物（如桥墩）带来水面壅高变化的影响，不能自动绘制成果图。为此本水面曲线可视化计算系统力图与实际工程需要相结合，将既有的计算程序发展到能考虑流量、糙率沿程变化及地形变化和挡水建筑物带来的水面壅高

3、等综合影响，并且自动绘制出河道水面纵横剖面计算成果图，提高了工作效率。2 数学模型及算法2.1 计算方法天然河道的水流一般是非恒定非均匀流。但经验表明，天然河道水力要素随时间的变化是很缓慢的。因此，我们可以近似地认为天然河道的水流运动在一定时间内是恒定非均匀流，这种近似满足中小型天然河道规划设计阶段对水力计算的要求。计算恒定流一般采用伯诺里方程，通过二分法试算求解，即： (1)式(1)中：符号的下表1、2分别代表同一河段的下游断面和上游断面； 为沿程水头损失，；为局部水头损失，； 为局部阻力系数，对于收缩河段采用0，对于扩散河段，可以根据实际情况选取，如果出现两个试算水位均小于下游水位的情况，

4、表示局部阻力系数选用不当，计算机将每次自动给此系数增加“-0.2”，再进行计算,直到得出合理结果。2.2 二分法原理采用二分法试算求解，将伯诺里方程变为： (2)式(2)中：H为上游断面的试算水位，则当Z0时，H即为其解。试算时，将H的范围取为高于河底而低于断面最高点，如果河水超过断面最高点时，断面需要加高。 二分法的具体作法是：先计算河道上半部分，即初设水位为断面最高点，得到Z值，并令Z0=Z，次设水位为河底与断面最高点之间的中点，又得Z值，用Z0Z作判别式，如Z0Z<0，表示水面在此区间，如Z0Z>0，表示水面在另一区间，根据判别结果，对新区间继续二分和继续判别，逐步

5、缩小搜索区，直至试算的两个水位差在允许误差范围内为止，从而得出该断面的水位。2.3 河道横断面计算即计算过水面积A和湿周。计算方法为将过水断面分成若干梯形，分别计算每个梯形的面积和湿周后，叠加而成。3 可视化设计界面及计算机绘图功能的开发3.1 界面程序编制要点本系统采用VB6.0编制了可视化操作界面。首先输入工程名称、起推水位和所计算的横断面个数，选择流量、糙率特征以及是否具有挡水建筑物，如界面I。（略）其次输入各断面的河底高程、所测地形点数、各断面距上一断面的距离、动能修正系数、局部阻力系数、流量和糙率，如界面II。（略）接着根据提供的河道横断面示意图，输入各断面地形点的坐标(Xi，Yi，

6、Xi)，其中Yi为该点的高程，Xi和Xi分别为到中心点0的距离，这样充分考虑了河道地形变化带来的影响，如界面III。（略）其中在界面III中提供了CheckBox控件用来选择该断面是否具有挡水建筑物，选中该控件，将自动跳出挡水建筑物基本资料的填写，通过输入桥墩个数和桥墩直径，从而考虑挡水建筑物壅水带来的影响，如界面。（略）最后为计算控制界面，首先检查输入的数据是否正确（可以参考各横断面图，如果有误直接点击“修改数据”按钮），最后进行计算，并且自动绘制出河道水面曲线成果图，如界面。（略）为了使五个界面构成有机整体并与核心计算程序之间建立联系，在界面程序中设置四个模块，在第一个模块中定义了一系列的

7、公用类型（Public）变量，用来将当前界面中输入的数据赋给该模块中的相应变量并保存起来，使在操作过程中已经输入的数据不会丢失。另外三个模块为ConnectCAD()、AddLine和AddText，其功能在接下来的计算机绘图功能的开发中将加以详述，这样给设计者在进行试算时带来了很大的便利。3.2 计算机绘图功能的开发计算机自动绘图功能的开发，将减轻设计人员的工作强度，提高设计效率。本设计系统采用过程ConnectCAD()与AutoCAD实现连接，用AcadApp ActiveDocument. ModelSpace的AddLine、AddText、AddDimAligned属性分别实现绘制

8、直线、添加文字和进行尺寸标注，使系统能够快速准确的绘制出天然河道水面线成果图。其中，过程ConnectCAD用于连接AutoCAD，程序代码如下：Public ACADAPP As AcadApplication '检查运行AUTOCADPublic Sub 连接autocad()On Error Resume NextSet ACADAPP = GetObject(, "AUTOCAD.APPLICATION")If Err Then Err.Clear Set ACADAPP = CreateObject("AUTOCAD.APPLICATION&qu

9、ot;) If Err Then MsgBox ("不能运行AUTOCAD,请检查是否安装了AUTOCAD") Exit Sub End IfEnd IfACADAPP.Visible = True : ACADAPP.Documents.AddEnd SubAddLine属性用于绘制直线，程序代码如下：Public Sub DrawLine(X1 As Variant, Y1 As Variant, X2 As Variant, Y2 As Variant, Lcolor As Integer) Dim Lineobj As AcadLine : Dim Pstart(0

10、 To 2) As DoubleDim Pend(0 To 2) As Double Pstart(0) = X1 : Pstart(1) = Y1 : Pend(0) = X2 : Pend(1) = Y2 Set Lineobj = ACADAPP.ActiveDocument.ModelSpace.AddLine(Pstart, Pend) Lineobj.Color = LcolorZoomExtentsEnd Sub AddText属性用于文本编辑，程序代码如下：Public Sub WriteFont(P_Iinsert As Variant, Text_String As Str

11、ing, Text_Roation As Single, SIZE As Variant) 'CAD中写字Dim Styobj As AcadTextStyleDim Typeface As String:Dim Bold As BooleanDim Italic As Boolean : Dim CharSet As LongDim PitchandFamily As Long:Dim Textobj As AcadTextDim Insertp(0 To 2) As Double:Dim Height As DoubleSet Styobj = ACADAPP.ActiveDocu

12、ment.TextStyles.Add("M_Font1")Typeface = "宋体" : Italic = True : CharSet = 1PitchandFamily = 1 Or 16Styobj.SetFont Typeface, Bold, Italic, CharSet, PitchandFamilyACADAPP.ActiveDocument.ActiveTextStyle = StyobjInsertp(0) = P_Iinsert(0) : Insertp(1) = P_Iinsert(1) : Height = SIZESet

13、 Textobj = ACADAPP.ActiveDocument.ModelSpace.AddText(Text_String, Insertp, Height)Textobj.Rotation = Text_RoationTextobj.UpdateEnd SubAddDimAligned属性用于尺寸标注，程序代码如下：Private Sub Command1_Click()Dim dimObj As AcadDimAlignedDim point1(0 to 2) As Double : Dim point2(0 to 2) As Double : Dim location(0 to 2

14、) As DoubleSet dimObj=AcadApp.ActiveDocument.ModelSpace.AddDimAligned(point1, point2,location)End Sub4 结语经比较，本系统计算结果用广东省梅州市水电设计院的手算结果进行了验证，比较符合工程的实际情况。此外，在进行水面曲线计算时，糙率和局部阻力系数对计算成果的影响较大，在有资料的河段，可以根据水文资料来推求河段的平均糙率，在无资料的地区，可以结合实际情况选用。河道的分段原则上要求该段河道内的形状、糙率、底坡的变化不要太大，若不能满足时，计算断面可适当加密。 参 考 文 献1 徐正凡主编.水力学M

15、：北京：高等教育出版社，19862 (日)吉川秀夫著.水力学M：北京：北京冶金工业出版社，19863 刘润生主编.水力学M：上海：上海交通大学出版社，19825 张俊华编. 河道整治及堤防管理M：河南：黄河水利出版社，19986 张晋西. Visual Basic与AutoCAD二次开发M.北京：清华大学出版社，2002THE DEVELOPMENT AND RESEARCH TO VISUAL CALCULATE SYSTEM OF MEDIUM AND LITTLE-SIZED RIVERS SURFACE HYDROGRAPHRen xiang , Liu Chuan-shun(Stat

16、e Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 430072 , China)Abstract As for river administration and flood protection, it is essential to get the rivers surface hydrograph and ascertain some rivers hydraulic factors. Therefore, this paper research some applicative methodologies about medium and little-sized rivers hydraulic factor. At the same time, introduce a river surface