偏微分方程的有限元方法67页PPT
- 格式:ppt
- 大小:6.11 MB
- 文档页数:67
16.901讲义笔记一维有限%首先,我们考虑•个比上一节稍微复杂点的问题; 豎二f(X),卫冲,V(O) = O.V(L)=O在这里,f(X)是)C的般函数,我们来看•个特别的情形:f(x)=x(L-x),此时,方程的梏确解如F:有限元方法利用加权残差的方法■其中:(1)设va)=£«Ma), v()()是我们对v(x)的近似,省为未知常数9 V|(x)是用户选择的歯数,即形状朗数:(2)定义N个加权残差LRj = p^(x)R(V)dx • j = l-> N to其中,RV)二器・f为绒差凹⑴足“用户”选择的加权函数,即权函数:(3)令加权残并为冬•町以确定⑷的值,即求耳使得对所fi 1=I->N, Rj=Oe令限元方法( )是加权残若法的一种,下血看看我们是如何用它来解决问题的。
一维有限元方法有限元方法(〉扌野个连续区域离散化-系列小单尤,这些单元与有限差分法()或有限体积法()产牛的网格完全相同,而佼之前两者主耍的优点在于:能够容易地把握单元的变化范囤。
对于我们讨论的一维问题,可以将区域(数轴〉离散化为如下图所示:这里,叫三单•元的个数。
我们还会用別下血i些定义:个三角划分;尽管令限元法对于一维,二维,三维甚至高细问题都是仃效的,们我们还是要谈及区域离散化的一种方浓,即三角划分。
4 T定义为第I个单元所在的区域。
对于_维问题,这表明,TS-个满足片心的X的集合。
接卜来耍确定的是毎个单兀该用什么样的函数,典型的函数形式就是用从一个单元到卜一个单兀保持解连续的多项式。
例如:一个线性有限元如卜團;i示:在毎个单元内的函数是线形的,在毎两个单元的交点处足连续的。
对于专门诜择的满足线件变化的形状函数,右估计残差时有一个很明显的问题:回忆前曲的内容,RV)二器一f,它在一个单冗里等于什么呢?因为函数是线性的,所以器=0,则有:R(V)=f ,即R(V)与无关。
冋时,满足线性变化的形状函数似乎也是一个好的近似,我们举-个例子来说明。