handout-10 板单元.ppt

* 定义广义应变为薄板中面的曲率或扭率： 薄板问题的有限元 内力 注意：这里 单元的节点位移向量为： 单元内挠度表示为插值形式： * 薄板问题的有限元 将以前各式代入，得到求解挠度的微分方程 平衡方程 或者写为 * 非协调矩形薄板单元的形函数推导 将单元节点坐标值和节点自由度代入上式，得到 解得 所以 * 非协调矩形薄板单元 上述矩形薄板单元虽不协调，但能通过分片试验，在单元不断细化时可收敛于精确解 不能推广至任意四边形 * 如何回避C1连续性要求？ 基于经典薄板理论的单元（以挠度为场函数；前已介绍） 基于保持Kirchhoff直法线假设的其他薄板变分原理的板单元，例如基于Hellinger-Reissner变分原理的混合板单元、基于修正H-R变分原理或修正余能原理的应力杂交板单元、基于单元内或单元边界上若干点（而不是处处）满足Kirchhoff直法线假设的离散Kirchhoff单元（如DKT）等。 离散Kirchhoff（DK）单元的思想：采用独立的场量对横向位移和中面法线转角进行描述，但通过在选定位置处施加剪切应变为零（即Kirchhoff直法线假设）的条件（场量显式耦合） 考虑横向剪切变形的Mindlin板单元 目前，平板单元大体上分为3类： * Mindlin板问题的有限元 Mindlin板单元（可看作8节点实体元派生的）考虑剪切变形，适用于描述中厚板 Mindlin板单元易于满足协调性要求Mindlin理论放弃了Kirchhoff薄板理论中关于中面法线变形后仍然垂直于变形后的中曲面的假设，认为变形前中面的法线，在变形后仍然为直线，但不一定垂直变形后的中曲面。根据Mindlin理论，板内任一点的位移可描述为 * Mindlin板问题的有限元 节点自由度向量为 挠度和转角各自独立插值 插值函数只需满足C0连续性 Mindlin板单元对中面挠度和中面法线转角进行独立插值： 即 * Mindlin板问题的有限元 其中， 对于各向同性线弹性材料，Mindlin板的力矩、剪力与曲率、中面转角关系式是 kt 代表了剪切变形的等效厚度; ??反映了剪应力沿??方向的非线 性变化，可以认为??=5/6 即 * Mindlin板问题的有限元 单元内任一点的曲率和横向剪切应变为 * Mindlin板问题的有限元 单元应变能为 其中，单元刚度阵为 单元刚度阵还可表示为弯曲刚度和横向剪切刚度之和： * * * * * * * * * * * * * * * * * * * *