有限单元法课后习题全部的答案王勖成.pdf

习题1.2： 在用有限元法求解时，边界条件总是满足的，控制方程的不完全匹配，会产生误差。题中所 给出的近似函数：φ a =+a x +a x2 +a x3 ，应该满足边界条件，对于情况（1），代入边 0 1 2 3 2 1−a L −a L 界条件可得a0 0, a3 1 3 2 ，从而 L 3 3 3 x 2 x x φ a (x =− ) +a (x − ) + （1） 1 2 2 3 L L L 3 x 上式中的最后一项 前面没有待定系数，这是由于使用了在x=L 处φ=1 的强制边界条件。 L3 从物理意义上说，相当于给定边界条件的解为齐次方程的通解加一个特解的缘故。将（1） 式代入教材（1.2.26）式，得到残量： x 6x 6x R (x ) a (=−6 ) +a (2 − ) + +Q(x) 1 2 2 3 L L L 不同的求解方法，如配点法、子域法和伽辽金法，只是残量在某种意义上某个区域加权积分 为零。 配点法强制残量R(x)在有限个点严格为零，点的个数取决于未知数个数，这里为2 ，通常取 所选的点在域内均匀分布，则取x=L/3 和x=2L/3 处，R(x)=0,这样得到 L 2L R ( ) 0, R ( ) 0 ，从而可以解出待定系数a , a 。带入（1）式可以得到φ 。 1 2 3 3 配点法仅考虑了有限个点的局部特性，子域法则要求在有限个子域Ωi 内残量的积分 ∫Ω R (x )dx 0 为零，子域的个数仍然取决于未知函数个数，通常选取各子域的并集为整个 i 待求区域，一般情况可以选择各子域大小相同，但对于某些局部变化较复杂的区域，可以缩 小 子 域 的 大 小 ， 使 得 子 域 分 布 更 合 理 。 例 如 取 子 域 为 Ω {x | 0 =≤x ≤L / 2},Ω {x |L / 2 =≤x ≤L } ，则利用 R (x )dx 0, R (x )dx 0 ， 1 2 ∫Ω ∫Ω 1 2 可以求出待定系数a , a 。 1 2 伽辽金法作为加权余量法的特殊形式，权函数选择为插值函数N , N