1 引言 板料激光沖擊成形是近年來出現的一種柔性精密無模成形新方法，它能通過高能激光束對板料的逐點逐次沖擊，積小變形成宏觀的大變形。激光沖擊成形是在激光沖擊強化基礎上發展而來的，用于金屬板料塑性成形的一種新技術，它是利用激光誘導的沖擊波壓力作為板料塑性成形的變形力，從而實現金屬板料的宏觀塑性變形。與激光沖擊強化相比，約束的解除為板料的成形提供了空間。圖1為實驗裝置示意圖。金屬板料激光沖擊波成形涉及激光與材料的相互作用、沖擊波傳播及其對材料的加載機制、動態塑性成形理論以及沖壓成形工藝等眾多學科。受到激光器等硬件條件的制約，目前國內能開展激光沖擊成形試驗的研究機構尚且不多，并且試驗費用也是不小的開支，而沖擊成形又是一個極其復雜的過程，存在很多不確定性因素。以上原因大大限制了對于激光沖擊成形技術的研究應用。在這種情形下，采用計算機數值模擬技術來獲得板料變形與激波加載之間的規律不失為一種行之有效的方法。因此本章將借助于有限元分析軟件對金屬板料激光沖擊成形過程進行模擬，為將來建立變形量與激光能量等參數之間的關系以及大面積沖擊成形的理論研究和軟件編制提供依據。 針對本人前期所做的Al-Mg-Sc材料的激光沖擊成形實驗，利用有限元模擬軟件MSC. Marc對此過程進行模擬，以便于將模擬結果和試驗結果進行比較，優化激光參數。 2 金屬板料激光沖擊成形模擬過程 激光沖擊成形模擬前要處理的是:建立模型，將所求問題圖形化:指定單元類型;給出材料數據;對模型進行加載;設置邊界條件;選擇適當的求解器并設置輸出要求;求解;查看計算結果。 2.1 建立模型 激光沖擊成形所用的模型為約束模型，即激光輻照能量轉換體生成等離子體并發生爆炸產生強沖擊波壓力作用于金屬板料表面，使金屬板料產生宏觀塑性變形，這一過程比較復雜，涉及能量轉換體如何吸收激光能量、等離子體怎樣產生和爆炸、沖擊壓力波的傳播以及金屬板料產生宏觀塑性變形等等。由于金屬板料的成形與作用在其上的沖擊波的壓力和作用時間直接相關，其中，沖擊波的壓力是與激光參數、能量吸收層、約束層等因素決定，因此在對金屬板料激光沖擊成形進行模擬時，不再對沖擊波產生和傳播進行模擬，而是根據有關沖擊波理論和已有的研究，將激光誘導的沖擊波簡化為作用在沖擊區域內隨時間變化壓力載荷，然后將其作為載荷直接作用在金屬板料表面。圖2為激光沖擊模型示意圖。 圖2 激光沖擊模型示意圖2.2 網格劃分及單元選取 對于激光沖擊成形過程的有限元模擬來說，其建模過程不是很復雜，直接使用Marc軟件的前處理器MSC. Mentat進行網格劃分。由于在激光沖擊成形過程中所用的金屬板料為薄板，其厚度遠小于另外兩維尺寸，因此在模擬過程中考慮采用殼體單元。另外，考慮到激光沖擊成形本身的特點，非線性(包括幾何非線性和材料非線性)、大變形、而且材料的應變率非常高，而且殼體單元的節點數應該與網格相對應，對這些因素綜合考慮后，在選擇單元類型時選用S4R。單元大小為0.05mm*0.05mm，單元類型選取140號殼單元，它是一種四節點厚殼單元，以整體位移和旋轉為自由度，坐標值、平移和轉動都采用雙線性插值法。這種殼幾何形狀采用MITC4殼的幾何，采用了和物理穩定方法結合的ANS(假設自然應變)法，構建的單元就是縮減積分單元，這種平面單元沒有對翹曲進行任何人工修整。這些方法提高了一個積分點單元的精度而沒有影響計算速度，可進行大的非線性分析。 2.3 邊界條件定義 激光沖擊成形過程中邊界條件的處理對整個分析過程來說是比較重要的一環。這在多點，多次沖擊的時候尤為重要。實際在激光沖擊成形過程中，邊緣板料沿板料面內的自由度沒有限制，壓邊圈的存在只是限制了板厚方向的自由度，由于單次激光沖擊下板料的變形量在mm級，經過試驗研究發現板料變形后的尺寸在徑向方向很小，這樣在模擬的時候就可以對邊界條件進行簡化處理，既可以將壓邊圈壓住的部分自由度全部限制住.這對最后的變形結果影響不大。板料固定在凹模基座上，板料在其平面內的轉動、移動都沒有受到限制，因此在進行模擬過程中對邊界條件進行處理時限制除板料平面內的其他自由度如圖3所示為激光沖擊成形模型力與位移的邊界條件定義。 圖3 激光沖擊邊界條件定義2.4載荷工況定義 在建立金屬板料激光沖擊成形模型時，激光被直接簡化成隨時間變化的壓力載荷。在激光沖擊成形過程中主要是對沖擊載荷和加載時間定義。但是在模擬過程中，確定壓力載荷值大? ⒃睪勺饔檬奔湟約凹釉胤絞蕉閱Ｄ飩峁己苤匾? 2.4.1 激光沖擊峰值載荷及加載方式選取 激光沖擊峰值載荷通過在線測量最為準確。通過壓力測量方法得到峰值壓力，它是利用VISAR系統得到的，其主要原理基于氫氣探針探測受壓箔片背面產生的多普勒位移。在L.U.L.I實驗室的0. 6ns的測量壓力是利用PVDF壓電傳感器和EMV系統得到的，和PVDF測量不同的是VISAR和EMV測量包括自由速度測量。根據自由速度值，利用經典的蘭金雨貢紐關系(Rankine Hugoniot)(公式(1)和公式(2))可得到等離子體的壓力水平：峰值壓力的計算公式 P= 式中各個系數參都是已知的，可以估算出峰值壓力的大小，然后根據激光脈沖的作用時間和實驗記錄的脈沖信號，將誘導的沖擊波按照激光脈沖信號的分布，采取分段逐次逼近的辦法來實現激光脈沖的加載。這里需要注意的是沖擊波在整個作用時間內，是隨著時間的變化而變化的，先是上升，隨后衰減，呈準高斯分布形式。可以借助Fortran語言，來編寫得到的載荷曲線如圖5所示，激光脈寬20ns.激光沖擊波作用時間為60ns。2.4.2激光沖擊波時空分布 從現已發表的關于激光沖擊成形有限元仿真的文章來看，都未考慮到激光沖擊波的空間分布，很顯然這種簡化導致的誤差是不可忽視的。美國哥倫比亞大學的Y. Lawrence Yao在R. Fabbro模型的基礎上給出了激光沖擊強化中壓力沿光斑徑向的變化公式(3)，如下式中，R2l是激光光斑的半徑，r為離激光束中心的距離。 從式(3)可以看出，沖擊載荷是關于距離沖擊中心半徑的函數，考慮到激光脈沖呈高斯分布，將整個載荷作用區域(即光斑直徑內的區域)離散為多個帶狀區域，各區域之間的載荷按高斯分布加載，而帶狀區域之內將載荷視為均布載荷。 顯然，帶寬越小越接近于實際的加載狀況.當帶寬區域趨于0時，在整個光斑直徑內載荷呈高斯分布，與實際情況完全吻合。綜合考慮計算精度、計算效率和加載難易程度，將半徑為3. 5mm的加載區域以0. 5mm的間隔離散為7個環形帶，計算每個帶內的平均載荷。計算過程見附錄。圖6顯示了按高斯分布離散后的加載的帶狀云圖，可以看出，距激光沖擊中心越近，載荷越大，呈高斯分布。 在模擬過程中，還要考慮的是激光誘導的沖擊波壓力的作用時間，在已有的理論和實驗研究中，激光沖擊誘導的沖擊波壓力的作用時間一直是討論的主題。這是因為激光脈沖作用的時間很短，只有幾十納秒，而且在整個激光脈沖的作用過程中，上升沿時間和下降沿時間又不同，激光脈沖誘導產生的沖擊波究竟在板料表面上作用多少時間，有許多學者進行了研究。R.Fabbro等人的研究認為，激光誘導的沖擊波的作用時間大約為激光脈寬的2-3倍由于脈沖寬度為22ns，洪聽在其研究中也發現沖擊波的時間寬度是激光脈寬的3倍左右甚至更高。因此仿真時加載時間設置為70ns。在Marc中，可以用一張表(Table)來控制加載歷程。 2.5 材料參數定義 利用有限元分析軟件進行材料的成形仿真時，本構關系的建立是至關重要的一步。模擬中材料參數的選擇必須能夠反映真實情況，在金屬板料激光沖擊成形過程中，金屬板料發生塑性變形，而且應變率極高，同時材料還會出現硬化效應。因此在進行模擬時，材料參數的設定既要能夠反映板料的變形特點，又要考慮到高應變率下材料的動態屈服應力和硬化效應對變形的影響。 在用Marc做彈塑性分析時.彈性變形階段由彈性模量控制，所以本構關系中只考慮塑性變形的部分。粘塑性本構方程在一定范圍內較好的表達了材料的粘—塑性本構特性，但是對于實際問題，最好通過材料的動態性能試驗得到具體的粘一塑性本構方程，這樣使計算具有較高的精度。但由干激光沖擊成形的應變速率極高(達106~107s-1)，目前見到所做的試驗都在l04s-1量級以下，遠達不到要求的應變速率值。受試驗條件的限制，通過插值法，在有限元仿真時，本章使用應變率103~104s-1的本構模型，一樣能得出在激光沖擊成形高應變率下的數值模擬結果趨勢。Marc中自帶Johson-Cook模型，而且模型中各個參數我們已經通過靜態下拉伸實驗以及動態壓縮實驗求得，它們分別為A=314. 6MPa B=223, C= 0. 0051, n= 0.1220,假設材料為各項同性。 2.6 接觸 考慮到激光沖擊的實際情況，板料放在凹模上，再用壓板將板料壓緊。可以用定義接觸體的方法來很好地模擬實際沖擊時的約束情況，并且建立了接觸表來控制接觸體之間的接觸類型。如圖7所示，定義了三個接觸體:sheet為板料，定義為變形體;holddown為壓板，定義為剛體，和板料之間的接觸關系為touching;mold為凹模，定義為剛體，和板料之間的接觸關系為glue(粘合)。3激光沖擊成形數值模擬結果 3.1 激光沖擊成形過程 利用MSC. Marc的后處理器MSC. Mentat進行后處理，將分析生成的結果以圖形、曲線和動畫的形式展現出來，便于和實驗結果進行比較。以下都是單點單次激光沖擊成形數值模擬結果，圖8以帶狀云圖(CONTOUR BAND)形式顯示了激光沖擊板料變形過程，厚度為1. 2mm的AI-Mg-Sc合金板在激光沖擊載荷下變形量逐漸的變化。在Marc中變形過程設置成一百步，這里選取了四步來說明其過程的變化。圖9顯示了激光沖擊板料的最終變形量為1.50lmm。Y向位移沿路徑的變化曲線如圖10所示。3. 2激光能量與變形量的關系 為研究激光能量與變形量的關系，我們選取一組激光能量數據來模擬變形量隨激光能量變化的過程。圖11為激光能量分別為40. 030J, 36. 023J,30. 690J,25. 375J的激光對厚度為1.2mm的A1-Mg-Sc進行沖擊后的變形模擬圖，可以看出隨著激光能量的增大，板料變形量增大。表1為模擬結果和實驗結果的對比，從表中數據可以看出模擬結果均比實驗數據略有差距。這是因為在模擬中對實際情況進行了簡化處理.其中涉及到上面所述的激光沖擊波的作用時間和載荷大小以及加載方式等因素因此出現誤差，但從變化趨勢來看和實驗結果基本符合。但是從模擬的結果來看，都大于實驗結果，這可能是因為激光能量輸出有一定的誤差，實際中輸出的激光能量大于所顯示的讀數。圖12為不同激光能量下板料最大變形量模擬值與實驗值比較，當激光能量大于30)時，變形量增大幅度減小，這和實驗結果吻合。 圖12 不同激光能量下板料最大變形量模擬值與實驗值比較3. 3約束孔徑與變形量的關系 分別模擬約束孔徑分別為20, 25, 30mm激光沖擊板料變形，從模擬的結果可以得出:當激光能量大于30J時，板料最大變形量隨著約束孔徑的增大而增大，當激光能量小于30J時.板料最大變形量隨著約束孔徑的減小而增大。圖13為激光能量大于30J時，不同約束孔徑下板料最大變形量模擬值與實驗值比較，表明模擬結果與實驗結果變化趨勢是吻合的。表2為激光能量大于30J，約束孔徑對板料變形量的影響。 圖13 不同約束孔徑下板料最大變形量模擬值與實驗值比較3. 4等效塑性應力應變分布 圖14和圖15分別顯示了變形區域的等效MIS-ES應力和等效塑性應變。從等效MISES應力和呈“ω”狀的等效塑性應變曲線上可以看出，在距沖擊中心9mm(即凹模邊緣)附近和沖擊中心處，應力和應變值都達到了極值點，這些區域屬于易破裂區域。另外等效MISES應力還存在兩個峰值點，在趾沖擊中心3mm附近(即激光光斑邊緣)，應力值也比較大。 在其他參數相同，只是激光能量不同的情況下，等效應力隨著激光能量的增加而增加，但由于采用圓形光斑，在光斑周長上產生的應力波向中心傳播并在光斑中心聚合，產生一個巨大的脈沖，這個脈沖將會消除光斑中心點附近的殘余應力，在光斑中心產生應力空洞現象，所以在光斑的中心等效應力反而小(圖14)。4結論 本文對激光沖擊成形進行了計算機模擬分析，主要在這幾個方面進行了研究: 在對激光沖擊成形進行理論分析的基礎上，建立了有限元分析模型，并對其中的幾個關鍵步驟，如模型的建立、單元類型選? ⑼窕趾筒牧喜問難≡瘛⒊寤鞴討脅牧系撓不約凹す獬寤韃ǖ拇砑捌浼釉鼗坪褪笨輾植嫉確矯娼辛私檣堋? 對板料激光單點沖擊成形進行了模擬，并將模擬結果與實驗結果進行對比.通過對不同條件下的實驗結果進行數值模擬分析和修正，為預測板料的變形程度，優化加工工藝參數，為編制激光沖擊成形中板料變形的控制軟件提供理論依據。 分析了板料內等效應力和殘余應力的分布情況，在距沖擊中心9mm(即凹模邊緣)附近和沖擊中心處，屬于易破裂區域，距沖擊中心3mm附近(即激光光斑邊緣)，應力值也比較大。