FE-Ch04.1-5等参元及数值积分.ppt

　　 基于最小位能原来基础上建立的位移有限元，其解答具有下限性质。即有限元的计算模型具有较实际结构偏大的整体刚度。选取缩减积分方案将使有限元计算模型的刚度有所降低，因此可能有助于提高计算精度。 另外，这种缩减积分方案对于泛函中包含罚函数的情况也常常是必须的，用以保证和罚函数相应的矩阵的奇异性（见相应教程），否则将可能导致完全歪曲了的结果。 　　 保证结构总刚度矩阵[K]是非奇异的　　求解系统方程[K]{a}＝{P}，要求方程有解则必须系数矩阵的逆矩阵[K]-1是存在的，即在引入强迫边界条件后K必须是非奇异的。系数矩阵K非奇异的条件是|K|≠0，或称[K]是满秩的。所谓满秩就是系数矩阵中独立的行（列）数。如果K是N阶方阵，则要求它的秩为N，数值积分应该保证K的满秩，即K的所有行（列）的系数都是相互独立的。否则将使求解失败。 二、关于矩阵的秩，有下面两个基本原则： 矩阵相乘的秩原则为 　　若 B=UAV 则 秩B≤min（秩U，秩A，秩V） （2.1） 　　秩B就是矩阵B的秩，它必然小于最多等于相乘矩阵中秩的最小者。 矩阵相加的秩原则 　若 C=A+B 则 秩C≤秩A+秩B (2.2) 即矩阵和的秩必然小于最多等于矩阵秩的和。 　　　现在再来考察单元刚度矩阵的计算公式: 　　 （2.3） 其中弹性矩阵[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]必然是奇异的。 实际计算中需要检查上述矩阵K奇异的必要条件是否满足，这在采用缩减积分方案计算K时是必须检查的，因为采用精确积分方案计算矩阵K时，矩阵K非奇异的必要条件是恒被满足的。这是由于任何有别于刚体运动的位移模式加于一个真实的结构系统，必将在系统内产生大于零的应变能。如果采用精确积分方案进行K的计算，实际上就是对于给定离散方案下的应变能aTKa进行精确的计算。因为aTKa0，所以K必然是正定的，也即K是非奇异的。 而在采用减缩积分时，情况就不一样，系统的应变能并未被精确计算。对于某些除刚体运动而外的位移模式，在减缩积分点的应变可能正好等于零。如果结构不能约束这种位移模式，则该点的独立关系式实际上就不能对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是奇异的。 图4.1 给定位移模式下的8结点单元 三、完全积分 * * * * * * * * * * * * * 当单元具有规则形状时,所用的高斯积分点可以对单元刚度矩阵中的多项式进行精确地积分.对六面体和四面体单元而言,所谓“规则积分”是指单元的边相交成直角并且任何边上的节点是边中点.线性单元要完全积分,则在每一方向要有两个积分点.相应的,线性平面单元要有2×2个积分点,线性空间单元要有2×2×2个积分点.同理,二次单元要完全积分,则在每一方向要有三个积分点.相应的,二次平面单元要有3×3=9个积分点,二次空间单元要有3×3×3=27个积分点. 二次平面单元有3×3个积分点 线性平面单元要有2×2个积分点 已经证明对于梁的弯曲问题,完全积分容易导致所谓”锁定(locking)”现象.以下为典型的悬臂梁算例,在自由端受集中力作用.分别采用1×6, 2×12, 4×12, 8×24的网格,用四边形4节点和六面体8节点等参元求解. 在下表4-1列出计算结果与自由端的理论位移为3.09mm的比值.可看出,粗糙的网格的计算结果最差,但是8×24的网格所得结果也仅是精确值的56%.特别注意,当采用完全积分时在厚度方向布置的单元数对计算结果影响不大,这是由于剪力自锁(shear locking)引起的，对所有完全积分一阶实体单元均存在此问题． 为清楚起见,画过积分点的虚线. 很明显,上面直线的长度增加,说明1