FEM03-变分原理基础cai.pptVIP

有限元方法与应用变分原理基础 本节要点 有限元列式技术 加权残量法 配点法 最小二乘法 迦辽金法 变分原理 虚位移原理 虚应力原理 最小位能原理 最小余能原理 1 加权残量法与变分原理简介 有限元法不直接从问题的微分方程和相应的定解条件出发，而从与其等效的积分形式出发去求解原问题的近似解 加权残量法 加权残量法是等效积分的一般形式，适用于普通微分方程的求解 变分原理 若原问题具有某些特定的性质，则其等效积分形式可归结为某个泛函的变分，相应的近似解法转化为求泛函的驻值 2 微分方程的等效积分形式 物理问题 边界条件 域Ω内 边界Γ上 未知函数 微分方程组 Ω内 Γ上 (1) (2) 2 微分方程的等效积分形式 例1：一维弹性静力问题 ΓP上 Γu上 Ω内 且 精确解 特例 (3) 将以此为例说明典型列式技术的应用 2 微分方程的等效积分形式 例2：一维稳态热传导问题 Γq上 ΓT上 Ω内 2 微分方程的等效积分形式 物理问题(1)的等效积分形式 边界条件(2)的等效积分形式 原问题等效积分形式 Ω内 Γ上 (4) (5) (6) 3 加权残量法 未知函数近似解 残差（残量） Ω内 Γ上 (7) (8) 其中 试探函数（形函数） 待定系数 3 加权残量法 原问题近似等效积分形式 Ω内 Γ上 (9) 或 (9)式为通过选择待定系数，强迫残量在某种平均意义上等于零 采用使残量的加权积分为零来求得微分方程近似解的方法称为加权残量法 若将(9)式进行分部积分，则可获得等效积分的“弱”形式， 其可降低对近似函数连续性的要求 3 加权残量法 加权残量法求特例问题的一般列式 且 精确解 取试探函数 试探函数满足边界条件，但不满足微分方程 残量加权积分为零，即 由于试探函数满足边界条件，故边界条件不出现在等效积分表达式里面 3 加权残量法 加权残量法求特例问题的一般列式 且 精确解 残量为 残量为 3 加权残量法－配点法 权函数的选取 配点法求特例问题 (10) 3 加权残量法－配点法 配点法求特例问题 3 加权残量法－最小二乘法 权函数的选取 最小二乘法求特例问题 (11) 一项近似解： 3 加权残量法－最小二乘法 最小二乘法求特例问题 两项近似解： 3 加权残量法－迦辽金法 权函数的选取 迦辽金法求特例问题 (12) 所得到的系数矩阵是对称的，因此在用加权残量法建立有限元列式时几乎都采用迦辽金法。 一项近似解： 3 加权残量法－迦辽金法 迦辽金法求特例问题 两项近似解： 3 加权残量法－讨论 加权残量法的特点 加权残量法可用于广泛的方程类型 选择不同的权函数，可以产生不同的加权残量法，总体上迦辽金法的精度较配点法和最小二乘法要好 若近似函数取自完全的函数系列，并满足连续性要求，则当试探函数项数不断增加时，近似解可趋近于精确解 加权残量法的不足 解的收敛性没有严格的理论证明 近似解通常不具有明确的上下界性质 4 变分原理 内积 转化后的内积 边界条件 4 变分原理 泛函的构造（线性自伴随微分方程） 等效迦辽金提法 Ω内 Γ上 (13) Ω内 Γ上 利用自伴随线性算子的性质可得 则原问题的变分原理为 4 变分原理－虚位移原理 弹性力学平衡方程及力边界条件 其中 应力 体积力 微分算子 应变 面力 4 变分原理－虚位移原理 等效积分形式 虚位移原理（即弱形式，由分部积分获得） 虚位移原理的力学意义：若力系是平衡的，则它们在虚位移上所作之功的总和为零。反之，若力系在虚位移（及虚应变）上所做之功等于零，则它们一定满足平衡。即虚位移原理表述了力系平衡的充分必要条件。 虚位移原理不仅可用于线弹性问题，而且可以用于非线性弹性及弹塑性等非线性问题。 (14) 4 变分原理－虚应力原理 弹性力学几何方程及位移边界条件 其中 位移 微分算子 应变 指定位移 4 变分原理－虚应力原理 等效积分形式 虚应力原理（即弱形式，由分部积分获得） 虚应力原理的力学意义：若位移是协调的，则虚应力和虚边界约束反力在它们上面所作之功的总和为零。反之，如果上述虚力系在它们上面所作之功的和为零，则它们一定是满足协调的。即虚应力原理表述了力系平衡的充分必要条件。 虚应力原理可用于线弹性问题以及非线性弹性及弹塑性等不同力学问题。 (15) 4 变分原理－最小位能原理 (16) 系统总位能 弹性变形位能 外力位能 最小位能原理的力学意义：在所有区域内连续可导的并在边界上满足给定位移边界条件的可能位移中，真实位移使系统总位能最小。 可用于线弹性问题的求解。 其中 最小位能原理 由虚位移原理 4 变分原理－最小余能原理 (17) 系统总余能 弹性余能 外力余能 最小余能原理的力学意义：在所有区