有限单元法课件第五章 空间问题的有限元法.pptx

第一章 绪论1第二章 有限元法的基本原理2第三章 轴对称问题的有限元解法3第四章 杆件系统的有限元法4第五章 空间问题的有限元法5第五章 空间问题有限元法第一节 引言 当分析结构的形状,尺寸和边界条件不具备某种特殊性时,这种结构便属于空间问题。这时结构内的任一点i具有三个位移分量,即和六个应变分量以及六个应力分量式中,的表达式见式(2-6)以上15个分量都是空间坐标(x,y,z)的函数. 空间问题的规模要比平面问题和轴对称问题要大得多,且网格划分较困难,自动分网需要占用较长的时间,所以分析这类结构时要充分利用求解问题的特点(如对称性,相似性和重复性等)尽量减小有限元模型规模。有些结构可以先近似为轴对称问题或平面问题进行粗算,然后再按空间问题求解。 工程中的绝大部分结构属于空间问题,这类结构的几何模型为三维实体模型,网格划分时可采用四面体,五面体或六面体形式的空间实体单元。第二节 空间问题有限元法一、结构离散 在空间结构的离散中,常用的单元是四节点四面体单元,单元形状如图所示.本节将讨论这种单元的有限元法。常应变四节点四面体单元二、单元分析1.位移函数 图示的四面体单元,其节点为四面体的4个顶点,节点坐标分别为i(xi,yi,zi), j(xj,yj,zj), m(xm,ym,zm), p(xp,yp,zp)。每个节点有3个位移分量,即3个自由度,一个单元共有12个自由度。单元的位移列阵为假设单元的位移函数为线性函数 把4个节点的坐标值和位移分量分别代入式(5-5),联立求解可以求得待定系数,代回到式(5-5),经整理得(5-6)因此,单元位移可写成式中称为形函数矩阵。其中I是三阶单位矩阵;Ni,Nj,Nm,Np称为形函数,表达式分别为(5-8)式(5-8)中得各个系数分别为(5-9)V是四面体ijmp的体积,即 为了保证V的非负性,单元的4个节点必须按照图5-1所示的顺序进行编号,即在右手坐标系中,以右手螺旋按i-j-m的转向转动时是向p的方向前进为编号顺序。 由于四面体单元的位移函数是线性函数,可以证明,单元体任意三角形界面变形后仍保持为平面,且该平面由面上的3个节点位移决定。因此在变形过程中,两个相邻单元在交界面上的位移是一致的,即两个相邻单元的位移在交界面上连续,所以常应变四面体是协调元。2.单元应变矩阵 知道单元内各点的位移后,就可确定单元内任一点的应变,把式(5-6)代入几何方程式(2-2)得(5-10)式中(5-11) 式(5-11)中各矩阵元素的表达式见式(5-9)。显然应变矩阵 是常量矩阵,因此这类单元是常应变单元。3.单元应力矩阵由物理方程(2-5)得式中,的表达式见式(2-6);称为应力矩阵,它也是常数矩阵。由式(2-6)和式(5-11)得到(5-13)式中4.单元刚度矩阵由虚位移原理可以得到单元刚度矩阵的表达式为写成分块矩阵形式为其中任一分块矩阵 的计算公式为 式中式中, , ,的表达式见式(5-9),的表达式见式(5-13),可见单元刚阵是由节点坐标和材料弹性常数确定的,它是一个常数矩阵。三、总刚矩阵集成 求出每一个单元的刚度矩阵后,可按照第二章介绍的总刚矩阵的集成方法得到结构的总刚矩阵[K],即四、载荷移置 作用在结构的各类载荷同样需要移置为等效的节点载荷,移置方法与平面问题相同,下面给出移置结果。1.集中力的移置设作用在单元上的集中力为由 移置后产生的等效节点载荷根据第二章介绍的方法,上式各矩阵元素的表达式为式中,是形函数在载荷作用点处的值。2.表面力的移置设作用在单元上的表面力为由 移置后产生的等效节点载荷为则有式中,是形函数,其表达式见式(5-8)。3.体积力的移置设作用在单元上的体力为由 移置后产生的等效节点载荷为则有因此,由各类载荷移置产生的单元总的等效节点载荷列阵为由此可求出整个结构的节点载荷列阵,即五、约束处理和求解线性方程组 经约束处理消除[K]的奇异性后,即可求解结构刚度方程,求出每个节点的位移,进而求出单元内任一点的位移以及各单元的应变和应力。