有限单元法基础教学课件U8.pptxVIP

第8章 板壳结构力学问题 8.1 引言 8.2 薄板单元 1 薄板弯曲基本理论基本假设中面（z=0）上的点面内位移为零受横向力和弯曲载荷作用的平板垂直于中面的法线上的点具有相同的竖向位移垂直于中面方向的正应力可以忽略不计变形前正交于板中面的直线变形后仍正交于中面且长度不变板中面外的其它点的面内位移直法线假设薄板中任意一点的位移薄板的变形特征板中任意一点的应变横向剪切应变为零广义应变（曲率和扭率表达）面内应力分量薄板微元体及应力分布广义应力应变关系（本构关系）板的抗弯刚度板的抗弯刚度矩阵平衡方程基本方程几何方程本构方程平衡方程板弯曲问题的三种边界条件位移边界n:法向；s:切向外力边界混合边界条件或板的边界条件弹性薄板的弯曲应变能板弯曲时为平面应力状态，其应力和应变关系为系统的势能泛函2 基于薄板理论的平板单元结点变量插值函数绕y轴的转角：绕x的转角：离散化：w矩形4结点单元对板的中面进行离散，中面以外其它点的位移由中面挠度计算得到。结点位移：单元位移插值模式：代入结点位移和坐标：利用坐标变换xo和yo为单元中心点的整体坐标仅作用有竖向均布载荷q收敛性检查a）完备性检查z 向刚体位移刚体位移：绕 x 轴刚体转动绕 y 轴刚体转动常应变：满足完备性要求。b）连续性检查w 连续；非协调单元一般不连续。讨论：a）为非协调单元。b）能通过分片试验，能收敛于精确解，但收敛不是单调的。c）不能推广到一般的四边形单元。3三角形3结点单元3结点，9个参数每个结点自由度：2若用完备三次项，应包含10项1用面积坐标构造位移函数w 协调；可以证明：当c=1/2时，满足常应变要求；不协调。非协调单元8.3 厚板单元 1 厚板弯曲基本理论基本假设变形前正交于板中面的直线变形后仍为直线，其它三个假设与薄板相同。中面位移：厚板界面的转动关系中面外其它点位移中面法线的转角板中任意一点的应变层间剪切应变曲率和扭率（广义应变）板的面内应力状态为平面应力状态厚板微元体及应力分布平衡方程抗弯刚度抛物线分布通过修正的剪应变能等于按实际剪应力和剪应变计算得到的应变能相等，可得到矩形截面的k=5/6。k为截面剪切校正因子广义应力应变关系（本构方程）平衡方程板的应变能弯曲应变能剪切应变能板的势能泛函第一项：弯曲应变能其中挠度和转角均为独立变量，由此构造的有限元格式适用于较厚的平板弯曲问题。该类单元的位移和转动各自独立插值，称为 Mindlin板单元。剪切应变能第二项：后几项：外力势板很薄时，忽略层间剪切变形可以理解为剪切刚度罚函数法第二项可以理解为利用罚函数法引入薄板约束条件2 Mindlin板单元插值：C0连续函数，二维C0型插值函数均可。由驻值条件剪切锁死和零能模式当板很薄时， 为避免剪切锁死，Ks需奇异。 采用减缩积分可能致使K奇异，使解答中包含除刚体运动之外的且对变形能无贡献的变形模式——零能模式a）Midlin板单元保证K非奇异的必要条件：高斯积分点数；Me：单元数；：弯曲应变和剪切应变分量数；N：系统独立自由度数。b）保证Ks奇异性的充分条件积分方案的选择分别采用不同阶次的积分方案。应用中的处理实际应用和通用程序中，普遍采用减缩积分方案。零能位移模式的出现：板较厚且单元较少，边界约束较少时容易出现。较厚板：精确积分，不会出现零能模式和锁死。应用中薄板分析时，推荐采用4结点和9结点Lagrange单元。8.4 平面壳体单元1 弹性薄壳基本假设（Kirchhoff-Love假设）直法线假设 在变形前正交于壳体中面的直线，变形后仍然为正交于中面的直法线段，且保持长度不变。切平面应力假设正应力比应力分量小得多，相对而言可以忽略，即：模拟一般三维空间壳结构的平面壳体单元用矩形单元组成曲壳用三角形单元组成曲壳用三角形和矩形组合代替曲壳壳体与平板的异同a）相同点a）厚度方向的尺度比其它两个方向小得多b）均使用Kirchhoff假设b）不同点平板：u，v 位移[平面应力（薄膜应力状态）]与弯曲状态不耦合。壳体：u，v，w 同时发生，弯曲状态与薄膜状态相互耦合。1 局部坐标系下的平板单元平面壳元＝平面应力单元＋平板弯曲单元结点位移：结点力：1）平面应力状态（薄膜应力状态）三角形平面薄壳单元的结点位移和结点力2）平板弯曲状态在局部坐标系下在整体坐标系下2 向整体坐标的转换局部坐标系下整体坐标系下平面壳单元局部坐标系和整体坐标系1）结点位移变换L：正交矩阵2）单元结点位移变换3）单元刚度矩阵变换4）单元载荷向量的变换组集刚度矩阵后特殊情况的处理：若汇交于结点 i 的各单元在同一平面内，处理方法：（1）在局部坐标系建立该点的平衡方程，删去对应的方程（2）在此点上给出任意的刚度系数，不会影响计算结果。8.5 超参数壳体单元 1 几何形状三维实体单元用于壳体结构时