矩形管挤压过程数值模拟.pdfVIP

矩形管挤压过程数值模拟 1 1 2 1 唐妍 , 高锦张 , 贾俐俐 ，蒋松 1 东南大学材料科学与工程学院，江苏南京（211189 ） 2 南京交通职业技术学院，江苏南京（211188 ） E-mail ：tangyanmm@ 摘 要：本文以Al6063 矩形管为例，分别利用有限元法和有限体积法对矩形管平面分流模挤 压过程进行数值模拟分析，比较了2 种模拟方法针对矩形管挤压过程的模拟结果，得出有限 体积法可以更有效模拟矩形管挤压过程的结论。同时分析矩形管挤压过程的金属流动规律、 应力应变以及温度分布情况，为铝材矩形管分流模的结构设计和参数优化提供了分析依据。 关键词：矩形管；有限体积法；有限元法；数值模拟 中图分类号：TG376.9 1. 引 言 铝合金型材挤压加工过程伴随着很大的塑性变形，既有材料非线性（应力、应变之间的 非线性），又有几何非线性（应变、位移之间的非线性），加上复杂的边界条件，金属流动 规律十分复杂。如仍依靠传统的试模方法来开发挤压模具，将无法适应现在的市场要求，因 此利用数值模拟方法研究铝型材挤压机理和金属流动规律，来为挤压工艺和模具设计提供理 论指导，就显得尤为重要。 目前用于数值模拟塑性成形过程的方法主要有有限元法（FEM ）和有限体积法（FVM ）。 FEM 能很好地处理变形体与模具的接触问题，又具有较高的模拟精度，因此FEM 广泛地应 用于金属塑性成形的数值模拟之中[1~4] 。FVM 最初主要用于流体力学的数值模拟，近年来 被逐渐用于金属塑性成形的数值模拟研究[5~8] 。通过研究发现目前这两种方法大多被用于模 拟针对实心型材的平模挤压过程[9~11]，针对矩形管等空心型材的平面分流模挤压过程的模拟 比较研究还不是很多。因此，利用FVM 和FEM 模拟研究平面分流模挤压过程，对于合理 设计挤压模具结构和提高型材质量具有十分重要的意义。 本文以Al6063 矩形管为例，分别采用FEM 和FVM 对其挤压过程进行模拟，比较了两 种模拟方法针对矩形管挤压过程模拟的可行性，同时分析了矩形管挤压过程的变形规律。 2. FEM 和FVM 基本思想 FEM 基本思想是将材料的整个求解域离散为一系列相互连接的单元体（如三角形、四 边形、四面体或六面体等单元），通过求解各单元和节点在每一加载步上的物理场量，来描 述材料在成形过程中的变形情况。由于有限元采用 Lagrange 方法，每一个单元代表着材料 的一个局部求解域，各节点代表着材料的一个金属质点，所以，在模拟过程中，节点随着金 属的流动而一起运动，单元的形状也随之发生改变。于是，在大变形塑性成形的有限元数值 模拟中，网格经常发生畸变。为了让模拟能够顺利进行下去，必须对畸变后的网格进行重新 划分。 FVM 基本思想是将材料流动所要经历的空间用Euler 网格（如六面体网格）进行离散， Euler 网格的单元和中心格点在空间固定不动，材料只是从一个单元流到另一个单元，并且 材料的质量、动量和能量也随之从一个单元流到另一个单元，此过程必须满足质量守恒、动 量守恒、能量守恒、本构方程、状态方程和热平衡方程等控制方程[12] 。材料的变形行为则 通过均布在材料中的大量质点来表述。由于Euler 网格并不随材料一起运动，它在物体的成 - 1 - 形过程中不会发生畸变，因此 FVM 可以避开网格重划问题。但是在模拟过程中，FVM 网 格需要事先占据材料所要经过的所有空间。对于壁厚较薄的型材截面，就需要将网格划分很 细，因此采用FVM 模拟会占用大量计算机内存。 3. FEM 建模及模拟结果分析 在对铝型材进行模拟挤压之前, 首先应建立几何模型。模拟中所挤型材尺寸是 60 ×