有限元分析第七第一部分.doc

四阶问题（板的弯曲） 在常见的工程结构中，板或板梁结构较为普遍。在有限元分析中板可分为薄板和中厚板，在对薄板分析时采用了克希霍夫（Kirchhoff）假设。板中面上任一点（x, y）允许有三个位移分量，其中面内位移u、v构成一平面应力问题（二阶问题）。横向位移则构成一个四阶问题（弯曲问题）。对于线性问题（小挠度），这两个问题之间没有耦合。可以分别进行研究，再将结果迭加。平面应力问题在本章以前进行了研究。本章着重讨论板的弯曲问题。在对中厚板进行分析时，重点将介绍位移和转角各自独立插值的板单元，这种板单元考虑了板的剪切变形。而且，弯曲问题可以降阶为二阶问题来描述。（对插值函数要求CO连续） §7-1薄板小挠度弯曲的基本方程 　设板中面的横向位移（挠度）为w(x, y)　　　　 几何关系 中面法线绕x, y轴的转角 曲率 其中　　　　　又称为扭率 弹性关系 弯矩 其中Mxy又称为扭矩。t为板的厚度 平衡方程 设板在单位面积上受到的横向载荷为q(x,y)，则有 将（7-1-3）代入则得到以位移w为基本未知量的平衡方程 这是一个四阶椭圆型方程──双调和方程。 边界条件 设边界的切线方向为s。外法线方向为n则可将边界条件分为四类 挠度w的边界条件 要求在某一段边界（例如，Г1上）满足 转角边界条件 要求在某一段边界（例如Г2上）满足 弯矩边界条件 要求在某一段边界（例如Г3上）满足 剪力边界条件 要求在某一段边界（例如Г４上）满足 其中，关于w和　　 的边界条件为强制边界条件。 §7-2有限元解法 广义解 在所有满足强制性边界条件(即关于w,　 的边界条件)和协调条件（ 　　 在Ω上连续）的可能位移中使总势能泛函 取驻值的位移w称为板弯曲问题在Ritz意义下的广义解。 为了保证πP存在，必须要求 w∈H2(Ω)。故H2(Ω)为板弯曲问题（属于四阶问题）的容许空间。 收敛条件 如果试探函数取分片插值形式。我们将得到板弯曲问题的有限元形式的解，一般为近似解，为了保证有限元解能收敛到真实解。分片（单元）假定的位移场应满足以下条件 （１）　在单元内　　　　　　连续； （２）　包括足够的刚体位移模式。这种模式有三种：沿z方向的刚体平移和x、y轴的两个刚体型转动。为了保证这一点，假定的单元位移场中应包含x、y的完全一次多项式： （３）能够描述任何一种常曲率状态。为实现这一点。假定的单元位移场w中还应包含x、y的三个二次项 从（２）和（３）可知：w应至少包括x、y的完全二次多项式 (4) 协调条件 四阶问题要求穿过单元边界时　　　　　 连续。但如果沿边界的局部坐标系n－s考察。若穿过单元边界时w连续，则　　一定连续。故协调条件更恰当的提法应是：穿过单元边界时w（位移）和 （转角）连续。 上述四个条件为有限元解收敛到真实解的充分条件，其中条件（1）~（3）为必要条件。不满足条件（４）的单元，只有能够通过分片检验时才能保证收敛性。 为了同时保证位移和转角的协调性，一般采用Hermite型插值。这样至少可以保证节点处的协调性。既便如此。我们将会看到：实现协调性仍然是一件困难的事。下面介绍几种典型的单元。以说明构造平板单元的方法。 §7-3十二自由度矩形元（ＡＣＭ元） 单元和节点参数 单元：边与x、y轴平行的矩形。 取矩形的四个角点为节点。 取 为节点参数。 单元位移场 取 其中　　　　　（广义坐标）可由１２个节点参数唯一定出。 收敛性分析 上述位移场为x、y的四次多项式。完全到x、y的三次项。故收敛条件（1）~（3）可以满足。下面分析协调性。以２-３边为例。 沿２-３边： x = 常数。位移w是y的三次多项式。可以完全被节点２、３处的四个节点参数 所决定。故沿２－３边位移w是协调的。沿２、３边转角 是y的三次函数，不能仅由节点２、３处剩下的两个节点参数 所决定。故沿２－３　　不协调。对其它各边可得到类似的结论。 如果认为　　　　　　　　可以决定出沿边２－3线性变化的、协调的　，那么 在节点处为零，且仍然包括y的二次项和三次项，而其中的偶次项对通过分片检验是不利的，但是由于单元是形状十分规则的矩形，仍然可以通过分片检验；在广义解w∈H3(Ω)的情况下可以在理论上证明有限元解的收敛性。证明过程较为繁琐。这里从略。 ＡＣＭ元是出现较早的一种单元，虽然它不满足转角协调条件，但可以保证收敛性，加上它的单元刚度矩阵的元素可以用精确的解析表达式求得。这种单元在自编程序中得到广泛应用。