有限单元法讲座.pptVIP

B：按静力等效原理移置单元载荷。 一、体力的移置： y 0 x i j m 　　如果任意三角形单元ijk的重心c上受有自重 则按刚体静力等效原理可把W直接移置到i，j，m三个节点上而组成 ： 　　如果单元的重心c受有惯性力Pu作用，且 ，则Pu移置到i，j，m节点上的等效结点力为： 式中： ——旋转角速度 ——是单元重心处x坐标。 二、面力的移置 y 0 x i j m 已知在ij边受有面力q，则移置到i、j 结点上的等效节点力为： y 0 x i j m 当某一边上有三角形分布的面力时， 可由刚体静力等效直接写出 三、集中力的移置 如集中力G做用于其一边界上如图： 先将G分解为 ， 后，分别按线段的 比例把 和 分别移置到i，j两点上。 即： 0 x i j m y 下面我们通过一个只有四个方程的简单例子来说明。 总刚方程 　　在前面讨论整体刚度矩阵时，已经提到，整体刚度矩阵 的奇异性可以通过考虑边界约束条件来排除弹性体的刚体 位移，以达到求解的目的。 第六节 边界条件的处理 假定该系统中节点位移u1 和u3分别已知为 当引入这些节点的已知位移之后，总刚方程就变成为 然后，就用这组维数不变的方程来求解所有的节点位移。显然，其解答仍为原方程的解答。 u1 = ?1 ， u3 = ?3 方案(1)---置零法 事实上，该方程组的第一个方程为 方案(2)---乘大数法 假定该系统中节点位移u1 和u3分别已知为 当引入这些节点的已知位移之后，总刚方程就变成为 u1 = ?1 ， u3 = ?3 方案(3)---方程缩减法 假定该系统中节点位移u1 和u3分别已知为 当引入这些节点的已知位移之后，总刚方程就变成为 u1 = 0 ， u3 = 0 求位移 求反力 　　可以证明，对于一个给定的位移模式，其刚度系数的数值比精确值要大。所以，在给定的载荷之下，有限元计算模型的变形将比实际结构的变形小。因而，当单元网格分得越来越细时，位移的近似解将由下方收敛于精确解，即得到真实解的下界。 第七节 收敛准则 　　对于一个数值计算方法，一般总是希望随着网格的逐步细分所得到的解答能够收敛于问题的精确解。根据前面的分析，我们知道，在有限元分析中，一旦确定了单元的形状之后，位移模式的选择将是非常关键的。由于载荷的移置、应力矩阵和刚度矩阵的建立等等，都依赖于单元的位移模式，所以，如果所选择的位移模式与真实的位移分布有很大的差别，那么就很难获得良好的数值解。 　　⑵ 位移模式必须能包含单元的常应变。每个单元的应变一般都是包含着两个部分：一部分是与该单元中各点的坐标位置有关的应变（即所谓各点的变应变）；另一部分是与位置坐标无关的应变（即所谓的常应变）。从物理意义上看， 　　为了保证解答的收敛性，要求位移模式必须满足以下三个条件，即 　　⑴ 位移模式必须包含单元的刚体位移。也就是说，当节点位移是由某个刚体位移所引起时，弹性体内将不会产生应变。所以，位移模式不但要具有描述单元本身形变的能力，而且还要具有描述由于其它单元形变而通过节点位移引起单元刚体位移的能力。? 刚体位移: ⑶ 位移模式在单元内要连续、且在相邻单元之间的位 移必须协调。当选择多项式来构成位移模式时，单元内的 连续性要求总是得到满足的，单元间的位移协调性，就是 要求单元之间既不会出现开裂也不会出现重叠的现象。通 常，当单元交界面上的位移取决于该交界面上节点的位移 时，就可以保证位移的协调性。 当单元尺寸无限缩小时，每个单元中的应变应该趋于常量。 因此，在位移模式中必须包含有这些常应变，否则就不可 能使数值解收敛于正确解。很显然，在前面的位移模式 (b)中，与?2、?3、?5、?6 有关的线性项就是提供单元中的 常应变的。 在有限单元法中，把能够满足条件1和2的单元，称为完备单元；满足条件3的单元，叫做协调单元或保续单元。前面讨论过的三角形单元和矩形单元，均能同时满足上述三个条件，因此都属于完备的协调单元。在某些梁、板及壳体分析中，要使单元满足条件3比较困难，所以实践中有时也出现一些只满足条件1和2的单元，其收敛性往往也能够令人满意特别是放松条件3的单元，即完备而不协调的单元，已获得了很多成功的应用。不协调单元一般不象协调单元那样刚硬（即比较柔软），因此有可能会比协调单元收敛得快。 　　在选择多项式作为单元的位移模式时，其阶次的确定，要考虑解答的收敛性，即单元的完备性和协调性要求。实践证明，这两项确实是所要考虑的重要因素，但并不是唯一的 * 湖南