第6章形函数要点.pdf
《第6章形函数要点.pdf》由会员分享,可在线阅读,更多相关《第6章形函数要点.pdf(14页珍藏版)》请在三一文库上搜索。
1、第六章单元形函数的构造. 1 6.1 形函数构造的一般原理. 1 6.2 形函数的性质. 7 6.3 用面积坐标表达的形函数. 8 6.4 有限元的收敛准则. 10 6.5 等效结点载荷列阵 11 6.5.1 单元载荷的移置. 11 6.5.2 结构整体载荷列阵的形成 . 11 6.5.3 载荷移置与静力等效关系 12 习题 . 14 第六章单元形函数的构造 在有限单元法的基本理论中,形函数是一个十分重要的概念,它不仅可以用作单元的内插函数, 把单元内任一点的位移用结点位移表示,而且可作为加权余量法中的加权函数,可以处理外载荷, 将分布力等效为结点上的集中力和力矩,此外,它可用于后续的等参数单
2、元的坐标变换等。 根据形函数的思想,首先将单元的位移场函数表示为多项式的形式,然后利用结点条件将多项 式中的待定参数表示成场函数的结点值和单元几何参数的函数,从而将场函数表示成结点值插值形 式的表达式。在本节中,重点讨论几种典型单元的形函数插值函数的构造方式,它们具有一定的规 律。然后以平面三角形单元为例,讨论了形函数的性质,在此基础上分析了有限元的收敛准则。 6.1 形函数构造的一般原理 单元的类型和形状决定于结构总体求解域的几何特点、问题类型和求解精度。根据单元形状,可 分为一维、二维、三维单元。单元插值形函数主要取决于单元的形状、结点类型和单元的结点数目。 结点的类型可以是只包含场函数的
3、结点值,也可能还包含场函数导数的结点值。是否需要场函数导 数的结点值作为结点变量一般取决于单元边界上的连续性要求,如果边界上只要求函数值保持连续, 称为 C0 型单元,若要求函数值及其一阶导数值都保持连续,则是C1 型单元。 在有限元中, 单元插值形函数均采用不同阶次的幂函数多项式形式。对于 C0 型单元, 单元内的 未知场函数的线性变化仅用角(端)结点的参数来表示。结点参数只包含场函数的结点值。而对于 C1 型单元,结点参数中包含场函数及其一阶导数的结点值。与此相对应,形函数可分为Lagrange 型(不需要函数在结点上的斜率或曲率)和Hermite 型(需要形函数在结点上的斜率或曲率)两大
4、 类,而形函数的幂次则是指所采用的多项式的幂次,可能具有一次、二次、三次、或更高次等。 另外,有限元形函数 N是坐标 x、y、z的函数,而结点位移不是x、y、z的函数,因此静力学 中的位移对坐标微分时,只对形函数N作用,而在动力学中位移对时间t 微分时,只对结点位移向 量作用。 (1)一维一次两结点单元 i x j x ij AIEL, 图 3-8 一维一次两结点单元模型 设位移函数u(x)沿 x 轴呈线性变化 , 即xaaxu 21 )( (5.90) 写成向量形式为 2 1 1)( a a xxu(5.91) 设两个结点的坐标为 ji xx ,;两结点的位移分别为 ji uu ,,可以代入
5、上式并解出 21,a a,得 j i j i u u x x a a 1 2 1 1 1 (5.92) 位移函数u(x)记作形函数与结点参数乘积的形式 j i j i u u x x xxu 1 1 1 1)(5.93) 得到形函数为 ij i ij j jiij j ij i xx xx xx xx NNxxxx x xx x xN 1 1 1 1 1 1 1 (5.94) 在自然坐标系内进行定义,则可得到形函数的标准化形式 2 1 2 1 ji NNN(5.95) 其中,自然坐标的变换公式为1,1, 2 21 LLL。 图 3-9 一维一次两结点单元的局部坐标表达 (2)二维一次三结点单元
6、(平面三角形单元) 在总体坐标系统下,任一点的某一方向的位移是 123 ( , )u x yaa xa y(5.96) 设三个结点的坐标是 kkjjii yxyxyx,, kji uuu,为三个结点在某方向上的位移,具有如下关 系 k j i kk jj ii u u u yx yx yx a a a a a a yxu 1 3 2 1 3 2 1 1 1 1 1(5.97) 得到形函数矩阵如下式 1 1 1 1 1 kk jj ii yx yx yx yxN(5.98) 1 i x1 j x L 1 1 L1 2 L 上述推导可用如下MATLAB程序实现: clear v=sym(1, x,
7、y) m=sym(1,x1,y1;1,x2,y2;1,x3,y3) mm=inv(m) N=v*mm simplify(factor(N) (3)三维一次四结点单元(三维四面体单元) 在总体坐标系统下,任一点的某一方向的位移是 zayaxaaxu 4321 )(5.99) 按相似的方法可以得到 e k j i eee kkk jjj iii e k j i u u u u zyx zyx zyx zyx zyx u u u u N a a a a zyxu 1 4 3 2 1 1 1 1 1 11(5.100) 形函数矩阵如下式 1 1 1 1 1 kk jj ii yx yx yx yxN(
8、5.101) (4)一维二次三结点单元(高次单元) 图 3-10 一维二次三结点单元模型 设位移函数为 3 2 1 22 321 1 a a a xxxaxaau(5.102) 用结点位移 kji uuu,代入并求解 T aaa 321 , 3 2 1 2 2 2 1 1 1 a a a xx xx xx u u u k i k jj i k j i (5.103) 得到 k j i jkik ji kjij ki kiji kj k j i k jj i u u u xxxx xxxx xxxx xxxx xxxx xxxx u u u xx xx xx xxu k i 1 2 2 2 2
9、1 1 1 1 (5.104) 上式等号右端第一项矩阵即为形函数。 (5)一维三次四结点单元(Lagrange 型) ijkl i j L k 图 3-11 一维三次四结点单元模型 位移函数为三次方程 4 3 2 1 32 1 a a a a xxxu(5.105) 需要四个结点参数才能唯一地确定其中的常系数。这四个结点可以分别取两个端点和两个三分点。 类似地,可以得到如下形函数方程 , 1 1 1 1 1 1 32 32 32 32 32 lkji l k j i ll kk jjj i NNNNN u u u u xxx xxx xxx xxx xxxu l k ii (5.106) 其中
10、形函数中的各元素为 likiji lkj i xxxxxx xxxxxx N, ljkjij lki j xxxxxx xxxxxx N, lkjkik lji k xxxxxx xxxxxx N, kljlil kji l xxxxxx xxxxxx N. (5.107) (6)一维三次二结点单元(Hermite 型) (平面梁单元) ii u , jj u , x 图 3-12 一维三次二结点单元 这类单元的位移函数为 4 3 2 1 32 1 a a a a xxxu(5.108) 对应的转角方程为 4 3 2 1 2 3210 a a a a xx dx du (5.109) 用结点参
11、数 T jiji uu代入求解 4321 aaaa,即 j i j i jj ii jjj i jj ii jjj i j i j i u u xx xx xxx xxx a a a a a a a a xx xx xxx xxx u u iiii 1 2 2 32 32 4 3 2 1 4 3 2 1 2 2 32 32 3210 3210 1 1 3210 3210 1 1 (5.110) 得到 j i j i jj ii jjj i u u xx xx xxx xxx xxxu ii 1 2 2 32 32 32 3210 3210 1 1 1 jiujui NNNNN(5.111) 其
12、中形函数矩阵中各元素为 3 2 32 ji jij ui xx xxxxx N, 3 2 32 ji iji uj xx xxxxx N 2 2 ji ji i xx xxxx N , 2 2 ji ji j xx xxxx N(5.112) 上述结果可用MATLAB程序进行验证: clear x=sym(x); j=0:3; v=x.j % v=1 x x2 x3; m=sym(1,x1,x12,x13;1,x2,x22,x23;0,1,2*x1,3*x12;0,1,2*x2,3*x22) mm=inv(m) N=v*mm; simplify(factor(N) (7)二维一次四结点单元(平
13、面四边形单元或矩形单元) 用形函数表达的位移方程如下 lkji l k j i llll kkkk jjjj iiii NNNN u u u u yxyx yxyx yxyx yxyx xyyx a a a a xyyxu 1 4 3 2 1 1 1 1 1 11 (5.113) 其中形函数矩阵的元素为 )( )( 2121 22 yyxx yyxx Ni, i=1,2,3,4 (5.114) 对于平面四边形单元和矩形单元,可用局部坐标系统很好地加以解释。局部坐标的范围定义为 -1+1,四个结点的值固定。局部坐标系下的形函数为 4 )1)(1( i N(5.115) -1,-1 1,-1 -1
14、,1 1,1 图 3-13 二维一次四结点单元 (8)三维一次八结点单元 在三维一次单元形函数中,函数值沿三坐标轴(x、y、z 轴)呈线性变化。假设位移函数沿各坐 标轴的线性变化),(zyxuu可写成 xyzayzaxzaxyazayaxaau 87654321 (5.116) 假设在 i 结点的位移值为ui,并将数值代入上式,其他各结点(j,k,l,m,n,p,q) 亦类推,共有8 个式 子,其中第1 式如下 iiiiiiiiiiii zyxazyazxayxazayaxaau 87654321 (5.117) 可是以求得系数解 q p n m l k j i qqqqqqqqqqqq pp
15、pppppppppp nnnnnnnnnnnn mmmmmmmmmmmm llllllllllll kkkkkkkkkkkk jjjjjjjjjjjj iiiiiiiiiiii u u u u u u u u zyxzyzxyxzyx zyxzyzxyxzyx zyxzyzxyxzyx zyxzyzxyxzyx zyxzyzxyxzyx zyxzyzxyxzyx zyxzyzxyxzyx zyxzyzxyxzyx a a a a a a a a 1 8 7 6 5 4 3 2 1 1 1 1 1 1 1 1 1 (5.118) 则有 q p n m l k j i qqqqqqqqqqqq p
16、ppppppppppp nnnnnnnnnnnn mmmmmmmmmmmm llllllllllll kkkkkkkkkkkk jjjjjjjjjjjj iiiiiiiiiiii u u u u u u u u zyxzyzxyxzyx zyxzyzxyxzyx zyxzyzxyxzyx zyxzyzxyxzyx zyxzyzxyxzyx zyxzyzxyxzyx zyxzyzxyxzyx zyxzyzxyxzyx xyzyzxzxyzyxu 1 1 1 1 1 1 1 1 1 1(5.119) 最后得到形函数的表达式为 1 1 1 1 1 1 1 1 1 1 qqqqqqqqqqqq ppp
17、ppppppppp nnnnnnnnnnnn mmmmmmmmmmmm llllllllllll kkkkkkkkkkkk jjjjjjjjjjjj iiiiiiiiiiii zyxzyzxyxzyx zyxzyzxyxzyx zyxzyzxyxzyx zyxzyzxyxzyx zyxzyzxyxzyx zyxzyzxyxzyx zyxzyzxyxzyx zyxzyzxyxzyx xyzyzxzxyzyxN(5.120) (9)帕斯卡三角形 上述各种位移函数的构造有一定的规律,可以根据所谓的帕斯卡三角形加以确定,同时,这样 制定的位移模式,还能够满足有限元的收敛性要求。以下是几种典型情况。 一
18、维两结点单元的情况: 图 3-14 一维两结点单元的变量组成 1 x 2 x 3 x 4 x 5 x 一维三结点单元的情况: 1 x 2 x 3 x 4 x 5 x 图 3-15 一维三结点单元的变量组成 二维高阶单元的情况: 1 22 3223 432234 54322345 xy xxyy xx yxyy xx yx yxyy xx yx yx yxyy 图 3-16 二维高阶单元的变量组成 三维四结点单元的情况: 图 3-17 三维四结点单元的变量组成 6.2 形函数的性质 下面以平面三角形单元为例讨论形函数的一些性质。平面三角形单元的形函数为 Nab xc y iiii 1 2 , (
19、i =1, 2 , 3)(a) 其中, 33 22 11 1 1 1 2 yx yx yx ,为三角形单元的面积, iii cba,为与结点坐标有关的系数,它们分别等 于2公式中的行列式的有关代数余子式,即a1 、b1、 c1,a2、b2、 c2和 a3、b3、c3分别是 行列式2中的第一行、第二行和第三行各元素的代数余子式。 常数项 线性项 二次项 三次项 四次项 五次项 x 2 x 3 x 4 x 5 x z 2 z 3 z 4 z 5 y 5 z y 2 y 3 y 4 y 对于任意一个行列式, 其任一行(或列)的元素与其相应的代数余子式的乘积之和等于行列式 的值,而任一行(或列)的元素
20、与其他行(或列)对应元素的代数余子式乘积之和为零。因此有: 第一,形函数在各单元结点上的值,具有“本点是1、它点为零”的性质,即在单元结点1 上, 满足 1 2 1 , 11111111 ycxbayxN(b) 在结点 2、3 上,有 0 2 1 , 21211221 ycxbayxN(c) 0 2 1 , 31311331 ycxbayxN(d) 类似地有 1,0,0, 0,1,0, 333223113 332222112 yxNyxNyxN yxNyxNyxN (e) 第二,在单元的任一结点上,三个形函数之和等于1,即 1 2 1 2 1 , 33321321 333222111 321
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 函数 要点
链接地址:https://www.31doc.com/p-5211171.html