第四章 空间问题有限单元法 有限单元法与程序设计 教学课件.ppt

第五章 板壳问题有限单元法 三、三角形薄板单元 1. 位移插值函数： 需要有10个系数，在直角坐标下问题很难解决，较好的方案是设： 三角形板单元有9个自由度， 而一个完全三次式： 第五章 板壳问题有限单元法 三、三角形薄板单元 1. 位移插值函数： 但当三角形两边分别平行两坐标轴时，确定广义坐标的系数矩阵奇异 利用面积坐标三个分量不相互独立的特性，可解决该问题： 第五章 板壳问题有限单元法 三、三角形薄板单元 1. 位移插值函数： 利用结点的位移参数条件可确定w中的广义坐标，得到： 形函数矩阵 第五章 板壳问题有限单元法 三、三角形薄板单元 1. 位移插值函数： 形函数 试证明9自由度三角形薄板单元为非协调元 第五章 板壳问题有限单元法 三、三角形薄板单元 2. 单元分析 试推导9自由度三角形薄板单元的 1）应变矩阵 2）应力矩阵 3）单元刚度矩阵 y x L L 例：四边简支正方形薄板 A B D C y x L L A B D C 受均布荷载q及中心集中荷载P两种工况作用，分别用矩形单元和三角形单元计算最大挠度 第五章 板壳问题有限单元法 四、用矩形薄板单元进行薄壳分析 1. 薄壳结构受力特点： 1） 薄壳的概念： 两个曲面所限定的物体，如果曲面之间的距离比物体的其他尺寸小，就称为壳体。 t/R 0.05－薄壳 2） 薄壳的受力特点（相对于薄板)： a) 在荷载作用下除产生弯曲变形外，还存在中面变形。 b) 弯曲变形和中面变形同时发生，且相互耦合。 第五章 板壳问题有限单元法 一、薄板弯曲基本假定和基本方程 2. 薄板弯曲问题的基本假定－克希霍夫假定： 1）中面法线变形后既不伸长也不缩短； 2）板中面法线变形前是直线，变形后仍保持直线，且与变形后的中面保持垂直； 3）中面各点没有平行于中面的位移。 假定（1）与梁弯曲问题的互不挤压假定相似 ?z=0 即： w=w(x,y) 所以： 第五章 板壳问题有限单元法 一、薄板弯曲基本假定和基本方程 2. 薄板弯曲问题基本假定： 假定（2）与梁弯曲问题的平面假定相似，即剪切应变： ?zx=?zy=0 即： 有： 利用：w=w(x,y) 第五章 板壳问题有限单元法 一、薄板弯曲基本假定和基本方程 2. 薄板弯曲问题基本假定： 所以： 再使用假定（3），得： f1(x,y)=0， f2(x,y)=0 第五章 板壳问题有限单元法 一、薄板弯曲基本假定和基本方程 3. 薄板弯曲问题的应变： ?x=Xxz ?y＝Xyz ?xy＝2Xxyz ?z = ?yz = ?zx ＝0 六个应变分量中，根据假定，已知： 其余三个分量： 第五章 板壳问题有限单元法 一、薄板弯曲基本假定和基本方程 3. 薄板弯曲问题的应变： 曲率： 扭率： 薄板的形变分量 第五章 板壳问题有限单元法 一、薄板弯曲基本假定和基本方程 4. 薄板弯曲问题的应力： ( ?x、?y、?xy）通过平面问题的物理方程由应变求出 ( ?z、?zx、?zy）则必须由三个平衡微分方程求解给出 需注意： 应力分量（?z、?zx、?zy）尽管相对面内应力分量（?x、?y、?xy）很小， 它们对应的应变分量?z、?zx、?zy可略去不计，但它们本身由于是平衡所必须的而不能忽略不计。 第五章 板壳问题有限单元法 一、薄板弯曲基本假定和基本方程 4. 薄板弯曲问题的应力： 应力分量（?x、?y、?xy）： 特点： 均沿厚度呈线性分布，在中面处为零， 在板的上、下板面达到最大。 第五章 板壳问题有限单元法 一、薄板弯曲基本假定和基本方程 4. 薄板弯曲问题的应力： 应力分量（ ?z、 ?zx、?zy ）： 考虑薄板上、下板面的边界条件 解得横向剪应力，为 第五章 板壳问题有限单元法 一、薄板弯曲基本假定和基本方程 4. 薄板弯曲问题的应力： 特点： 横向剪应力?zx、?zy沿板厚度方向呈抛物线分布， 在板的上、下板面为零，在板中面最大。 第五章 板壳问题有限单元法 一、薄板弯曲基本假定和基本方程 4. 薄板弯曲问题的应力： 将z方向所有力作用等效移置到板面上，板上、下表面的边界条件变成 利用z方向的平衡条件求?z 利用边界条件可解得： 第五章 板壳问题有限单元法 一、薄板弯曲基本假定和基本方程 4. 薄板弯曲问题的应力： 特点： ?z沿板厚度方向呈三次方变化 最大值发生在板面为q，最小值在板底为0。 第五章 板壳问题有限单元法 一、薄板弯曲基本假定和基本方程 5. 薄板弯曲的平衡微分方程： 上式中，利用板下面的边界条件