有限元法基础董纪伟有限元第十章讲解.pptVIP

* 第 10 章 平板弯曲问题的有限单元法 （简介） 主 要 内 容 10.1 引 言 10.2 基于薄板理论的非协调板单元 10.3 基于薄板理论的协调板单元 10.4 Mindlin 板单元（位移和转动各自独立插值的板单元） 10.5 基于离散 Kirchhoff 理论(DKT)的板单元 10.6 应力杂交板单元 10.5 小 结 10.1 引 言 一、平板问题的特点 中面(xy 面) x y z 1. 几何特点 —— 薄板 h b 2. 载荷 —— 厚板 —— 薄膜 —— 垂直于板面作用 3. 截面上的内力 ——弯矩； ——扭矩； ——剪力； —— 薄膜力 (当板面内不受力时， = 0 ) (相对弯矩、扭矩很小 ) x y z 4. 截面上的应力 ——弯曲应力； ——扭转剪应力； ——剪力引起的剪应力 —— 垂直于板面方向的正应力 5. 变形描述与结点位移 —— 垂直于板面方向（z方向）的位移 挠度（w）: 板截面转角（?x、?y） 1 2 3 4 结点位移： 二、弹性薄板弯曲变形的基本假定 —— G. Kirchhoff 假定 1. 垂直于中面方向的线应变 ?z 忽略不计，即取： 2. 应力分量： 远小于应力分量： 即： 薄板的物理方程： 或： 3. 薄板中面内的各点都没有平行于中面的位移，即 1 2 3 4 x y z 由 得： 中面(xy 面) —— 中性平面 变形前后，平面在压 xy 面上的投影形状不变。 将： 对变量 z 积分，有 结论： 平板的全部应力、应变都可用板平面内的挠度 w(x,y) 表示 三、弹性薄板的基本方程 1. 几何方程 其中： 板弯曲的曲率矩阵 2. 物理方程 中面(xy 面) x y z h b 板的广义力： 由板内应力沿z方向线性分布的性质，有： 中面(xy 面) h b x y z 其中： —— 称为板的弯曲刚度（抗弯刚度） ——板的物理方程 将Mx、My、Mxy代入，得 3. 平衡微分方程 中面(xy 面) x y z h b 4. 边界条件 （1）在边界S1上，给定位移: 对固支边界有: （2）在边界S2上，给定位移 和力矩 ，即 （n为边界的外法线方向） 4. 边界条件 （1）在边界S1上，给定位移: 对固支边界有: （2）在边界S2上，给定位移 和力矩 ，即 （n为边界的外法线方向） 对简支边界有: （3）在边界S3上，给定力矩 和横向载荷 ，即 其中：S为边界的切线方向；Qn 为边界截面上单位长度的横向剪力。 （3）在边界S3上，给定力矩 和横向载荷 ，即 其中：S为边界的切线方向；Qn 为边界截面上单位长度的横向剪力。 对自由边界有: 四、弹性薄板的位能泛函 ?P 说明: ?P 中w导数的最高次数为2。 五、板壳有限单元法的步骤 （1）划分单元； （2）构造单元位移模式（插值函数）； （3）单元分析； （4）整体分析； 引入位移边界条件 （5）求解方程，结果处理。 说明： （1） ?P 中w导数的最高次数为2，由收敛性准则，在单元交界上，必须保证w及其一阶导数 的连续性。 即，保证w具有C1类连续性。 （2） 二维、三维问题的C1类连续函数构造非常困难。 （3） 三类平板单元： （a） 基于经典薄板理论的板单元—。——非协调元 （b） 基于 Kirchhoff 假设和广义变分原理的板单元 如： Hellinger-Reissner 广义变分原理——混合板单元； 修正的Hellinger-Reissner 变分（余能）原理——应力杂交元； 离散 Kirchhoff 假设板单元 （c） 基于 考虑横向剪切变形的 Mindlin 平板理论的板单元 Mindlin 平板理论——厚板理论 10.2 基于薄板理论的非协调板单元 10.2.1 矩形单元 1 2 3 4 x y z 1. 结点位移 每个结点三个自由度： 单元结点位移列阵： （10.2.1) （10.2.2) 2a 2b 2. 单元位移模式 （10.2.3) 1 2 3 4 x y z 2a 2b 刚体位移 常应变 线性应变 不完全的四次项 为确定待定系数： 将4个结点的坐标代入w及导数的表达式， 得到方程组： （10.2.4) 将上述方程组表示成矩阵形式： （10.2.5) C ——为12×12阶系数矩阵； 1 2 3 4 x y z 2a 2b （10.2.3) （10.2.5) 由此可解出： （10.2.6) 将其代入式（10.2.3）, 并用知阵表示： （10.2.7） 式中： （10.2.8） 其中： —— 单元中心的坐标 3. 应变矩阵（广义）