[理学]有限单元法课后习题全部答案王勖成

东南大学土木工程学院考博专业课复习资料

习题1.2：

在用有限元法求解时，边界条件总是满足的，控制方程的不完全匹配，会产生误差。题中所，代入边给出的近似函数：φ=a0+a1x+a2x+a3x，应该满足边界条件，对于情况（1）

2

3

1 a1L a2L2

界条件可得，从而 =a00,=a3

3L

x3x3x32

φ=a1(x 2)+a2(x )+3 （1）

LLL

x3

上式中的最后一项3前面没有待定系数，这是由于使用了在x=L处φ=1的强制边界条件。

L

从物理意义上说，相当于给定边界条件的解为齐次方程的通解加一个特解的缘故。将（1）式代入教材（1.2.26）式，得到残量：

R(x)=a1( 6

x6x6x)+a(2 )+3+Q(x) 22LLL

不同的求解方法，如配点法、子域法和伽辽金法，只是残量在某种意义上某个区域加权积分

为零。

配点法强制残量R(x)在有限个点严格为零，点的个数取决于未知数个数，这里为2，通常取所选的点在域内均匀分布，则取x=L/3和x=2L/3处，R(x)=0,这样得到

L2L

R()0,R()0，从而可以解出待定系数a1,a2。带入（1）式可以得到φ。

33

配点法仅考虑了有限个点的局部特性，子域法则要求在有限个子域 i内残量的积分

∫

i

R(x)dx=0为零，子域的个数仍然取决于未知函数个数，通常选取各子域的并集为整个

待求区域，一般情况可以选择各子域大小相同，但对于某些局部变化较复杂的区域，可以缩小子域的大小，使得子域分布更合理。例如取子域为

={x|0≤x≤L/2}, ={x|L/2≤x≤L}，则利用==12∫R(x)dx0,∫R(x)dx0，

1

2

可以求出待定系数a1,a2。

伽辽金法作为加权余量法的特殊形式，权函数选择为插值函数N1,N2，

x3x32这里N1(x)x 2,N2(x)x ，这样，利用∫Ni(x)R(x)dx=0,=i1,2可以求出待

LL

定系数a1,a2。

对于其余边界条件情况可依此类推。

东南大学土木工程学院考博专业课复习资料

练习题1.4，注意近似函数要满足边界条件，从而可知截面及坐标系如图所示：

，很多同学把积分区域弄错了，也有不少同学计

算错误。这里，由于边界为零，采用泛函及其弱形式得到的积分结果是相同的。最终计算得到：a1=4608/(13π4), a2=-512/(15π4), a3=-1536/(85π4)。

练习题1.5，泛函的欧拉方程基本没太多问题，泛函为零得到边界条件：

dw dw dw

0 dx2δ dx dx3δw =

0

1.5 如有一问题的泛函为=Π(w)

23

L

∫

L

EI d2w 2kw2

+qw dx，其中E, I, k是常数，q 2 +

dx22

是给定函数，w是未知函数，试导出原问题的微分方程和边界条件.

=δΠ(w)

∫

L

d2w d2w

++EIkwwqwδδδ 2 2 dx

dxdx

L

∫

L

L d2w d2δw d2wd(δw)d3wd(δw)

∫EI3EI dxEI2dx2 2 0dxdxdxdxdxdx 0

Ld2wd(δw)d3wd4w

=EI2 EI3δw+∫EI4δwdx

0dxdx0dxdx0

L

L

=δΠ(w)

∫

L

d2w d2w

δw dx EI 2 δ 2 +kwδw+q=

dxdx dwd(δw)dw

EI δw23

dxdx0dx0

2

L

3

L

∫

L

d4w

++EIkwq δwdx4

dx

+EI

d4w

微分方程：EI+kw+q=0 4

dx

d2wd2wd3wd3w

边界条件：0， 0 2233

dxx0dxxLdxx0dxxL

