有限元与数值法-讲稿3.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 残数（内部） 残数（边界） 考虑微分方程和边条 权函数 ( ) x W 加权残数法（Weighted Residual Method） * 此处，一个方程，n个未知数（C1~Cn） 加权残数法（Weighted Residual Method） * 选 n 个权函数 Wj （j=1~n） j=1~n n 个方程 求得C1~Cn 加权残数法（Weighted Residual Method） * 近似解构造方法 基函数系　　选择原则 连续性 线性无关 正交 完备 典型的基函数系 多项式 三角级数 梁振动振形 柱稳定函数 B-样条函数 通常取近似解为基函数的线性组合----基函数的选择方法 * 域内残数法 选取的基函数满足边界条件但不满足微分方程 边界残数法 选取的基函数满足域内微分方程但不满足边界条件 混合残数法 选取的基函数域内微分方程和边界条件都不满足 按基函数的性质进行分类 * 1.子域法 强迫余量在n个子域 的积分为零 n个方程，求得 C1~Cn 取 子域上近似 按权函数的性质进行分类 * 2.配点法 取 j 个方程 当子域法中，令面积→0，退化为配点法 * （最小二乘法的残数方程） (*) 对应每一点误差的平方和最小，即接近真解。 3.最小二乘法（Least Square Method） * 一次矩 二次矩 n 次矩 R的 j 次矩 4.矩法 伽辽金方程 把基函数作为权函数： 5.伽辽金法（Galerkin Method） 误差与解函数空间“正交” * 以上方法的比较 以上方法都将原问题转化为代数方程组的求解 Aa=c 配点法、子域法得到的是非对称的系数矩阵A； 最小二乘法、Galerkin法得到的是对称的系数矩阵A 最小二乘法易于产生病态矩阵A；并且不能通过分部积分法降低被积函数的微分阶次，因此要求单元间函数的充分的连续性 * 设 n=2 时 0 1 x 例题 即 余量为 * 1.子域法求解： x 0 1 0.5 例题 解方程，得到 不对称的系数矩阵 * 2.配点法： x 0 0.75 0.25 1 3.最小二乘法： 解方程，得到 解方程，得到 其中 该方程显然有对称的系数矩阵 * 4. Galerkin法： 经与精确解比较，Galerkin法结果具有较高精度 解方程，得到 * 例题：一维稳态热传导问题 取近似解为 加权余量法格式： * Galerkin法：取 ，代入上式中，得到 即 显然， 其中 * 结果的比较 * WRM推导虚功原理 三维弹性固体的平衡方程和边界条件： 取权函数为 ，则加权残数方程（等效积分形式）为 * WRM推导虚功原理 分部积分给出 引用虚位移（微小位移）与虚应变的关系及力的边界条件，上式可写为 此即虚功原理 弱形式 * WRM推导虚功原理 注解： 这里假定虚位移在域内连续可导，否则不能通过分部积分建立等效积分的弱形式 这里假定虚位移满足位移边界条件，否则外力虚功项中还应包括位移边界上约束反力的虚功 推导虚功原理的过程中，没有涉及本构关系，所以虚功原理可以用于非线性弹性及弹塑性等非线性问题 虚功原理表述了平衡条件 这里给出的虚功原理是基于小变形理论的，因此不能直接用于基于大变形理论的力学问题（对于大变形问题需要采用恰当的应力和应变度量） * 练习 推导下列方程的弱形式： 解： * 作业：平面应力问题的解 提示：u,v 用多项式做近似展开 求：位移和应力。用最小二乘法 参数 a=20mm b= 10mm * 作业解答 * * * * * * * * * * 有限元与数值方法第三讲微分方程的等价积分形式 授课教师：刘书田 Tel Email：stliu@dlut.edu.cn 教室：综合教学楼 351 时间：2013年3月29日：8：00—10：20 * 微分方程的等价积分形式 描述物理问题的三种数学工具： １ 微分方程边值问题（强式形式） 　　　　　　相应的数值方法是差分方法 ２ 能量原理（例如，最小总势能原理，最小余能原理） 　　　　　　相应的数值方法是有限元法 ３ 微分方程的弱形式 　　　　　　相应的数值方法是加权残数法、有限元法 　　　　　　和虚功原理可比较 　　　　　　系统如果没有能量也可使用 * 微分方程的算子形式 在域内： 边