有限单元法的基本知识和地震波传播正演模拟的应用.pptVIP

地球物理数值计算方法 地球物理数值计算方法 Numerical Methods in Geophysics 地球物理数值计算方法 第六章 有限单元方法 有限单元法 有限单元法简介 有限单元法简介 有限单元法简介 有限单元法简介 有限单元法简介 有限单元法简介 线性代数基础 线性代数基础 线性代数基础 线性代数基础 线性代数基础 线性代数基础 线性代数基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法基础 有限单元法-1D单元 有限单元法-1D单元 有限单元法-1D单元 有限单元法-1D单元 有限单元法-1D单元 有限单元法-1D单元 有限单元法-2D单元 有限单元法-2D单元 有限单元法-2D单元 有限单元法-2D单元 有限单元法-2D单元 有限单元法-2D单元 有限单元法-2D单元 有限单元法-2D单元 有限单元法-2D单元 有限单元法-2D单元 有限单元法-2D单元 有限单元法-2D单元 有限单元法-例子 有限单元法-例子 有限单元法-例子 有限单元法-例子 有限单元法-例子 有限单元法-例子 有限单元法-例子 有限单元法 声波方程： 对于特殊选取的基函数，系数值为离散节点上的函数值u。 求解时间相关项：时间微分采用有限差分近似： 得到tk+1时刻的解为： 我们知道如何计算刚度矩阵A，如何计算质量矩阵M？ 声波方程： M为： 基函数为： 一维问题： 刚度矩阵(stiffness matrix)： 方程中Aij的计算： 一维问题： 刚度矩阵(stiffness matrix)： 方程中Aij的计算： 一维问题： 边界条件和源项： 原始方程为： 如果假设： 最后变换为以下方程： 加上边界条件后，设： 其中u(0)和u(1)是区域[0，1]边界处的值。 一维问题： 边界条件和源项： 原始方程为： 加上边界条件后，设： 最后变换为以下方程： 可以写为矩阵形式： 一维问题： 边界条件和源项： 用图形形象表示（系统通过修改后的源项来考虑边界条件）： 边界条件 源项 边界条件 一维问题： 数值算例（等间距离散点）： 偏微分方程： 区域:[0,1];nx=100; dx=1/(nx-1);f(x)=δ(1/2); 边界条件：u(0)=u(1)=0 Matlab 有限单元程序 Matlab 有限差分程序 一维问题： 数值算例（等间距离散点）： 偏微分方程： 区域:[0,1];nx=100; dx=1/(nx-1);f(x)=δ(1/2); 边界条件：u(0)=u(1)=0 有限单元Matlab程序 计算结果（蓝色） 有限差分Matlab程序 计算结果（红色） 一维问题： 数值算例（等间距离散点）、边界条件不为零： 偏微分方程： 区域:[0,1];nx=100; dx=1/(nx-1);f(x)=δ(1/2); 边界条件：u(0)=0.15 u(1)=0.05 有限单元Matlab程序 计算结果（蓝色） 有限差分Matlab程序 计算结果（红色） 一维问题： 不等间距离散点时的刚度矩阵： 一维问题： 数值算例（不等间距离散点）： 偏微分方程： 区域:[0,1];nx=100; dx=1/(nx-1);f(x)=δ(1/2); 边界条件：u(0)=u0 u(1)=u1 刚度矩阵的计算程序 一维问题： 数值算例（不等间距离散点）： 偏微分方程： 区域:[0,1];nx=100; dx=1/(nx-1);f(x)=δ(1/2); 边界条件：u(0)=0.15 u(1)=0.05 总结： 1、有限单元方法的分析中，我们用一组正交基函数的线性组合来给出定义在区域D上的函数的近似值，线性组合的系数相当于某些离散节点上的函数值。 2、对某些偏微分方程组，离散节点上的值可以通过求解线性方程组获得，求解线性方程组得过程包括矩阵（有时是稀疏矩阵）的求逆。 3、边界条件在变换以后的方程中自然满足，这是有限单元方法的优势之一（和有限差分法相比） 坐标变换： 我们希望用一组基函数的线性组合来近似给出区