[2018年最新整理]FE-Ch-3单元与插值函数的构造.pptVIP

一.Lagrange插值多项式： 1.n个结点构造n-1次Lagrange插值多项式 注： 1)结点i的插值函数,2)?i为第个i结点坐标,3)?为自然坐标 即： 为结点当n-1次插值函数。 i=1,2…n 2. 的性质 i=1,2…n 1) n-1次插值函数,共有n个 2) 3) 3.构造一维单元插值函数： a. Lagrange线性插值 (n=2) 或记为： b. 二次Lagrange插值(n=3) 二.二维、三维Lagrange单元族： 1.二维情况 单元结点： ?方向n+1个点 n阶插值函数 ?方向m+1个点 m阶插值函数 1)插值函数：（一般情况） 很明显，有： 可证明： 常用当有：一次单元、二次单元 2).一次单元4结点单元（矩形） 双线性Lagrange 这里插值函数以结点编号，而不是两个方向。原因是为了在计算机上实现方便。 3).二次单元 9结点矩形单元 (记作： ) 插值函数： 角结点： 边中结点： 内部中结点： 2.Lagrange单元族的特点： 1)Lagrange单元族、插值函数构造方便 2)但存在一些问题：以二维为例， a)内部结点较多 (n-1)(m-1)个 对于二次以上的单元,就减少内部结点： 1.8结点四边形二维单元 右下图表示边中结点5插值函数的变化规律： ?方向二次Lagrange插值， ?方向一次Lagrange插值， i=5,6,7,8 显然满足， 2)角结点插值函数：双线性插值函数构成 以 为例，显然： a) b) 为了满足 ，需要修正上述函数 可取： 同理： 角结点插值函数也可以按照以下划线法构成:以(1±ξ)和(1±η)的乘积与通过其相邻的两边中点的直线方程相乘构成,以保证在两个邻边和对边满足插值函数为零. 3)插值函数性质： 4)位移模式 单元内： 较Lagrange单元族少 项，结点减少了（自由度） 但完全多项式次数不减。 边界上： 所以单元间界面上位移协调。 角结点插值函数可以表示成一双线性函数和用适当的分数乘以相邻的两个边界上的边内节点函数之和: 5） 边结点插值函数：如果只有一个边中点,可以直接表示成ξ(或η)方向二次和η(或ξ)方向一次式的乘积,以保证在两个邻边和对边满足插值函数为零. 对于p次单元,一般边内节点可以采用沿主方向(例如设ξ方向为p次,插值函数表示成ξ方向p次和η方向含待定常数的一次式的乘积.角结点插值函数可以表示成一双线性函数和用适当的分数乘以相邻的两个边界上的边内节点函数之和. 6 5 2 3 1 (-1,1) (1,1) (-1,-1) ξ η (1,-1) 4 (-1/3,-1) (1/3,-1) 例如:构造一个六节点的过渡单元, 以便进行3次—1次过渡变化. 要求在底边放置两个边内节点,边结点插值函数可以直接表示成: 代入条件: N5(-1/3,-1)=1, N5(1/3,-1)=0; N6(-1/3,-1)=0, N6(1/3,-1)=1; 可以求得: A=9/32,B=3,C=9/32,D=-3 得到: 6 5 2 3 1 (-1,1) (1,1) (-1,-1) ξ η (1,-1) 4 (-1/3,-1) (1/3,-1) 其中: 6 5 2 3 1 (-1,1) (1,1) (-1,-1) ξ η (1,-1) 4 (-1/3,-1) (1/3,-1) 代入条件: N1(-1/3,-1)=0, N1(1/3,-1)=0; N2(-1/3,-1)=0, N2(1/3,-1)=0; 可以求得: E=2/3, F=1/3, G=1/3, H=2/3 得到: 6 5 2 3 1 (-1,1) (1,1) (-1,-1) ξ η (1,-1) 4 (-1/3,-1) (1/3,-1) 右图表示沿边1562,插值函数的变化规律： N5 N6 例如:构造一个8节点的过渡单元, 边内节点均匀布置。求：边内节点5，7和角结点1的插值函数表达式: η ξ 1 2 3 4 5 6 7 8 ξ η 由教材知道，边内点5的插值函数可以直接表示成: η ξ 1 2 3 4 5