導流模又稱前室模，是生產薄壁復雜實心型材常用的擠壓模具。采用導流模具生產型材時，坯料首先通過導流孔預成形，獲得與型材相似的幾何形狀，然后再進行二次變形，擠壓出各種斷面形狀的型材。導流孔不僅增大了坯料與型材的幾何相似性，有效地控制了金屬流動，減少了產品的扭擰和彎曲變形，而且改善了模具的受力條件，提高了模具壽命，所以其設計是否合理是導流模具設計的關鍵。傳統模具設計主要依賴設計者的經驗，經過反復試模修模，最后生產出合格零件，既費時又費力。利用數值仿真的方法不僅周期短、成本低，而且可以獲得材料的應力、應變、溫度、速度等現場難以測量的場量，對于掌握擠壓過程中材料在模具中的流動規律進而提高型材質量非常重要。為此，國內外許多學者采用數值模擬的方法對擠壓成形過程進行了大量的研究，并取得了較大的進展。

目前關于擠壓過程的數值模擬研究大多采用拉格朗日法和歐拉法。采用拉格朗日法分析大變形的擠壓問題時，單元容易產生畸變，需要頻繁的網格重劃，導致體積損失過大，嚴重影響了計算精度；而采用歐拉法模擬時，需要利用復雜的數學映射來描述自由表面的運動狀況并且對材料所有可能流過的區域劃分網格，占用大量計算機內存，特別是分析薄壁、空心、復雜斷面型材時計算時間令人難以接受。ALE(Arbitrary Lagrangian Eulerian, ALE) 算法吸收了拉格朗日和歐拉算法的優點，通過引入一個獨立于物質構型和空間構型的參考構型，使得計算網格能夠任意移動，從而避免了單元畸變以及自由界面追蹤所帶來的困難。本文采用基于ALE算法的商業化軟件HyperXtrude對一異形鋁型材的擠壓過程進行數值模擬，以導流孔形狀為影響因素，采用Box-Behnken設計確定試驗方案建立分析模型，利用最小二乘法對模擬結果進行擬合，得到型材截面速度均方差的響應曲面公式，并以此為優化目標采用粒子群算法對導流孔形狀進行優化。

2 優化模型的相關理論

2.1 ALE算法基本理論

ALE算法以物質運動的ALE模型為理論基礎，不同于拉格朗日或歐拉描述模型，ALE模型另外引進了一個獨立于物質構型(ΩX)和空間構型(Ωx)的參考構型(Ωξ)。計算網格的劃分是在參考構型中進行的，獨立于物體和空間運動的，可以根據需要自由選擇。根據任意單元體中的質量、動量和能量守恒定律，獲得參考坐標系下的控制方程：
    質量守恒方程：

 

    動量守恒方程：

 

    能量守恒方程： 

  

式中ρ為材料的密度，wi、wj為物質點在ξ空間中的位置矢量對時間的導數，vi為物質點在空間中的速度矢量，Tji為定義在參考構型下的第一類皮奧拉-克希荷夫應力張量，fi為作用于物體中單位質量的體力，e為物體單位質量中的內能。

在采用ALE算法求解具體問題時，需要引入正確的材料本構關系并對模型施加適當的邊界條件，然后對整個控制方程進行求解。

2.2 響應曲面模型和Box-Behnken試驗設計

響應曲面法(response surface methodology, RSM)，也稱回歸設計，是采用多元二次回歸方程來擬合因素和響應值之間的函數關系，通過對回歸方程的分析來尋求最優工藝參數，解決多變量問題的一種統計方法。RSM法與其他數據統計方法相比，不僅考慮了自變量之間的交互作用，提高了擬合精度，而且還可以運用圖形技術將二者之間的函數關系顯示出來，使結果更加直觀。本文選擇二階響應曲面方程，其模型可表示為：

  

式中，xi為設計變量，ε為殘余誤差，βo、βi、βii、βpi均為待定系數。

Box-Behnken試驗設計是一種基于三水平的二階試驗設計方法，可以評價指標和因素之間的非線性關系，是RSM常用的試驗設計方法之一。Box-Behnken設計的每個因素只需要三個水平，與其他方法(如星點設計法)相比所需的試驗次數較少，效率更高，且所有的影響因素不會同時處于高水平，所有的試驗點都落在安全區域內，因此本文選取Box-Behnken試驗設計方法對導流孔的形狀進行優化設計。

2.3 PSO算法基本原理

粒子群優化算法（particle swarm optimization, PSO）是Kennedy和Eberhart于1995年提出的一種以鳥類覓食行為為基礎的全局優化算法，因其原理簡單，算法實現容易，運行速度快，已被廣泛應用于機械、化工、生物醫學等領域。該算法模仿鳥類的群體覓食過程，將待優化問題的解看作為搜索空間中的一個粒子（鳥），解的優劣程度由適應函數決定。各粒子通過不斷追蹤自身所經過的最佳位置和整個種群所經過的最佳位置這兩個極值來更新自己的速度和位置，最終到達全局最優解所在的位置。

3 模擬試驗方案的確定

為了獲得最優的導流孔形狀，進而更好地控制材料在模具中的流動，提高型材質量，在保持其他工藝參數不變的情況下，選取導流孔不同部位的寬度作為優化變量，如圖1所示，并結合實際情況確定各個變量的范圍如表1所示。以模具出口型材截面上的金屬流速均方差(SDV)為優化目標，建立優化模型如下：

式中，n為所考慮截面內的節點個數；vi為出口截面上第i個節點的流速；v為出口截面上所有節點的平均流速。

表1 設計變量及其取值范圍

  

 
圖1 設計變量的位置

選取6063鋁合金作為模擬材料，擠壓筒預熱溫度為450℃，坯料加熱溫度為480℃，擠壓比為31.2，模具與坯料間的對流換熱系數為3000 W/m2·℃，擠壓速度為2mm/s，模具工作帶處設為庫侖摩擦，摩擦系數為0.3，其他部位設置為粘性摩擦。將材料流經區域分為棒料、導流孔、工作帶和型材四個部分進行網格劃分，并建立分析模型如圖2所示，其中工作帶和型材部分采用三棱柱網格，導流孔和棒料部分采用四面體網格。將上述4個設計變量在各自范圍內取3個水平，采用Box-Behnken設計確定試驗方案，對所設計29組試驗方案進行數值模擬，并計算出相應的型材出口流速均方差，如表2所示（由于數據較多，在此只列出部分結果）。

 

圖2 分析模型
表2 Box-Behnken試驗設計表及計算結果

   

 3.2 型材壁厚