非圓形工件落料模具設(shè)計(jì)(復(fù)合模-UG）【說明書+CAD+UG】,說明書+CAD+UG,非圓形工件落料模具設(shè)計(jì)(復(fù)合模-UG）【說明書+CAD+UG】,圓形,工件,模具設(shè)計(jì),復(fù)合,ug,說明書,仿單,cad 畢業(yè)設(shè)計(jì) 文獻(xiàn)翻譯 院（系）名稱 工學(xué)院機(jī)械系 專業(yè)名稱 機(jī)械設(shè)計(jì)制造及其自動化 學(xué)生姓名 周凱 指導(dǎo)教師 薛東彬 2012年 03 月 10 日 黃河科技學(xué)院畢業(yè)設(shè)計(jì)（文獻(xiàn)翻譯） 第24頁 沖壓半徑和角度對管向外卷曲過程的影響 摘要：跟據(jù)更新的拉格朗日理論的公式，本研究采用彈塑性有限元法和擴(kuò)展增量的測定方法，最小半徑法，包括元素的產(chǎn)生，節(jié)點(diǎn)接觸或分離的工具，最大應(yīng)變和極限旋轉(zhuǎn)增量。使用改性的庫侖摩擦法建立有限元方法的計(jì)算機(jī)代碼。不同半徑和角度的錐形沖頭是用于模擬成的硬銅及銅合金管的兩端。影響的各種要素，包括沖床半頂角α和半徑R的比值，管的平均直徑，管壁厚度比，力學(xué)性能，并對潤滑管向外卷曲，進(jìn)行調(diào)查。模擬結(jié)果顯示，當(dāng)沖床入口的彎曲半徑滿足的條件ρ≦ρc，卷曲加工出現(xiàn)在管的末端。另一方面，如果ρ≧ρc，管端遭受過擴(kuò)口?？勺儲裞稱為臨界彎曲半徑。ρc的值隨著α的值增加而增加。此外，研究結(jié)果還表明，ρc既不與管的材料，也不與潤滑劑有關(guān)。 關(guān)鍵詞：彈塑性；有限元;半頂角;向外卷曲過程。 1介紹 形成金屬管的凸邊的過程中，經(jīng)常被用于連接兩個(gè)管件，連接和鎖管的一部分，其組成部分，和流體管道連接管的末端相結(jié)合的互補(bǔ)與加固管道。它是一種常見的與管兩端相關(guān)的工業(yè)技術(shù)。攻絲機(jī)擴(kuò)口，法蘭管端進(jìn)程一般分為分為擴(kuò)口，加強(qiáng)擴(kuò)口，卷曲成型。本文所討論的管卷曲變形卷曲管的結(jié)束了一圈向外開放。管卷曲變形的模擬是一項(xiàng)復(fù)雜而艱巨的任務(wù)，因?yàn)樽冃芜^程是高度非線性的。由于非線性變形的特點(diǎn)是： 1．大位移，旋轉(zhuǎn)，變形過程中金屬變形。 2．金屬材料遇到大變形的非線性材料變形行為。 3．由摩擦產(chǎn)生的金屬和工具之間的界面，和他們的接觸條件的非線性邊界條件。 這些特點(diǎn)，使有限元方法，使用最廣泛的金屬工藝分析。改善工藝和提高工業(yè)生產(chǎn)率，本研究開發(fā)的彈塑性有限元計(jì)算機(jī)代碼，使用選擇性降低集成（SRI）的模擬方法，它采用一個(gè)矩形的四個(gè)節(jié)點(diǎn)內(nèi)四個(gè)整點(diǎn)元素。目標(biāo)是到模擬的管卷曲的變形過程。 Nadai在1943年研究管收口。理論歸納為弧形殼體理論的延伸。Nadai假設(shè)摩擦系數(shù)為常數(shù)，而忽略了殼的有效應(yīng)力變化的存在?？唆?shù)呛蜏丈M(jìn)行了一系列管減噪實(shí)驗(yàn)，建立管減噪的各種限制，并評估各種參數(shù)的影響。Manabe and Nishimura兩個(gè)錐形燃燒管和管減噪過程進(jìn)行了一系列實(shí)驗(yàn)，調(diào)查研究的參數(shù)包括不同角度的錐形打孔，潤滑，材料成形負(fù)荷和應(yīng)力應(yīng)變分布的各種參數(shù)的影響管壁厚度。唐和小林提出了剛塑性有限元理論，并開發(fā)出一種計(jì)算機(jī)代碼來模擬105毫米AISI1018鋼冷收口。黃等。模擬冷收口過程中，通過彈塑性和剛塑性有限元法的方法，并比較模擬結(jié)果與實(shí)驗(yàn)數(shù)據(jù)來驗(yàn)證剛塑性有限元分析的準(zhǔn)確性。北澤等人。錐形沖頭，與假設(shè)的銅及黃銅材料進(jìn)行實(shí)驗(yàn)，管材料剛性完美的塑料，他們探討了以前的沖頭圓弧半徑，角度，管壁厚，管卷曲變形過程的影響。北澤采用碳鋼，銅，黃銅的管材管卷曲變形試驗(yàn)。當(dāng)在管端前緣擴(kuò)口的能量增量是大于卷曲邊緣的能量增量。管底料的向外卷曲，是在不屈的情況下形成。另一方面，如果管端材料附著因?yàn)椴磺腻F形打孔表面和失敗卷曲或形成擴(kuò)口。能源規(guī)則被用來誘導(dǎo)管形成的比較和驗(yàn)證的目的時(shí)的變形能，從而確立了向外卷曲的標(biāo)準(zhǔn)。 在這項(xiàng)研究中，在其他研究中使用相同的材料常數(shù)和尺寸都采用向外卷曲過程的模擬結(jié)果相比，報(bào)告的結(jié)果驗(yàn)證了彈塑性有限元計(jì)算機(jī)代碼的準(zhǔn)確性。 2基本理論的說明 2.1剛度方程 采用更新的拉格朗日公式（超低頻）金屬成形過程中的增量變形的應(yīng)用框架（體積成形，板料成形）是描述塑性流動規(guī)則的增量特性的最切實(shí)可行的辦法。在每個(gè)變形階段的當(dāng)前配置中的ULF用于評估一個(gè)小的時(shí)間間隔t變形，一階理論的精度要求是一致的，作為參考狀態(tài)。 更新的拉格朗日方程的虛擬工作速率方程寫成 是基爾霍夫應(yīng)力的Jaumann率，是柯西應(yīng)力，是變形率，這是笛卡爾坐標(biāo)，是速度，的名義牽引速度，是速度梯度，是固定的空間直角坐標(biāo)，V和材料的體積和表面上牽引規(guī)定。 圖1.幾何初始向外卷曲變形管和一個(gè)特定階段的邊界條件。單位毫米;沖床半徑R，半頂點(diǎn)沖角α，沖床入口彎曲半徑ρ。 J2流動規(guī)律的模型采用的彈塑性行為的金屬片，是應(yīng)變硬化率，是有效應(yīng)力，E是楊氏模量，是泊松比，是的偏量。取1為塑性狀態(tài)和0為彈性狀態(tài)或卸載。 據(jù)推測，速度{V}，在離散化元素中的分布 其中[N]是形函數(shù)矩陣和表示節(jié)點(diǎn)的速度。變形率和速度梯度寫成 其中[B]和[E]中，分別代表應(yīng)變率的速度矩陣和速度梯度速度矩陣。將式（4）和（5）代入式（1），得到元素的剛度矩陣。 由于虛功率式的原則。本構(gòu)關(guān)系是非線性方程。差餉，他們可以更換任何單調(diào)增加的措施方面定義的工具，如位移增量，增量。 Alpha有限元的標(biāo)準(zhǔn)程序，以形成整體的整體剛度矩陣， 在其中 在這些方程中，[K]是全球切線剛度矩陣，元素彈塑性本構(gòu)矩陣，表示位移增量，表示在規(guī)定的結(jié)點(diǎn)力增量。和被定義為應(yīng)力校正矩陣，由于在電流應(yīng)力變形的任何階段。 圖2.（a）變形的幾何形狀，（b）節(jié)點(diǎn)速度分布在12個(gè)不同的形成階段管向外卷曲硬銅管，N= 0.05; T0 =0.8毫米;=60擄; r =3.4毫米。 2.2斯里蘭卡方案 由于四邊形元素充分整合計(jì)劃的實(shí)現(xiàn)，導(dǎo)致過度的約束效應(yīng)時(shí)，材料近不可壓縮成型過程中的彈塑性情況[13]，休斯提出應(yīng)變率的速度矩陣分解為擴(kuò)張矩陣和偏差矩陣，即 矩陣中和由傳統(tǒng)的四點(diǎn)一體化集成。當(dāng)材料變形到幾乎不可彈塑性狀態(tài)，修改后的擴(kuò)張矩陣是由一個(gè)點(diǎn)的一體化集成，即必須更換擴(kuò)張矩陣，即 是修改后的應(yīng)變率速度矩陣。將式（8）代入式（7），修改后的應(yīng)變率速度矩陣 明顯的，速度梯度-速度矩陣，取而代之的是修改后的速度梯度–速度矩陣 3數(shù)值分析 在分析模型管卷曲過程的是軸對稱。因此，只有右邊一半的中心軸是考慮過的。部分的有限元是由計(jì)算機(jī)自動處理。由于是從彎曲，并在管的末端卷曲的急劇變形，這一節(jié)需要更精細(xì)的單元劃分，以獲得精確的計(jì)算結(jié)果。如圖左邊一半。1表示零件和模具開始的尺寸。表1給出了詳細(xì)的尺寸。在局部坐標(biāo)系，軸1表示的切線方向，管材料和工具之間的聯(lián)系而軸n為同一聯(lián)系的正常方向。常數(shù)（X，Y）坐標(biāo)和地方坐標(biāo)（l，n）描述結(jié)點(diǎn)力，位移和元素的應(yīng)力和應(yīng)變。 圖3.正如圖2，但R= 3.8毫米。 表1 沖頭的角度和半徑 沖頭頂點(diǎn) 角度α（度） 沖頭半徑R（毫米 60 2.0 2.4 2.8 3.1 3.4 3.8 4.1 4.4 4.8 5.1 65 2.4 2.8 3.1 3.4 3.8 4.1 4.4 4.8 5.5 5.4 70 3.1 3.4 3.8 4.1 4.4 4.8 5.1 5.4 5.8 6.1 75 3.4 3.8 4.1 4.4 4.8 5.1 5.4 5.8 6.1 6.4 80 4.8 5.4 5.8 6.1 6.4 6.8 7.2 7.6 8.0 8.4 85 5.8 6.1 6.4 6.8 7.4 9.0 9.4 9.8 10.2 10.6 11.0 表2.所用材料的機(jī)械性能 材料 外徑×厚度 真空熱處理 n F(Mpa) 屈服應(yīng)力(Mpa) 銅 25.4×0.8 原樣 0.09 380 220 500℃ 1 h 0.53 630 26 原樣 0.05 450 280 300℃ 1 h 0.09 450 260 400℃ 1 h 0.46 610 50 600℃ 1 h 0.50 640 42 70/3黃銅 25.4±0.8 原樣 0.18 730 280 應(yīng)力 ；應(yīng)變 表2給出了我們模擬的物質(zhì)條件。管外半徑保持在25.4毫米不變;但也有三種不同的管壁厚度值，即，0.4，0.6和0.8毫米，在實(shí)驗(yàn)和計(jì)算。硬銅管和黃銅管的泊松比和楊氏模量分別是0.33和110740兆帕，0.34和96500兆帕。 3.1邊界條件 如圖右手半管。1是指管卷曲變形在某一階段的變形形狀。邊界條件包括以下三個(gè)部分： 1.FG和BC段的邊界： 這里是節(jié)點(diǎn)的切向摩擦力增量，是普通的力量增量。假定涉及到摩擦材料，工具的接觸面積，和不等于零。，表示工具的配置文件的正常方向的位移增量，從規(guī)定的沖頭位移增量中確定。 2.邊界上的CD，DE，EF，和GA部分： 上述情況反映了這個(gè)邊界上的節(jié)點(diǎn)都是自由的。 3.AB段的邊界： AB部分是固定在管底的邊界。在此位置沿Y軸方向的節(jié)點(diǎn)位移增量設(shè)置為零，而節(jié)點(diǎn)B是完全固定的沒有任何動靜。 在管卷曲過程中，邊界將會改變。因此，有必要檢查FG和BC沿邊界部分的變形階段，在每個(gè)接觸節(jié)點(diǎn)的法向力。如果達(dá)到零，則節(jié)點(diǎn)將成為自由節(jié)點(diǎn)和邊界條件（1）（2）轉(zhuǎn)移。同時(shí)，GA和EF沿管段的自由節(jié)點(diǎn)還檢查了計(jì)算。如果該節(jié)點(diǎn)接觸到?jīng)_頭，自由邊界條件的約束條件（1）改變。 3.2 彈塑性和接觸問題的處理 當(dāng)解釋以前明確暗示的邊界條件時(shí)，接觸條件應(yīng)維持在一個(gè)增量變形過程。為了滿足這個(gè)要求，山田等人提出了γ最小值方法。[ 14]采用擴(kuò)展處理了彈塑性和接觸問題[ 15]。每個(gè)加載增量的最小值均被控制，有以下六個(gè)標(biāo)準(zhǔn)。 彈塑性狀態(tài)。元素的應(yīng)力大于屈服應(yīng)力時(shí)，R1的計(jì)算方法[15]因此，以確定達(dá)到一樣的表面屈服應(yīng)力。 最大應(yīng)變增量。獲得R2默認(rèn)的最大應(yīng)變增量ψ主應(yīng)變增量比dε，即，限制增長到這樣的規(guī)模，一階理論是有效的措施。 最大旋轉(zhuǎn)增量。默認(rèn)的最大旋轉(zhuǎn)增量β的旋轉(zhuǎn)增量，即，限制增長到這樣的規(guī)模，一階理論是有效的措施。 滲透條件。形成的收益時(shí)，管的自由節(jié)點(diǎn)可能滲透的工具。比[16]的計(jì)算方法等，剛剛接觸到的工具的自由節(jié)點(diǎn)。 分離條件。當(dāng)形成的收益，接觸節(jié)點(diǎn)可以脫離接觸面。 [ 16 ]計(jì)算每個(gè)接觸節(jié)點(diǎn)，這種正常的組成部分，節(jié)點(diǎn)力為零。 滑粘摩擦條件。修正的庫侖摩擦定律提供了選擇的接觸狀態(tài)，即滑動或粘狀態(tài)。這種狀態(tài)檢查每個(gè)接觸節(jié)點(diǎn)由以下條件： （a）如果，還有節(jié)點(diǎn)處于滑動狀態(tài) (b)如果還有則節(jié)點(diǎn)處于固定狀態(tài) 一個(gè)方向滑動節(jié)點(diǎn)是相反的方向前漸進(jìn)的步驟，使接觸節(jié)點(diǎn)粘貼節(jié)點(diǎn)在下一個(gè)增量步。然后在這里得到比R，，產(chǎn)生變化的摩擦狀態(tài)的基體，在是一個(gè)很小的公差。 最大應(yīng)變增量Ψ，最大的旋轉(zhuǎn)增量β這里使用的常數(shù)是0.002o和0.5o，分別。這些常量都被證明是有效的一階理論。。此外，一個(gè)小的檢查程序的滲透和分離條件是允許的。 3.3 卸載過程 卸載后的回彈現(xiàn)象是管形成過程中具有重要意義。假設(shè)是固定的，對管底的節(jié)點(diǎn)卸載程序執(zhí)行。所有的元素都將復(fù)位彈性。節(jié)點(diǎn)來接觸到的工具的力量在扭轉(zhuǎn)成為管規(guī)定力的邊界條件，即 同時(shí)，核查的滲透，摩擦，分離條件被排除在仿真程序。 四．結(jié)果與討論 圖2顯示沖角α=60°，圓弧半徑R=3.4毫米銅管卷曲成形仿真結(jié)果。圖2（a）表示的幾何形狀的變形，在最后的形狀是卸載后的最終形狀的一部分。圖2（b）顯示了變形過程中的節(jié)點(diǎn)的速度分布。最后圖顯示卸貨后節(jié)點(diǎn)的速度分布。這些數(shù)字表明，管端材料處于起步階段順利進(jìn)入沿打孔表面。管底料后，逐漸彎曲向外卷曲，造成所謂向外卷曲成形。圖3顯示的R=3.8毫米的情況下的模擬結(jié)果。管端材料還顯示所謂的燃燒形成合格后，電弧感應(yīng)部分和圓錐面形成。 圖4顯示了α=75°，R=6.1毫米卷曲銅管的模擬的結(jié)果。相比之下，圖。5顯示了擴(kuò)口中的R= 6.4毫米的情況下形成的模擬結(jié)果。 圖4 如圖2條件，而α=75o；R=6.1毫米 圖5 如圖2條件，而α=75o；R=6.4毫米 圖6顯示了應(yīng)變分布相應(yīng)的那些圖2和圖3在沖床行程12毫米。圖7顯示了應(yīng)變分布相應(yīng)的那些圖4和圖5在沖床進(jìn)步14毫米。變形主要是由于拉在圓周方向彎曲，在正北方向和厚度方向的剪切。最大直徑的管端擴(kuò)口成形。因此，其應(yīng)變值在圓周方向遠(yuǎn)高于向外卷曲。減少管壁厚度也更重要的擴(kuò)口成形。因此，在相同的條件下，在管端擴(kuò)口比在其他情況下更容易斷裂。 圖6計(jì)算硬銅管應(yīng)變分布下的卷曲和擴(kuò)口過程比較 n=0.05;to=0.8毫米；α=60o 圖7 條件如圖6，而α=75o 圖8和圖9顯示仿真結(jié)果硬銅和黃銅，分別。同管壁厚，管半徑，潤滑條件，和半頂沖床是用來在模擬。然而，不同的價(jià)值觀的沖頭圓弧半徑的使用，從而產(chǎn)生不同類型的形成，即，擴(kuò)口和卷曲。就是說，如果圓弧半徑，，在沖頭入口大于臨界彎曲半徑，ρ，然后管端的材料形成沿電弧燃燒沖壓部分和圓錐面。反之，如果，管底材料葉片穿孔表面形成卷曲。這一數(shù)字表明臨界彎曲半徑隨半頂沖床增加。比較結(jié)果表明，我們的數(shù)值模擬是一致的實(shí)驗(yàn)數(shù)據(jù)報(bào)告[ 12]。所謂臨界彎曲半徑，ρ可表示如下： 圖8 影響壁厚硬銅管彎曲半徑和沖角之間至關(guān)重要的關(guān)系 n=0.05 圖9 條件如圖8 黃銅管n=0.05 這里 R c=沖頭半徑相應(yīng)的ρ丙 R fmin=最小值所需的耀斑沖壓半徑 R cmax=最大值所需的卷曲沖壓半徑 t0=壁厚管 圖10和圖11表明之間的相關(guān)性無量綱臨界彎曲半徑（ρ丙）和半頂角的2種不同壁厚徑比的情況下，努力按銅管，分別。價(jià)值ρ可以定義如下： 如果ρ是無關(guān)的壁厚徑比，那么幾何相似定律的存在臨界彎曲半徑。事實(shí)判斷出來的曲線顯示在圖10和11是幾乎相同的，我們一定會存在的幾何相似的臨界彎曲半徑管材料。 圖10 圖11 圖12和圖13表明之間的關(guān)系，最大半徑和半頂角的時(shí)候，在數(shù)值模擬的壁厚度值為0.4，0.6，和0.8毫米的情況下硬銅管和銅管分別發(fā)生卷曲。當(dāng)雙方的半徑和相應(yīng)的半頂角的沖落在曲線的條件，將導(dǎo)致擴(kuò)口。另一方面，當(dāng)兩個(gè)值低于曲線的條件，將導(dǎo)致冰壺。這些數(shù)字表明，最大半徑增加的半頂角的擴(kuò)大。 圖12 圖13 圖14顯示了相關(guān)材料的加工硬化指數(shù)（n）和臨界彎曲半徑。無論變化值之間的0.03和0.53，由于不同的管材料，數(shù)值模擬的結(jié)果表明，值的臨界彎曲半徑保持不變，從而形成一條直線。這些調(diào)查結(jié)果是幾乎相同的實(shí)驗(yàn)數(shù)據(jù)表明[ 12 ]，這意味著，臨界彎曲半徑的時(shí)候，卷曲是無關(guān)的加工硬化指數(shù)。 圖14 圖15顯示了影響潤滑的關(guān)系臨界彎曲半徑和半頂角的沖頭。臨界彎曲半徑的增加而增加半頂角的沖頭的情況μ=0.05和μ=0.20。此外，結(jié)果顯示，有一個(gè)微不足道的差異程度臨界彎曲半徑之間的潤滑條件。 圖15黃銅管依賴的臨界彎曲半徑對潤滑條件。 五．總結(jié) 從更新的拉格朗日公式，帶有一個(gè)圓錐形的沖壓工藝模擬管卷曲的彈塑性有限元計(jì)算機(jī)代碼開發(fā)。高非線性的過程是一個(gè)漸進(jìn)的方式考慮與一個(gè)最小半徑技術(shù)是通過限制每個(gè)增量步的線性關(guān)系的大小。 摩擦改性庫侖法開發(fā)是一個(gè)連續(xù)函數(shù)。如前所述，此功能可以處理工具和金屬界面的滑動和粘性很難與普通連續(xù)摩擦模型描述的現(xiàn)象。斯里蘭卡有限元是四節(jié)點(diǎn)四邊形單元中的4個(gè)集成點(diǎn)。然后再加上與大變形有限元分析的SRI進(jìn)一步成功地應(yīng)用于錐形管卷曲形成工藝的分析。根據(jù)上述分析和討論，可以得出以下結(jié)論： 1、一個(gè)關(guān)鍵的彎曲半徑，存在。如果拱半徑ρ在沖頭入口大于臨界彎曲半徑，，然后管端的材料在沖壓拱感應(yīng)部分和錐形面形成擴(kuò)口。另一方面，如果，管底材料從表面和結(jié)果卷曲。彎曲半徑的臨界值隨沖頭的半頂角增加而增加。 2、臨界彎曲半徑，ρ，遵循幾何相似定律。ρ值隨沖床的半頂角的增加而增加。我們也知道，ρ值是不相關(guān)的管材料，潤滑條件或壁厚直徑半徑比。 3、擴(kuò)口圓周應(yīng)變大于卷邊的。擴(kuò)口厚度也有比卷曲減少。因此，擴(kuò)口更容易發(fā)生管端斷裂。相比之下，卷曲的產(chǎn)品可以保證處理的安全，并加強(qiáng)管材料的結(jié)束部分。 4、模具的形狀表示為一個(gè)數(shù)值函數(shù)。因此，在這項(xiàng)研究中開發(fā)的有限元模型，可用于在連續(xù)形狀的任何工具，來模擬擴(kuò)口或卷曲過程中的所有類型。 參考文獻(xiàn) 1、A. Nadai，《塑性應(yīng)力狀態(tài)在彎曲外殼：形成鋼制外殼，機(jī)頭高爆外殼鍛造所需的力量》，在美國機(jī)械工程師學(xué)會年度會議，紐約，美國， 1943.11.23-12.1。 2、A. K. Cruden 和 J. F. Thompson《反向擠塑罐最終的關(guān)閉者》，第511卷，國家工程實(shí)驗(yàn)室，，瑞典，格拉斯哥，1972。 3、K. Manabe and H. Nishimura，《成形載荷和成形極限，錐管收口–收口與擴(kuò)口管研究雜志Ⅰ》，日本社會對技術(shù)的可塑性，23（255），335 -–342頁，1982–4。 4、K. Manabe and H. Nishimura，《實(shí)驗(yàn)應(yīng)力應(yīng)變分布在錐管收口–收口與擴(kuò)口管的研究Ⅱ》雜志，日本社會對技術(shù)的可塑性，23（255），335 -–342頁，1982–4。 5、K. Manabe and H. Nishimura，《應(yīng)力分析和載荷在錐形收口管–收口與擴(kuò)口管的研究Ⅲ》，日本社會的塑性技術(shù)學(xué)報(bào)，23（258），650 -–657頁，1982–7。 6、K. Manabe and H. Nishimura，《錐形沖頭在擴(kuò)口管形成負(fù)載–收口與擴(kuò)口管的研究》雜志五“，日本社會對技術(shù)的可塑性，24（264），47 -52頁 1983–1。 7、K. Manabe and H. Nishimura，《擴(kuò)口管與錐形沖頭的應(yīng)力和應(yīng)變分布–收口與擴(kuò)口管研究》雜志六，日本社會對技術(shù)的可塑性，24（264），pp.276–282，1983–3。 8、M. C. Tang and S. Kobayashi，《收口過程的有限元方法的調(diào)查，第1部分：在室溫下的收口（冷收口）》，《交易應(yīng)用工程學(xué)報(bào)》，104， 305 -–311頁，1982。 9、Y. M. Huang, Y. H. Lu and M. C. Chen，《分析過程中使用的彈塑性和塑性冷去噪方法》，《材料加工技術(shù)雜志》，30， 351 -–380頁，1992。 10、K. Kitazawa and M. Kobayashi 《在管端卷曲變形機(jī)制的實(shí)驗(yàn)研究 - 殼的卷曲–Ⅱ》雜志，日本社會對技術(shù)的可塑性，28（323），pp.1267–1274，1987–12。 11、K. Kitazawa, M. Kobayashi and S. Yamashita，《管端卷曲的基礎(chǔ)能源分析 - 錐狀殼卷曲理論分析Ⅰ》，29（331），845 -–850頁，1988–8。 12、K. Kitazawa，《標(biāo)準(zhǔn)向外卷曲管》，《交易應(yīng)用工程學(xué)報(bào)》，115， 466 -–471頁 1993。 13、J. C. Nagtegaal, D. M. Parks and J. R. Rice,，《準(zhǔn)確的數(shù)值模擬，在全塑范圍內(nèi)的有限元解決方案》，《交易中的計(jì)算方法應(yīng)用于機(jī)械工程》，4， 153 -–177頁，1974。 14、Y. Yamada, N. Yoshimura and T. Sakurai，《塑性應(yīng)力應(yīng)變矩陣及其應(yīng)用的解決彈塑性問題的有限元方法》，《國際機(jī)械科學(xué)學(xué)報(bào)》，10， 343 -–354頁，1968。 15、A. Makinouchi，《有限元建模用于加工鈑金》，在K. Mattiasson，A. Samuelsson, R. D. Wood and O. C. Zienkiewicz， numiform'86程序會議，哥德堡， 327–332頁，1986。 16、Y. M. Huang, Y. H. Lu and A. Makinouchi, 《彈塑性有限元分析V形薄板彎曲》，《材料加工技術(shù)雜志》，35， 129 -–150頁，1992