有限元方法 第五章 平面三角形单元.pptVIP

节点的布置是与单元的划分互相联系的。通常，集中载荷的作用点、分布载荷强度的突变点，分布载荷与自由边界的分界点、支承点等都应该取为节点。并且，当物体是由不二.节点的选择及单元的划分图3-11第64页,共89页，2024年2月25日，星期天节点的多少及其分布的疏密程度（即单元的大小），一般要根据所要求的计算精度等方面来综合考虑。从计算结果的精度上讲，当然是单元越小越好，但计算所需要的时间也要大大增加。另外，在微机上进行有限元分析时，还要考虑计算机的容量。因此，在保证计算精度的前提下，应力求采用较少的单元。为了减少单元(a)(b)图3-12同的材料组成时，厚度不同或材料不同的部分，也应该划分为不同的单元。第65页,共89页，2024年2月25日，星期天在进行节点编号时，应该注意要尽量使同一单元的相邻节点的号码差尽可能地小，以便最大限度地缩小刚度矩阵的带宽，节省存储、提高计算效率。如前所述，平面问题的半带宽为B=2(d+1)，在划分单元时，对于应力变化梯度较大的部位单元可小一些，而在应力变化比较平缓的区域可以划分得粗一些。还有一点值得注意的是，单元各边的长度不要相差太大，以免出现过大的计算误差或出现病态矩阵。例如，图3-12所示的(a)、(b)两种单元划分，虽然都是同样的四个节点，但(a)的划分方式显然要比(b)的方式好。三.节点的编号第66页,共89页，2024年2月25日，星期天若采取带宽压缩存储，则整体刚度矩阵的存储量N最多为N=2nB=4n(d+1)其中d为相邻节点的最大差值，n为节点总数。例如在图3-13中，(a)与(b)的单元划分相同，且节点总数都等于14，但两者的节点编号方式却完全不同。(a)是按长边进行编号，d=7，N=488；而(b)是按短边进行编号，d=2，N=168。显然(b)的编号方式可比(a)的编号方式节省280个存储单元。(a)(b)图3-13第67页,共89页，2024年2月25日，星期天四.单元节点i、j、m的次序在前面章节中，我们曾指出，为了在计算中保证单元的面积?不会出现负值，节点i、j、m的编号次序必须是逆时针方向。事实上，节点i、j、m的编号次序是可以任意安排的，只要在计算刚度矩阵的各元素时，对?取绝对值，即可得到正确的计算结果。在实际计算时，应该注意所选有限元分析软件的使用要求。五.边界条件的处理及整体刚度矩阵的修正在前面讨论整体刚度矩阵时，已经提到，整体刚度矩阵的奇异性可以提高考虑边界约束条件来排除弹性体的刚体位移，以达到求解的目的。第68页,共89页，2024年2月25日，星期天一般情况下，求解的问题，其边界往往已有一点的位移约束条件，本身已排除了刚体运动的可能性。否则的话，就必须适当指定某些节点的位移值，以避免出现刚体位移。这里介绍两种比较简单的引入已知节点位移的方法，这两种方法都可保持原[K]矩阵的稀疏、带状和对称等特性。下面我们来实际考察一个只有四个方程的简单例子。⒈保持方程组为2n×2n系统，仅对[K]和{R}进行修正。例如，若指定节点i在方向y的位移为vi，则令[K]中的元素k2i,2i为1，而第2i行和第2i列的其余元素都为零。{R}中的第2i个元素则用位移vi的已知值代入，{R}中的其它各行元素均减去已知节点位移的指定值和原来[K]中该行的相应列元素的乘积。第69页,共89页，2024年2月25日，星期天假定该系统中节点位移u1和u2分别被指定为当引入这些节点的已知位移之后，方程(a)就变成然后，就用这组维数不变的方程来求解所有的节点位移。显然，其解答仍为原方程(a)的解答。u1=?1，u2=?2第70页,共89页，2024年2月25日，星期天⒉将[K]中与指定的节点位移有关的主对角元素乘上一个大数，如1015，同时将{R}中的对应元素换成指定的节点位移值与该大数的乘积。实际上，这种方法就是使[K]中相应行的修正项远大于非修正项。若把此方法用于上面的例子，则方程(a)就变成事实上，该方程组的第一个方程为第71页,共89页，2024年2月25日，星期天例图5-16所示为一平面应力问题离散化以后的结构图，其中图（a）为离散化后的总体结构，图（b）为单元1，2，3，4的结构，图（c）为单元3的结构。用有限元法计算节点位移、单元应变及单元应力（为简便起见，取泊松比