} { } { ](
[ } { ]
[K ∆u =− Kfr ∆u − ∆ut （3.46）
其中{ u∆ 與} {∆ut}分別為為節點位移向量與模具位移向量，將（3.46）
式重新整理可得有限元素方程式之最後形式如下所示：
} ]{
[ } { ) ] [ ] [
( K + K fr ∆u = Kfr ∆ut （3.47）
而[K fr]摩擦剛性矩陣(其矩陣大小為3× ）僅存在於與模具接觸時之3 節點自由度。
3.6 廣義 rmin法之增量步驟的計算
每一次時間增量的過程中，板金成形是以其增量開始之座標為參 考基準，其採用 ULF 理論並結合 Yamada 等所提之廣義之 rmin方法來 決定每一次計算步驟的增量值的大小，下列六種情況下 r 值中之最小 者定為 rmin。
(1) 元素最大容許應變增量 (2) 元素最大容許旋轉增量 (3) 元素最大容許等效應力增量
(4) 材料降伏之判定 (5) 自由節點與模具接觸 (6) 接觸節點與模具分離
應用rmin方法的程序及增量大小的計算包含下列步驟:
1. 在增量開始前，先給予一個假設模具位移增量∆a；
2. 相對於前述模具的位移增量∆a，利用剛性方程式解出其虛擬 解∆u；
3. 經由前述(1)~(6)的r值，找出最小的值即為rmin的值；
4. 將所得的解 u∆ 乘上rmin係數(0＜rmin＜1)作為其加權，故此增 量的大小將變為rmin∗ u∆ ，此將用來改善成形的變形狀態、總 位移、總應力及每一個元素降伏極限的數值。
由於本研究採增量理論，每一步驟與材料之幾何形狀及元素所處 狀態有關，故以rmin決定該步驟之實際增量。在第一次的虛擬增量計 算後，經由上述各種增量限制之判定，每一步驟僅允許一元素或節點 改變其原有狀態，並選擇最小之rmin去決定該步驟之實際增量。
圖 3-1 退化殼元素之座標
圖 3-2 殼元素之自由度
圖 3-3 接觸彈簧元素
四、結果與討論
本計畫為整合型計畫之子計畫一，主要以開發金屬板材成形三維 解析 CAE 軟體有關前處理(pre- processor)的部份為主。其步驟如下：
1. 首先使用 I-DEAS 套裝軟體對料片進行元素分割、節點編號，並 將其節點座標值、元素號碼等數據檔輸出並轉換成金屬板材成形 三維解析 CAE 軟體能接受之輸入格式。
2. 建立工具(包括衝頭、模具與壓料板)之有限元素三維模式，並轉換 成金屬板材成形三維解析 CAE 軟體能接受之輸入格式。
3. 撰寫 CAE 軟體有關邊界條件與材料機械性質輸入部份之程式，以 建立完整的 CAE 軟體有關之前處理部份，並與其他子計畫整合成 完整的金屬板材成形三維解析 CAE 軟體，並將其應用於金屬板材 突緣引伸成形之分析。
今年為本程式完成階段，本子計畫負責將程式系統化，並配合子 計畫二之平行運算，已使本程式具有「操作簡單」、「執行速度快」與
「精準度高」等優點。將本計畫所推導之前處理(pre- processor)的部份 程式化，並與其他子計畫之前後處理程式進行整合後，以金屬板材帽 型引伸加工之解析為例進行分析。在模具方面有限元素網格分割是採 用三角形元素，其目的是解決模具之幾何形狀問題，以避免網格分割 時在模具表面同時產生三角形和四邊形元素，增加分析程式的複雜程
度。而工件之有限元素網格分割則是採用四邊形元素。上述模具與工 件之幾何圖形皆由套裝軟體 I-DEAS 建立，並透過 Simulation 模組產 生有限元素網格資料，再轉出 Universal 檔案，以作為本整合計畫程 式之前處理輸入資料，經分析後再經由 I-DEAS 軟體作後處理。以下 為金屬板材帽型引伸加工為例，進行分析：