FE-Ch04.1-5等参元与数值积分汇.pptVIP

　　 基于最小位能原来基础上建立的位移有限元，其解答具有下限性质。即有限元的计算模型具有较实际结构偏大的整体刚度。选取缩减积分方案将使有限元计算模型的刚度有所降低，因此可能有助于提高计算精度。 另外，这种缩减积分方案对于泛函中包含罚函数的情况也常常是必须的，用以保证和罚函数相应的矩阵的奇异性（见相应教程），否则将可能导致完全歪曲了的结果。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 　　 保证结构总刚度矩阵[K]是非奇异的　　求解系统方程[K]{a}＝{P}，要求方程有解则必须系数矩阵的逆矩阵[K]-1是存在的，即在引入强迫边界条件后K必须是非奇异的。系数矩阵K非奇异的条件是|K|≠0，或称[K]是满秩的。所谓满秩就是系数矩阵中独立的行（列）数。如果K是N阶方阵，则要求它的秩为N，数值积分应该保证K的满秩，即K的所有行（列）的系数都是相互独立的。否则将使求解失败。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 二、关于矩阵的秩，有下面两个基本原则： 矩阵相乘的秩原则为 　　若 B=UAV 则 秩B≤min（秩U，秩A，秩V） （2.1） 　　秩B就是矩阵B的秩，它必然小于最多等于相乘矩阵中秩的最小者。 矩阵相加的秩原则 　若 C=A+B 则 秩C≤秩A+秩B (2.2) 即矩阵和的秩必然小于最多等于矩阵秩的和。 　　　现在再来考察单元刚度矩阵的计算公式: 　　 （2.3） Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 其中弹性矩阵[D]是d×d方阵，秩D=d。d是应变分量数(或独立关系数)。对于二维平面问题，d＝3，轴对称问题d＝4，三维问题d＝6。应变矩阵B是d×nf矩阵，nf是单元的结点自由度数。在一般情况下dnj, 所以秩B＝d。根据矩阵的秩的基本规则（2.1）式及（2.2），可以得到结论：秩Ke≤ng·d，ng是高斯积分点数。如果系统的单元数为M，再一次利用（2.2）式可得秩　　　　　　　　　　　K≤M·ng·d （2.4） 因此，刚度矩阵K非奇异的必要条件是　　　　　　　　　　　M·ng·d≥N 　（2.5）N是系统的独立自由度数，也就是刚度矩阵K的阶数。 （2.5）式表明：假如未知场变量a的元素数目超过全部积分点可能提供的独立关系数，则矩阵[K]必然是奇异的。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 实际计算中需要检查上述矩阵K奇异的必要条件是否满足，这在采用缩减积分方案计算K时是必须检查的，因为采用精确积分方案计算矩阵K时，矩阵K非奇异的必要条件是恒被满足的。这是由于任何有别于刚体运动的位移模式加于一个真实的结构系统，必将在系统内产生大于零的应变能。如果采用精确积分方案进行K的计算，实际上就是对于给定离散方案下的应变能aTKa进行精确的计算。因为aTKa0，所以K必然是正定的，也即K是非奇异的。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 而在采用减缩积分时，情况就不一样，系统的应变能并未被精确计算。对于某些除刚体运动而外的位移模式，在减缩积分点的应变可能正好等于零。如果结构不能约束这种位移模式，则该点的独立关系式实际上就不能对K的秩作出贡献，因而最后导致K非奇异的必要条件（2.5）式不能满足。图4.1所示就是一个例证。对于一个8结点单元的情况，在给予刚体位移约束以后，独立自由度数N＝2×8－3＝13。如采用2×2缩减积分方案计算刚度矩阵K，则其秩最大可能M·ng·d＝1×4×3＝1213。显然此时矩阵K非奇异的必要条件（2.5）式未被满足，因此K是奇异的。 Evaluation only. Created with Aspos