第五章薄板弯曲.pptVIP

单元的内力 如已解出板结构的全部节点位移{δ}，则对任意的e单元都可以找出相应的单元节点位移{δ} e ，再应用应变矩阵[B]和薄板弯曲的弹性矩阵[D]，即可得到单元的内力 上式中[D] [B]= [S]为薄板弯曲应力矩阵，为3X12的长方矩阵。 第二十九页，共四十一页，2022年，8月28日 第一页，共四十一页，2022年，8月28日 5.1 薄板的弯曲变形 如h以表示板厚，以l表示其他方向的尺寸，当h/l15时，可认为是薄板。 板内厚度中点构成的平面称中面。 板件一般常驻有垂直于中面的载荷（横向载荷），在载荷作用下，板面发生弯曲，中面由平面变为曲面，称为挠曲面。 第二页，共四十一页，2022年，8月28日 以未变形的中面为xy坐标面，中面各点沿z轴的横向位移以w表示，称为挠度，如图5-1所示。 一般挠度为中面各点坐标的函数，即 w=w(x,y) 称为挠曲面方程。 第三页，共四十一页，2022年，8月28日 薄板弯曲时，板内各点的应变为 其中z为点到中面的距离 为挠曲面沿方向的曲率 为扭曲率 第四页，共四十一页，2022年，8月28日 当板弯曲挠度很小时，曲率、扭曲率与挠曲面的关系为 板挠曲面的曲率、扭曲率表示出板弯曲变形的程度，这3个分量也可合称为曲率。可用列阵表示为 第五页，共四十一页，2022年，8月28日 (5.1) 第六页，共四十一页，2022年，8月28日 因此，板内的应变可用列阵表示为 (5.2) 第七页，共四十一页，2022年，8月28日 应力与应变的关系为 其中[Dp]即平面应力问题的弹性系数矩阵 板的中面处z=0，有 即中曲面内没有面内应变，也没有面内应力。 (5.3) 第八页，共四十一页，2022年，8月28日 板内各点应变与其z坐标呈正比关系。 应力与z坐标也成正比，沿板厚度方向线性变化。 正应力σx,σy在板的横截面上将合成为弯矩，剪应力将合成为扭矩。分别表示如下： 第九页，共四十一页，2022年，8月28日 第十页，共四十一页，2022年，8月28日 弯矩Mx、 My与扭矩Mxy 3项为薄板弯曲的内力，合在一起用列阵表示为 第十一页，共四十一页，2022年，8月28日 将式(5.3)代入上式，并完成积分有 其中 [D]为薄板弯曲的弹性系数矩阵，它与平面应力问题弹性系数矩阵相似，只是多一个系数h3/12。 第十二页，共四十一页，2022年，8月28日 薄板弯曲的弹性应变能为 其中V为板的体积域。 将式(5.2)及(5.3)代入上式，并沿厚度方向积分，可得 第十三页，共四十一页，2022年，8月28日 其中S为板中面的面积域，[D]为薄板弯曲的弹性系数矩阵。 由上式可见，薄板弯曲变形时，单位面积中面的弹性应变能为其曲率的二次型。 板弯曲的曲率是其挠度w的二阶导数，因而薄板弯曲的弹性应变能为包括w二阶导数的二次泛函数。 (5.6) 第十四页，共四十一页，2022年，8月28日 5.2 四节点的矩形薄板单元 对于薄板弯曲，可以只研究其中面的变形； 对于矩形板单元，可以只研究一个矩形平面，但是，此单元上一点实际上代表着一个长度为板厚的法线段。 按基本假设，此法线段长度不变，其位移应包括中心点的挠度w和法线绕x、y轴的转角θx和θy 。 因而板单元任一节点i应有3个位移分量。 第十五页，共四十一页，2022年，8月28日 图5-3为矩形板单元，规定位移的正方向： w i沿z轴方向； 转角θxi和θyi绕x、y轴按右手螺旋规定正方向。 第十六页，共四十一页，2022年，8月28日 节点位移 按直法线假设，以小挠度变形下，法线的转角可由挠曲面的斜率表示。 i节点的3项位移可用列阵表示为 第十七页，共四十一页，2022年，8月28日 板单元的每个节点有3项独立位移，即有3个自由度，4个节点共有12个自由度。 如e单元4个节点的编号为k、l、m、n，则此单元全部节点位移可以列阵表示为 为分析方便，此顺序是按节点分组排列的。如按节点分块，上述节点位移应表示为 第十八页，共四十一页，2022年，8月28日 形状函数 取矩形单元的对称轴为x、y轴，可假定单元内挠度具有如下的多项式形式 其中a1、 a2… a12为待定系数。 (5.7a) 第十九页，共四十一页，2022年，8月28日 12个待定系数对应于单元的12个自由度。 前3项为常数项及线性项，反映出中面平板无弯曲的刚体位移。 3个二次项经二阶微分后给出常曲率，反映出中面变形的3种常应变形式。 因此，前6项满足了单元的完备性要求。 第二十页，共四十一页，2022年，8月28日 含有完全的三次多项式，其四次项是不完全的，此种近似的挠度函数具有三次多项式的精度。 不完全四次项的两项是对称的，这使单元对x及y轴具有同等的变形能力；当坐标轴转90o时，单元不会表现出不同的