第七篇 薄板弯曲问题有限元法.ppt

式中前三项反映刚体位移，次三项对应于常应变。二次项只取了后三项是为了用结点位移表示参数 时考虑计算上的方便。同理，三次项不取前三项，剩下六个，只能挑选三个或进行某种线性组合。为了考虑每项面积坐标对称地出现，作出了如上的最简单可行的线性组合。组合未取“+”号，是由于 ，所以使用加号的最简单线性组合是不合宜的。 为了将位移模式写成标准形式，就需要求得形函数。为了方便起见，求形函数的工作可以分成两步进行。第一步是选取w、 、 作为结点自由度，求得对于它们的形函数，在这里把L3 = 1－L1－L2看作是L1和L2的函数。第二步是利用关系式 将第一步中所用的结点自由度变换成(7-38)式所指定的结点位移，再通过合并整理就很容易地得到形函数。式中b1 = y2 – y3，c1 = x3 – x1。对于b2、c2的值可以用下标轮换定出。 (7-26) 现在来决定参数 到 。将三角形单元的三个结点的面积坐标代入(7-25)式，立即得到 ， ， 。利用(7-25)式计算 和 ，得到 将结点的面积坐标代入上式，得六个方程如下 式中 表示对Li的偏导数在j点的值。从上式解得 将上式代入(7-25)式，并归并 和 前的各项，就可以得到对应于它们的形函数。 利用(7-26)式，将 和 变换为 和 ，于是得到相应于 和 的形函数 和 的计算公式。 最后得如下形式的形函数为 (d) 用下标轮换可得结点2和3的形函数，将上式写成矩阵形式 其中 (7-29) 而 是一个的系数矩阵，它是 式中 它们分别是结点1、2和3的面积坐标。 位移模式可写成如下的标准形式 可以验证在相邻单元间的挠度是连续的，但它的法向斜率仍不连续。事实上，在任何一条边上，挠度可表达成边线方向s的三次式，并且不包含与此边相对结点的结点位移（因为与它相对应的形函数在此边上等于零）。也就是说，一条边上的挠度可以由端部两个结点处的w和 所完全决定，而对于 则不然。因此这个三角形单元是一个完备的非协调单元。 在推导刚度矩阵和弯矩公式时，要计算 、 和 由于形函数是用面积坐标表示的，因此必须写出两个坐标系中的偏导数之间的关系。仍然取L1、L2作为独立坐标，而L3 = 1－L1－L2作为L1和L2的函数，根据复合函数的求导数规则，并且利用坐标变换公式，可以得到下列两个关系式 。 二．三角形单元的刚度矩阵 式中 是三角形的面积，而 将w的标准式代入几何方程，得单元应变列阵 式中 记号[ N ] i,11等表示[ N ] i对于L1的两次偏导数等。将(7-29)式代入上式，得到 由于[ L ]是面积坐标的三次函数，它的三个对L1和L2的二阶偏导数将是面积坐标的一次函数。把这三个偏导数算出后，容易把上式写成如下形式 其中 (7-38) 而 可以将[ C ] [ Ai ]乘出并记作矩阵[ Gi ]，于是 式中 (7-40) 若将单元刚度矩阵写成如下形式 其中子矩阵 若令 (7-41) 于是得 式中D是平板的弯曲刚度。 注意到(7-38)式，上式的积分是容易计算的。实际上 (7-43) (7-42) 式中 它的积分是 而 把(7-40)和(7-43)式代入(7-42)式并把它展开，可得单元刚度矩阵的子矩阵 式中[ H ]的所有元素，按(7-41)式展开 进行计算。 如果平板单元受有分布横向载荷q的作用，于是等效结点力是 这里以使用(d)式所示的形函数较为方便。 三．三角形单元的等效结点力和内力矩的计算 当q = q0是常量时，将(d)式代入上式，并积分得 对于内力矩阵，按下式计算 式中 (7-49) 把上式代入(7-49)式并进行一系列运算之后，可以得到内力矩列阵 式中 对于任意形状平板，它的边界条件可能是指定为沿着曲线边界切线方向的弯矩Ms或转角 和沿着法线方向的扭矩Mn或转角 ，这里所用的