航空航天结构有限元法_6789.pdfVIP

第 6 章 变分原理与有限元法 6.1 微分方程的变分解法 连续介质的场问题，如应力场、温度场、电磁场等，在数学上可用偏微分方程或微分方 程组及其相应的边界条件和初始条件来描述。其定态问题，常常称为边值问题。它们的解是 在由已知边界条件所定义的区域中寻求的。但是，许多问题由于边界条件比较复杂，直接从 微分方程求精确解比较困难，甚至不可能。因此，有限元法在求解场问题上得到了广泛的应 用。 这里，我们将有限元法解释为求解改进的变分问题和加权余量问题的一种近似方法，以 便把有限元分析广泛地应用到求解场问题的领域中。 由于求解微分方程的边值问题和变分法（泛函求极值方法）具有等价性，故对复杂的微 分方程连同它的边界条件（自然边界条件），首先转化为求泛函的极值问题，然后采用有限 元的离散方法求泛函的极值，从而建立有限元方程求解。很多工程物理问题，它们变量之间 的关系，既可以用微分形式表达，也能找到它们的积分形式（泛函形式），其共同规律就是 能量积分。特别是这些问题在定常和平衡状态下，它们的势能最小，只是形式不同而已。因 此，求定常和平衡状态下的场量分布，都可以利用变分法求解。 利用经典的变分直接法求解偏微分方程，对于复杂的边界条件在选择近似函数时遇到困 难，于是将整个求解域进行剖分（分片插值），使变分解法又前进一步，这就是根据变分原 理发展而来的有限元法。它在求解椭圆型微分方程方面得到了广泛的应用。这里，我们首先 简要地介绍一下变分法的基本原理。 6.1.1 泛函极值求解与欧拉方程 如图 6.1，要求点A 到点 B 之间距离最短的一条曲线，这就是一个简单的变分问题。列 出A 、B 两点之间曲线的关系式为： 2 2 2 dL dx =+dy （6-1 ） 解出L 即得 x1 ′2 L ⎡y x ⎤ 1=+y dx （6-2 ） ⎣ ( )⎦ ∫x 0 式中：L 是 A 到 B 的曲线长度。由式（6-2 ）可以看出，y x 不同，则L ⎡y x ⎤也不同， ( ) ⎣ ( )⎦ y A(x , y ) 0 0 x o B (x , y ) 1 1 P (x, y ) y (x ) B (x , y ) 1 1 A(x , y ) 0 0