有限元发展历史.docxVIP

有限元方法历史简介 数学有限元方法 （ FEM ）是用来求偏微分方程式 （ PDE）的近似解， 也求积分方程式，例如热传输方程式。 求解方法是基于完全取消微分方程式 （稳态问题）， 或把偏微分方程式（PDE ）译成等效的常微分方程式，然后采用像有限差等标准的技术求解。 在解偏微分方程式时，主要的挑战是创建近似研究的方程式，但数字稳定，这意味着在输入数据和中间计算都不会聚集错误，并造成无意义的输出结果。有许多这么做的方法， 它们都有各自的优缺点。 对于求解复杂域 （像汽车和油管道） 偏微分方程式，或当希望在全部范围精确变化时， 有限元方法是好的选择。 例如，在模拟地球气候模式时，在土地和完全 开放的海域之上有着准确的预测是非常重要的， 采用有限元方法， 这个要求是可以做得到的。 历史 有限元方法起源于需要解决市政工程和航空工程方面复杂的弹性结构分析问题。它的 开发可以追溯到 A.Hrennikoff （ 1941）和 R.Courant （ 1942）的工作。虽然这些先驱者使用这些方法， 并且引人注目的不同， 但他们都共享一个基本的特性： 把连续域的网格离散化进入一组离散的子域里。 Hrennikoff 的工作是采用格子使域离散，而与之类似，为了求解起源 于汽缸扭转的问题的二阶椭圆的偏微分方程式（ PDEs）， Richard Courant 的方法是把域划 分成有限的三角形子域。 对于由 Rayleigh ，Ritz 和 Galerkin 开发的偏微分方程式 （PDEs）， Richard Courant 的贡献是改进， 绘制了大量的早期结果。 针对机身和结构分析的有限元方法 的开发最早开始于 1950 年代中期，并且用于市政工程的有限元方法许多是 1960 年代在伯克 利开始启动（见伯克利早期有限元研究）。在 1973 年 Strang 和 Fix 出版的《有限元方法的 分析》 里，提供的方法采用了严格的数学基础，并且已经在广泛变化的工程学科， 即电磁和 流体力学里，针对物理系统的数字建模，归纳成为应用数学的分枝。 在结构力学里，有限元方法的开发常常是基于能量理论，即虚功原理或最小总潜能原理，对于结构工程师来说，早就强烈要求提供综合的，直觉的和物理的依据。 技术讨论 我们将从可以推断的普通方法里取二个简单问题来举例说明有限元方法。我们假设读者是熟悉微积分学和线性代数。我们将采用一维空间 式中 f 是假设的，而 u 是 x 的未知函数，并且 u〞是与 x 有关的 u 的二阶导数。二维空间取样问题是狄利克雷问题 式中 Ω 是在（x，y）平面内连接开区域， 那些边界是 “和谐的 ”（即平滑流形或多边形） ，并且 uxx 和 uyy，分别表示与 x 和 y 有关的二阶导数。 通过计算不定积分，可以 “直接 ”求解问题 P1。然而，只有当只有一个维度空间时，才 使用这个方法求解边界值问题，并且不推广到更高空间的问题，或像 u+u”=f问题。出于这 个原因，我们将针对 P1 开发有限元方法，并且略微叙述它对 P2 的广义性。 我们的解释将发生在二个步骤里，反映出二个本质的步骤，第一步必须采用有限元方 法（ FEM ）求助于求解边界值问题（ BVP ）。在第一步，在它的弱或变分形式上重新描述初始的边界值问题（ BVP ）。通常这一步几乎不需要作计算，只是在纸上手工进行转换。第二步是离散化，在有限的维度空间里，把弱形式离散化。在这个第二步之后，对于大的，但是 有限空间的线性问题，我们有具体的公式，那些解将近似解答初始的边界值问题（ BVP ）。 然后就在计算机里执行这个有限的空间问题。 变分公式化 第一步是把 P1 和 P2 转换为它们的变分公式。如果 u 求解 P1，那么对于任何平滑函数v，我们有 反过来，如果对于假设的 u，⑴控制每个平滑函数 v(t) ，那么一步就可以显示这个 u 将 求解 P1。（证据是非平凡的，并且采用 Sobolev 空间） 通过在⑴的右边采用部分积分法，我们获得 式中我们已经做了另外的假设 v(0)=v(1)=0 。 4.1 存在的证据概要和解的唯一性 我们可以定义 是有界变分的 (0，1)的函数， 在 x=0 和 x=1 是 0。这样的函数 是“一次可微分的 ”，并且它产生出相对称的双线性图 Φ，然后把定义的内积 转换 成为 Hibert 空间（详细的证据是非平凡的） 。在另一方面， 左侧 也是内积， 这次在 Lp 空间 L 2(0， 1)。针对 Hibert 空间的 Riesz 表示法则显示有一个唯一的 u 解⑵和因 此的 P1。 4.2 P2 的变分形式 如果我们采用 Green 的理论做部分积分， 我们看到如果 u 求解出 P2，那么对于任何 v： 式中 表示梯度，并且 ·表示二维