数值模拟与仿真基础——5 空间问题的有限元法.pptVIP

数值模拟与仿真技术基础 麻向军 华南理工大学 聚合物新型成型装备国家工程研究中心 聚合物成型加工工程教育部重点实验室 5、空间问题的有限元法 5.1 4节点四面体单元位移函数 5.2 单元刚度方程 5.3 载荷移置与等效节点载荷 5.1、4节点四面体单元的位移函数 1、位移函数形式 12个节点位移只能确定12个多项式的系数，所以空间问题的 4 结点四面体单元的位移函数如下： 2、待定系数 β1,β2, …β12 的求解 由克来姆法则可知：当 6V?0，上述方程有唯一解： 为了描述方便，引入系数 3、位移函数的插值函数形式 写成矩阵形式： 4、位移函数及其性质 当节点位移一定时，单元形态完全决定于Ｎi ,Ｎj ,Ｎm。因而形函数就具有如下的性质： 　　性质1、形函数 Ｎi 在节点 i 处的值为１，而在其它两个节点 (j , m) 处的值为零。 1、位移函数： 2、应变矩阵——用单元节点位移表示单元应变 　　Ｂ —— 单元的应变矩阵，其物理意义反映了单元任一点的应变与单元位移之间的关系。 　　对于一个给定的单元，节点坐标一定，系数 bi , ci , di也随之确定， V 也为常数，因此，应变矩阵为常量矩阵。 　　由弹性理论中关于平面问题的物理方程可知，当不考虑变温影响时，单元中任一点的应力 ? 为:　　 4、四面体常应变单元刚度矩阵的显式： 5.3 载荷移置与等效节点载荷 (1)、集中力 P 的移置公式： 设① 单元 (i, j, m) 中任意一点 (x, y) 作用集中载荷 推导： ① 根据单元内位移与节点位移关系 (2)、体力 g 的移置公式 设：单元 (i, j, m) 上作用有体力 (3)、表面力 q 的移置公式 设： 单元 (i, j, m) 的边 jm 上有表面力 3、载荷移置举例 　　以单元自重 (或作用在单元形心处的集中力) 为例。 　　设一个均质等厚的三角形单元 ijm，其厚度为 t，面积为A，材料密度为 ρ，则单元的自重为：W=t?A? ρ , 且其作用在单元形心 c 处。 (1) 按照集中力处理：重力 P=[0 0 －ρV]作用在形心处 (2)、按照体积力的移植公式 谢谢大家！ z y O x p(xp ,yp,zp) i(xi ,yi, zi) j(xj ,yj, zi) m(xm ,ym, zp) uj vj wj ui vi wi um vm wm up vp wp 将：节点位移 节点坐标 代入(2-2)式，得到： (i,j,m,p)表示下标轮换，即 假定 性质2、在单元任一节点处，三个形函数之和等于1。 5.2 4节点四面体单元刚度方程 Ｄ —— 弹性矩阵，反映了单元材料方面的特性。 S —— 称为单元的应力矩阵，其物理意义是反映了单元中任一点的应力与节点位移之间的关系。 　　对于4节点四面体单元， D、B 为常量矩阵， S 也为常量矩阵，这种常应变单元，也是一种常应力单元。 Krs——子刚度矩阵 单刚——决定于单元的形状，大小，方向和材料性质 　　——无关单元平移或坐标轴改变 ④ 各节点相应的虚位移为： ③ 发生微小位移时，集中力作用点相应的虚位移为： ② 各节点上的等效节点载荷向量为： 　　　 ② 根据能量等效原则： ② 积分在整个单元上有 推导： ① 在单元中取微元体dxdydz，将其体积力 gdxdydz 作为集中力，应用集中力的移置公式　 上述公式适用于任何单元及任意坐标方向。 推导：将单元面积 tds 上的表面力 qtds 当作集中力 j Rjx Rjy m Rmx Rmy i Rix Riy O x y * *