第七章平板弯曲问题的有限元分解.ppt

关于平板壳元 平板壳元是平面应力单元和平板弯曲单元的组合。平板弯曲单元稍加扩充就可以应用于壳体分析。 然而用折板代替壳体，网格需要合理的密度才能得到满足实际要求的计算精度。 采用曲面壳元能够更好地反映壳体的真实几何形状，通常可以得到比平板壳元更好的结果。 曲面壳元 * 曲面壳元 基于薄壳理论的曲面壳元需要构造具有C1连续性同时满足完备性要求的插值函数是非常困难的。 基于Mindlin壳理论，构造位移和转动独立插值的曲面壳元，只要求满足C0连续性，自然要容易的多，而且这类壳单元经过适当的处理也可以进行薄壳分析（注意剪切锁死）。 下面介绍一类从三维实体单元退化而来的超参数壳元（最简单的4节点单元） 。 * 退化的Mindlin超参壳元 1 2 4 3 在三维实体中引入壳体理论假设，将壳体上下表面上一对节点的6个自由度退化为壳体中面上一个节点的5个自由度。 * 节点坐标： 节点中面法线单位向量： 坐标插值函数的构造 * 如果平板单元受有分布横向载荷q的作用，于是等效结点力是 这里以使用(d)式所示的形函数较为方便。 3．三角形单元的等效结点力和内力矩的计算 * 当q = q0是常量时，将(d)式代入上式，并积分得 * 对于内力矩阵，按下式计算 式中 (7-43) * 把上式代入(7-43)式并进行一系列运算之后，可以得到内力矩列阵 式中 对于任意形状平板，它的边界条件可能是指定为沿着曲线边界切线方向的弯矩Ms或转角 和沿着法线方向的扭矩Mn或转角 ，这里所用的记号与一般平板理论书中的记号恰好相反。这里取边界的外法线方向n为正方向，而使用右手坐标系定出切线s的正方向。于是，有下列关系式存在 * 式中c、s是外法线n方向对于x和y轴的方向余弦。 利用上式可以对刚度矩阵作某种变动，使之达到结点位移和结点力的变换，从而直接利用曲线边界上的边界条件。这种变动可以在单元刚度矩阵中进行，也可以在整体刚度矩阵中进行。 * 上述矩形单元和三角形薄板单元都是非协调元，且都会产生在单元公共边上的法线转角 （或法向斜率 ）不连续问题。为了实现薄板单元的协调性（协调板元或保续板元），完善板元的计算理论和提高计算精度，人们做了大量研究，提出了很多方法。其中一种方法是增加单元结点自由度数目，例如在矩形板元的每个结点上增加一个扭率 ，使单元变成16个自由度的协调板单元；在三角形板单元中每个结点上增加 、 、 做为结点自由度，且在每边中点取其法线斜率 做为结点自由度、而构造出21个自由度的协调三角形板单元。 * 7.3 考虑横向剪切变形影响的平板弯曲单元 另一种方法，是将三角形（或四边形）板单元划分为3个子三角形，把每个三角形各边中点的法线斜率 做为自由度，每个子三角形都有12个自由度，使各子三角形板元之间协调，然后再根据原三角形板元内部的连续性和限制条件，利用“凝聚法”消去内自由度，从而构造出协调的三角形板元 （12个外自由度，3个内自由度）或协调的四边形板元。 * 第一种方法的明显缺点：①在实际应用时涉及到高阶导数的自由度的边界条件，难以处理；②它不是一种普遍适用的方法，在某些板问题中，曲率或扭率在结点上不一定连续，在有限元计算中强令其连续，当然会使所得结果不可能收敛到精确解。 第二种方法在SAP-5程序中得到应用。但该法的计算公式比较复杂，消去每个单元的内自由度的凝聚过程所耗计算时间比较长。 在板弯曲问题的有限元法中，构造协调元的困难是单元公共边上的法向转角 的连续性难以满足，如果考虑板横向剪切变形的影响，放弃经典薄板理论中的中面法线n-n始终保持为直线的假设，就可以绕开这个困难，而使板问题的有限元分析前进一步。 * 一般来说，板变形前的中面法线n-n，在变形后将变成一条曲线，但可以近似地用一条直线n2-n2表示，n2-n2线已不是经典薄板理论中变形后的中面法线n’-n’。此时，n2-n2线绕x轴和y轴的转角仍用 和 表示，但 。 经典薄板理论中的其它两个假设，仍然有效。 * 根据上述假定，板内任意点的3 个位移分量具有如下形式 把上式代入几何方程式，可以得到平板应变列阵是 (7-45) (7-44) * 通过应力应变关系，可以得到应力列阵如下 式中弹性矩阵是 而它的子矩阵 (a) (7-46) (7-47) * 由