有限元4-薄板弯曲问题研究.ppt

有限单元法 第4章  弹性薄板弯曲问题的有限元法 板壳结构 平板上所承受的荷载分类 大挠度问题与小挠度问题 薄板与中厚板 当最大挠度w远小于t 时, 称为小挠度问题（or刚性板） 当最大挠度w与t相差不大时,称为大挠度问题（or柔性板） 4.1 基本理论 一、基本假定 １、略去垂直于中面的法向应力。（ ），并以中面上沿Z方向的挠度w代表板的挠度。 ２、变形前垂直中面的任意直线，变形后仍保持为垂直中面的直线。（─法向假定, ） ３、板弯曲时，中面不产生平面内应力。(─中面中性层假定， ) 上述假定常称为薄板小挠度问题假定（or 柯克霍夫Kirchhoof假定）。符合上述假定的平板即为刚性板。 二、基本方程 以上述假定为基础，板分析中常用挠度w(x,y)作为基本未知量。 １、几何方程（应变─挠度关系） ①弹性曲面沿x, y 方向的倾角 设A点挠度为w, 则沿x方向倾角(绕y轴) ③ z平面的应变分量和曲、扭率 基本假定，由于 , 故板内任意点的应变与平面问题相同: 上式中的 , , 为曲面在x,y方向的曲、扭率。 这三个参数称作弹性曲面的弯扭变形分量，它们完全确定了薄板内各点的应变分量。记为： ２、物理方程（应力─挠度关系） 由于忽略σz 对变形的影响, 因此z平面的应力─应变关系具有与平面问题相同的形式: 板与梁的比较 4、板的平衡方程 在经典解析法中，w(x, y)常设为三角级数形式。例如，四边简支矩形板的w(x, y)设为: (纳维尔解) 4.2 有限元分析方法 一、矩形单元的典型形式 将图示矩形薄板沿x,y方向划分成若干小矩形(常取等分)。从中取出一小矩形(单元)，共有四个结点，此时不能象在平面问题中一样，将结点视为“铰”，而是“刚性的”，即每个结点有三个位移分量： 即结点i的位移 符号重新定义是为了有限元表示的方便，由此得单元结位移向量 二、 位移模式（函数） １、位移模式的选取 插值多项式取为： ２、位移模式的检验 （三个基本要求： 刚体位移，常应变，尽可能的边界协调） ① 前三项含单元的刚体位移状态： 第一项α1 与坐标x, y无关, 表示z方向的挠度是─常量, 刚体移动 ③ 能保证相邻单元在公共边界上挠度的连续性 以单元的１～２边界为例，在此边界上y=-b，为常量，代入位移模式 (4-2-1)，可知边界上的挠度w是x的三次函数，合并整理后可得： ④ 不能保证相邻单元在公共边界上法线转角的连续性 仍以1-2边界为例，将y=-b代入后，此时 三、形函数和形函数矩阵 分别将单元结点1, 2, 3, 4的坐标值代入(4-2-1) ，并事先求出 ,便可得 到各结点的位移值。一共可得12个关于αi 的方程组，联立求解可得αi 。 如果采用局部坐标（ξ，η）并用单元节点位移值来表示位移函数，则 式中形函数: 四、单元的几何矩阵[B]和内力矩阵[S] 1．几何矩阵[B] 由前可知 ,将（4-2-2）的w代入,得到几何矩阵： 2.内力矩阵[S] 由基本方程（4-1-5） 可得到： 五、单元刚度矩阵 由一般公式得: 六、荷载等效变换 由荷载等效变换的一般公式可得 八、计算例题 例题1 计算图示四边固定方板。 方板的边长为l，厚度为t，弹性模型量为E，波松比μ=0.3，全板承受均布法向荷载ｑ，求薄板中的挠度和内力。 此时，单元的a, b是:a=b=l/4 计算节点荷载: 由前面的均布荷载计算公式得： 解出 4.3 薄板有限元程序设计 一、总框图 根据弯曲板有限元分析方法的解题过程，可写出其总框图如下： 二、子框图 １、单元坐标结点编号及单刚形式 为了取挠度向下为正，又能与前述坐标系统统一，特将前述坐标前翻180° 程序中单刚数组为 DK(20, 20), 子程序: Subroutine DG(A, B, E, T, U)为其形成单刚的子程序。 ２、自动形成单元编号信息（ 单元信息数组：[IB]）。 ３、结点定位向量。 ４、形成荷载列向量。（ a. 结点力； b. 非结点力(只考虑均布力)） ５、总刚，Subroutine ZG(M, N, LD, A, B, E, T,U) ６、解方程。 FJZG( ), HUD( ) ７、算单元力。 Subrout