廖敦明《有限差分法基础》第9章 应力场模拟.pptVIP

应力应变模拟 箱体铸钢件 实物图（箭头方向出现热裂缺陷） * 槽板铸件应力场数值模拟 Return * 浇注后30秒时铸件Y向应变分布 * 浇注后360秒时铸件等效应力分布图 * Engineering application instance A shell box of ZL101 (ZAlSi7Mg) * 应力应变模拟 铝合金铸件 * Stress simulating results (a) whole castings (b) cutaway view Stresses distribution of X-axis direction * Equivalent stresses and strain Distribution of equivalent stresses Distribution of equivalent strains * Real casting The crack of the practice casting * 应力应变模拟 槽板铸件 * 应力应变模拟 大型轧辊铸件 不同时刻轧辊长度方向温度场剖切图 * 应力应变模拟 大型轧辊铸件 不同时刻轧辊长度方向应力分布剖切图 * 应力应变模拟 大型轧辊铸件 轧辊实物图及裂纹 * * 谢 谢！ 华铸软件WWW.INTECAST.COM 027华铸软件 — 21世纪铸造企业信息化整体解决方案 InteCAST Software – 21st Century Foundry Solutions 热弹塑性理论及力学模型分析 15个基本待求变量 15个基本方程 塑性力学特点 屈服准则（屈服条件） 基于增量理论的热弹塑性－蠕变本构方程 增量迭加 * 15个基本待求变量 3个位移分量： 6个应变分量： 6个应力分量： * 一点的应力状态分布图 A B D C 微小单元的受力图 15个基本方程 6个几何方程： 3个平衡方程: 6个物理方程: * 简化力学模型 双线性强化弹塑性力学模型 * 屈服准则（屈服条件） Von Mises 屈服准则： * 基于增量理论的弹塑性—蠕变本构方程 其中 —— 弹性应变增量 —— 塑性应变增量 —— 热应变增量 —— 蠕变应变增量 * 增量迭加 * 求解应力场问题的方法 位移法： 取单元结点的位移为基本的未知量，将 各方程中的应力和应变都用位移表示。 这样，首先解出三个位移未知量u、v、w， 再由它们去求应变和应力。 应力法： * 以位移变量表示的平衡方程组的导出 * 最终待求解的线性代数方程组： 与FEM相比，FDM没有必要存储系数矩阵。可以大大节省内存， 避免巨型稀疏矩阵的存储等复杂问题，节省开发工作量。 * 迭代方法的选择 采用 SOR方法， 即逐次超松弛迭代法 超松弛因子（1ω2），一般取 ω =1.6 * 栅形应力框热应力场模拟与残余应力校核 栅形应力框主要尺寸 应力框铸件的FD模型 * 浇注后180秒、1650秒两时刻铸件的温度分布 * 栅形铸件粗、细杆中心点模拟冷却曲线 * 约束条件 * 应力场计算结果 浇注后180秒、1650秒时 沿直杆方向上的应力分量分布图 * 主要部位中间杆方向应力变化曲线 * 应力应变模拟 应力框 * θ形应力框热应力场模拟与残余应力校核 θ形应力框主要尺寸 应力框三维几何实体 * 浇注后5042秒时的色温图 * 浇注后5042秒时沿中间杆方向应力分布 * 浇注后5042秒时Mises等效应力分布 * 主要部位中间杆方向应力变化曲线 * 选取的试样 栅形试样 θ 形试样 * 铸件主要部位实测所得的冷却曲线 * 浇注所得实际铸件 * * 切断的试样照片 5 6 1 2 4 3 θ试件切断形貌及应变片位置 左图为上表面，右图是下表面 栅形试件断口形貌 * 减速箱铸钢件的热应力场模拟 减速箱铸件几何模型 * 浇注后50秒时箱体等效应变分布图 * * 第9章 应力场数值模拟 intecast@163.com 华铸软件 — 21世纪铸造企业信息化整体解决方案 InteCAST Software – 21st Century Foundry Sol