东南大学土木工程学院考博专业课复习资料

分强制边界和自然边界。

补充题

试作加权余量发的最小二乘配点法，并给出所得到的求解方程系数矩阵的特点分析。 （最小二乘配点法思路是，利用使求解域内所选各点处误差平方的总和为最少的条件，去建立求解试函数系数的方程。配点法是强迫余量误差在所选点上为0，最小二乘配点法则是余量在所选点上的误差，满足平方和最小。） 解：近似函数为u(x)=Ni(x)ai，不失一般性

余量为：R(x)=A(u) f(x)=A(Ni(x)ai) f(x) 最小二乘配点法取权函数

wjA(Niai)δ(x xk)其中j=1,...,n; k=1,...,m 且m≥n aj

加权余量要求

∫

wjRd =0

=jRd ∫ w

T

i(x)ai)δ(x xk)[A(Ni(x)ai) f(x)]d ∫ ajA(N

=∫AT(Nj(x))δ(x xk)[A(Ni(x)ai) f(x)]d

=

=

∑{A

k=1m

m

T

(Nj(xk))[A(Ni(xk)ai) f(xk)]

m

}

(Nj)A(Ni)ai ∑AT(Nj)f

=k1=k1

T

∑A

=Ka-P (写成矩阵形式)

因此，kij=

=A(N)A(N)∑

T

j

i

k=1

m

kji， 系数矩阵对称，且无需积分。

复习题1.7

自然边界条件强制边界条件的区别何在?为什么这样命名？对于一个给定的微分方程，如何区分这两类边界条件？

自然边界条件与强制边界条件，二者都是针对边值条件来说的。边值条件一般有三类边界条件。第一类：狄里克莱(Dirichlet)条件；第二类，诺依曼(Neumann)条件；第三类，前两者的混合条件，也叫洛平(Robin)条件

在选择近似函数时，已经事先满足的边界条件为强制边界条件。而自然边界条件则是在将等效形式化为弱形式时包含在边界积分场上的边界条件。

对于2m阶微分算子，含0到m-1阶导数的边界条件称为强制边界条件，近似函数应该事先满足。含m到2m-1阶导数的边界条件称为自然边界条件，近似函数不必事先满足。对于给定的微分方程，判断其阶次，再依据边界所含导数阶数可区分两类边界。

思考题1.8 泛函在什么条件下有极值？了解泛函是否有极值的意义何在？

=δΠ0且或δ2Π>

<0， 泛函极值性对于判断解的近似性质有意义，利用它可以对解的

东南大学土木工程学院考博专业课复习资料

上下界做出估计。

思考题1.9 什么是里兹法？通过它建立的求解方法有什么特点？里兹方法收敛性的定义是什么？收敛条件是什么？ 里兹法：在某一函数空间寻找试探函数，利用加权值的独立变分性将该函数的驻值问题转化为该函数关于权值的极值问题。其特点是：试探函数是全域的，解的精度依赖于试探函数的选取，其收敛性有明确的结论。

收敛性意义：当在∞维空间中选取试探函数，当试探函数的数目趋于∞时，利用里兹法得到的近视解将收敛于精确解。

收敛条件：1完备性，2试探函数满足Cm 1连续性

思考题1.0 里兹法的优缺点？举例说明

优点：理论简单，收敛性有严格的理论基础，得到的求解方程的系数矩阵是对称的，在场函数事先满足强制边界条件情况下，解具有上下界性质。 缺点：当求解域的形状很不规则时候，里兹法所要找的试探函数难以满足全部的强制边界条件，这样会降低精度。另外，由于其是基于变分原理，对于没有等价泛函的问题无法处理。

习题1.6 两端简支弹性基础上的梁受均不载荷。

=Π(w)

∫

L

EI d2w 2kw2 …… 此处隐藏：8494字，全部文档内容请下载后查看。喜欢就下载吧 ……