FE-Ch05.1-3有限元应用中及实际考虑.ppt

例如，对一个实际物体按照3维问题进行离散处理 当实际物体太大，可以考虑先离散成子结构(Substructures)，用超级单元处理 一.子结构法的理论根据： 1. 最小位能原理： 不失一般性，考虑仅有一级子结构， 上述位能泛函可以写成： 转化矩阵： 表示了单元内结点编号和子结构内结点编号的关系。 表示了子结构内结点编号和整体结构内结点编号的关系。 首先可以对子结构集成方程： 2.子结构集成方程的讨论 将子结构中自由度 分离为 和 相应的方程可写为： i) 与内部自由度有关的项 已集成完毕。 内部结点与子结构以外的任何单元不相联, 所以相应的项不受其它结构的影响. ii) 尚未集成完毕 iii) 由前面讨论过的Gauss消去法的特点，可以 边集成边消元。 a) 元素集成完毕的行可以作为消元行。 b) 被消元行元素可以没有集成完毕，集成与 消元可以同时进行，这里我们可以先用 的有关行进行消元。 我们叫做内部自由度的凝聚 二.内部自由度的凝聚 由第一组方程，可解出： 代入第二组方程得： 或： ( * ) 1. 对总体方程的贡献仅外部自由度有关项 因为内部自由度已经消元完毕，没有必要 再在总刚中出现，因此大大降低了总刚的 自由度数。 由总体方程可解出 得到每一个子结构的 2. 由(*)式可以计算 ，并计算各单元应力。 （注：需要对每个子结构保留一些信息 外存） 3. 程序实现不采用 求逆的方法，而是采用 消元 回代方法。 消元： 将内部自由度有关各行作为消元行。 b) 、 正是集成总体方程所需的项。 c) 为回代所需项，要保留 外存。 d) 总体方程求解完后， i) 对每个子结构取出 ii) 由外存读入 iii) 由下往上回代得 《注》：适用范围，整体分析、分叉结构… §5.４ 非协调元与分片试验 问题提出：经典的矩形单元，具有双线形位移模式。位移模式是具有完全一次式，非完全二次式。然而，有限元的计算精度取决于插值函数的完全项的阶数．因此，导致了经典的矩形单元产生其固有的缺陷，如：剪切锁定和奇异能量模式． 一. Wilson平面非协调单元 Wilson平面非协调单元的思路：引入内部位移自由度，使位移模式具有完全二次式的形式 。然而，内部位移自由度并不和实际的节点相联系，因此可以理解为增加了虚拟的内部节点．这种节点不能保证两个单元的界面位移连续。因此，称为“非协调单元” 在矩形双线性单元基础上，增加一位移修正项： 显然在四个结点处修正项等于零，因此它只影响单元内部位移，可见[α]e是单元内部位移参数。所求得的结点位移将仍直接是实际问题的近似值。 Wilson平面非协调单元的构造 根据这一位移模式，可得应变矩阵为 经虚位移原理做单元列式，可得 双线性 修正项 相关项 相关项 结点力 双线性等效载荷 修正等效载荷 各项的具体计算公式见下页。 双线性 修正项等效 相关项 修正项 双线性等效 由于[α]e是单元内部位移参数，与单元集合体（整体）的其他单元无关，因此在集成之前可做如下消去(称:静态凝聚) 单元刚度静态凝聚 单元等效荷载 由此可见做了上述消去后，单元刚度方程形式上和双线性单元完全一样。 协调元是肯定收敛的，非协调元收敛性不能保证。下一章将介绍如何检查非协调元的收敛性。这里只给出结论： 1) 当单元网格为矩形或平行四边形时，分析将试收敛的。 2) 任意四边形时，以?=? =0点的[J]作为单元的各点的[J]，这样计算的修正单元刚度也能使分析收敛。 二. 三维非协调元 1.位移模式 以三维8结点，6面体为例 其中： 其位移模式包含： 一次项已完全了，二次项缺： 选择什么样的附加插值函数为好？ 其好处是在各角结点， 均为0，这样保证了： 在结点处位移是协调的。 此处： 称为单元内附加自由度。 协调性：在边界面上，位移的分布是双向二次 分布，但仅有4个边界结点，而内部附 加自由度是与其它单元无关的。 在各单元内值不变。 所以，在边界上位移不协调，此类单元称为非协调元。 2.Wilson非协调元的有限元一般格式： 代入? 可得： 近似取为零 可以得到： 类似子结构法中凝聚内部自由度 可得： 内部自由度与其它单元无关，不出现在总体 方程，不增加总刚的自由度数。 * Chap 5 有限单元法应用中的若干实际考