理论教育 确定单元基函数-偏微分方程数值解法

确定单元基函数-偏微分方程数值解法

时间:2023-10-26 理论教育 版权反馈
【摘要】:任取e单元作为讨论基函数选取的典型单元。全域的基函数可看做所有单元基函数的总和,每个结点相应有一个基函数。在本例题中,将式进行叠加求和,即可获得总体基函数。相应结点序号为n的总体基函数可表示为上述表达式适用于n=2,3,4的内部结点所对应的基函数。图3.8.2 一维线性基函数

确定单元基函数-偏微分方程数值解法

任取e单元作为讨论基函数选取的典型单元。记e单元中的基函数为

{Nie}

记单元中的近似函数为

ue=uieNieuie是待定系数)

构造单元基函数{Nie}应使uie正好是相应结点上的函数值,且每个结点分别对应一个基函数。若单元e中有I个结点,则应有I个单元基函数Niei=1,2,…,I)。本例题中,线段单元有两个结点,单元e有两个自由度,基函数有两个,即

Niei=1,2)

根据拉格朗日基函数需满足的条件

Niepjijij=1,2)

其中,pj是单元结点序号为j的结点。于是得到本例中,基函数的插值条件为

978-7-111-44528-9-Chapter03-221.jpg

按上述原则构造的基函数,保证了近似函数中的待定系数正好是函数在结点上的值。我们根据单元中的结点数目、结点坐标以及插值条件,并假定基函数是多项式函数,就可完全确定Nie。本例中,单元中有2个结点,将选取2个基函数N1eN2e,假定基函数是多项式函数,则根据插值的条件,它必定是线性多项式

Nie=aie+biexi=1,2)

应满足

Niexje=aie+biexjeijij=1,2)

其中,xje是单元e中第j号结点的坐标。不难解得

978-7-111-44528-9-Chapter03-222.jpg

于是得单元e中近似函数的表达式为

978-7-111-44528-9-Chapter03-223.jpg

其中,Δh=x2e-x1e。(www.daowen.com)

引进无量纲自然坐标

978-7-111-44528-9-Chapter03-224.jpg

则有

978-7-111-44528-9-Chapter03-225.jpg

在自然坐标中,单元中序号为1,2的结点的坐标值分别为ξ=0,ξ=1,将其代入式(3.8.6)中,可以看到,单元基函数的函数值在相应结点上的值为1,其余结点0,单元中近似函数的系数正是相应结点上的函数值。

全域的基函数可看做所有单元基函数的总和,每个结点相应有一个基函数。由于一个结点可能分属于若干单元,因此这个结点的基函数将是它所分属的单元中相应基函数的叠加。若全区域中的基函数记为

{Nn}(n=1,2,…,N

则根据unineuie的表示方法,有

978-7-111-44528-9-Chapter03-226.jpg

全域的近似函数可以表示为

978-7-111-44528-9-Chapter03-227.jpg

又根据uieineun的转换关系,上式可写成

978-7-111-44528-9-Chapter03-228.jpg

可见,全区域的近似函数是各个单元中的近似函数的叠加。在本例题中,将式(3.8.5)进行叠加求和,即可获得总体基函数。相应结点序号为n的总体基函数可表示为

978-7-111-44528-9-Chapter03-229.jpg

上述表达式适用于n=2,3,4的内部结点所对应的基函数。当n=1,5时,边界结点所对应的基函数只需将无定义的边界外区域去掉即可。内部结点所对应的基函数如图3.8.2所示。它是“尖顶形”的函数,在相应结点上的值为1,其他结点的值都为0,相邻结点间呈线性变化。

978-7-111-44528-9-Chapter03-230.jpg

图3.8.2 一维线性基函数

免责声明:以上内容源自网络,版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。

我要反馈