薄板理论的板单元.ppt

其中： 可通过轮换下标：1——2——3 得到。 3. 单元特性矩阵 （1） 单元应变矩阵 —— 广义应变矩阵 其中： （2） 单元刚度矩阵 （3）单元等效结点载荷列阵： 注：需将上述积分转化为面积坐标下的积分。 4. 单元的收敛性 （1） 完备性 由位移模式的构造过程可知，上述位移函数w能描述刚体位移、常应变，所以，满足完备性要求。 （2） 协调性 相邻单元挠度 w 的连续性 单元边界上，w为三次变化函数，需由4个参数确定。单元每个节点存在： 1 2 3 x y z w 每边2个节点，具有4个上述参数，可唯一确定w 的函数，即 w 满足连续性要求。 相邻单元法向导数 的连续性 因为，w为三次变化函数，所以 为二次变化函数，需由3个 参数才能唯一确定。 而单元边界上，只存在 2个法向导数量，如： 1、 2边界上，只存在法向导数： 即单元边界上，法向导数不协调。 1 2 3 x y z w 3节点三角形板单元 —— 非协调单元 对这种单元，当单元尺寸不断缩小时，不能保证其解收敛于精确解。其收敛性保障需通过分片试验。 单元形状对收敛性的影响： Irons等证明， 若单元网格是由三组 等间距直线产生的，如 4×4, 4×4A 可通过单元的分片试验，因而解收敛。 对网格 4×4B 虽然也收敛，但位移值误差达1.5%。而不规则网格误差更大一些。 结果： （1）沿板中心线的挠度； （2）沿板中心线的 Mx 。 例：不同支承条件的方板承受分布载荷或集中载荷。 10.3 应力杂交元介绍 早期的应力杂交元主要用于薄板的有限元分析，避免构造C1类协调单元的困难。基于下面的修正的余能原理： 此变分原理要求单元应力场事先满足域内的平衡方程，单元边界位移事先满足给定的位移约束。 假定 其中 是单元节点位移向量, N, M分别是单元边界和单元内应力的插值函数矩阵。 b是单元内的应力参数, 对平面问题, 单元应力场可假设为： 上面的假设保证了单元内平衡方程的满足。 另一方面，因为位移只是在单元的边界上假设，即使是要求C1类连续，也不再困难。 将假设单元应力和位移代入修正的泛函，得到： 单元边界上的力也可以用向量表示为 其中: j为单元边界方向余弦组成的矩阵. 其中: 注意到应力参数b定义在每个单元内部，可在单元一级消去： * 第 10 章 平板弯曲问题的有限单元法 （简介） 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 节点位移： 二、弹性薄板弯曲变形的基本假定 ——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上，给定位移: 对固支边界有: （n为边界的外法线方向） 对简支边界有: （3）在边界S3上，给定力矩 和横向载荷 ，即 其中：S为边界的切线方向；Qn 为边界截面上单位长度的横向剪力。 对自由边界有: （2）在边界S2上，给定位移 和力矩 ，即 四、弹性薄板的位能泛函 ?P 说明：?P 中w导数的最高次数为2。 五、板壳有限单元法的步骤 （1）划分单元 （2）构造单元位移模式（插值函数） （3）单元分析 （4）整体分析 引入位移边界条件 （5）求解方程，结果处理。 说明： （1） ?P 中w导数的最高次数为2，由