基于有限体积法的铝型材挤压过程数值模拟.docxVIP

基于有限体积法的铝型材挤压过程数值模拟 铸造模拟网格划分 由于其质量轻、强度高、加工和应用性能好、焊接、耐腐蚀性好、造型美观、断面多样性等优点，广泛应用于民用装饰、包装、汽车制造、航空航天等领域。国内外学者采用理论解析、物理和数值模拟等方法对铝型材挤压成形过程进行了研究。很多学者采用有限元法对各种铝型材挤压过程进行了不同层次的模拟。国内外也出现了不少商业化的金属体积成形有限元仿真软件, 这些商品化软件在铝型材挤压成形理论研究方面发挥了很大作用。但是相对于一般的体积成形工艺, 铝型材挤压属于三维流动大变形问题, 其挤压比通常为50～90;金属流动剧烈, 规律复杂。用有限元法模拟铝型材挤压, 不可避免地会碰到网格再划分问题, 而且坯料尺寸和型材的厚度相差太大, 这给网格再划分带来极大的困难, 即使是商品化程度很高的有限元分析软件, 在金属进入工作带后, 也常常出现网格不能生成和计算精度严重下降的现象, 从而导致模拟过程中断。因此众多的铝型材有限元分析研究, 都只局限在金属刚流入工作带的这一阶段。 由于铝型材挤压产品的质量主要取决于金属在模具型腔内的流动是否平衡, 制品前沿的金属流速是否均匀, 即在铝型材挤压过程中, 更关注的是制品断面单位时间内的金属流动量, 因此, 完全可以采用流体力学中的Euler方法来描述铝型材挤压过程中金属的流动, 从而避免有限元方法和Lagrange难以处理且又无法回避的三维网格重划分等问题。基于有限体积法, 日本锻造协会在MSC/DYTRAN的基础上开发了商品化分析软件MSC/Superforge。国内的学者利用Superforge对铝型材挤压成形进行了模拟。但是Superforge并不是专门针对铝型材挤压模拟的软件, 并且从很多学者的应用经验来看, 用Superforge模拟大变形的铝型材挤压存在一些问题, 如:网格过粗会造成模拟效果粗糙, 增加网格数可解决此问题, 但网格数目增加又会造成计算量和运算时间呈三次方增长, 对计算机内存和CPU要求过高, 其实用性受到严重限制。因此, 研究专门针对铝型材挤压过程的有限体积法数值模拟关键技术和开发相应的数值模拟软件, 对于铝型材挤压工艺和模具开发, 以及提高挤压产品质量具有十分重要的意义。 1 euter网格简介 有限体积法基本思想是将材料流动的空间用Euler网格 (如六面体网格) 进行离散, Euler网格是一个固定在空间的参考框架, Euler网格的单元和中心格点在空间固定不动, 材料只是从一个单元流到另一个单元, 材料的质量、动量和能量也随之从一个单元流到另一个单元, 在此过程中, 必须满足质量守恒、动量守恒、能量守恒、状态方程等控制方程。 1 质量常数守固定方程 ?ρ?t+?(ρvi)?xi=0(1)?ρ?t+?(ρvi)?xi=0(1) 2 pij-sij自由基vivvivvivvivvivvivvivvivvivvivvivviv2.2.2.2 ?(ρvi)?t+?(ρvivj+pδij-Sij)?xj=0(2)?(ρvi)?t+?(ρvivj+pδij?Sij)?xj=0(2) 3 [viviij]t+e0.[vi.5.3.3.3.3.3.3.3.3.3.34.5.2.35.5.2.35.2.35.2.34.5.2.35.2.35.2.35.2.35.2.35.2.35.3.4.3.4.3.4.3.4.3.4.3.4.3.4.3.3.3.3.3.3.3.3.3.3.3.3.3.4.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3 ?(ρE)?t+?(ρviE)?xi=?[vi(Sij-pδij)]?xj(3)?(ρE)?t+?(ρviE)?xi=?[vi(Sij?pδij)]?xj(3) 4 t+vi规 ?(ρΦ)?t+?(ρΦvi)?xi=0(4)?(ρΦ)?t+?(ρΦvi)?xi=0(4) 式中ρ——材料密度 v——速度矢量 p——静水压力 δij——克氏符号 Sij——应力张量 E——内能 Φ——任意状态变量,