《机械结构有限元分析及强度设计》_第7章.pptVIP

图7.5SHELL181单元的几何特性、节点的位置及坐标系返回7.1薄板的弯曲变形7.1.1克希霍夫假设在弹性力学里，把两个平行面和垂直于这两个平行面的柱面或棱柱面所围成的物体称为平板，简称为板。两个板面之间的距离t称为板的厚度，而平分厚度t的平面称为板的中间平面，简称中面。如果板的厚度t远小于中面的最小尺寸b，且满足，则该板就称为薄板，否则就称为厚板。下一页返回7.1薄板的弯曲变形

薄板承受的载荷一般可分解为两个分量，一个是作用在薄板的中面之内的纵向载荷，另一个是垂直于中面的所谓的横向载荷。对于纵向载荷，可以认为它们是沿厚度方向均匀分布的，所以，由它们引起的应力、应变和位移都可以按平面应力问题进行计算。而横向载荷将会使薄板产生弯曲，所以，由它引起的应力、应变和位移可以按薄板的弯曲问题进行计算。本章只讨论由于横向载荷使薄板发生小挠度弯曲所引起的应力、应变和位移问题。当薄板弯曲时，其中面所弯成的曲面称为薄板的弹性曲面，而中面内各点在垂直于中面方向的位移称为挠度。对于薄板小挠度弯曲问题，通常采用克希霍夫假设，其内容包括以下4个方面。上一页下一页返回7.1薄板的弯曲变形

1.厚度不变假设由于垂直中面方向的正应变（即应变分量εz）极其微小，所以可以忽略不计。取εz=0，因而有：因此，板内各点的挠度w与z坐标无关，只是x、y的函数。2.直线假设在薄板弯曲变形前垂直于板中面的直线，在簿板弯曲变形后仍为直线，且垂直于弯曲后的中面。这说明在平行于中面的面上没有剪应变，即：上一页下一页返回7.1薄板的弯曲变形

3.正应力假设中面上的正应力远小于其他应力分量的假设：平行于中面的各层相互不挤压、不拉伸，则沿z向的正应力可忽略，即：4.小挠度假设薄板中面内的各点都没有平行于中面的位移，即：上一页下一页返回7.1薄板的弯曲变形

因故有因此，中面的任意一部分虽然弯曲成弹性曲面的一部分，但其在xy面上的投影形状却保持不变。上一页返回7.2矩形薄板单元分析7.2.1矩形单元的位移模式将平板中面用一系列矩形单元划分，则得到一个离散的系统以代替原来的平板，要使各单元至少在节点上存在挠度及其斜率的连续性，必须把挠度及其在x和y方向的一阶偏导数指定为节点位移（或称广义位移）。通常将节点i的位移列阵写成以下的形式：下一页返回7.2矩形薄板单元分析

则与之对应的节点力列阵可以表示为：它们的符号规定：对于挠度w和与之对应的节点力W，以沿z轴的正方向为正；对于转角θx、θy和与之对应的节点力矩MθxMθy，则以右手定则标出的矢量沿坐标轴的正方向为正。上一页下一页返回7.2矩形薄板单元分析

对于矩形单元，可在平面问题中引入一个自然坐标系oξη来研究其单元特性，如图7.3所示。由于矩形单元的每个节点有三个位移分量，一个单元有4个节点，即共有12个节点位移分量，因此，我们要选取含有12个参数的多项式作为位移模式，即：上一页下一页返回7.2矩形薄板单元分析

最后两项的选取是使单元在边界上有三次式的形式。按照式（7.20）可以算出转角，即：上一页下一页返回7.2矩形薄板单元分析

将矩形单元的4个节点坐标(ξi,ηi)分别代入式（7.20），就可以得到用12个参数来表示的节点位移分量的联立方程组，求解这12个方程，从中解出a1～a12，再代入式（7.21），经归纳并整理后就可以改写成如下的形式：或者写成标准形式，即：上一页下一页返回7.2矩形薄板单元分析

其中如果把形函数写成通式，即：上一页下一页返回7.2矩形薄板单元分析

于是有：其中，上一页下一页返回7.2矩形薄板单元分析

由式（7.9）可以看出，整个薄板的位移完全由平面在z方向上的挠度w决定，而在中面各点不产生x和y方向上的位移。因此，薄板可能产生的刚性位移就只有沿z方向的平动以及绕x轴和y轴的转动，而z轴方向的旋转是没有的。在位移模式式（7.20）中，其前三项反映了薄板单元的这三个刚体位移。再由式（7.10）可以看出，板内各点的应变完全由挠度w的三个二阶导数决定。如果这三个二阶导数不随坐标变化而变化，则可将其描述为平板单元的一个常应变状态，式（7.20）中的第四、五、六三个二次项反映了这个常应变状态（或称常曲率状态）。因此，总是能够保证存在一组节点位移，其可