1、o6.1单元节点编号与带宽存储o6.2形状函数矩阵与刚度矩阵的性质o6.3边界条件的处理与支反力计算o6.4单元刚度矩阵的缩聚o6.5为以函数构造与收敛性要求o6.6C0型单元与C1型单元o6.7单元的拼片试验o6.8有限元分析数值解的精度与性质o6.9单元应力的计算结果的误差与平均处理o6.10控制误差和提高精度的h方法和p方法6.1单元节点编号与带宽存储o计算机进行有限元分析时，需要存储所有单元和节点信息，随着所求解问题自由度的增大，计算规模的增大，整体刚度矩阵的规模非常巨大。o由于整体刚度矩阵中显现出相邻单元之间的关联性，因此矩阵中的大部分数据都为零，反映非零数据的一个指标就是带宽。o由
2、于刚度矩阵是对称的，可以看出，若节点的自由度数目为m，则每一个单元在整体刚度矩阵的半带宽为 di=（第i个单元中节点编号的最大差值+1）m d=max（di）（i=1，2n）其中n为整个结构系统的单元数。显然 对于二维问题，m=2 对于三维问题，m=36.2形状函数矩阵形状函数矩阵与刚度矩阵的性质的性质o以一维杆单元为例，杆单元的位移场为12()1ijeeiuijujxxu xaa xuullN uN uiujuNNN形函数矩阵形函数矩阵1、左端发生单位位移，右端固定2、右端发生单位位移，左端固定3、发生刚体位移6.2形状函数矩阵与刚度矩阵的性质刚度矩阵的性质o仍然以一维杆单元为例，它的刚度方
3、程为 1、考虑单元左端发生单位位移，右端固定情况 2、考虑单元右端发生单位位移，左端固定情况 3、考察刚体位移111211212222kkupkkupo性质1：单元刚度矩阵的对角元素kii表示要使单元的第i个节点产生单位位移，而其它的节点位移为0时，需要在i点施加的节点力。o性质2：单元刚度矩阵的对角元素kij（ij）表示要使单元的第j个节点产生单位位移，而其它的节点位移为0时，需要在i点施加的节点力。o性质3：单元刚度矩阵是对称的。这可以由功的互等定理得到。对于线弹性体，力所做的功跟加载次序无关，这可以利用上面的性质1和2得到。11121121222110kkkkkk 11121221222
4、201kkkkkk 第一种加载状态第二种加载状态第一种加载状态下的外力在第二种加载状态下移动相应位移做的功为1121210+1=kkk第二种加载状态下的外力在第一种加载状态下移动相应位移做的功为1222121+0=kkk根据功的互等定理，可以得到结论：刚度矩阵刚度矩阵是对称对称的。o性质4：单元刚度矩阵是半正定的。o性质5：单元刚度矩阵是奇异的。o性质6：单元刚度矩阵的任意行或列代表一个平衡力系，当节点位移全部为线位移时，任意行或列的代数和应该为0。同样，由单元刚度矩阵所组装的整体刚度矩阵也有以下性质：1）对称性 2）奇异性 3）半正定性 4）稀疏性 5）非零元素呈现带状分布112211 11
5、2121121 122122121 121211102211.221.2nntijijijnnnnnnnnnUq Kqk u uk uk u uk u uk uk u uk u uk uk u uk u u6.3边界条件的处理与支反力的计算o位移边界条件在大多数情况下有两种类型。1、零位移边界条件 2、给定具体数值的位移边界条件 根据上述两类边界条件，刚度方程的求解有以下几种方法：1、直接法 2、置“1”法 3、乘大数法 4、罚函数法o直接法直接法abbabbbbK qPK qPaaabaababbbbKKqPKKqP0aq aaabbababbbbK uK qPK uK qPaaabaaba
6、bbbbKKqPKKqPaqu-1bbbbbaqKPK u1、既可以处理零约束，又可以处理非零约束的情况。2、处理过程直观。3、待求矩阵的规模变小（维数变小），适合于手工处理。4、矩阵的节点编号及排序改变，不利于计算机的规范化处理。o置置“1”法法bbbbK qP1000abbbbqKqP0aq 1、只能处理零约束情况。2、待求矩阵的规模不变，不需重新排列，适合于计算机处理。3、保持整体刚度矩阵的对称性，利于计算机的规范化处理。abbabbbbK qPK qPaaabaababbbbKKqPKKqP0aq 直接法o乘大数法乘大数法aaabaaababbbbMKKqMKuKKqP1、既可以处理零
7、约束，又可以处理非零约束的情况。2、待求矩阵的规模不变，不需重新排列。3、保持整体刚度矩阵的对称性，利于计算机的规范化处理。abbabbbbK qPK qPaaabaababbbbKKqPKKqP0aq 直接法aaaabbaababbbbMKqK qMKuK uK qPo罚函数法罚函数法 罚函数法的最大好处是可以直接求出位移边界上的支反力。支反力的计算：支反力的计算：除了罚函数法能够求出支反力以外，其它的方法都需要求解一定的方程得到。6.4单元刚度矩阵的缩聚o采用高次位移函数的单元也常被称为高阶单元。对于高次单元来说，除了几何端点以外，其余的那些节点可能与其它的单元不发生关系，当中间的节点与其
8、它单元无关时，我们称作是内部节点。而其余的节点是外部节点。既然内部节点与其他单元无关，那么在组成整体刚度之前，就可以把他们消去，也就是把内部节点的位移用外部节点的位移来表示。o以一维三节点杆单元为例1 12233()u xN uN uN u111213112122232231323333kkkupkkkupkkkupo以一维三节点杆单元为例111213112122232231323333kkkupkkkupkkkupaaabaababbbbKKqPKKqPaaaabbabaabbbbK qK qPK qK qP11aaaabbbbbaaabbbbbaaK qK KPK qPqKPK q其中a代
9、表的是外部节点外部节点，b代表的是内部节点内部节点。6.5位移函数构造与收敛性要求o单元中的位移模式一般采用设有待定系数的有限多项式作为近似函数，优先多项式的选取原则应该考虑以下几个方面：1、待定系数是由节点位移条件确定的，因此它的个数应该与节点位移DOF个数相等。2、在选取多项式时，必须选择常数项和完备的一次项。单元位移模式中的常数项和一次项可以反映单元的刚体位移合唱应变的特性。这是因为当划分的单元数趋于无穷时，即单元缩小趋于一点，此时单元应变趋于常数。3、选择多项式应该由低到高，尽量选取完全多项式以提高单元的精度。o因此，在构造一个单元的位移函数时，应该参考由多项式函数构成的Pascal三
10、角形和上述原则进行函数项次的选取与构造。o收敛性问题收敛性问题 在有限元分析中，当节点数目或单元插值函数的项数趋于无穷大时，即单元尺寸趋于零时，最后的解答如果能够无线的逼近准确解，那么这样的位移函数或形函数是逼近于真实的，这就称为收敛。为使有限元分析的解答收敛，位移函数必须满足一些收敛准则，这些准则都经过过严密的理论验证。主要包括以下三个方面。o收敛性准则收敛性准则 定义：当单元尺寸趋于零时，有限元的解趋于真实解。准则准则1：完备性准则（针对单元内部）。如果在势能泛函中所出现的位移函数的最高阶导数是m阶，则有限元解答收敛性的条件之一是选取单元内的位移场函数至少是m阶完全多项式。准则准则2：协调
11、性准则（针对单元之间）。如果在势能泛函中位移函数所出现的最高阶导数是m阶，那么位移函数在单元交界面上必须具有直至m-1阶的连续导数，即Cm-1连续性。o以一般的梁问题为例 21dd2llUWEIz vxp x v xx 从上式可以看出，所出现的物理量是关于位移的最高阶导数为2，因此假定形状函数的时候，形函数至少应该包含完整的二次多项式。由准则2可知，位移函数为C1连续，即在单元之间的位移函数至少要求一阶导数连续。o以一般的平面问题为例1d2ddpxxyyzzxyxyyzyzzxzxxyzxyzsUWb ub vb wp up vp w A 从上式可以看出，所出现的物理量是关于位移的最高阶导数为
12、1，因此假定形状函数的时候，形函数至少应该包含完整的一次多项式。由准则2可知，位移函数为C0连续，即在单元之间的位移函数要求零阶导数连续。亦即函数的本身连续，而其一阶导数可以不连续。o协调元与非协调元协调元与非协调元 当单元的位移函数满足完备性要求时，称单元是完备的（一般都比较容易满足），当单元的位移函数满足协调性条件时，称单元是协调的（在单元与单元之间的公共边界上对于高阶连续性要求较难满足）。当单元的位移函数即完备又协调时，则有限元分析的解答是收敛的，即当单元尺寸趋于零时，有限元分析的解答趋于真实解。我们称这种单元为协调单元。o一般情况下，当泛函中的导数高于一阶时，则要求许可函数在单元交界面
13、上具有C1或更高的连续性，这时构造单元的插值函数往往比较困难。如果在单元之间的交界面上位移或导数不连续，将在交界面上引起无限大的变形，这时必须产生附加应变能，而我们建立泛函时，并没有考虑这种情况。因此，基于最小势能原理得到的有限元分析解答就不可能收敛于正确解。在某些情况下，可以放松对协调性的要求，只要这种单元能够通过拼片试验，有限元分析的解答仍然可以收敛于正确的解答。这样的单元称为非协调性单元。6.6C0和C1型单元oC0型单元 在泛函中位移函数的最高阶导数为1，在交界面上具有0阶的连续导数，即节点上仅仅要求位移连续。杆单元、平面问题单元、空间问题单元等6.6C0和C1型单元oC1型单元 在泛
14、函中位移函数的最高阶导数为2，在交界面上具有1阶的连续导数，即节点上除要求位移连续外，还要求1阶导数连续。梁单元、板单元、壳单元等6.7单元的拼片试验o由于非协调单元之间的位移不能保证位移协调，可以通过拼片试验来考证是否能描述常应变和刚体位移，若能通过拼片试验，则解得收敛性就能得到保证。o如图所示的单元状况，其中至少一个节点被单元所完全包围，若节点i完全被单元所包围，节点i的平衡方程为10meeijjieK uPo对于非协调单元，需要考察它的收敛性，即考察它是否具有常应变的能力，因此，我们设计这样一个试验（拼片试验拼片试验）：o当对单元片中的各个节点赋予对应于常应变状态的位移和载荷值时，核对对
15、i点平衡方程的正确性，如果能够满足，也就是单元满足常应变要求，因此当单元尺寸不断减小时，有限元解能够收敛于真正解。o以平面问题为例012012jjjjjjuaa xa yvbb xb y 由片面问题的平衡方程可知，当单元内的应变或应力都为常数时，则对应的体积力为零。对应于图中的i点，它的边界力也为零，因此 。所以此时，通过拼片试验的前提是，当赋予各节点以上位移模式的位移时，i点的平衡方程变为0eiP 10meijjeK u 即必须在节点i施加附加约束，该约束力所作的功等于单元交界面上位移不协调引起的附加应变能。o仍以平面问题为例012012jjjjjjuaa xa yvbb xb y 由片面问
16、题的平衡方程可知，当单元内的应变或应力都为常数时，则对应的体积力为零。对应于图中的i点，它的边界力也为零，因此 。所以i节点以外节点有以上位移模式的位移时，对于i点的平衡方程0eiP 如果求解上式得到的位移值和常应变状态下的位移相一致，则认为通过拼片试验。否则认为不能通过拼片试验。10meijjeK u6.8有限元数值解的精度与性质o求解精度估计 以平面问题为例，单元的位移场可以展开成以下形式.iiiuvuuxyxy 如果单元尺寸为h，则上式中的x和y都是h量级，若单元的位移函数采用p阶完全多项式，即它能逼近上述泰勒级数的前p阶多项式，那么位移解u的误差将是O（hp+1）量级。3节点3角形单元
17、(p次多项式)：量级位移应变应变能误差收敛速度误差误差h/1量级:O（h2）O（h2）O（h）O（h2）h/2量级:O（h2/4）。h/3量级:O（h2/9）。h/4量级:O（h2/16）。o这里讨论的都是仅仅局限于网格的离散误差，即当一个连续的求解域被离散成有限个子域，由单元的试函数来逼近整体的域的场函数所引起的误差。另外，实际误差还应该包括计算机的数值运算误差。21224/4iiiiO huuuuO h21143iiiuuu精确解与不同网格计算结果之间的关系o有限元分析的下限性质 有限元是把结构无限多的自由度简化为有限多的自由度，结构的刚度被夸大了，即使是用无限多个自由度来描述，也必然使得
18、原系统刚度增加，变得更加刚硬，即刚度矩阵的总体数值变大，由刚度方程知，计算出的位移结果偏小。由于位移函数的收敛性准则包含完备性和协调性两方面的要求，而完备性要求比较容易满足，而协调性则较难满足，因此这往往是研究的重点。o位移解的下限性质是基于协调单元单调收敛的前提得到的，在有些情况下，使用非协调单元也可以得到工程上的满意解答，有时甚至更好，这是由于位移不协调引所造成的误差与其它误差相抵消的缘故。6.9单元应力计算结果的误差和平均o应力结果的误差性质 对于弹性问题，其三大变量 对于一个具体问题，成了求2关于的极值问题。它是一个误差泛函。truetruetrueuuu1ddd2pTTTtrueij
19、itrueitruestrueUWb up uA trueuu o可见，对于求近似解极值近似解极值的问题 从力学上看，是求位移变分引起的总势能为极小值的问题。从数学上看，是求应变差和应力差在弹性矩阵加权意义下的最小二乘问题。因此，应变和应力的近似解的性质，是在加权残值最小二乘意义上对真实应变和真实应力的逼近。o高斯点上的应力性质高斯积分点上的应力和应变的近似解将具有比其它位置高得多的精度，这可以从图中看出。o公共节点上的应力平均绕节点直接平均法绕节点加权平均法，可以按体积或面积加权平均二单元平均法 11reeklkleiir 1eereeAA6.10控制误差和提高精度的h方法和p方法oh方法：
20、不改变各单元基函数，只通过逐步加密单元使计算结果向正确解逼近。它往往采用比较简单的单元。一般可以将误差控制在510%范围内。o其收敛性比p方法差，但是由于不用高阶多项式位移模式，数值稳定性和可靠性都好。op方法：保持网格固定剖分不变，增加单元上基底函数的阶次，从而改善计算精度。o实际证明：p方法收敛性好，由于使用高次多项式，会出现数值不稳定现象，另外，计算机容量和速度的限制，多项式的阶次不能太高，尤其在振动和稳定性问题求解高阶特征值时，这两个方法都不能得到令人满意的结果，这是多项式插值本身的局限性造成的。or方法：不改变单元类型和数目，只通过移动节点来减小误差，因而，单元的自由度数目不变。o自适应方法：它利用反馈原理，利用上一部的计算结果来修正有限元模型，其计算量小，计算精度却显著提高它是一种需要多次计算的方法，即多进程的方法，可以分别和h方法、p方法以及h-p方法相结合。自适应方法由误差指示算子来监控，而收敛程度则由误差估计算子来表征。它是一种能够自动调节其算法以改进求解过程的数值方法，它包括以下很多技术：误差估计 自适应网格改进 非线性问题中载荷增量的自适应选取 瞬态问题中时间步长的自适应调整等
