偏微分方程的有限元法.pptx

单击此处编辑母版标题样式,Harbin Institute of Technology Yangkun kyang,#,/43,计,算,物,理,学,单击此处编辑母版文本样式,第二级,Harbin Institute of Technology Yangkun kyang,1,/104,第五章 偏微分方程的有限元法,5.1,泛函与变分原理,5.2,基于变分原理的有限元法,5.3,matlab,有限元法工具箱,Harbin Institute of Technology Yangkun kyang,2,/104,第五章 偏微分方程的有限元法,有限元法（,FEA,，,Finite Element Analysis,，,FEM,）,有限元法的基本思想是用,较简单的问题,代替,复杂问题,，然后再对简单问题进行求解的数值计算方法。,有限元法将求解域看成是由许多被称为,有限元,的小的互连子域组成，对每一单元假定一个较简单的近似解，然后推导求解这个域总的满足条件，从而得到问题的解。这个解不是准确解，而是近似解。有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的数值计算方法。,有限元法于上世纪,50,年代首先在力学领域,-,飞机结构的静、动态特性分析中得到应用，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。,有限元法主要用于求解,拉普拉斯方程,和,泊松方程,所描述的各类物理场中。,Harbin Institute of Technology Yangkun kyang,3,/104,第五章 偏微分方程的有限元法,有限元法,-,变分原理,基于变分原理的有限元法是逼近论、偏微分方程、变分与泛函分析的巧妙结合。,基于变分原理的有限元法以,变分原理,为基础，把所要求解的,微分方程,定解问题，首先转化为相应的,变分问题,，即,泛函求极值,问题；它将求解域看成是由许多称为有限元的小的互连子域组成，然后利用,剖分插值,，对每一单元假定一个合适的,(,较简单的）近似解，把离散化的,变分问题,转化为普通多元函数的,极值问题,，然后推导求解这个域总的满足条件,(,边界条件），即最终归结为一组多元的,代数方程组,，求解代数方程组，就得到待求边值问题的数值解。,Harbin Institute of Technology Yangkun kyang,4,/104,第五章 偏微分方程的有限元法,有限元法,-,加权余数法,自从,1969,年以来，某些学者在流体力学中应用加权余数法中的迦辽金法或最小二乘法等同样获得了有限元方程，因而有限元法可应用于以任何微分方程所描述的各类物理场中，而不再要求这类物理场和泛函的极值问题有所联系。,加权余数法的核心思想是：近似解与解析解相比会存在误差,R,，但是可以通过一个准则使,R,尽量小，求解这个等式，就可以得到待定常数的值，也就得到了近似解。,Harbin Institute of Technology Yangkun kyang,5,/104,第五章 偏微分方程的有限元法,有限元法特点,有限元法的物理意义直观明确，理论完整可靠。因为,变分原理,描述了支配物理现象的物理学中的,最小作用原理,（如力学中的最小势能原理）。,优异的解题能力。有限元法对,边界几何形状复杂,以及,媒质物理性质变异,等复杂物理问题求解上，有突出优点：,不受几何形状和媒质分布的复杂程度限制。,不必单独处理第二、三类边界条件。,离散点配置比较随意，通过控制有限单元剖分密度和单元插值函数的选取，可以充分保证所需的数值计算精度。,Harbin Institute of Technology Yangkun kyang,6,/104,5.1,泛函与变分原理,数学上，通常自变量与因变量间的关系称为函数，而泛函则是函数集合的函数，也就是,函数的函数,，即自变量为函数，而不是变量。,5.1.1,泛函的定义,泛函通常是指一种定义域为函数，而值域为实数的“函数”。,设,C,是函数的集合，,B,是实数集合。如果对,C,中的任一元素,y(x),，在,B,中都有一个元素,J,与之对应，则称,J,为,y(x),的泛函，记为,Jy(x),。,Harbin Institute of Technology Yangkun kyang,7,/104,5.1,泛函与变分原理,例,5.1.1,质点在重力作用下，沿一条光滑的从,A,点到,B,点的曲线运动，如图所示。求下落时间最短的曲线。,曲线上任一小段线元长度为：,A,B,x,y,O,x,0,x,1,捷线问题,Harbin Institute of Technology Yangkun kyang,8,/104,5.1,泛函与变分原理,线元处的质点速度为,A,B,x,y,O,x,0,x,1,ds,线元下落时间为,从,A,点到,B,点的下落时间为,Harbin Institute of Technology Yangkun kyang,9,/104,5.1,泛函与变分原理,5.1.2,函数的变分,设,y(x),是泛函,J,定义域内任一函数，如果,y(x),变化为新函数,Y(x),，且,Y(x),属于泛函,J,的定义域，则,Y(x),与,y(x),之差为函数,y(x),的变分。,变分,y,是,x,的函数，它不同于函数的增量,y,。,性质：函数求导与求变分可以交换次序,Harbin Institute of Technology Yangkun kyang,10,/104,5.1,泛函与变分原理,5.1.3,泛函的变分,定义,最简泛函,F(x,y,y),称为泛函的“核函数”,泛函的变分,最简泛函,:,核函数只包含自变量,x,、未知函数,y(x),以及导数,y(x),Harbin Institute of Technology Yangkun kyang,11,/104,5.1,泛函与变分原理,利用二元函数的泰勒展开,Harbin Institute of Technology Yangkun kyang,12,/104,5.1,泛函与变分原理,其中,分别称为泛函的,一阶变分,和,二阶变分,。,Harbin Institute of Technology Yangkun kyang,13,/104,5.1,泛函与变分原理,泛函取极值的必要条件：,一阶变分为零,性质：对于最简泛函，变分运算可以与积分、微分运算交换次序,Harbin Institute of Technology Yangkun kyang,14,/104,5.1,泛函与变分原理,5.1.4,泛函的极值问题,泛函的一阶变分,利用,1,泛函的极值问题的间接,解法,转化为微分方程：欧拉方程,Harbin Institute of Technology Yangkun kyang,15,/104,5.1,泛函与变分原理,对于驻定问题，,两边界固定,这就是最简泛函的,欧拉方程,，等价于泛函取极值的必要条件。,把,变分问题转化微分方程的定解问题（边值问题）来求解,。,Harbin Institute of Technology Yangkun kyang,16,/104,5.1,泛函与变分原理,对于例,5.1.1,求下落时间最短的轨迹,利用最简泛函的,欧拉方程,。,Harbin Institute of Technology Yangkun kyang,17,/104,5.1,泛函与变分原理,代入欧拉方程,Harbin Institute of Technology Yangkun kyang,18,/104,5.1,泛函与变分原理,变换得到,进一步化简得到,积分,Harbin Institute of Technology Yangkun kyang,19,/104,5.1,泛函与变分原理,做变量替换,得,而,Harbin Institute of Technology Yangkun kyang,20,/104,5.1,泛函与变分原理,对上式积分得到,这样就得到了下落时间最短曲线的参数方程,式中常数,c,1,和,c,2,由始末两点位置确定,练习：画出经过,(0,0),和,(1,1),的下落时间最短曲线。,连接两个点上凹的唯一一段旋轮线,Harbin Institute of Technology Yangkun kyang,21,/104,5.1,泛函与变分原理,2,泛函的极值问题的直接解法,基本做法：,瑞利,-,里兹,(Rayleigh-Ritz),法,(1),选定一组具有相对完备性的基函数，构造一个线性组合的近似函数,(2),将含有,n,个待定系数的构造函数作为近似的极值函数，代入泛函,Harbin Institute of Technology Yangkun kyang,22,/104,5.1,泛函与变分原理,(3),为了求泛函的极值，按照多元函数取极值的必要条件,(4),求解以上方程组，求出 就可以得到极值函数的近似解,(5),再将含有,n+1,个待定系数的函数,作为近似极值函数，重复,(2)(4),，就可以得到极值函数,新的,近似解。如果连续两次所得到的结果接近，就认为最后得到的函数就是极值函数的近似解。,Harbin Institute of Technology Yangkun kyang,23,/104,5.1,泛函与变分原理,例,5.1.2,求下列泛函的极值函数。,解：为了满足边界条件，取基函数为,近似函数为,Harbin Institute of Technology Yangkun kyang,24,/104,5.1,泛函与变分原理,当,n=1,时,代入泛函,取极值,Harbin Institute of Technology Yangkun kyang,25,/104,5.1,泛函与变分原理,计算得到,近似函数,同理,n=2,时,利用欧拉方程，得到的精确解,Harbin Institute of Technology Yangkun kyang,26,/104,5.1,泛函与变分原理,Harbin Institute of Technology Yangkun kyang,27,/104,5.1,泛函与变分原理,泛函的极值问题可以通过变分运算产生一个微分方程和相应的边界条件，即欧拉方程，其解对应于最简泛函的极值函数。也就是,泛函的极值问题,可以等价为,在一定边界条件下求解微分方程问题。,变分原理,通过求解一个相应的泛函的极小函数而得到偏微分方程边值问题的解。,有限元法,正是,里兹法,与,有限差分法,相结合的成果，它取长补短地在理论上以变分为基础，在具体方法构造上又利用了有限差分法网格离散化处理的思想。,Harbin Institute of Technology Yangkun kyang,28,/104,5.1,泛函与变分原理,20,世纪,60,年代初首次提出结构力学计算有限元概念的克拉夫（,Clough,）教授形象地将其描绘为：“,有限元法,=Rayleigh Ritz,法,分片函数,”。,有限元法是,Rayleigh Ritz,法的一种局部化情况。不同于求解满足整个定义域边界条件的允许函数的,Rayleigh Ritz,法（往往是困难的），有限元法将函数定义在简单几何形状（如二维问题中的三角形或任意四边形）的单元域上（分片函数），且不考虑整个定义域的复杂边界条件，这是有限元法优于其它近似方法的原因之一。,Harbin Institute of Technology Yangkun kyang,29,/104,5.2,基于变分原理的有限元法,对于具有不同物理性质和数学模型的问题，有限元法的基本做法是相同的，只是具体公式推导和运算求解不同。,有限元法基本做法,首先把待求的偏微分方程边值问题转化为等价的变分问题。,然后通过有限单元剖分的离散处理，构造一个分片解析的有限元子空间。,通过构造近似函数，把变分问题近似地转化为有限元子空间中的多元函数极值问题，由此直接利用,Rayleigh Ritz,法探求变分问题的近似解（极值函数解），以此作为所求边值问题的近似解。,Harbin Institute of Technology Yangkun kyang,30,/104,5.2,基于变分原理的有限元法,有限元法具体求解步骤,建立积分方程,根据变分原理或方程余量与权函数正交化原理，建立与微分方程初边值问题等价的积分表达式，这是有限元法的出发点。,区域单元剖分,根据求解区域的形状及实际问题的物理特点，将区域剖分为若干相互连接、不重叠的单元。区域单元划分是采用有限元方法的前期准备工作，这部分工作量比较大，除了给计算单元和节点进行编号和确定相互之间的关系之外，还要表示节点的位置坐标，同时还需要列出自然边界和本质边界的节点序号和相应的边界值。,Harbin Institute of Technology Yangkun kyang,31,/104,5.2,基于变分原理的有限元法,确定单元基函数,根据单元中节点数目及对近似解精度的要求，选择满足一定插值条件的插值函数作为单元基函数。有限元方法中的基函数是在单元中选取的，由于各单元 具有规则的几何形状，在选取基函数时可遵循一定的法则。,单元分析,将各个单元中的求解函数用单元基函数的线性组合表达式进行逼近；再将 近似函数代入积分方程，并对单元区域进行积分，可获得含有待定系数,(,即单元中各节点 的参数值,),的代数方程组，称为单元有限元方程。,总体合成,在得出单元有限元方程之后，将区域中所有单元有限元方程按一定法则进行累加，形成总体有限元方程。,Harbin Institute of Technology Yangkun kyang,32,/104,5.2,基于变分原理的有限元法,边界条件的处理,一般边界条件有三种形式，对于第二类边界条件，一般在积分表达式中可自动得到满足。对于第一类边界条件和第三类边界条件，需按一定法则对总体有限元方程进行修正满足。,解有限元方程,根据边界条件修正的总体有限元方程组，是含所有待定未知量的封闭方程组，采用适当的数值计算方法求解，可求得各节点的函数值。,Harbin Institute of Technology Yangkun kyang,33,/104,5.2,基于变分原理的有限元法,有限元分析可分成三个阶段：,前置处理、计算求解和后置处理。,前置处理是建立有限元模型，完成单元网格划分；后置处理则是采集处理分析结果，使用户能简便提取信息，了解计算结果。,Harbin Institute of Technology Yangkun kyang,34,/104,5.2,基于变分原理的有限元法,1.,求解区域离散,离散单元基本要求：,各单元只能在顶点处相交。,不同单元在边界处相连，既不能相互分离又不能相互重叠。,各单元节点编号循序应一致，一律按,逆时针方向,，从最小节点号开始。同一单元节点编号相差不能太悬殊，对多区域的编号，按区域连续编号。,把求解区域分割成有限个单元体的集合。单元体形状原则上是任意的，一般取有规则形体。,有限元法计算步骤,Harbin Institute of Technology Yangkun kyang,35,/104,5.2,基于变分原理的有限元法,三角单元是经常使用的单元剖分方法，剖分时应注意几下几点：,三角形不能重叠。,不能把一个三角形的顶点取为相邻三角形的边上。,剖分的三角形应该避免钝角。,三角形不可过于狭长，最长边一般不大于最短边的,3,倍。,三角形三边之比尽量接近,1,。,不能把一个三角形跨越不同的介质。,每个三角形最多只有一个边在边界上。,三角形单元面积越小，计算精度越高,Harbin Institute of Technology Yangkun kyang,36,/104,5.2,基于变分原理的有限元法,把求解区域划分,m,个三角形有限单元，共有,n,个节点,在有限单元,e(j,k,l),上进行分片线性插值，插值函数为,2.,选择近似函数,Harbin Institute of Technology Yangkun kyang,37,/104,5.2,基于变分原理的有限元法,在单元节点上,求解以上方程组可以得到,3.,求解单元形函数,Harbin Institute of Technology Yangkun kyang,38,/104,5.2,基于变分原理的有限元法,同理可以求出,Harbin Institute of Technology Yangkun kyang,39,/104,5.2,基于变分原理的有限元法,Harbin Institute of Technology Yangkun kyang,40,/104,5.2,基于变分原理的有限元法,则插值函数可以写为,单元形函数（基函数）,Harbin Institute of Technology Yangkun kyang,41,/104,5.2,基于变分原理的有限元法,三角元,e,插值函数可以改写为矩阵形式,Harbin Institute of Technology Yangkun kyang,42,/104,5.2,基于变分原理的有限元法,下面以泊松方程为例讨论有限元解法,所对应的泛函为,4.,建立单元特征式,难点,：寻找与微分方程对应的泛函,Harbin Institute of Technology Yangkun kyang,43,/104,5.2,基于变分原理的有限元法,在第,e,个三角元的泛函,由于,Harbin Institute of Technology Yangkun kyang,44,/104,5.2,基于变分原理的有限元法,改写为矩阵形式,Harbin Institute of Technology Yangkun kyang,45,/104,5.2,基于变分原理的有限元法,其中,同理,Harbin Institute of Technology Yangkun kyang,46,/104,5.2,基于变分原理的有限元法,三角元,e,的泛函,其中,Harbin Institute of Technology Yangkun kyang,47,/104,5.2,基于变分原理的有限元法,改写,K,e,到所有,n,个节点，即把扩充部分添零，以方便总体矩阵的处理,其中,Harbin Institute of Technology Yangkun kyang,48,/104,5.2,基于变分原理的有限元法,求解区域上的总体泛函,其中,变分问题被离散化的多元二次函数的极值问题,5.,建立系统有限元方程,Harbin Institute of Technology Yangkun kyang,49,/104,5.2,基于变分原理的有限元法,根据多元函数极值理论,得到第,i,点有限元方程,即,求解上述有限元方程（线性代数方程组），就可以得到节点上的函数值。,Harbin Institute of Technology Yangkun kyang,50,/104,5.2,基于变分原理的有限元法,获得有限元方程之后，就可以选择各种方法求解相应的代数方程组，常用方法有高斯消去法、列元素消去法、迭代法等等。,在变分问题中第二类、第三类边界条件已经自然包含在泛函达到极值的要求中，不必单独处理，称为,自然满足的边界条件,，只需考虑,第一类强加边界条件,，强加边界条件的处理方法因代数方程组的解法而异。,6.,有限元方程求解与边界条件处理,Harbin Institute of Technology Yangkun kyang,51,/104,5.2,基于变分原理的有限元法,迭代法求解,：凡是遇到边界节点所对应的方程均不迭代，节点值始终保持给定值，不必单独处理边界。,直接法求解,：,节点,m,为边界，函数值,u,m,=u,0,，处理方法为，把对角元素的特征元素设置为,1,，即,k,mm,=1,然后把,m,行与,m,列的其它元素全部设置为,0,，方程的,等式右边,改为给定的函数值,u,0,，其它元素则要减去该节点处理前对应的,m,列的特征系数,k,im,与,u,0,的乘积。,Harbin Institute of Technology Yangkun kyang,52,/104,5.2,基于变分原理的有限元法,例,5.2.1,一个边长为,1,的二维正方形静电场域，电位函数为,(x,y),，边界条件如图所示，试用有限元法确定二维静电场域的电位分布。,解：该二维静电场域的电位函数,(x,y),可以用下列第一类边界条件的偏微分方程描述：,Harbin Institute of Technology Yangkun kyang,53,/104,5.2,基于变分原理的有限元法,按照右图进行三角形单元剖分，单元编号按照从左到右，从下到上的顺序编号。,节点编号：,1(0,0)2(0,1)3(0.5,0.5)4(1,0)5(1,1),三角形单元编号：,e(j,k,l),单元内顶点按逆时针编号,1(1,3,2)2(1,4,3)3(2,3,5)4(3,4,5),第一类边界条件,Harbin Institute of Technology Yangkun kyang,54,/104,5.2,基于变分原理的有限元法,对于三角元,1(1,3,2),Harbin Institute of Technology Yangkun kyang,55,/104,5.2,基于变分原理的有限元法,扩展到全部节点,Harbin Institute of Technology Yangkun kyang,56,/104,5.2,基于变分原理的有限元法,同样，对于三角元,2(1,4,3),Harbin Institute of Technology Yangkun kyang,57,/104,5.2,基于变分原理的有限元法,扩展到全部节点,Harbin Institute of Technology Yangkun kyang,58,/104,5.2,基于变分原理的有限元法,同样，对于三角元,3(2,3,5),Harbin Institute of Technology Yangkun kyang,59,/104,5.2,基于变分原理的有限元法,扩展到全部节点,Harbin Institute of Technology Yangkun kyang,60,/104,5.2,基于变分原理的有限元法,同样，对于三角元,4(3,4,5),Harbin Institute of Technology Yangkun kyang,61,/104,5.2,基于变分原理的有限元法,扩展到全部节点,Harbin Institute of Technology Yangkun kyang,62,/104,5.2,基于变分原理的有限元法,全部节点,K,Harbin Institute of Technology Yangkun kyang,63,/104,5.2,基于变分原理的有限元法,有限元方程,可以采用迭代法和直接解法，求解此线性代数方程组。,Harbin Institute of Technology Yangkun kyang,64,/104,5.2,基于变分原理的有限元法,迭代法：迭代公式为,代入初值,Harbin Institute of Technology Yangkun kyang,65,/104,5.2,基于变分原理的有限元法,直接解法：需要处理第一类边界条件,Harbin Institute of Technology Yangkun kyang,66,/104,5.2,基于变分原理的有限元法,差分法,得到差分递推公式,Harbin Institute of Technology Yangkun kyang,67,/104,5.2,基于变分原理的有限元法,差分网格,利用边界条件,Harbin Institute of Technology Yangkun kyang,68,/104,5.3 matlab,有限元法工具箱,大型通用有限元商业软件,国外软件,ANSYS,、,ADINA,、,ABAQUS,、,MSC,等,国内软件,FEPG,、,JFEX,、,KMAS,等,Matlab,偏微分方程工具箱（,PDE Toolbox,）提供了利用,有限元法、,图形界面,求解偏微分方程的计算环境。,PDEtool,有较大的局限性，可以求解特殊,PDE,问题，比如只能求解二阶,PDE,问题，并且不能解决偏微分方程组，但是它提供了,GUI,界面，从繁杂的编程中解脱出来了，同时还可以通过,File-Save As,直接生成,M,代码。,Harbin Institute of Technology Yangkun kyang,69,/104,5.3 matlab,有限元法工具箱,PDE Toolbox,求解偏微分方程类型,1,椭圆型方程（,Elliptic,）,2,抛物线型方程（,Parabolic,）,3,双曲型方程（,Hyperbolic,）,Harbin Institute of Technology Yangkun kyang,70,/104,5.3 matlab,有限元法工具箱,4,特征值方程（,Eigenmodes,）,上述微分方程中,c,、,a,、,d,、,f,在,椭圆型方程,中可以为,函数,，但在其它方程中必须为,常数,。,Harbin Institute of Technology Yangkun kyang,71,/104,5.3 matlab,有限元法工具箱,PDE Toolbox,边界条件,1,狄里赫利条件（,Didchlet,）,2,诺依曼条件（,Neumann,）,n,为边界上的单位外法线矢量，,h,、,r,、,q,、,g,可以为函数,Harbin Institute of Technology Yangkun kyang,72,/104,5.3 matlab,有限元法工具箱,PDE Toolbox,启动,1,启动,2,界面,Harbin Institute of Technology Yangkun kyang,73,/104,5.3 matlab,有限元法工具箱,PDE Toolbox,菜单,Options,打开或关闭栅格,调整栅格大小,打开或关闭捕捉栅格功能,绘图轴的坐标范围,打开或关闭绘图方轴,关闭帮助信息,图形缩放,选择应用模式,重新显示图形,Harbin Institute of Technology Yangkun kyang,74,/104,5.3 matlab,有限元法工具箱,PDE Toolbox,菜单,Draw,进入绘图模式,对角点绘矩形,固定中心绘矩形,矩形对角点绘椭圆,固定中心绘椭圆,绘多边形,旋转已选图形,将几何描述矩阵输出到主工作空间,Harbin Institute of Technology Yangkun kyang,75,/104,5.3 matlab,有限元法工具箱,PDE Toolbox,菜单,Boundary,进入边界模式,对已选边界输入条件,显示边界区域标识开关,显示子区域标识开关,删除已选的子域边界,删除所有 的子域边界,将分解几何矩阵、边界条件矩阵输出到主工作空间,Harbin Institute of Technology Yangkun kyang,76,/104,5.3 matlab,有限元法工具箱,PDE Toolbox,菜单,PDE,进入偏微分方程模式,显示子区域标识开关,调整,PDE,参数和类型,将,PDE,参数输出到主工作空间,Harbin Institute of Technology Yangkun kyang,77,/104,5.3 matlab,有限元法工具箱,PDE Toolbox,菜单,Mesh,输入网格模式,初始化三角形网格,加密当前三角形网格,优化网格,退回上一步,用数字化的颜色显示网格质量，大于,0.6,可接受,显示网格节点标识,显示三角形网格标识,修改网格生成参数,输出网格矩阵到主工作空间,Harbin Institute of Technology Yangkun kyang,78,/104,5.3 matlab,有限元法工具箱,PDE Toolbox,菜单,Solve,对已经定义的偏微分方程求解,调整解,PDE,的参数,输出解到主工作空间,Plot,显示图形解,绘图参数设置,输出动画,Harbin Institute of Technology Yangkun kyang,79,/104,5.3 matlab,有限元法工具箱,PDE Toolbox,求解步骤,求解区域设置,应用模式设置,输入边界条件,微分方程参数设定,网格剖分,初值和误差设置,解方程,图形解显示参数设置,File-Save As,直接生成,M,代码,Harbin Institute of Technology Yangkun kyang,80,/104,5.3 matlab,有限元法工具箱,例,5.3.1,如图带有矩形孔（,0.1,*,0.8,）的金属板（,1,*,1.6,），金属板左侧保持在,100,，右侧热量可以向环境定常流动，上下侧及内孔保持绝热，初始温度为,0,。求,t=0.1,、,0.3,、,0.5,、,1.5s,时金属板温度分布,解：此问题可以表示为如下定解问题,Harbin Institute of Technology Yangkun kyang,81,/104,5.3 matlab,有限元法工具箱,求解区域设置,提示符输入,pdetool,1,选择画矩形,2,画矩形,3,双击矩形，弹出对话框，,输入准确矩形参数,Harbin Institute of Technology Yangkun kyang,82,/104,5.3 matlab,有限元法工具箱,同样画出矩形孔，利用两个矩形运算得到求解区域,图形运算,Harbin Institute of Technology Yangkun kyang,83,/104,5.3 matlab,有限元法工具箱,2.,应用模式设置,Harbin Institute of Technology Yangkun kyang,84,/104,5.3 matlab,有限元法工具箱,3,输入边界条件,1,点击，显示边界,2,双击边界，弹出边界条件窗口,h=1,r=100,Harbin Institute of Technology Yangkun kyang,85,/104,5.3 matlab,有限元法工具箱,输入每个边的边界条件,红色：,Dirichlet,蓝色：,Neumann,Harbin Institute of Technology Yangkun kyang,86,/104,5.3 matlab,有限元法工具箱,4.,微分方程参数设定,1,点击，设置方程，弹出窗口,2,抛物线型,d=1 c=1 a=0 f=0,Harbin Institute of Technology Yangkun kyang,87,/104,5.3 matlab,有限元法工具箱,5.,网格剖分,点击，网格剖分,点击，加密网格,Harbin Institute of Technology Yangkun kyang,88,/104,5.3 matlab,有限元法工具箱,6.,初值和误差设置,单击,Solve,菜单中的,Paramenters,Harbin Institute of Technology Yangkun kyang,89,/104,5.3 matlab,有限元法工具箱,7.,解方程,点击，解方程,Harbin Institute of Technology Yangkun kyang,90,/104,5.3 matlab,有限元法工具箱,8.,图形解显示参数设置,Harbin Institute of Technology Yangkun kyang,91,/104,5.3 matlab,有限元法工具箱,不同时刻温度分布,Harbin Institute of Technology Yangkun kyang,92,/104,5.3 matlab,有限元法工具箱,不同时刻温度分布,Harbin Institute of Technology Yangkun kyang,93,/104,5.3 matlab,有限元法工具箱,MATLAB,除了提供有限元工具箱求解二阶偏微分方程之外，还提供了,pdepe,函数，可以直接求解,一维椭圆型和抛物线型偏微分方程或偏微分方程组,。,sol=,pdepe,(m,pdefun,pdeic,pdebc,xmesh,tspan),sol,：,是一个三维数组，,sol(:,:,i),表示,第,i,个微分方程,u,i,的解,。,xmesh,tspan,：,空间和时间离散向量,Harbin Institute of Technology Yangkun kyang,94,/104,5.3 matlab,有限元法工具箱,pdefun,：,PDE,描述函数,c,f,s=pdefun(x,t,u,du),m,x,t,就是对应于,偏,微分方程,中相关参数，,du,是,u,的一阶导数，由给定的输入变量即可表示出,c,f,s,这三个,输出,函数,。,Harbin Institute of Technology Yangkun kyang,95,/104,5.3 matlab,有限元法工具箱,pdebc,：,PDE,的边界条件描述函数,pa,qa,pb,qb=pdebc(xa,ua,xb,ub,t),其中,a,表示,下,边界，,b,表示,上,边界,。,pdeic,：是,PDE,的,初始,条件,描述函数,u,0,=pdeic(x),Harbin Institute of Technology Yangkun kyang,96,/104,5.3 matlab,有限元法工具箱,例,5.3.2,试求解下列的偏微分方程组,初始条件：,边界条件：,Harbin Institute of Technology Yangkun kyang,97,/104,5.3 matlab,有限元法工具箱,对照偏微分方程的标准形式，则方程组可以改写为,可见,m=0,Harbin Institute of Technology Yangkun kyang,98,/104,5.3 matlab,有限元法工具箱,Harbin Institute of Technology Yangkun kyang,99,/104,5.3 matlab,有限元法工具箱,下边界,上边界,边界条件,Harbin Institute of Technology Yangkun kyang,100,/104,5.3 matlab,有限元法工具箱,Harbin Institute of Technology Yangkun kyang,101,/104,5.3 matlab,有限元法工具箱,初始条件,Harbin Institute of Technology Yangkun kyang,102,/104,5.3 matlab,有限元法工具箱,Harbin Institute of Technology Yangkun kyang,103,/104,5.3 matlab,有限元法工具箱,Harbin Institute of Technology Yangkun kyang,104,/104,5.3 matlab,有限元法工具箱,Harbin Institute of Technology Yangkun kyang,105,/104,上机,5,利用有限元工具箱计算二维热传导问题,Harbin Institute of Technology Yangkun kyang,106,/104,上机,5,利用有限元工具箱，求解下列双曲线型偏微分方程,求解域,s,为,边界条件：构成求解域的边界值都为,5,初始条件：,