计力第2章-4.pptVIP

4．温度变化引起的等效结点载荷 由于温度变化，体内将产生应变，可被作为初应变来处理。 假如温度改变为T，产生的初应变是 上式中代入（2．6．7）式的Bi和表1．2中的弹性矩阵D以及（2．6．7）式的?0可得 这样就得到（2．6．34）的近似表达式 2．6．4 单元刚度矩阵的精确积分 由计算单元刚度矩阵的（2．6．12）式可以得到应变矩阵B不是常量阵时，与（2．6 15）式相应的分块矩阵的计算公式为 把应变矩阵分成两部分 将（2．6．37）式代入（2．6．36）式可以得到 因此（2．6．39）式中第二，三两项积分为零，即 在 2 ?2的分块矩阵中只需修正一个元素。 考虑到（2．6．40）式，修正项 经运算、整理并采用下列符号 （2．6．46）式可合并写作 因此如果单元某个结点的r 为零时，只要将其对应的对数项除去即可。 （2）当单元的两个结点的径向坐标相等时，例如rj= rm ，见图2．23， 然而，在此情况下，梯形jljmml的面积为零，因此对于j－m部分无需进行积分，只需令Ajm=Bjm=0即可。 以上两种情况在网格划分时通常都会出现。综合以上讨论，可以将（2．6．45）式作如下修改即可消除计算中的奇异性： （3）对于有两个结点在对称轴上的单元，例如结点i和j，则有rj=ri=0，此时除了分块矩阵外，其他分块矩阵的积分均发散，但是此时应有ui=uj=0，引入此位移条件修正刚度矩阵是将第2i－1行和2i－1列以及 2j-1行和 2j-1列划掉。其结果就是令在此情况下单元刚度矩阵的修正项 如果稍作分析以上结论是不难理解的。因为3结点三角形单元的径向位移u可表示成 既然 其中任一分块krs (3?3)由下式计算 2．7．2 六面体单元 如同在平面问题中采用矩形单元一样，在空间问题中也可以采用六面体单元。 图2．27（b）表示的是20结点六面体单元。 除了角点外，每边中点也是结点。 单元共有60个自由度。 位移函数可选取 2．8 小 结 本章以弹性力学静力分析问题为例讨论了有限单元法的基本原理和表达格式，其基本点可以归纳如下： 1．建立有限单元法的原理和选择单元的类型。 本章所采用的最小位能原理和基于它的以结点位移为基本未知量的位移元是固体力学有限单元法中应用最为广泛，也是最成熟的一种选择。 现行的有限单元法的通用程序几乎无例外地都以位移元作为它最主要的、甚至是唯一的单元形式。 2．建立有限元的计算模型。 这里包括单元形式的选择、有限元网格的划分和边界条件的设置。 本章着重讨论了平面问题3结点的三角形单元，它是在有限单元法中最早采用、而且至今仍有广泛应用的一种单元。通过它可以更好地了解有限元的概念以及表达格式的形成。 本章所引入的平面或轴对称的高阶三角形单元和矩形单元以及空间的四面体单元和长方体单元实际上都可以看成是平面3结点三角形单元的推广 至于平面或空间的任意形状单元的插值函数和特性矩阵的形成将在第3、4章中详细讨论。关于有限元网格划分和边界条件设置的一般原则通过本章的讨论可以有个基本的认识，在第5章中将给予进一步的讨论。 3．建立有限元的特性矩阵和求解方程。 单元特性矩阵包括插值函数矩阵、应变矩阵，因最终要建立的是单元刚度矩阵和载荷向量，并用以形成有限单元法的求解方程。 本章给出的广义坐标有限单元法的表达格式具有普遍意义。在以后各章的讨论中将可以看出，对于静力学问题如采用现今更为广泛应用的等参元，所不同的只是一些具体计算方法，而基本格式是一致的。 在2．6．2节中我们把Br、Bs中的r 和z用单元形心处的 和 来近似，得到了单元 刚度矩阵的近似积分。 这一节将讨论它的精确积分。 其中 是用单元形心的 和 代入Br中 得到的应变矩阵中的不变部分，就是计算近似单元刚度矩阵时用的应变矩阵； 是应变矩阵的变化部分 上式中 ，D都是常量阵，可以提在 积分号的外面； 或 的积分包含下面关系 于是（2．6．39）式可以写作 其中 就是（2．6．16）式中给出的 近似单元刚度矩阵， 是它的修正部分 为 式中I1、I 2、I 3代表积分 上述积分需在三角形ijm上进行，具体积分时，可以在三个梯形（见图2．22）i1imm1，m1mjj1以及i1ijj1上进行，由前面二个梯形面积积分和中减去第三个梯形面积上的积分。 上述积分需在三角形ijm上进行，具体积分时，可以在三个梯形（见图2．22）i1imm1，m1mjj1以及i1ijj1上进行，由前面二个梯形面积积分和中减去第三个梯形面积上的积分。 （2．6．44）式的积分是 （2．6．46）式中 表示对括号中的各项 i，j，m 指标轮换后再求和。 采用（2．6．43）及（2．6．47）式计算修正的刚度系数时，有两