哈尔滨工业大学 04 结构力学——平面有限元1.pptVIP

结构力学II 代入结点座标得 引人 则有 广义座标法和试凑法得到的形函数完全相同，但试凑法要方便一些。 ? ? 单元位移函数 ? 单元协调性 位移函数中包含常数项和一次项，故满足完备性条件；边界上线性变化且由结点确定，故满足位移连续。所以矩形双线性单元是收敛的。 ? 5-3-2 应变和应力矩阵 应变矩阵 ? 应变矩阵 ? 如在 ij 边 3. 在三角形单元内任一边上的形函数只与该边的两端点坐标有关且是线性的，而与其它结点坐标无关。 2 1 i、j边的直线方程为 在i、j边上形函数 根据性质1 ? 注意应用 则有 2 1 ? 在 ij 边的位移 由此性质可知：单元边界是线性变化的。相邻单元在公共结点有相同位移，所以能保证相邻单元的位移协调。 ? 5-2-2 用面积座标建立单元位移场 为以后方便理解和确定复杂单元内任意一点处形函数的值，在此引入面积坐标的概念 在单元内任一点 P(x , y) 可用三个数（Li，Lj，Lk）来确定 ? 因为 所以 可见面积坐标就是形函数 ? 由面积坐标的定义不难发现，平行于某边直线（ik）上的所有点具有相同的坐标（Lj），并且该坐标就等于“该直线到至（ik）边的距离”与“点 j 到（ik）边的距离”之比。 ? 三个结点的面积坐标分别为： 不难验证面积坐标与直角坐标的变换关系 ? 5-2-3 基于势能原理的单元分析 微分算子矩阵 ? 几何矩阵只与结点坐标有关，它的元素都是常数，所以单元内每点的应变也必为常数，所以三结点三角元称常应变三角元. ? 应变势能 外力势能 体系势能 ? 由最小势能原理 ? 荷载列阵 （1）单元上任意点的集中荷载 直接加在相应位移方向上 非结点荷载等效到结点上 所以 ? （2）单元上体积力 所以 ? 结果分析 单元应力 ? ?x ? xy ? yx ?y ? n ? 应力圆半径 平均应力 最大主应力 最小主应力 主方向角 ? 2? ?1 ?2 ?y ?m ?x ? ? O ?xy ? ? 5-2-4 几点结论 一.结点的选择和单元划分 1.集中力作用点、分布力突变点、支承点应选作结点。 2.不同厚度、不同材料的部分不应划在同一个单元。 3.应力变化大处单元应密集一些。结点的多少与疏密要考虑计算机的容量和计算精度。 有限元分析应注意的问题和结果整理 ? 4.单元边界的边长之比应尽可能靠近1。 不宜 5.相邻单元的尺寸尽可能接近。 6.结点所连接的单元个数尽可能一致。 适宜 适宜 不宜 ? 二.结点编码 尽可能使相关结点的结点编码差值最小，以缩小刚度矩阵的带宽。从而提高储存和运算效率 总刚半带宽=(相关结点最大差值+1)*结点位移数 ? 总刚半带宽=(相关结点最大差值+1)*结点位移数 总刚半带宽=16总刚需占用的存贮空间为: 16?14?2=448 总刚半带宽=6总刚需占用的存贮空间为: 6?14?2=168 1 2 3 4 5 6 7 8 9 10 11 12 13 14 1 3 5 7 9 11 13 2 4 6 8 10 12 14 ? 三.充分利用结构的对称性 FP FP FP FP FP ? 四.应力结果的整理 位移的计算结果一般比应力、内力结果精度高。位移达到满意结果，由几何方程求应变，再由物理方程求应力，结果的精度较差。上述三角形单元为常应力，矩形单元应力线性变化，而工程问题的应力是比较复杂的。为更好地反应实际应力情况，需要对计算结果进行整理。常用处理方法有两种:绕结点平均法和两单元平均法。 ? 以交于同一结点各单元此结点处某应力分量的代数平均值，作为此结点该实际应力的近似值。对于边界处的结点，由内结点结果的外得到。 结点4的应力由结点1、2、3的应力外插得到 1、绕结点平均法（相关单元应力平均值作为结点应力值） A B C D E F 1 2 3 4 ? 2、两单元平均法（相邻单元应力平均值作为边界应力值） 三角形单元时，以两相邻单元应力平均值作为边中点的应力近似值。矩形单元时，以两相邻单元公共边两端结点四个应力的平均值作为边中点的应力近似值。对于边界处的结点，同样由内结点结果的外插得到。 A B C D E F 1 2 3 4 5 6 7 8 ? 5-2-5 计算实例 例题1 Page 104 ? 5-2-6 收敛准则 ? 希望随着离散网格的逐步细分，所得到的解答能收敛于问题的精确解。 ? 单元形状确定后，位移模式的选择非常关键（因为荷载的位置、应力矩阵、刚度矩阵的建立都依赖于单元的