四节点矩形单元有限元分析.ppt

平面问题有限元分析四节点矩形单元 天津大学建筑工程学院 tianjinuniversity 本节内容提要 1 分析提高有限元法求解精度的途径2 简要回顾三节点三角形单元有限元分析过程3 全面介绍四节点矩形单元有限元分析过程4 总结 分析提高有限元求解精度的途径 tianjinuniversity 一 三节点三角形单元的缺点 三节点三角形单元精度低 收敛慢 由于单元内应力和应变均为常量 故在单元内不能很好地反映应力和应变的变化 该单元只有三个节点 单元自由度少 单元位移插值函数 位移模式 只能是线性函数 描述单元内位移变化的能力差 分析提高有限元求解精度的途径 tianjinuniversity 二 提高有限元求解精度的途径 第一个途径是对某一种特定类型的单元采用网格划分加密 依靠单元的收敛性提高求解精度 第二个途径是对一定的单元网格和单元尺寸 采用高精度单元来提高求解精度 分析提高有限元求解精度的途径 tianjinuniversity 三 建立高精度单元的原理和途径 原理 提高单元位移插值函数多项式的阶次 从而提单元拟合局部区域位移 应力变化的能力 途径 增加单元的节点数目 对于平面有限元问题 除三节点三角形单元外 还可以考虑六节点三角形单元和四节点矩形单元 tianjinuniversity 三节点三角形单元有限元分析过程 设位移函数 求位移函数中的未知量 代入函数中 整理可得形函数 性质 几何方程求解应变 几何矩阵 物理方程求解应力 弹性矩阵应力矩阵 运用虚功原理求解 由合成 方法 建立节点荷载列阵 方法 处理位移约束条件 方法 组成 tianjinuniversity 四节点矩形单元有限元分析过程 一 四节点矩形单元位移函数 单元节点编号为k l m n 逆时针 单元节点位移列阵为 设位移函数为 或写为 tianjinuniversity 四节点矩形单元有限元分析过程 二 求解位移函数中的未知系数 将节点坐标 代入函数中 并写成矩阵形式 解上述方程组可得 tianjinuniversity 四节点矩形单元有限元分析过程 三 将所求值代入位移函数中 tianjinuniversity 四节点矩形单元有限元分析过程 四 整理位移函数可得形函数 展开上式可得 四节点矩形单元有限元分析过程 其中 形函数为 tianjinuniversity 四节点矩形单元有限元分析过程 单元位移插值函数可以由单元形状函数与节点位移值的乘积表示 即可以表示为 tianjinuniversity 四节点矩形单元有限元分析过程 由此可见 位移插值函数完全由形函数决定 因此抛开节点位移 只讨论形函数的性质 就可以了解单元的变形性质 例如 四节点矩形单元 若则由可得 因此可以看出 单元变形完全由形函数决定 tianjinuniversity 四节点矩形单元有限元分析过程 另外 可以验证形函数另外两个性质 1 同理对于其余三个形函数 tianjinuniversity 四节点矩形单元有限元分析过程 tianjinuniversity 2 即在单元内任意一点处的形函数之和等于1 四节点矩形单元有限元分析过程 tianjinuniversity 五 几何方程求解应变 将位移插值函数代入几何方程中 形函数矩阵经过微分算子矩阵作用后得到3 8几何矩阵 四节点矩形单元有限元分析过程 tianjinuniversity 六 物理方程求解应力 由平面问题物理方程可得 其中 因此 应力矩阵为 四节点矩形单元有限元分析过程 tianjinuniversity 结论 对于平面四节点矩形单元 其单元上的应力 应变不再是常数 而是在一定程度上呈线性变化 即 方向的正应力和正应变随坐标线性变化 方向的正应力和正应变随坐标线性变化 剪应力沿坐标和坐标均成线性变化 因此 若在弹性体中采用相同数目的节点时 矩形单元的精度要比常应变三角形单元的精度高 四节点矩形单元有限元分析过程 tianjinuniversity 七 运用虚功原理求解 由虚功原理 节点力在节点的虚位移上所做的虚功应等于单元内部应力在虚应变上所做的虚功 即内力虚功 外力虚功 也即 将 和 代入上式 可得 由此 可得 四节点矩形单元有限元分析过程 tianjinuniversity 四节点矩形单元有限元分析过程 tianjinuniversity 引入无量纲坐标 四节点矩形单元有限元分析过程 tianjinuniversity 由几何方程可得单元应变场表达式 可记为 四节点矩形单元有限元分析过程 tianjinuniversity 几何矩阵可表示成分块形式 其中 四节点矩形单元有限元分析过程 tianjinuniversity 由应力与应变关系 可得单元应力场表达式 应力矩阵可表示成分块形式 其中 对于平面应变问题 四节点矩形单元有限元分析过程 tianjinuniversity 其中 即 单元刚度矩阵 四节点矩形单元有限元分析过程 tianjinuniversity 其中 四节点矩形单元有限元分析过程 tianjinuniversity 八 由合成 刚度集成法 首先求出各单元的贡献矩阵 然后将它们叠加形成整体刚度矩阵 但是由于编程时需先将各单元的贡献矩阵储存起来 而贡献矩阵的阶数与整体刚度矩阵阶数相同 因此占用非常大空间 不利于节约空间资源 单元定位数组法 将单元的节点位移编码按照节点顺序排成一行形成一维数组 利用各单元的定位数组 采用 边定位 边累加 的方法 四节点矩形单元有限元分析过程 tianjinuniversity 九 建立节点荷载列阵 节点荷载列阵的组成 其中 为节点荷载 为等效节点荷载 可按照虚功等效原则求解 即将单元内的荷载移置到节点上后 应当与原荷载所作虚功等效 集中力 分布体力 均质等厚单元自重 分布面力 均布侧压 x方向均布荷载 x方向三角形荷载 四节点矩形单元有限元分析过程 tianjinuniversity 十 处理位移约束条件 1 降阶法 若第r个自由度方向位移分量为0 则将整体刚度矩阵第r行 第r列划掉 后一行上移 右一列左移 这样总刚减少一阶 未知数减少一个 例 四节点矩形单元有限元分析过程 tianjinuniversity 2 对角元素置1法 例 已知位移边界条件 可为零 四节点矩形单元有限元分析过程 tianjinuniversity 3 对角元素乘大数法 例 已知位移边界条件 可为零 略去小量有即 总结 tianjinuniversity 四节点矩形单元采用双线性位移插值函数 应力和应变沿坐标轴呈线性变化 因而求解精度比三节点三角形单元高 由于位移模式在单元边界上线性变化 内不能很好地反映应力和应变的变化 并且根据单元公