水力学中常用的基本计算方法.doc

蒀袀螃膀薂蚃肂腿节衿羈膈莄蚁袄膈蒇袇螀芇蕿蚀肈芆艿蒃羄芅蒁蚈羀芄薃薁袆芃芃螆螂节莅蕿肁节蒇螅羇莁薀薇袃莀艿螃蝿荿莂薆肈莈薄螁肄莇蚆蚄羀莇莆衿袆羃蒈蚂螁羂薁袈肀羁芀蚁羆肀莃袆袂肀蒅虿螈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薁螅袅肅芁薈螁肄莃螄聿膄蒆薇羅膃薈螂袁膂芈薅袇膁蒀袀螃膀薂蚃肂腿节衿羈膈莄蚁袄膈蒇袇螀芇蕿蚀肈芆艿蒃羄芅蒁蚈羀芄薃薁袆芃芃螆螂节莅蕿肁节蒇螅羇莁薀薇袃莀艿螃蝿荿莂薆肈莈薄螁肄莇蚆蚄羀莇莆衿袆羃蒈蚂螁羂薁袈肀羁芀蚁羆肀莃袆袂肀蒅虿螈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薁螅袅肅芁薈螁肄莃螄聿膄蒆薇羅膃薈螂袁膂芈薅袇膁蒀袀螃膀薂蚃肂腿节衿羈膈莄蚁袄膈蒇袇螀芇蕿蚀肈芆艿蒃羄芅蒁蚈羀芄薃薁袆芃芃螆螂节莅蕿肁节蒇螅羇莁薀薇袃莀艿螃蝿荿莂薆肈莈薄螁肄莇蚆蚄羀莇莆衿袆羃蒈蚂螁羂薁袈肀羁芀蚁羆肀莃袆袂肀蒅虿螈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薁螅袅肅芁薈螁肄莃螄聿膄蒆薇羅膃薈螂袁膂芈薅袇膁蒀袀螃膀薂蚃肂腿节衿羈膈莄蚁袄膈蒇袇螀芇蕿蚀肈芆艿蒃羄芅蒁蚈羀芄薃薁袆 水力学中常用的基本计算方法水力学中经常会遇到一些高次方程，微分方程的求解问题。多年来，求解复杂高次方程的基本方法便是试算法，或查图表法，对于简单的微分方程尚可以用积分求解，而边界条件较为复杂的微分方程的求解就存在着较大的困难，但随着计算数学的发展及计算机的广泛使用，一门新的水力学分支计算水力学应运而生，但用计算机解决水力学问题，还需要了解一些一般的计算方法。在水力学课程中常用的有以下几种，现分述于后。 一、高次方程式的求解方法：（一）二分法1、二分法的基本内容：在区间X1，X2上有一单调连续函数F（x）=0，则可绘出F（x）关系曲线。如果在两端点处函数值异号即F（x）F（x2）0，（见图（一），则方程F（x）=0，在区间X1，X2之间有实根存在，其根的范围大致如下：取1若F（x2）F（x3）0， 则解X1，X3 2若F（x2）F（x3）0， 则解X3，X23若F（x2）F（x3）=0， 则解=X3 对情况1，可以令x2=x3，重复计算。 对情况2，可以令x1=x3，重复计算。 当规定误差之后，只要|x1-x2|，则x1(或x2)就是方程F(x)=0的根。 显然，二分法的理论依据就是高等数学中的连续函数介值定理。它的优点是思路清晰，计算简单，其收敛速度与公比为的等比级数相同；它的局限性在于只能求实根，而不能求重根。 2、二分法的程序框图(以求解明渠均匀流正常水深为例)最后必须说明，二分法要求x2值必须足够大，要保证F1F20，否则计算得不到正确结果。为了避免x2值不够大，产生计算错误，在程序中加入了判别条件F1F20。也可以给定xJ及步长x，让计算机选择x2(x2=x1+x)。(二)牛顿法，1、牛顿法的基本内容：设有连续函数F(x)=0，则可以绘出F(x)x关系曲线，选取初值xo，过点（xoF(xo)）作一切线，其斜率为辅F（xo），切线与x轴的交点是x1，则有：再过(x1，F(x1)作切线，如此类推得到牛顿法的一个迭代序列：xn+l=xn-F(xn)/F(xn)，令xn=xn+1，重复计算，直至满足给定的精度要求，即|xn+1-xn|，从而得到方程F(x)=0的根。牛顿法具有平方收敛速度，比较快，但计算工作量大，每次运算除计算函数值外，还要计算微商值。对于牛顿法来讲，只要F(x)在零点附近存在连续的二阶微商，是F(x)的一重零点，且初值xo充分接近于，那么牛顿迭代就一定收敛。2、牛顿迭代法的程序图框可以参看临界水深计算程序。(三)迭代法xn+1=(xn)xn=xn+11、迭代法的基本内容：设有一方程F（x）=0，可以写成等价形式x=(x)，再做迭代式xn+1=(xn),只要给定初始值xo，用式子 进行反复运算，就可以得到了个序列xn。如果序列xn收敛，并假设= xn+1=（xn）=()也就是说是方程F(x)=0的根。一般可以根据问题的性质人为地规定精度，认为当|xn+1-xn|时，xn+1或xn就是方程的根。y=(x)y=x从上面可以看出，迭代法的实质就是把方程F(x)=0的根问题转化为求迭代序列xn的极限问题，迭代过程的几何意义就是求两曲线的交点问题，交点处的横坐标正好就是方程的解。但是，用计算机进行迭代计算，最关键的问题就是要构造出正确的迭代计算公式；如果构造的迭代计算公式不适当，或者初值给的不合适，则都可能导致计算不收敛或者发生计算错误。例如，泄水建筑物下游收缩水深的计算，其基本数学公式是： 若进行迭代计算选用如下公式： （hc1为初值）则只要初值hc10，Eo，则计算结果均收敛于正确解。如若选用迭代公式，。只要初值hc10，则计算均收敛，但收敛于下游一个较大的淹没水深，并非正确解，发生计算错误。在非棱柱体明渠水面线计算中，同样存在这一问题。可见迭代公式选择的重要性。 迭代公式的收敛性可以用下述定理描述：把方程F(x)=0改写成x=(x)，如果(x)满足lipschitz条件：即对任意的x1和x2，都有|（x1）-（x2）|L|x1-x2|，其中L为一与x1和x2均无关的正常数(简称李氏常数)；若L1，则迭代收敛，收敛的准确解为，其收敛速度为：|xn-|高次方程的解往往不是一个，因此构造出迭代公式，并判定其收敛之后，尚要根据问题的背景及实际情况对收敛解进行分析讨论，确认为是真解方可，否则，尚需修改初值或重新构造迭代公式。2、迭代法的程序框图可以参看水跃水力计算程序。本程序集中迭代法应用较多，主要原因是迭代法思路简单明了，收敛较快，程序也短。二、插值法和数据拟合法插值法和数据拟合法都是根据给定的数据表(如水位流量资料)，设法寻找一个解析形式的函数(x)，来近似地代替这些数据间的函数关系f(x)。但插值法和数据拟合法之间又有不同之处。插值法要求求得的插值函数y=(x)在插值结点阵字库x1上都满足yi=(xi)，即要求插值函数曲线通过所有的插值结点(xi，yi)，但一般的实验数据中总是有观测误差存在，因此曲线y=(x)通过所有的实验点会使曲线保留全部的观测误差的影响，这是我们所不希塑的。数据拟合法则克服了这点不足，它不要求曲线(x)通过所有的点(xi，yi)，仅要求曲线和数据之间存在着较高的相关系数即可。关于插值法和数据拟合法的具体内容可参看有关方面的书籍。本程序集中，插值法及数据拟合法均有应用，例如，实例二中，上游水位与过水断面面积之间的关系就是应用分段线性插值得到它的分段函数表达式的；再如，较多的查图或者查表计算问题，我们都是应用数据拟合法求出它的函数式，应用的程序是熊运章教授的EDPS程序，分析得出的拟合曲线其相关系数均在0.99以上。三、数值积分法1、数值积分法的基本内容：在区间a、b上定义的黎曼可积函数f(x)，由于f(x)可能以表格形式给出，也可能f(x)的原函数无法用初等函数描述，因此直接积分存在着困难，但我们可以通过插值先构造一个多项式P(x)去逼近被积函数f(x)，并以P(x)在区间a、b上积分去代替f(x)在区间a、b上的积分，即，其中P(x)是比较容易积分的代数多项式。最常用的P(x)函数有如下两种形式：(1)P(x)为一条过两端点的直线，其函数形式是：，从右图可以看出，这实质上是用直角梯形面积去代替边梯形的面积，故称为梯形求积公式(2)P(x)为一条过两端点及其中点的抛物线，它的函数形式是： 从右图可以看出，这实质上是以抛物线围成的曲边梯形面积代替原f(x)所围成的曲边梯形面积，故称之为抛物线求积公式，也称辛浦生公式。为了提高精度，常将区间a，b分为足够多的N等分，然后求各个小部分的和。2、数积积分的程序框图可以看棱柱体明渠水面线计算程序。四、直接差分法直接差分法是以偏差商代替偏导数，把基本的微分方程化为差分方程，在自变量域st平面上建立差分网格，根据问题的初始条件及边界条件，求得各网格结点上的近似解。直接差分法的差分网格形式有好几种，不尽相同，各种形式所要求的步长也不一样，再者，差分方程还存在着相客性，收敛性，稳定性问题，比较复杂，读者可参考有关方面的专著。直接差分法的程序例子可以看清华大学水力学下册。五、有限单元法有限单元法的基本思想是将微分方程的定解问题(包括边值和初值问题)通过变分途径化为求泛函数的极值问题；然后将求解区域分割为许多小单元，通过区域剖分和分片插值，把求泛函数极值的问题又转化为代数方程组的求解问题，解此代数方程组就可以得到问题的解答。关于代数方程组的解法，对于线性代数方程组常见的有高斯消去法，LU分解法，迭代法等，对于非线性方程组常见的有牛顿迭代法等。多有现成程序可供利用。有限单元法的程序例子可以参看清华水力学下册。 芇芆螀螆莆莈薂肄莅蒁螈羀莄薃薁羆莃莃袆袂莂蒅虿膁莁薇袄肇莁蚀蚇羃莀荿袃衿肆蒂蚆螅肅薄袁肃肅芃蚄聿肄蒆罿羅肃薈螂袁肂蚀薅膀肁莀螀肆肀蒂薃羂腿薅蝿袈腿芄薂螄膈莇螇膃膇蕿薀聿膆蚁袅羅膅莁蚈袁膄蒃袄螆膃薆蚆肅芃芅袂羁节莈蚅袇芁蒀袀螃芀蚂蚃膂艿莂薆肈芈蒄螁羄芇薆薄袀芇芆螀螆莆莈薂肄莅蒁螈羀莄薃薁羆莃莃袆袂莂蒅虿膁莁薇袄肇莁蚀蚇羃莀荿袃衿肆蒂蚆螅肅薄袁肃肅芃蚄聿肄蒆罿羅肃薈螂袁肂蚀薅膀肁莀螀肆肀蒂薃羂腿薅蝿袈腿芄薂螄膈莇螇膃膇蕿薀聿膆蚁袅羅膅莁蚈袁膄