平面问题有限元解法(公式推导讲解).ppt

* (f) 其中，Fe 是单元的结点力， k称为单元劲度列阵； 对三角形板单元，节点力为： (e) （5）应用虚功方程，由单元的结点应力 求出单元的结点力。假设把单元和节点切开，对右图中的i节点： 节点对单元的作用力为节点力， 作用于单元上。 三角形单元离散化结构分析步骤(续） （4）应用物理方程，由单元的结点位移 求出单元的应力： (d) 其中，S称为应力转换矩阵； Fe是作用于单元的外力，此外，单元内部还作用有应力。根据虚功方程，可以将单元的节点力Fe用应力来表示，从而得到节点力的公式： * （7）列出各结点的平衡方程，组成整个结构的平衡方程组。由于节点i受有环绕节点的单元移置而来的节点载荷 和节点力 因而i节点的平衡方程为： （i=1,2,…,n） (h) 三角形单元离散化结构分析步骤(续） （6）应用虚功方程，将单元中的外力载荷向结点移置，化为结点载荷（即求出单元的节点载荷）： (g) 将（f）代入（h），整理得： （j) 其中，K称为整体刚度矩阵，FL是整体结点载荷列阵，δ是整体结点位移列阵。 在上述求解步骤中，（2）至（6）是针对每个单元进行的，称为单元分析； （7）是针对整个结构进行的称为整体分析。 对三角形三个结点i,j,m结点，位移函数应当等于该节点的位移值，即： 三角形单元的位移模式 对每个单元，只要求得单元中的位移函数，就可以应用几何方程求得应变，再应用物理方程求得应力。有限单元法中常取结点位移为基本未知量，由单元的结点位移求出单元中的位移函数是首先必须解决的问题。 可以假定一个位移模式，来表示单元中的位移函数（即在单元中做出位移插值函数）。三角形单元中，可以假定位移分量只是坐标的线性函数，即假定： 6个方程解出α1-6，代入u，v式整理得： 其中： * 三角形单元的位移模式 Ni也可以该写成为： 其中系数ai,bi,ci是： 其中A就等于三角形ijm的面积： 按照解析几何学，在图示的坐标系中，为了得出的面积A不致成为负值，节点i，j，m的次序必须是逆时针转向的。 Ni，Nj，Nm这三个函数，表明了单元ijm的位移次形态（也就是位移在单元内的变化规律），因而称为形态函数，简称形函数。 * 三角形单元的位移模式 位移模式的表示式可用矩阵表示为： 简写为： 其中 是单元的节点位移列阵。 是形态函数矩阵或形函数矩阵。 有限单元法中，应力转换矩阵和劲度矩阵的建立以及载荷的移置等，都依赖于位移模式。 * 简写为： 其中矩阵B可写成分块形式： 其子矩阵为： 单元的应变列阵和应力列阵 利用几何方程和物理方程，求出单元中的应变和应力，用结点位移表示： 将位移函数(16)和(18)代入几何方程(6)，得出用结点位移表示单元应变。 * 将D表达式(9)和B表达式(27)代入上式，并写成分块形式，即得到平面应力问中的应力转换矩阵： 单元的应力列阵（续） 再将单元的应变式(26)代入物理方程(8)，得出用结点位移表示单元中应力的表达式。