加权余量法简介
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
(2)试函数应具有直到比消除余量的加权积分表达式中最高 阶导数低一阶的导数连续性。
(3)试函数应与问题的解析解或问题的特解相关联。若计算 问题具有对称性，应充分利用它。
5.1.2 基本方法概述
下面以内部法为例，介绍按权函数分类时加权余量的五种基本 方法。对内部法来说，消除余量的统一格式是:
V WIiRI dV 0 (i 1, 2,L , n)
子域法是令余量在一个子域上的总和为零。而配点法是使
余量在指定的n个点上等于零，这些点称为配点。此法的权
函数为:
WIi （P Pi）
Dirac(犹拉克) 函数，它的定义为:
0 (x xi )
x xi x Βιβλιοθήκη Baiduxi
b
(x
a
xi )dx
0 1
xi a , b xi a , b
P、Pi—分别代表求解域内任一点和配点。 由于此法只在配点上保证余量为零，因此不需要作积分计算，
5.1 加权余量法的基本概念 加权余量法(Method of Weighted Residuals)或称加权残值
法或加权残数法，是一种直接从所需求解的微分方程及边 界条件出发，寻求边值问题近似解的数学方法。
早在20世纪30年代就在数学领域得到应用，随着计算机 的发展，它受到了国内外学者的普遍重视，得到了迅速的 发展。自1982年召开“全国加权残数法学术会议”后，我
前面所介绍的各种有限元法都是将整个求解域进行离 散，以单元为分析对象，设法建立单元的变量场，然后应 用能量原理或广义变分原理导出单元列式及整体分析方法， 从而解出用于建立单元变量场的基本未知量，进而求得其 他所需的物理量。
为解决工程实际计算还有一些其他数值方法，如加权余 量法、边界元法、样条有限元法、半解析法等，它们在计算 力学中形成了自己独特的理论和方法，内容也非常丰富，已 有大量文献资料和专著。本章只能对加权余量的基本概念、 方法和基本思路等作一简单介绍，为深入研究或进一步学习 打下必要的基础。
显然，混合法对于试函数的选取最方便，但在相同精度 条件下，工作量最大。对内部法和边界法必须使基函数事先 满足一定条件，这对复杂结构分析往往有一定困难，但试函 数一经建立，其工作量较小。
无论采用何种方法，在建立试函数时均应注意以下几点：
（1）试函数应由完备函数集的子集构成。已被采用过的试函 数有幂级数、三角级数、样条函数、贝赛尔函数、切比雪 夫和勒让德多项式等等。
所以是最简单的加权余量法
3．最小二乘法(Least Square Method) 本法通过使在整个求解域上余量的平方和取极小来建立消除余 量的条件。
若记余量平方和为I(C)，即 I (C) V RI2dV V RIT RIdV
则极值条件为:
I (C) C
2
( RI V C
)T
RI dV
0
由此可见，本法权函数为:
WIi
RI Ci
(i 1, 2,L , n)
4．伽辽金法(Galerkin Method) 本法是使余量与每一个基函数正交，也即以基函数作为权函数
WIi Ni (i 1, 2,L , n)
当试函数 u%包含整个完备函数集时，用本法必可求得精确解。
5．矩法(Method of Moment) 本法与伽辽金法相似，也是用完备函数集作权函数。
国加 权余量法在结构分析领域内的应用已从静力发展到动力、 稳定、材料非线性和几何非线性等各方面。
大量的结构分析问题，如杆系结构分析、二维及三维 弹性结构分析，以及板、壳应力分析等等，都可归结为在 一定的边界条件(或动力问题的初始条件)下求解微分方程 的解，我们称这些微分方程为问题的控制方程。
下面以加权余量法的数学模型和基本方法两个方面来介 绍加权余量法的基本概念。
f g
在V域内 在S边界上
(5.1.4)
显然 RI 、 RB反映了试函数与真实解之间的偏差，它们分别称
做内部和边界余量。
若在域V内引入内部权函数 WI ，在边界S上引入边界权函数 WB 则可建立n个消除余量的条件，一般可表示为：
V WIi RI dV S WBi RBdS 0 (i 1, 2,L , n)
V WIiRI dV 0 (i 1, 2,L , n)
(5.1.6)
2．边界法 试函数满足控制WI方i 程，也即 RI L(u%) f 0 此时消除余量的条件为:
SWBi RBdS 0 (i 1, 2,L , n)
(5.1.7)
3．混合法 试函数不满足控制方程和边界条件，此时用式(5.1.5)来消除余量。
(5.1.5)
不同的权函数 WIi 和 WBi 反映了不同的消除余量的准则。从上 式可以得到求解待定系数矩阵C的代数方程组。一经解得待定 系数，由式(5.1.3)即可得所需求解边值问题的近似解。
由于试函数 u%的不同，余量 RI 和 RB 可有如下三种情况， 依此加权余量法可分为： 1．内部法 试函数满足边界条件，也即 RB B(u%) g 0 此时消除余量的条件成为:
式中：
n
u% Ci Ni NC i 1
(5.1.3)
Ci —— 待定系数，也可称为广义坐标；
Ni ——取自完备函数集的线性无关的基函数。
由于 u% 一 般只是待求函数u的近似解，因此将式(5.1.3)
代入式(5.1.1)和式(5.1.2)后将得不到满足，若记：
RI RB
L(u%) B(u%)
1．子域法(Subdomain Method) 此法首先将求解域V划分成n个子域 Vi ，在每个子 域内令权函数 等于1，而在子域之外取权函数为零，也即:
1 WIi 0
(Vi内) (Vi外)
如果在各个子域里分别选取试函数，那么它的求解在形式上将 类似于有限元法。
2. 配点法(Collocation Method)
5.1.1 方法概述及按试函数分类
设问题的控制微分方程为:
在V域内
L(u) f 0
(5.1.1)
在S边界上
B(u) g 0
(5.1.2)
式中 :
L、B——分别为微分方程和边界条件中的微分算子；
f、g ——为与未知函数u无关的已知函数域值；
u——为问题待求的未知函数。
当利用加权余量法求近似解时，首先在求解域上建立一个试函数 u% ， 一般具有如下形式：
(3)试函数应与问题的解析解或问题的特解相关联。若计算 问题具有对称性，应充分利用它。
5.1.2 基本方法概述
下面以内部法为例，介绍按权函数分类时加权余量的五种基本 方法。对内部法来说，消除余量的统一格式是:
V WIiRI dV 0 (i 1, 2,L , n)
子域法是令余量在一个子域上的总和为零。而配点法是使
余量在指定的n个点上等于零，这些点称为配点。此法的权
函数为:
WIi （P Pi）
Dirac(犹拉克) 函数，它的定义为:
0 (x xi )
x xi x Βιβλιοθήκη Baiduxi
b
(x
a
xi )dx
0 1
xi a , b xi a , b
P、Pi—分别代表求解域内任一点和配点。 由于此法只在配点上保证余量为零，因此不需要作积分计算，
5.1 加权余量法的基本概念 加权余量法(Method of Weighted Residuals)或称加权残值
法或加权残数法，是一种直接从所需求解的微分方程及边 界条件出发，寻求边值问题近似解的数学方法。
早在20世纪30年代就在数学领域得到应用，随着计算机 的发展，它受到了国内外学者的普遍重视，得到了迅速的 发展。自1982年召开“全国加权残数法学术会议”后，我
前面所介绍的各种有限元法都是将整个求解域进行离 散，以单元为分析对象，设法建立单元的变量场，然后应 用能量原理或广义变分原理导出单元列式及整体分析方法， 从而解出用于建立单元变量场的基本未知量，进而求得其 他所需的物理量。
为解决工程实际计算还有一些其他数值方法，如加权余 量法、边界元法、样条有限元法、半解析法等，它们在计算 力学中形成了自己独特的理论和方法，内容也非常丰富，已 有大量文献资料和专著。本章只能对加权余量的基本概念、 方法和基本思路等作一简单介绍，为深入研究或进一步学习 打下必要的基础。
显然，混合法对于试函数的选取最方便，但在相同精度 条件下，工作量最大。对内部法和边界法必须使基函数事先 满足一定条件，这对复杂结构分析往往有一定困难，但试函 数一经建立，其工作量较小。
无论采用何种方法，在建立试函数时均应注意以下几点：
（1）试函数应由完备函数集的子集构成。已被采用过的试函 数有幂级数、三角级数、样条函数、贝赛尔函数、切比雪 夫和勒让德多项式等等。
所以是最简单的加权余量法
3．最小二乘法(Least Square Method) 本法通过使在整个求解域上余量的平方和取极小来建立消除余 量的条件。
若记余量平方和为I(C)，即 I (C) V RI2dV V RIT RIdV
则极值条件为:
I (C) C
2
( RI V C
)T
RI dV
0
由此可见，本法权函数为:
WIi
RI Ci
(i 1, 2,L , n)
4．伽辽金法(Galerkin Method) 本法是使余量与每一个基函数正交，也即以基函数作为权函数
WIi Ni (i 1, 2,L , n)
当试函数 u%包含整个完备函数集时，用本法必可求得精确解。
5．矩法(Method of Moment) 本法与伽辽金法相似，也是用完备函数集作权函数。
国加 权余量法在结构分析领域内的应用已从静力发展到动力、 稳定、材料非线性和几何非线性等各方面。
大量的结构分析问题，如杆系结构分析、二维及三维 弹性结构分析，以及板、壳应力分析等等，都可归结为在 一定的边界条件(或动力问题的初始条件)下求解微分方程 的解，我们称这些微分方程为问题的控制方程。
下面以加权余量法的数学模型和基本方法两个方面来介 绍加权余量法的基本概念。
f g
在V域内 在S边界上
(5.1.4)
显然 RI 、 RB反映了试函数与真实解之间的偏差，它们分别称
做内部和边界余量。
若在域V内引入内部权函数 WI ，在边界S上引入边界权函数 WB 则可建立n个消除余量的条件，一般可表示为：
V WIi RI dV S WBi RBdS 0 (i 1, 2,L , n)
V WIiRI dV 0 (i 1, 2,L , n)
(5.1.6)
2．边界法 试函数满足控制WI方i 程，也即 RI L(u%) f 0 此时消除余量的条件为:
SWBi RBdS 0 (i 1, 2,L , n)
(5.1.7)
3．混合法 试函数不满足控制方程和边界条件，此时用式(5.1.5)来消除余量。
(5.1.5)
不同的权函数 WIi 和 WBi 反映了不同的消除余量的准则。从上 式可以得到求解待定系数矩阵C的代数方程组。一经解得待定 系数，由式(5.1.3)即可得所需求解边值问题的近似解。
由于试函数 u%的不同，余量 RI 和 RB 可有如下三种情况， 依此加权余量法可分为： 1．内部法 试函数满足边界条件，也即 RB B(u%) g 0 此时消除余量的条件成为:
式中：
n
u% Ci Ni NC i 1
(5.1.3)
Ci —— 待定系数，也可称为广义坐标；
Ni ——取自完备函数集的线性无关的基函数。
由于 u% 一 般只是待求函数u的近似解，因此将式(5.1.3)
代入式(5.1.1)和式(5.1.2)后将得不到满足，若记：
RI RB
L(u%) B(u%)
1．子域法(Subdomain Method) 此法首先将求解域V划分成n个子域 Vi ，在每个子 域内令权函数 等于1，而在子域之外取权函数为零，也即:
1 WIi 0
(Vi内) (Vi外)
如果在各个子域里分别选取试函数，那么它的求解在形式上将 类似于有限元法。
2. 配点法(Collocation Method)
5.1.1 方法概述及按试函数分类
设问题的控制微分方程为:
在V域内
L(u) f 0
(5.1.1)
在S边界上
B(u) g 0
(5.1.2)
式中 :
L、B——分别为微分方程和边界条件中的微分算子；
f、g ——为与未知函数u无关的已知函数域值；
u——为问题待求的未知函数。
当利用加权余量法求近似解时，首先在求解域上建立一个试函数 u% ， 一般具有如下形式：