加权余量法简介43977

加权余量法简介43977

ID:24276614

大小:580.50 KB

页数:14页

时间:2018-11-13

加权余量法简介43977_第1页
加权余量法简介43977_第2页
加权余量法简介43977_第3页
加权余量法简介43977_第4页
加权余量法简介43977_第5页
资源描述:

《加权余量法简介43977》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库

1、加权余量法加权余量法的基本概念加权余量法(MethodofWeightedResiduals)或称加权残值法或加权残数法,是一种直接从所需求解的微分方程及边界条件出发,寻求边值问题近似解的数学方法。当n有限时,定解方程存在偏差(余量)。取权函数,强迫余量在某种平均意义上为零。采用使余量的加权积分为零的等效积分的“弱”形式。来求得微分方程近似解的方法称为加权余量法。设问题的控制微分方程为:在V域内在S边界上式中:L、B——分别为微分方程和边界条件中的微分算子;f、g——为与未知函数u无关的已知函数域值;u——为问题待求的未知函数。方法概述及按试函数

2、分类当利用加权余量法求近似解时,首先在求解域上建立一个试函数,一般具有如下形式:式中:——待定系数,也可称为广义坐标;——取自完备函数集的线性无关的基函数。由于一般只是待求函数u的近似解,因此将式代入控制方程后将得不到满足,若记:在V域内在S边界上显然反映了试函数与真实解之间的偏差,它们分别称做内部和边界余量。若在域V内引入内部权函数,在边界S上引入边界权函数则可建立n个消除余量的条件,一般可表示为:不同的权函数和反映了不同的消除余量的准则。从上式可以得到求解待定系数矩阵C的代数方程组。一经解得待定系数,即可得所需求解边值问题的近似解。由于试函数

3、的不同,余量和可有如下三种情况,依此加权余量法可分为:1.内部法试函数满足边界条件,也即此时消除余量的条件成为:2.边界法试函数满足控制方程,也即此时消除余量的条件为:3.混合法试函数不满足控制方程和边界条件,此时用下式来消除余量。显然,混合法对于试函数的选取最方便,但在相同精度条件下,工作量最大。对内部法和边界法必须使基函数事先满足一定条件,这对复杂结构分析往往有一定困难,但试函数一经建立,其工作量较小。无论采用何种方法,在建立试函数时均应注意以下几点:试函数应由完备函数集的子集构成。已被采用过的试函数有幂级数、三角级数、样条函数、贝赛尔函数、

4、切比雪夫和勒让德多项式等等。试函数应具有直到比消除余量的加权积分表达式中最高阶导数低一阶的导数连续性。试函数应与问题的解析解或问题的特解相关联。若计算问题具有对称性,应充分利用它。基本方法概述下面以内部法为例,介绍按权函数分类时加权余量的五种基本方法。对内部法来说,消除余量的统一格式是:1.子域法(SubdomainMethod)此法首先将求解域V划分成n个子域,在每个子域内令权函数等于1,而在子域之外取权函数为零,也即:如果在各个子域里分别选取试函数,那么它的求解在形式上将类似于有限元法2.配点法(CollocationMethod)子域法是令

5、余量在一个子域上的总和为零。而配点法是使余量在指定的n个点上等于零,这些点称为配点。此法的权函数为:P、Pi—分别代表求解域内任一点和配点。由于此法只在配点上保证余量为零,因此不需要作积分计算,所以是最简单的加权余量法3.最小二乘法(LeastSquareMethod)本法通过使在整个求解域上余量的平方和取极小来建立消除余量的条件。若记余量平方和为I(C),即则极值条件为:本法权函数为:4.伽辽金法(GalerkinMethod)本法是使余量与每一个基函数正交,也即以基函数作为权函数当试函数包含整个完备函数集时,用本法必可求得精确解。5.矩法(M

6、ethodofMoment)本法与伽辽金法相似,也是用完备函数集作权函数。但本法的权函数与伽辽金法又有区别,它与试函数无关。消除余量的条件是从零开始的各阶矩为零,因此对一维问题对二维问题其余类推这五种基本方法在待定系数足够多(称做高阶近似)时,其精度彼此相近。但对低阶近似(n较小)情况下,后三种的精度要高于前两种。基本方法举例为说明上述基本概念,以图所示等截面悬臂梁,受满跨均布荷载作用,求悬臂端B的竖向位移为例,说明基本方法的应用。图示梁的控制方程为:其边界条件为:若取试函数为:不难验证其满足边界条件,也即。而控制方程的内部余量为:子域法解由于试

7、函数仅一个待定常数,因此只需取一个子域(等于全域)即可,消除余量的条件为:由此可解得:代回原式可得:配点法解同上所述,只需选一个配点来建立消除余量的条件。若令:可得:若令:则得:可见不同的配点结果是不一样的。最小二乘法解此时消除余量的条件为:可得:伽辽金法解此时,消除余量的条件为:由此可得:本例各方法的精度比较本问题的精确解由梁位移计算可得为:由此可得,上述各方法对本例计算的误差依次为:-33.3%;1.75%(22.2%);13.9%;0.96%;-33.3%上面22.2%为式结果。

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。