有限元及数值方法-讲稿3-4.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 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年4月07日：8：00—10：20 * 基于积分方程的数值方法的基本思想 微分提法：真实解在任意点均满足微分方程 积分提法：对于所有可能的解（u(x))中，真实的解应满足下式 积分形式的近似解法： 在有限个可能的解中，真实解的近似解为使下式取极小的解。 * 微分方程的算子形式 在域内： 边界上： 其中，A，B1，B2为微分算子 微分方程的等价积分形式 对于满足微分方程及其边界条件的解 u ,上式显然是成立的； 如果对任意的函数v(x)，上式成立，则可以证明u是微分方程的解。 * 　显然，如果在区域Ω上， 几乎处处为零，则对任意的　　　　有　 或采用数学的语言描述 上式成立的条件是要求函数可积。 引理:如果对任意的 　 ，恒有 则 如果F(X)=0代表了微分方程，则上面定理和引理建立了微分方程和其积分形式之间的联系 (一) 预备知识 微分方程的等价积分形式 * 等价积分形式 若对任意函数列向量 有 则该积分表达式与微分表达式 完全等效。 同理，若对任意函数列向量 有 则该积分表达式与微分表达式 完全等效。 故称 为原微分方程 的等价积分形式。 * 等价积分形式可积的条件： 1. 单值且在域内和边界上可积分 若 A 的最高阶导数为n，则u 的n-1 阶导数 必须连续，即u 具有 连续性 等价积分方程对函数连续性的要求：函数是可积的。 被积函