Altair 2009 HyperWorks 技术大会论文集 基于响应曲面法和粒子群算法的铝型材挤压模具优化设计* 基于响应曲面法和粒子群算法的铝型材挤压模具优化设计* 陈浩 张存生 马新武 赵国群 陈浩 张存生 马新武 赵国群 山东大学模具工程技术研究中心 山东大学模具工程技术研究中心 - 1 -Altair 2009 HyperWorks 技术大会论文集 基于响应曲面法和粒子群算法的铝型材挤压模具 基于响应曲面法和粒子群算法的铝型材挤压模具 优化设计 优化设计 Optimal Design of Aluminum Extrusion Die Based on Response Surface Method and Particle Swarm Optimization 陈浩 张存生 马新武 赵国群 山东大学模具工程技术研究中心，济南，250061 摘 要摘 要：导流孔形状是导流模具设计的主要参数，对挤压产品的质量、模具寿命和挤压力有着非常重要的影响。 本文以模具出口型材截面的速度均方差为目标， 以导流孔的形状为设计变量，采用基于 ALE 算法的 HyperXtrude 软件，结合响应曲面法和粒子群算法建立了导流孔形状优化模型并对其进行非线性优化。 与优化前的结果相比， 优化后型材截面上速度分布更加均匀，型材质量得到了大幅度提高。 关键词关键词：铝型材，ALE 算法，响应曲面法，粒子群算法，HyperXtrude Abstract：Diversion hole shape is the main parameters of diversion die design. The quality of extrusion products, die life and extrusion pressure has a very significant impact. In this paper, the speed of mold export profile cross-section mean square deviation as the goal, with diversion hole shape as design variables, based ALE algorithm HyperXtrude software. Combined with response surface method and particle swarm optimization has established diversion hole shape optimization models and of its non-linear optimization. Compared with the results before and optimization, the optimized velocity profile more uniform cross-section, profiles quality has been greatly improved. Keyword: aluminum profile，ALE algorithm，response surface method，Particle Swarm Optimization，HyperXtrude 1 引言 1 引言 导流模又称前室模， 是生产薄壁复杂实心型材常用的挤压模具。 采用导流模具生产型材时，坯料首先通过导流孔预成形，获得与型材相似的几何形状，然后再进行二次变形，挤压出各种断面形状的型材。 导流孔不仅增大了坯料与型材的几何相似性， 有效地控制了金属流动，减少了产品的扭拧和弯曲变形，而且改善了模具的受力条件，提高了模具寿命，所以其设计是否合理是导流模具设计的关键【1】。传统模具设计主要依赖设计者的经验，经过反复试模修模， 最后生产出合格零件， 既费时又费力。 利用数值仿真的方法不仅周期短、 成本低，而且可以获得材料的应力、应变、温度、速度等现场难以测量的场量，对于掌握挤压过程中- 2 -Altair 2009 HyperWorks 技术大会论文集 材料在模具中的流动规律进而提高型材质量非常重要。 为此， 国内外许多学者采用数值模拟的方法对挤压成形过程进行了大量的研究2-6，并取得了较大的进展。 目前关于挤压过程的数值模拟研究大多采用拉格朗日法和欧拉法。 采用拉格朗日法分析大变形的挤压问题时，单元容易产生畸变，需要频繁的网格重划，导致体积损失过大，严重影响了计算精度； 而采用欧拉法模拟时， 需要利用复杂的数学映射来描述自由表面的运动状况并且对材料所有可能流过的区域划分网格， 占用大量计算机内存， 特别是分析薄壁、 空心、复杂断面型材时计算时间令人难以接受。ALE(Arbitrary Lagrangian Eulerian, ALE) 算法吸收了拉格朗日和欧拉算法的优点， 通过引入一个独立于物质构型和空间构型的参考构型， 使得计算网格能够任意移动，从而避免了单元畸变以及自由界面追踪所带来的困难7。本文采用基于ALE算法的商业化软件HyperXtrude对一异形铝型材的挤压过程进行数值模拟，以导流孔形状为影响因素，采用Box-Behnken设计确定试验方案建立分析模型，利用最小二乘法对模拟结果进行拟合， 得到型材截面速度均方差的响应曲面公式， 并以此为优化目标采用粒子群算法对导流孔形状进行优化。 2 优化模型的相关理论 2 优化模型的相关理论 2.1 ALE 算法基本理论算法基本理论 ALE 算法以物质运动的 ALE 模型为理论基础，不同于拉格朗日或欧拉描述模型，ALE模型另外引进了一个独立于物质构型(X)和空间构型(x)的参考构型()。 计算网格的划分是在参考构型中进行的，独立于物体和空间运动的，可以根据需要自由选择。根据任意单元体中的质量、动量和能量守恒定律，获得参考坐标系下的控制方程： 质量守恒方程： 0iiw t +=动量守恒方程： jiii ji jjTvvwft+=+能量守恒方程： ii iji ijvqeewTti+=式中为材料的密度，、为物质点在iwjw空间中的位置矢量对时间的导数，为物质点在空间中的速度矢量，为定义在参考构型下的第一类皮奥拉-克希荷夫应力张量，ivjiTif(3) (2) (1) - 3 -Altair 2009 HyperWorks 技术大会论文集 为作用于物体中单位质量的体力，为物体单位质量中的内能。 e在采用 ALE 算法求解具体问题时，需要引入正确的材料本构关系并对模型施加适当的边界条件，然后对整个控制方程进行求解8。 2.2 响应曲面模型和响应曲面模型和 Box-Behnken 试验设计试验设计 响应曲面法(response surface methodology, RSM)，也称回归设计，是采用多元二次回归方程来拟合因素和响应值之间的函数关系， 通过对回归方程的分析来寻求最优工艺参数， 解决多变量问题的一种统计方法。RSM 法与其他数据统计方法相比，不仅考虑了自变量之间的交互作用，提高了拟合精度，而且还可以运用图形技术将二者之间的函数关系显示出来，使结果更加直观。本文选择二阶响应曲面方程，其模型可表示为： 2 0 11nniiiiipipi iip iyxxx x=+ (4) 式中，为设计变量，ix为残余误差，0、i、ii、pi均为待定系数。 Box-Behnken 试验设计是一种基于三水平的二阶试验设计方法， 可以评价指标和因素之间的非线性关系，是 RSM 常用的试验设计方法之一。Box-Behnken 设计的每个因素只需要三个水平，与其他方法(如星点设计法)相比所需的试验次数较少，效率更高，且所有的影响因素不会同时处于高水平，所有的试验点都落在安全区域内，因此本文选取 Box-Behnken试验设计方法对导流孔的形状进行优化设计9。 2.3 PSO 算法基本原理算法基本原理 粒子群优化算法（particle swarm optimization, PSO）是 Kennedy 和 Eberhart 于 1995 年提出的一种以鸟类觅食行为为基础的全局优化算法，因其原理简单，算法实现容易，运行速度快，已被广泛应用于机械、化工、生物医学等领域。该算法模仿鸟类的群体觅食过程，将待优化问题的解看作为搜索空间中的一个粒子（鸟） ，解的优劣程度由适应函数决定。各粒子通过不断追踪自身所经过的最佳位置和整个种群所经过的最佳位置这两个极值来更新自己的速度和位置，最终到达全局最优解所在的位置10。 3 模拟试验方案的确定 3 模拟试验方案的确定 为了获得最优的导流孔形状，进而更好地控制材料在模具中的流动，提高型材质量，在保持其他工艺参数不变的情况下，选取导流孔不同部位的宽度作为优化变量，如图 1 所示，并结合实际情况确定各个变量的范围如表 1 所示。 以模具出口型材截面上的金属流速均方差(SDV)为优化目标，建立优化模型如下： 21()ni ivv SDVn= =(5) 式中，为所考虑截面内的节点个数；为出口截面上第 个节点的流速；niviv为出口截- 4 -Altair 2009 HyperWorks 技术大会论文集 面上所有节点的平均流速。 图 1 设计变量的位置设计变量 取值范围 导流室宽度 a(mm) 916 导流室宽度 b(mm) 612 导流室宽度 c(mm) 69 导流室宽度 d(mm) 916 表 1 设计变量及其取值范围 选取 6063 铝合金作为模拟材料，挤压筒预热温度为 450，坯料加热温度为 480，挤压比为 31.2，模具与坯料间的对流换热系数为 3000 W/m2，挤压速度为 2mm/s，模具工作带处设为库仑摩擦，摩擦系数为 0.3，其他部位设置为粘性摩擦。将材料流经区域分为棒料、导流孔、工作带和型材四个部分进行网格划分，并建立分析模型如图 2 所示，其中工作带和型材部分采用三棱柱网格， 导流孔和棒料部分采用四面体网格。 将上述 4 个设计变量在各自范围内取 3 个水平，采用 Box-Behnken 设计确定试验方案，对所设计 29 组试验方案进行数值模拟，并计算出相应的型材出口流速均方差，如表 2 所示（由于数据较多，在此只列出部分结果） 。 - 5 -Altair 2009 HyperWorks 技术大会论文集 工作带 和型材 导流室 棒 料 图 2 分析模型 表 2 Box-Behnken 试验设计表及计算结果 设计变量 序号 a(mm) b(mm) c(mm) d(mm) 型材截面速度均方差 1 12.5 6 7.5 16 5.411 2 12.5 9 7.5 12.5 2.568 3 12.5 6 6 12.5 2.819 4 12.5 12 7.5 9 5.634 5 12.5 9 6 9 3.572 25 16 9 7.5 9 4.603 26 12.5 12 9 12.5 5.346 27 16 9 9 12.5 3.451 28 16 9 7.5 16 2.114 29 12.5 9 7.5 12.5 2.568 4. 响应曲面的建立和粒子群算法寻优4. 响应曲面的建立和粒子群算法寻优 根据表 2 中的数据， 利用最小二乘法拟合响应曲面， 得到型材出口速度均方差与导流孔形状的响应面函数如下： 222243.151 1.9832.1942.2551.6870.1020.1120.1240.0580.0880.1760.1190.1360.1950.107SDVRabcdabacadbcbdcdabcd=+(6) 表 3 为本文模型的方差分析结果，模型较小的 P 值表明回归模型对响应值的影响极为- 6 -Altair 2009 H