等参数单元课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《等参数单元课件.ppt》由用户(三亚风情)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 参数 单元 课件
- 资源描述:
-
1、 等参元的基本思想是:首先导出关于局部坐标系(Local coordinate, 或Natural coordinate, 自然坐标系)的规整形状的单元(母单元)的高阶位移模式,然后利用形函数多项式进行坐标变换,得到关于整体坐标系(Global coordinate)的复杂形状的单元(子单元),其中子单元的位移函数插值节点数与其位置坐标变换的节点数相等,位移函数插值公式与位置坐标变换式都采用相同的形函数与节点参数,这样的单元称为等参元。 等参数单元(Isoparametric elements)简称等参元,是根据特定方法设定的一大类单元,不一定具有相同的几何形状。因为等参元具有规范的定义原理和
2、较强的适应复杂几何形状的能力。在有限元理论中占有重要的地位。采用等参元,一方面能够很好地适应曲线边界和曲面边界,准确地模拟结构形状;另一方面,等参元一般具有高阶位移模式,能够较好地反映结构的复杂应力分布情况,即使单元网格划分比较稀疏,也可以得到比较好的计算精度。 (1) (1) 局部坐标系下的位移模式局部坐标系下的位移模式根据形函数的定义,在局部坐标系中,建立起几何形状简单且规整称之为母单元。下面以四边形单元为例说明等参元的基本概念。111111母单元是平面中的22正方形, 坐标形函数定义原点在单位形心上。单元边界是四条直线:,函数阶次相适应。对于具有线性形函数的四边形单元,共有四个节点,如图
3、6-1所示。三个,如图6-2。的单元,如图6-1所示,为保证用的未知量在相邻单元之间的连续性, 单元节点数目应与形如果是二次函数的四边形单元,单元每边的节点数为(a)母单元 (b)子单元图6-1 线性矩形单元及其平面坐标变换 对于如图6-1所示的线性4节点四边形等参元(Bilinear element),它在局部坐标系形函数如下:44443333222211111111, 1,N代入节点坐标值41100iN(i=1, 2, 3, 4) (6.1)isoparametricii00 其中 , . 该形函数是定义成自然坐标下的归一化变量 的函数。,对于8节点二次四边单元,角点上的形函数为11141
4、0000iN (i=1, 2, 3, 4) (6.2a) 边中点上的形函数为Ni121120(i = 5, 7) (6.2b) Ni121120(i = 6, 8) (6.2c) 可见,线性四边形单元和8节点二次四边单元用局部坐标形函数表达的位移模式如下11,4,8nniiiiiiuNu vNvn (6.3)(2)等参坐标变换等参元需要用坐标变换把形状规整的母单元转换成具有曲线(面)边界的、形状复杂的单元。转换后的单元称为子单元。子单元在几何上可以适应实际结构的各种复杂外形。即可以采用各种形状复杂的子单元在整体坐标系中对实际结构进行划分。,zyx,子单元通过坐标变换映射成一个局部坐标系下的规整
5、的母单元。坐标变换是指在局部坐标和整体坐标关系。在这里,坐标变换关系利用形函数建立起来。例如,对于上述4节点线性四边形单元,有4112244141122441,iiiiiixNxNxNxNxyNyNyNyNy (6.4a)之间建立一一对应 (6.4b)对于上述8节点二次四边形单元,有8112288181122881,iiiiiixNxNxNxNxyNyNyNyNy ,iNiiyx ,其中,是用局部坐标表示的形函数,是节点i的整体坐标上式即 为平面坐标变换公式。 如图6-1和图6-2所示的二维单元的平面坐标变换,其中母单元是正方形,子单元变换成曲边四边形,且相邻子单元在公共边上的整体坐标是连续的
6、,且在公共节点上具有相同的坐标,即相邻单元是连续的。 (3) 两种坐标系的关系雅可比矩阵(Jacobian matrix) 局部坐标系和整体坐标系之间具有如下偏导数的关系。根据复合函数的求导法则,有yzNyyNxxNNiiii (6.5)(6.6) 上式写成矩阵形式 yNxNJyNxNyxyxNNiiiiii 其中,J称为雅可比矩阵(Jacobian matrix) 将式(6.4)的表达式代入雅可比矩阵式(6.7),求得雅可比矩阵。对于式(7.4a)表示的双线性四边形单元,有 4411, iiiiiiNNxyxy(6.7)所以 1xyJ(6.9)例6-1:对一个四边形单元,在整体坐标系下的4个
7、节点(1, 2, 3, 4)的坐标分 别是(a, b), (c, b), (c, d), (a, d),求其雅可比矩阵及其逆矩阵。 解:根据式(6.1)、(6.4a) 和式(6.7)可以得到雅可比矩阵各元素)1 ()1 ()1 ()1 (41432111xxxxxJ11yyxxJJ(6.8)因为 的表达式为4 , 3 , 2 , 1,iyxii 为4个节点在整体坐标系下的坐标,将各坐标值代入得到也可比矩阵,具体为0202c adbJ进一步,由式(6.8)可求得其逆矩阵为12020c ad bJ)1 ()1 ()1 ()1 (41432112yyyyyJ)1 ()1 ()1 ()1 (41432
8、121xxxxxJ)1 ()1 ()1 ()1 (41432121yyyyyJ就可以 以二维二次8节点单元为例,说明等参元分析的一般原理。8节点四边形等参元的位移模式为8181, ,iiiiiivNvuNu(6.10)其中,ui和vi是节点i的位移。ii,x y, 采用坐标变换使母单元的8个节点与等参元的8个节点整体坐标和局部坐标的变换式为yi)一一对应。整体坐标8181, ,iiiiiiyNyxNx(6.11)将上述等参元的位移模式代入弹性力学平面问题的几何方程,将会得到如下形式的、用应变矩阵B表示的单元应变分量计算式值(xi,128xeeyxyuxvyuvyxBBBB (6.12)为了求得
9、应变矩阵B,进行如下推导。由于形函数,iN是局部坐标的函数, 需要进行偏导数的变换 1iiiiNNxNNyJ(6.13)其元素根据坐标变换式确定,即88118811;iiiiiiiiiiiiNNxyxyNNxyxy(6.14) 将单元任一点的应变列阵代入平面问题的物理方程,得到单元应力列阵。再利用虚功原理,进一步推导出单元刚度矩阵1111d dd deTTt x yt kB DBB DB J 式中,t为单元厚度。注意上式中的积分限,在整体坐标系中的积分相应地转化成了局部坐标系下的积分,是积分限为-1到+1的定积分。 (6.15) (6.16) 在等参元分析中,单元外载荷的计算分析如下。设外载荷
10、包括信集中载荷、体积力、表面力等,可以写成如下形式11NNeeeeee RRFQPF Q P 其中,对于集中载荷,设单元任意点c作用有集中载荷TxyGGG,移置到单元有关节点上的等效节点载荷为(6.17) 1,2,8eee TiixiyicFFNiFG (6.18) 对于体积力,设单元上作用的体力为Txyp pP,移置到单元各节点载荷为11111,2,8TxeeeiixiyiiypPPN tdxdyNt J d dpi Pp 对于表面力,设单元某边界上作用的表面力为Txyq qq,则这条边上的等效载荷为xeee TiixiyiyqQQNtdsq Q (6.19)式中, 是单元作用有面力的边界域
11、, ds是边界域内的微段弧长。在上述分析的基础上,利用结构中所有等参元的单元刚度矩阵集成结构整体刚度矩阵。列写结构有限元方程、引入约束条件,进而进行结构整体分析。上的等效三个节点将原来在整体坐标系统下的三角形平面单元,通过坐标转换,可以转换到局部坐标系中去,设在整体坐标系下的节点分别为( , )( , )( , )iijjkkx y x y x y、,对应到局部坐标系 中的坐标分别为(-1,1)、(-1,-1)、图6-3 三角形单元在整体坐标系和局部坐标系下的对应变换设坐标转换方程式为 (1,-1),如图6-3所示。(6.21)将各点坐标值带入上式得到:1231231 1111 11110 0
12、 1111xxxa b cyyyd ef 由上式可解出:11231231 1111 111111111xxxxyyyy (6.22)所以11231231 111 1 10 0 11 1 1 1 1 1a b cx x xd e fy y y J(6.23)10011xabcydef (6.20)(6.25)对于三角形单元,坐标转换方程式分别对与 求偏微分有如下xaxbydye代入上式可得:deta bdxdyd dd e (6.26)平面三角形单元的应变能公式可以写成:1111121det2TTTTUdAJ dd B DBB DB(6.27)三角形单元的面积的积分表达式可以写成如下形式 111
13、11111det|xxdAdxdyd dd dyy J(6.24)具体表达式dd11T11eJDBBK(6.29)VUW0V设位能等于应变能与外力功之和,即,由最小势能原理, 可知 0UW上式可写成:(6.30)其中,K即为单元刚度矩阵,如下(6.31)式中为节点位移。将应变能U对求微分得到:1111det TUJ d d B DB(6.28)0eeRK 如图6-4所示的矩形单元,其边长分别为2a和2b,两边分别平行于x、y轴。若取该矩形的四个角点为节点,每个节点位移有两个分量,所以矩形单元共有8个自由度。采用平面三角形单元的分析方法,同样可以完成对这种单元的等参数有限元分析。这里引入一个局部
14、坐标系、,可以推出比较简洁的结果。 u2 (U2)u3 (U3)v3 (V3)v2 (V2)v1 (V1)v4 (V4)u1 (U1)u4 (U4)xoyo12432a2b图6-4矩形单元(1) 坐标变换与形函数 在图中,取矩形单元的形心为局部坐标系的原点,和轴分别与整体坐标轴x和y平行,两坐标系存在有以下的坐标变换关系byyaxx00(6.32)式中2)(2)(2)(2)(2)(2)(2)(2)(/142334124132043210yyyybxxxxayyyyyxxxxx(6.33)在局部坐标系中,节点i的坐标是(i , i ),其值分别为1。取局部坐标系中的位移模式为uv 1234567
15、8(6.34)将4个节点的局部坐标值代入上式,可以列出4个节点处的位移分量,得到一组关于8个未知参数1,2,8的方程组,由此可求得这8个未知参数,最终可以得到用形函数矩阵和节点位移表示的位移模式方程4141 ,iiiiiivNvuNu4/)1)(1(00iN(6.35)(6.36)(2) 单元刚度矩阵由弹性力学几何方程求得单元内任一点的应变分量为11111xyxyuuubaxvvvaybabuvuvuvabbayx (6.37) (6.38)1234eBBBB 将式代入(6.35)得式中000001011001411iiiiiiiiiNbbNaaabababNNabB (i=1, 2, 3,
16、4) (6.39) 可以得出用节点位移表示的单元应力,即1234e= DSSSS (6.40)式中iiS=D B (i=1,2,3,4) (6.41) 对于平面应力问题0000200111141111122iiiiiiibaEbaabab S (6.42)将单元刚度矩阵写成分块形式11121314212223243132333441424344ekkkkkkkkkkkkkkkkk (6.43)其中的子矩阵按下式进行计算TijijtdxdykB DB (6.44)设单元厚度t是常量。上式可以积分出如下显式表达式1111211111132324 11111112323Tijijijijijijij
展开阅读全文