05板壳弯曲问题的有限单元法.docVIP

第五章 板壳弯曲问题的有限单元法 本章将讨论弹性板壳弯曲问题的有限单元法。为便于理解，首先介绍薄板弯曲问题理论。 弹性理论将板厚和平面尺寸相比很小的板称为薄板(工程上则将的板就认为是薄板)。当薄板受有垂直于板面(亦即垂直于中面)的外力时，横截面上将发生弯曲和扭转，引起弯曲应力和扭转应力。薄板的中面成为一个弹性曲面。这种问题称为薄板的弯扭问题，但通常简称为薄板的弯曲问题。 对于薄板小挠度(挠度)问题，经一定简化，可认为它的变形完全由横挠度所确定。因此，可以取和它的若干阶导数作为节点参数建立平板单元。目前已提出了非常多的板单元。下面将着重介绍比较常用的部分协调的矩形单元和一种三角形单元。虽然它们都不是完全协调的，但所得计算结果表明，它们都具有良好的收敛性和精确度。在这之后将介绍一种考虑了横向剪切变形影响的平板单元，它适用于薄板和厚板，由于是完全协调单元，它具有较好的性能。 对于薄板大挠度问题，属于几何非线性范畴，将在以后章节专门讨论，本章不予介绍。 在平板弯曲问题的基础上，本章在最后还将介绍薄壳弯曲问题的有限单元法。 §5-1 薄板弯曲问题的基本原理 一、基本假设 薄板弯曲问题实际上是三维问题。为使问题简化，同时在一定程度上又能反映实际情况，并满足工程精度要求，分析时先作如下假设： (1)平行于中面的各层材料互不挤压；板的厚度没有变化。 (2)变形前板中面的法线，变形后仍保持为弹性曲面的法线。 (3)板中面内的各点，没有平行于中面的位移。 上述假设统称为克希霍夫(Kirchhoff)假设。在满足上述假设条件下，板中面挠度只是板平面内二个坐标的函数，而板内应力又都可以表示为的形式。 图5-1 薄板弯曲 取薄板的中面为xoy面，z轴垂直于中面，如图5-1所示，则第(1)条假设可以表示为 (5-1) 对z求导数为零的函数必与z无关，即有，也就是说，薄板中每一法线上所有各点都有相同的位移，也就是等于弹性曲面的挠度。 第(2)个假定表示：薄板的法线(z方向的线段)与x方向或者y方向的线段都保证互相垂直，没有剪应变，也就是,。于是由几何方程可见有 , 从而有 , 对z积分(注意和都不随z而变)，得 , (a) 其中和是任意函数。 第(3)个假定表示,。由式(a)可见，必须，。于是式(a)简化为 , (5-2) 公式(5-2)用挠度表示了位移分量，它具有明显的物理意义。 如图5-2所示，在薄板中面上有一点，变形后仅发生z方向位移(挠度，假设(3))。原来的中性平面现变为弹性曲面，经过A点的法线虽然位置发生了变化，但其仍和该曲面垂直(假定(2))，且长度保持不变(假定(1))。所以由图5-2也可直接得点的位移分量 (b) 式中，为中面法线绕y轴倾斜角度，在小挠度时有 (c) 从而，得到(5-2)式的第一式。同样，分析yoz面内的变形，可得(5-2)式的第二式。 二、应力应变表达式 1. 应变表达式 由上面分析知，弯曲薄板内仅剩下三个不为零的应变分量，根据式(5-2)，利用几何方程，可以将它们用挠度表示 (5-3) 在微小位移形变的情况下，和分别代表薄板弹性曲面在x方向和y方向曲率，而代表它在x和y方向的扭率。这三者完全确定了薄板内所有各点的应变分量，因此可以称为薄板的形变分量，而三者的总体可以称为薄板的形变。用矩阵表示，薄板的形变是 (5-4) 通过式(5-3)，薄板内各点的应变与薄板的形变之间的关系可以简单地表示为 (5-5) 2. 应力表达式 根据假定(1)，薄板内各层不相互挤压，即是忽略了的影响。板内各点的非零应力分量也仅剩下、、三个分量。由物理方程，有 (5-6) 式中 (5-7) 是薄板弯曲问题的弹性矩阵，它和平面应力问题的弹性矩阵完全相同。 由于薄板的形变和z无关，所以由式(5-6)可以看出，薄板内的应力沿板厚方向成线性分布。图5-3为板上微分六面体(dxdyt)上的应力分布情况。 图5-3 薄板上的应力分布 三、薄板上的内力 现在来分析式(5-6)及图5-3所示的应力分量将合成什么样的内力。从薄板上割取一个微分小六面体，它在x和y方向都是一个单位的宽度，如图5-4所示。在垂直于x轴的截面上，正应力由于和z成正比，它只可能合成一个力偶，而这个力偶的矩