《第三章加权残值法.doc》由会员分享，可在线阅读，更多相关《第三章加权残值法.doc（11页珍藏版）》请在金锄头文库上搜索。
1、第三章加权残值法第三章 加权残值法加权残值法（Method of weighted Residuals）定解问题的近似求解方法。优点：原理统一，简便，工作量少，计算精度较高.加权残值法的发展：基本思想在19世纪初就已提出;20世纪20年代，由毕卡（Picone）用来求解微分方程；克兰德（Crandall)将这一方法统一，并定义为加权残值法。国内：20世纪60年代期间，最早由钱令希教授介绍了多种加权残值方法并用于分析薄板力学问题。徐次达教授自60年代开始利用加权残值法求解固体力学问题。3.1 加权残值法的基本概念设某一具体的工程定解问题：Luf=0（在域V内） （3.1。1）Gug=0（在边界S
2、上） （3。1.2）这里，u为待求的未知函数，L和G分别为控制方程（在域V内）和边界条件（在边界S上）的微分算子。f和g分别是域内和边界上的已知项.一般地，定解问题（3。1.1）、（3。1.2）的精确解难以求得，从而求助于近似解，这里我们假设一个待求函数u的试函数: （3。1。3）其中Ci为待定系数,vi为试函数项。将(3.1。3）代入定解问题的两个微分方程中，一般不会精确满足，于是就出现了内部残值（Residuals）RV和边界残值RS，即： (3.1。4) （3.1。5）为了消除残值，选取内部权函数（Weighted function）WV和边界权函数WS，使得残值RV和RS分别与相应权函
3、数的乘积在域内和边界上的积分为零，即: （3.1.6) (3.1。7）据此，我们就可以得到关于待定系数Ci(i=1，2，N)的代数方程组，求得了Ci后,即确定了近似解（3。1.3).按试函数是否满足控制方程和边界条件,将加权残值法分为三类，即内部法、边界法和混合法。这三种方法各有自己的优点，当然也存在不足.（1）在内部法中，对于一般比较规则的边界，选取满足边界条件的试函数是比较容易的.并且，由于边界条件已经满足，所以计长工作量较少.但是对于复杂的边界，这一方法就很不方便。(2)在边界法中，由于基本控制方程已经满足，近似计算仅在边界上进行,因而计算工作量少,精度较高，不足的是，要事先求得不同问题
4、控制方程的泛定解，比较困难。（3）混合法的优点在于，对试函数要求不严，复杂的边界条件和复杂的控制方程都能适应,缺点是计算工作量较大。总之,对于复杂控制方程,简单边界问题，宜采用内部法；对简单控制方程，复杂边界，适合用边界法;对控制方程和边界条件都较复杂的问题，采用混合法较好。这三种方法中，内部法一般应用较多。3。2 加权残值法的基本方法根据权函数的形式分类,主要有以下五种方法：（1)最小二乘法（Least Square Method）最小二乘法的基本思想是选取一个试函数，使得在域V内的残值平方积分: (3.2.1)最小。为使J（Ci）最小,取极值条件：,（i=1，2，N） (3。2。2)即可得
5、到最小二乘法的基本方程：，（i=1,2，N） （3。2.3）可见，最小二乘法就是将权函数取作。式（3.2。3）将给出N个代数方程，用于求解N个待定系数Ci（i=1，2，N）.这个方法一般计算精度高，但运算较为繁琐。（2）配点法（Collocation Method）如果选用狄拉克函数（Dirac Delta Function）作为权函数，即： （3。2.4）就得到了配点法.其中，函数又称单位脉冲函数，其具有以下性质： (3。2.5a) （3.2。5b) （3。2.5c） （3。2。5d)于是，将权函数（3。2。4）代入（3.1。6）中，便可得配点法的基本方程为:，（i=1，2,N） （3。2.
6、6)对于高维问题，例如二维问题的配点法基本方程为：,(i=1，2，N） （3。2。7）由残值R在N个配点xi（或二维（xi，yi)）处为零。得到N个代数方程，从而求得待定系数Ci(i=1，2，N）。配点法是加权残值法中最简单的一种,只是其计算精度相对差一些。（3）子域法（Subdomain Method）划分的子域总数应等于待定系数Ci的总数。如果将待求问题的整个区域V按任意方式划分为N个子域Vi（i=1，2,N）,并定义此时的权函数为: （3.2.8)于是在每个子域Vi内可列出消除残值的方程为：，(i=1，2，N） (3。2。9）这里,N个子域共有N个方程，联立求解即得待定系数Ci(i=1,
7、2，N）。需要说明的是，每个子域的试函数的选取可以相同，也可以不同.若各子域的试函数互不相同时，则必须考虑各子域间的连接条件。(4）伽辽金法（Galerkin Method）伽辽金法是俄国工程师伽辽金提出的并以他的名字而命名的方法。伽辽金法中的权函数就是试函数中的基函数，即：Wi=vi，（i=1,2，N） （3。2.10），（i=1，2，N） (3.2。11)由残值方程和试函数中的每一个基函数正交这一性质,不仅保证了解的收敛性,还使得伽辽金法精度高而计算工作量又不算太大，所以该方法应用广泛。（5）矩量法（Method of Moment）当权函数选取为xi（i=0,1，N1)时，就得到了矩量法
8、的基本方程为：，（i=0,1，N1） （3。2。12)由上式不难求得待定系数Ci（i=1，2，N）。至此，我们根据所选取的权函数类型，介绍了五种基本方法.在实际应用中，这五种基本方法可以单独使用，也可以相互结合而产生新的近似方法。3.4 加权残值法在力学中的应用本节将给出一些加权残值法求解力学问题的例子，使同学们能够熟练掌握应用这一方法求解具体力学问题的过程。例3.4.1：梁的弯曲问题考虑一两端固定，均布载荷作用下的直梁(图3.4。1),梁的跨度为L，均布载荷为q，梁的抗弯刚度为EI.图3。4.1 均布载荷作用下的固定梁解：梁的挠曲线微分方程为： （3.4.1)梁在两固定端所满足的边界条件为:
9、 (3。4。2a，b）（前面已经介绍过，按照试函数的类型，可将加权残值法分为三类：内部法、边界法和混合法。）1、内部法选取梁的挠度试函数为： （3。4。3)因此，所选取的梁的挠度试函数满足边界条件（3.4。2a，b），从而内部残值表达式为： （3.4。4）（a）最小二乘法求解最小二乘法相对应的残值方程为： （3。4.5)从而两端固定梁的挠度近似解为： （3。4。6)该解已经是材料力学中的精确解了。（b）配点法求解配点法相对应的残值方程为:R（xi）=24EICq=0 （3。4.7)xi为任意坐标都使得:此结果与最小二乘法所得结果完全相同。同理，我们还可以采用子域法，伽辽金法及矩量法，均可获得同
10、样的结果。2、混合法考虑到梁所满足的挠曲线方程中最高含有四阶导数，而且四阶导数的值为常数,这里,我们假设梁的挠度试函数为： （3.4.8)注意到，微分方程以及四个边界条件可得到5个方程，可以确定出（3.4。8)节中的5个待定系数Ci（i=0，1，2，3,4）。将（3。4。8）代入到(3。4。1）和（3.4。2a,b）中，得: （3.4.9)于是，挠度曲线函数为：与前面的方法所得结果一致.例3。4。2:简支梁的弯曲问题解：梁的挠曲线微分方程为： （3。4。1)边界条件为: (3.4.10）于是，把问题化为微分方程（3.4。1）和边界条件（3.4。10）的边值问题进行求解。首先凭经验选取试函数（Trial Function），以下采用低阶近似求解(所谓低阶近似，是指试函数中只含一个或几个待定参变量）。一阶近似选取的试函数为： （3。4.11）二阶近似选取的试函数为： （3.4。12）和已满足边界条件（3。4。10）,但不满足控制微分方程（3.4。1），将和代入(3。4。1）,得到内部
