pde中微元分析法的主要思想_第五章偏微分方程的有限元法解析.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Harbin Institute of Technology Yangkun kyang@hit.edu.cn */104 5.3 matlab有限元法工具箱 pdebc： PDE 的边界条件描述函数 [pa,qa,pb,qb]=pdebc(xa,ua,xb,ub,t) 其中 a 表示下边界，b 表示上边界。 pdeic：是PDE 的初始条件描述函数 u0=pdeic(x) Harbin Institute of Technology Yangkun kyang@hit.edu.cn */104 5.3 matlab有限元法工具箱 例5.3.2 试求解下列的偏微分方程组 初始条件： 边界条件： Harbin Institute of Technology Yangkun kyang@hit.edu.cn */104 5.3 matlab有限元法工具箱 对照偏微分方程的标准形式，则方程组可以改写为 可见 m=0 Harbin Institute of Technology Yangkun kyang@hit.edu.cn */104 5.3 matlab有限元法工具箱 Harbin Institute of Technology Yangkun kyang@hit.edu.cn */104 5.3 matlab有限元法工具箱 下边界 上边界 边界条件 Harbin Institute of Technology Yangkun kyang@hit.edu.cn */104 5.3 matlab有限元法工具箱 Harbin Institute of Technology Yangkun kyang@hit.edu.cn */104 5.3 matlab有限元法工具箱 初始条件 Harbin Institute of Technology Yangkun kyang@hit.edu.cn */104 5.3 matlab有限元法工具箱 Harbin Institute of Technology Yangkun kyang@hit.edu.cn */104 5.3 matlab有限元法工具箱 Harbin Institute of Technology Yangkun kyang@hit.edu.cn */104 5.3 matlab有限元法工具箱 Harbin Institute of Technology Yangkun kyang@hit.edu.cn */104 上机 5 利用有限元工具箱计算二维热传导问题 Harbin Institute of Technology Yangkun kyang@hit.edu.cn */104 上机 5 利用有限元工具箱，求解下列双曲线型偏微分方程 求解域s为 边界条件：构成求解域的边界值都为5 初始条件： * * * * * * * * * * * * * * * 对于任意的积分区间，积分为零，则被积函数一定为零。 * * * * * Harbin Institute of Technology Yangkun kyang@hit.edu.cn */104 5.2 基于变分原理的有限元法 有限元方程 可以采用迭代法和直接解法，求解此线性代数方程组。 Harbin Institute of Technology Yangkun kyang@hit.edu.cn */104 5.2 基于变分原理的有限元法 迭代法：迭代公式为 代入初值 Harbin Institute of Technology Yangkun kyang@hit.edu.cn */104 5.2 基于变分原理的有限元法 直接解法：需要处理第一类边界条件 Harbin Institute of Technology Yangkun