[由RayleighRitz法到FEM.docVIP

由Rayleigh-Ritz法到FEM 1. 引言 Rayleigh-Ritz法中，有一定的局限性：选取位移试函数必须满足位移边界条件，对复杂边界问题，很难做到。此缺陷使得Rayleigh-Ritz法难以应用于复杂工程问题的求解。 有限元方法（Finite Element Method，FEM）可以有效解决Ritz法的缺陷。它将求解区域划分为若干“单元”（Element），通过单元上的插值来构造位移函数的近似解。这样构造出的位移近似解就可以很容易地满足复杂边界上的位移边界条件。所以，FEM可以看成是Ritz方法的推广。 FEM和Ritz方法的实现过程基本相同，即包含以下三个主要步骤： ① 构造位移近似解，其中包含待定系数——这一步是FEM与Ritz法的主要区别所在，也是FEM的关键； ② 应用最小势能原理，获得关于待定系数的线性方程组； ③ 求解线性方程组，得到待定系数，代入位移近似解即可得到位移解，进而可得到应变、应力。 下面就按上述步骤来实现平面问题的FEM求解过程。 图1 考虑图1所示的平面弹性体。由于形状复杂，无法按照Ritz法的要求给出在内部连续、在位移边界上满足位移边界条件的位移近似解。FEM的重要突破在于：将求解区域划分为若干个element。以element为单位可以很容易地构造位移近似解。 2. 单元剖分 FEM的第一步是：将求解区域（弹性体）划分为e