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

《板壳弯曲问题的有限单元法.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第6章薄板弯曲问题的有限单元法薄板弯曲问题的基本方程薄板弯曲问题的非协调矩形单元非协调三角形板单元薄板弯曲问题的协调元6.1薄板弯曲问题的基本方程1弹性薄板的基本假设（克希霍夫假设）无挤压薄板弯曲时，平行于中面的各层面之间无挤压。这意味着薄板弯曲后厚度保持不变，因此可取。显然挠度w只是x,y的函数：直法线变形前垂直于中面的直线段，变形后仍为直线，且仍然垂直于弯曲后的中面。这意味着yz和zx平面内的剪应变为零从而得：无侧移薄板中面内各点都没有平行于中面的侧向位移，即结合几何方程可知，中面内形变分量均为零,即从上述的附加假设出发，可以将位移u、v

2、用w表示。推导得这就是薄板弯曲问题的克希霍夫(Kirchhoff)假设，使用克希霍夫假设计算的板称为克希霍夫板。将用w表示的位移u,v代入几何方程这里，记为称为薄板的广义应变分量。薄板中的应力[D0]是平面应力问题的物理矩阵．薄板内力[D]是板的弯曲刚度矩阵．显然最大应力发生在薄板的上下表面2弹性薄板的几点简化应力分量的减少应变分量的减少位移之间有了附加关系应力应变关系的简化1薄板弯曲问题节点位移参数的选择采用克希霍夫假设后，薄板的变形状态完全由一个变量，即中面挠度w(x,y)来确定。然而，在有限元法中只取挠度本身作为节点位移参数是不够的。按

3、克希霍夫理论，薄板内部非中面上各点的位移(u,v,w)是用相应的中面点的挠度w(x,y)和该点处中面法线转角θx和θy来表示的(2式)。因而，为了保证板内位移(u,v,w)在整个求解区域内单值连续，除要求w在全域内单值连续外，还必须要求θx和θy在全域内也是单值连续的。这里6.2矩形薄板单元将只要求函数本身连续的问题称为C0问题，如弹性力学平面问题；将不但函数本身，还要求其一阶导数连续的问题称为C1问题，如薄板弯曲问题。如果将位移模式仍然取为多项式，要求在全域内位移及一阶导数连续，这等价于在单元边界上要保证位移及一阶导数连续，因此在单元结点上

4、必须保证位移及一阶导数连续，即应选取三个结点位移参数如果取四节点单元，则取位移函数为两个四次项的选取，保证了在单元边界上，即x=const，y=const时，位移是三次多项式。位移连续性问题。在ij边上，y=const，共有四个参数，可由ij边两端节点的位移参数唯一确定，因此在相邻单元的公共边界上，位移w及其切向导数是连续的。仍有四个参数，但是节点参数只有两个，无法唯一确定法向导数。也就是说，在两个相邻单元的公共边界上，位移模式w的法向导数并不相同。再来看法向导数。法向导数为由以上讨论和进一步的研究可以得出结论，仅规定位移w及其一阶导数作为节

5、点位移参数时，取位移模式为简单多项式，要保证单元边界上位移w的法向导数连续是不可能的，常称这样的单元为不完全协调元。不完全协调元的位移模式只满足了“收敛准则”的完备条件，而未满足协调条件。有关其收敛性的问题需要再讨论。但是计算实践表明，这里所给出的不完全协调四节点矩形单元的计算结果是收敛的．2位移模式将矩形薄板沿坐标方向划分为若干矩形单元，每个单元设有四个节点，每个节点有三个位移分量，即挠度w，绕y轴转角θx,绕x轴转角θy。即单元的节点位移为节点荷载为单元的节点荷载为取位移函数为在位移函数中，前三项包含了单元的刚体位移状态，二次项代表了单元

6、的均匀应变状态。可以证明，此位移模式能够保证相邻单元的公共边界上挠度w和转角的连续性。分别求出上式中对x,y的导数将单元四个节点的坐标分别代入前三式后，可得12个关于αi的方程组，求解后代回(7)式，令其中称[N]为形函数矩阵，第i个子矩阵为为节点的坐标值，则将形函数(9)代入(3)式，得出这里的[B]称为应变矩阵第i个子矩阵[Bi]为3势能泛函与有限元模式板的势能泛函可写成将(10)式代入(11)式得按最小势能原理将(12)式代入(13)式得记得出4不完全谐调元的分片检验前面说明，薄板不完全协调矩形单元的位移插值函数不能满足“收敛准则”所要

7、求的协调条件，但是计算结果表明是收敛的。如何判断此种不完全协调元计算结果的收敛性呢?埃恩斯提出“分片检验”的概念，并指出：位移插值函数能否通过“分片检验”，是判断不完全协调计算结果是否收敛的充分必要条件。“分片检验”的具体做法如下，任意取一个至少有一个内部节点的，由若干个单元组成的拼片，并且：在内部节点上既不允许有载荷，也不允许有约束。当把任何一种与常应变状态对应的节点位移或节点力加到该单元拼片的边界节点上时，用某种位移插值函数计算得到单元拼片内部的位移符合常应变状态的条件，则说该位移插位函数能够通过“分片检验”。经检验表明，前面介绍的不完全

8、协调矩形元能够通过“分片检验”，因而计算结果是收敛的。6.3三角形板单元三结点板单元，每个结点三个位移参数每个单元共有9个参数。如果位移函数取为多项式，则一个完备的