薄板弯呛旺问题的有限单元法.pptVIP

第6章 薄板弯曲问题的有限单元法 ;6.1 薄板弯曲问题的基本方程;从而得：;将用 w 表示的位移 u ,v 代入几何方程 这里，记 为;薄板中的应力;2 弹性薄板的几点简化 应力分量的减少 应变分量的减少 位移之间有了附加关系 应力应变关系的简化;1 薄板弯曲问题节点位移参数的选择 采用克希霍夫假设后， 薄板的变形状态完全由一个变量，即中面挠度 w(x, y) 来确定。然而，在有限元法中只取挠度本身作为节点位移参数是不够的。 按克希霍夫理论，薄板内部非中面上各点的位移 (u,v,w) 是用相应的中面点的挠度w(x, y)和该点处中面法线转角θx和θy来表示的(2式)。因而，为了保证板内位移 (u,v,w) 在整个求解区域内单值连续，除要求w在全域内单值连续外，还必须要求θx和θy在全域内也是单值连续的。这里;将只要求函数本身连续的问题称为C0问题，如弹性力学平面问题；将不但函数本身，还要求其一阶导数连续的问题称为C1问题，如薄板弯曲问题。;两个四次项的选取，保证了在单元边界上，即x=const，y=const时，位移是三次多项式。;位移连续性问题。 在 ij 边上，y=const，;由以上讨论和进一步的研究可以得出结论，仅规定位移 w及其一阶导数 作为节点位移参数时，取位移模式为简单多项式，要保证单元边界上位移 w的法向导数连续是不可能的，常称这样的单元为不完全协调元。不完全协调元的位移模式只满足了“收敛准则”的完备条件，而未满足协调条件。有关其收敛性的问题需要再讨论。但是计算实践表明，这里所给出的不完全协调四节点矩形单元的计算结果是收敛的． ;2 位移模式 将矩形薄板沿坐标方向划分为若干矩形单元，每个单元设有四个节点，每个节点有三个位移分量，即挠度 w，绕y轴转角θx, 绕 x 轴转角θy。即;单元的节点位移为 节点荷载为 单元的节点荷载为;取位移函数为;将单元四个节点的坐标分别代入前三式后，可得12个关于αi 的方程组，求解后代回(7)式，令;将形函数(9)代入(3)式，得出;3 势能泛函与有限元模式 板的势能泛函可写成;将(12)式代入(13)式得;4 不完全谐调元的分片检验 前面说明，薄板不完全协调???形单元的位移插值函数不能满足“收敛准则”所要求的协调条件，但是计算结果表明是收敛的。如何判断此种不完全协调元计算结果的收敛性呢?　 埃恩斯提出“分片检验”的概念，并指出：位移插值函数能否通过“分片检验”，是判断不完全协调计算结果是否收敛的充分必要条件。 “分片检验”的具体做法如下，任意取一个至少有一个内部节点的，由若干个单元组成的拼片，并且：在内部节点上既不允许有载荷，也不允许有约束。当把任何一种与常应变状态对应的节点位移或节点力加到该单元拼片的边界节点上时，用某种位移插值函数计算得到单元拼片内部的位移符合常应变状态的条件，则说该位移插位函数能够通过“分片检验”。 经检验表明，前面介绍的不完全协调矩形元能够通过 “ 分片检验 ”，因而计算结果是收敛的。 ;6.3 三角形板单元;面积坐标的一、二、三 次式分别为;取三角形板的位移插值函数为;按前述方法求出 ，位移函数表为;6.4 基于薄板理论的协调元