有限元与数值方法-讲稿2.ppt

* * * * * * * * * * * * * * * * * * * * A点需要满足的差分方程： 第二类边界条件： 转化为 （补充方程） 复杂区域形状（复杂边界）的处理 3 4 0 2 1 A 在A点展开 向后差分 * 一维非定常热传导混合问题的差分解 0 b t x t j x kh * 一维非定常热传导混合问题的差分解 可得到热传导方程的显式差分格式 采用差分格式 j+1 j k k-1 k+1 时间域和空间域同时进行差分法离散; 时间域的离散基于前一时刻的变化率 * 一维非定常热传导混合问题的差分解 特点 1.知道时刻j的温度Ｔ就可直接计算求得时刻j+1的温度 2. 当 r≤1/2 及满足相容性条件有条件的收敛和稳定 r值可从系数阵的特征值确定，对r的限制影响到时间步长和空间步长的匹配关系，空间步长小时时间步长也要小 初始条件 边界条件 j+1 j k k-1 k+1 * 采用差分格式 我们可以得到热传导方程的隐式差分格式 已知时刻 j-1 的温度Ｔ，为了计算时刻 j 的温度要解线性代数方程组 一般来说是无条件稳定的 初始条件 边界条件 t k k-1 k+1 j-1 j 一维非定常热传导混合问题的差分解 时间域和空间域同时进行差分法离散; 时间域的离散基于当前时刻的变化率 * 有限差分法的特点 有限差分法的基本概念非常简单。但实施会遇到困难： 均匀网格只适用于非常简单的矩形区域 即使非常简单的矩形网格，第二、第三类边界条件也需引入虚拟域外节点 复杂边界处理： 或是在边界上设置节点，导致复杂的边界差分格式 或是在域内采用非均匀的格式，差分格式的自动生成非常困难 可以想象三维问题的复杂,可以想象高价微分方程的复杂 一旦将系数矩阵输入计算机，线性方程组求解可以高度自动化，无须人为干预 如何适应各种复杂边界形状、边界条件，在计算机上自动形成系数矩阵是差分法的瓶颈 对于抛物型、双曲型方程,显式格式有稳定性问题,时间步长和空间格式需要协调；隐式格式虽然稳定,但效率往往很低 * 不规则网格上的差分格式 插值和差分格式 插值和差分在某种意义上可看作函数关系描写的两种方法: 离散到连续之间的正映射和逆映射. 离散描述 连续描述 插值 差分 这一关系使我们可以采用插值的方法来建立差分格式 设想已知函数在 0，h 的值为 f0，f1 ，则可以利用这两个点的值，插值得到一个线性函数： 这样,我们就得到了一阶导数的向前差分表达式。但是由这个线性函数无法求得二阶导数的差分格式 * 插值和差分格式 如果知道函数在三点的值： 利用这个办法也可以构造多元函数的偏导数的差分格式，但具体实施也许并不容易 * 用格林公式构造非规则网格的差分格式 Ａ １ ２ ３ ４ ５ 对于左边的不规则网格，假设场量 u 在单元内是 x 和 y 的线性函数，由 可得到近似等式： 线积分方向由右手定则给出 * 要求满足 在ΔA12内构造线性插值函数 假定原点放在A, 各点坐标为 u1 u2 u3 u4 u5 u0 ５ １ Ａ ２ ３ ４ Ω Г 用格林公式构造非规则网格的差分格式 * ５ １ Ａ ２ ３ ４ Ω Г 用格林公式构造非规则网格的差分格式 其中， 为三角形ΔA12的面积 故 * 由 用格林公式构造非规则网格的差分格式 ５ １ Ａ ２ ３ ４ Ω Г 按以上方法求出各边上的线积分后，可利用格林公式将二阶导数近似表示为不规则单元边界上的线积分形式： * * * * * * * * * * * * * * * * * * * * * 有限元与数值方法第二讲 授课教师：刘书田 Tel：座机电话号码； Email：stliu@dlut.edu.cn 教室：研究楼 351 时间：2013年3月22日：8：00 - 9:35; 9:45 -10：30 * 固体力学控制方程的基本形式（弹性力学为例） 一般方程 基本未知数（15）：应力（6）；应变（6）位移（3） 基本方程（15）： 平衡方程（3），本构方程（6），几何方程（6） 边界条件：应力边界条件；位移边界条件 位移解法 基本未知数：位移分量（3个） 基本方程（3）：用位移表示的平衡方程（3） 边界条件：应力边界条件；位移边界条件 应力解法 基本未知数（6）：应力（6） 基本方程：平衡方程（3），相容方程（3+3） 边界条件：应力边界条件；位移边界条件 * 基本控制方程：微分方程的边值问题 求响应量 平衡方程 几何方程 物理方程 （热弹性） （各项同性） * 求解方法及相对应的控制方程： 力法：以应力为基本位移量 平衡方程：变形协调方程 位移法：以位移为基本未知数 平衡方程：位移表示的平衡方程 对应原理：变分原理 微分方程的积分形式，泛函变分与基本方程的对应。建立各种